譚凱, 王琪, 丁開(kāi)華, 李恒, 鄒蓉,聶兆生, 王迪晉, 楊少敏, 喬學(xué)軍
1 中國(guó)地震局地震研究所,地震大地測(cè)量重點(diǎn)實(shí)驗(yàn)室, 武漢 430071 2 中國(guó)地質(zhì)大學(xué)(武漢),地球物理與空間信息學(xué)院,行星科學(xué)研究所, 武漢 430074 3 中國(guó)地質(zhì)大學(xué)(武漢),信息工程學(xué)院, 武漢 430074
?
近場(chǎng)位移數(shù)據(jù)約束的2013年蘆山地震破裂模型及其構(gòu)造意義
譚凱1, 王琪2*, 丁開(kāi)華3, 李恒1, 鄒蓉2,聶兆生1, 王迪晉1, 楊少敏1, 喬學(xué)軍1
1 中國(guó)地震局地震研究所,地震大地測(cè)量重點(diǎn)實(shí)驗(yàn)室, 武漢 430071 2 中國(guó)地質(zhì)大學(xué)(武漢),地球物理與空間信息學(xué)院,行星科學(xué)研究所, 武漢 430074 3 中國(guó)地質(zhì)大學(xué)(武漢),信息工程學(xué)院, 武漢 430074
以往的研究顯示了2013年蘆山MS7.0級(jí)地震發(fā)震斷層的隱伏逆沖斷層基本特征,但是破裂深部細(xì)節(jié)差異較大.本文以近場(chǎng)密集的同震形變數(shù)據(jù)約束蘆山地震破裂面幾何形狀及滑動(dòng)分布, 結(jié)果顯示蘆山地震破裂面具有鏟狀結(jié)構(gòu),上部16 km為43°~50°高角度斷層,深部16~25 km為小于27°的低角度斷層,破裂深度與重定位的余震分布深度一致.破裂分布模型清楚顯示上下兩個(gè)斷層上各有一個(gè)滑動(dòng)幅度大于0.5 m的峰值破裂區(qū),最大滑動(dòng)量1.5 m位于13 km深處.重定位的余震分布基本都落在最大滑動(dòng)量等值線外部庫(kù)侖應(yīng)力增加的區(qū)域.蘆山地震破裂面幾何形狀和滑動(dòng)分布特征與2008年汶川MS8.0級(jí)地震映秀—北川破裂相似,支持龍門(mén)山?jīng)_斷帶發(fā)育大規(guī)模的近水平滑脫層, 是青藏高原東緣地殼縮短增厚、龍門(mén)山擠壓隆升的重要證據(jù).
蘆山地震; 同震形變; 斷層幾何形狀; 破裂滑動(dòng)分布
龍門(mén)山推覆構(gòu)造帶是中國(guó)大陸南北地震構(gòu)造帶中段的重要組成部分,繼2008年汶川MS8.0特大地震后,緊接著2013年龍門(mén)山西南段爆發(fā)蘆山MS7.0級(jí)地震,野外地質(zhì)考察沒(méi)有觀測(cè)到明顯的地表破裂(徐錫偉等,2013).USGS公布的震源機(jī)制解顯示震源深度約12.3 km,傾角33°.趙博等(2013)利用GSN/IRIS臺(tái)網(wǎng)和國(guó)家臺(tái)網(wǎng)及四川省區(qū)域臺(tái)網(wǎng)的波形數(shù)據(jù)對(duì)主震及部分余震進(jìn)行了震源機(jī)制解反演,表明主震為逆沖地震,震源機(jī)制解斷層傾角為45°.王衛(wèi)民等(2013)反演的破裂過(guò)程給出的傾角為38.5°,最大滑動(dòng)量位于10~12 km深度,以逆沖為主向上擴(kuò)展.Jiang等(2014)用GPS連續(xù)站數(shù)據(jù)反演了震中位置、破裂產(chǎn)狀和滑動(dòng)分布,破裂傾角為43°,最大滑動(dòng)位于深度13 km 處.這些結(jié)果最大的差別在于破裂斷層的傾角差異較大,主要原因可能是,遠(yuǎn)場(chǎng)地震波的空間分辨率較低,而近場(chǎng)形變約束數(shù)據(jù)較少,所以一般都假設(shè)單傾角的矩形平面破裂,無(wú)法對(duì)破裂深部的幾何形狀進(jìn)行精細(xì)約束,不能判斷蘆山地震是否具有類似汶川地震的深部滑脫層破裂特征.
本文利用蘆山地震震區(qū)各類GPS大地控制網(wǎng)復(fù)測(cè)資料以及強(qiáng)震臺(tái)網(wǎng)記錄獲取近場(chǎng)50多個(gè)站點(diǎn)同震變形資料,并基于彈性半空間位錯(cuò)模式反演地表位移,對(duì)龍門(mén)山西南段發(fā)震斷層幾何、運(yùn)動(dòng)學(xué)特征進(jìn)行更為準(zhǔn)確的約束.我們的研究結(jié)果證實(shí)龍門(mén)山構(gòu)造帶西南段同樣具有斷坡-滑脫構(gòu)造.
龍門(mén)山構(gòu)造帶南起瀘定、天全,向北東延伸經(jīng)寶興、都江堰、江油、廣元進(jìn)入陜西勉縣一帶,全長(zhǎng)約500 km,寬40~50 km,由走向N45°E、傾向NW、傾角50°~75°的4條逆斷裂疊瓦狀組合而成,自北西往南東分別為龍門(mén)山后山斷裂、中央斷裂、前山斷裂和山前隱伏斷裂(鄧起東等,1994;陳國(guó)光等,2007).跨龍門(mén)山區(qū)的GPS觀測(cè)表明其現(xiàn)今地殼縮短速率在3 mm·a-1以內(nèi),右旋走滑兼逆沖運(yùn)動(dòng)性質(zhì)(張培震等,2008).
龍門(mén)山斷裂帶具有長(zhǎng)期活動(dòng)性.2008年前,龍門(mén)山斷裂帶中、小地震(震級(jí)M<7的地震)活動(dòng)頻繁,但未有發(fā)生過(guò)7級(jí)以上大地震的歷史記載(張勇等,2008).但是探槽開(kāi)挖揭示龍門(mén)山地區(qū)存在多次古地震事件(冉勇康等,2008;陳立春等,2009;史翔,2009;張?jiān)罉蚝屠詈}垼?010;王煥,2011).在過(guò)去的6000年里,龍門(mén)山地區(qū)曾發(fā)生過(guò)三次古地震事件(包括汶川地震),平均復(fù)發(fā)間隔不超過(guò)3000年(Ran et al.,2013),汶川5.12地震之前的事件(倒數(shù)第二個(gè))與汶川MS8.0地震地表變形規(guī)模相當(dāng)(冉勇康等,2008),測(cè)年數(shù)據(jù)表明該事件發(fā)生在2300~3300 aBP(Ran et al.,2010),倒數(shù)第三個(gè)事件的發(fā)生時(shí)間晚于5920~5730 aBP(Ran et al.,2013).另有學(xué)者指出,在接近6000 aBP的一次古地震事件中,中央斷裂帶和前山斷裂帶曾同時(shí)活動(dòng)(江娃利等,2009).
2008至2013年不到5年間,龍門(mén)山斷裂中北段和南段相續(xù)爆發(fā)汶川MS8.0級(jí)和蘆山MS7.0級(jí)地震,在周邊引起可觀的地表變形,為研究龍門(mén)山構(gòu)造帶結(jié)構(gòu)、變形特征和地震活動(dòng)性提供了豐富的地殼形變信息.
2008年5月12日,四川汶川MS8.0級(jí)地震在地表錯(cuò)開(kāi)了龍門(mén)山中央斷裂和前山斷裂.地質(zhì)野外考察表明,灌縣—江油斷裂的地表破裂帶長(zhǎng)約72 km,純逆沖性質(zhì).映秀—北川斷裂的地表破裂帶長(zhǎng)約240 km,基本呈直線展布,從映秀的逆沖逐漸過(guò)渡到青川的走滑運(yùn)動(dòng)(徐錫偉等,2008).GPS、三角測(cè)量等多源數(shù)據(jù)獲得的近場(chǎng)密集同震形變顯示,從映秀到馬爾康,同震形變由最大5.3 m 減低到0.2 m,地表破裂帶160 km 以遠(yuǎn)變形依然醒目(楊少敏等,2012).譚凱等(2013)以同震形變?yōu)榧s束,用格網(wǎng)搜索法結(jié)合最小二乘法論證了破裂上部?jī)A角較高,而下部?jī)A角較為平緩的特征,表明龍門(mén)山中北段深部滑脫層的存在是青藏高原擠壓隆升的證據(jù).
2013年蘆山地震后,“中國(guó)大陸構(gòu)造環(huán)境監(jiān)測(cè)網(wǎng)絡(luò)”(簡(jiǎn)稱陸態(tài)網(wǎng)絡(luò))(李強(qiáng)等,2012)、“四川GPS連續(xù)觀測(cè)網(wǎng)絡(luò)”(廖華等,2013)以及在蘆山周邊加密布設(shè)的GPS連續(xù)站(Jiang et al.,2014),為同震形變計(jì)算提供了可靠的數(shù)據(jù)(圖1),GPS數(shù)據(jù)用GIPSY 軟件處理,數(shù)據(jù)處理方法見(jiàn)楊少敏等(2012)文獻(xiàn),水平位移中誤差在3 mm 以內(nèi).為增強(qiáng)同震位移資料約束,我們對(duì)蘆山震區(qū)的“陸態(tài)網(wǎng)絡(luò)”區(qū)域站和國(guó)家GPS大地控制網(wǎng)C級(jí)點(diǎn)約20多點(diǎn)進(jìn)行復(fù)測(cè),復(fù)測(cè)工作于震后一年內(nèi)完成,每站總計(jì)觀測(cè)12~36 h.在綜合考慮震間和震后形變的基礎(chǔ)上,獲得流動(dòng)觀測(cè)站同震位移場(chǎng),其水平位移中誤差在5 mm以內(nèi).另外,采用類似的零線校正方法(Iwan et al.,1985),對(duì)近場(chǎng)8臺(tái)強(qiáng)震儀記錄的加速度波形曲線進(jìn)行二次積分獲得同震靜態(tài)形變.與附近的GPS位移比較,兩者的大小和方向基本一致,估計(jì)強(qiáng)震儀水平位移中誤差約10 mm 左右.后續(xù)反演中,GPS連續(xù)站和流動(dòng)站數(shù)據(jù)用實(shí)際中誤差定其權(quán)重,強(qiáng)震儀水平位移按2倍中誤差降權(quán)處理.強(qiáng)震儀位移垂直精度較差,舍棄不用.總共新增約30個(gè)近場(chǎng)站點(diǎn)的同震靜態(tài)水平位移觀測(cè)值.
同震形變主要集中在邛崍與天全間方圓80~100 km的范圍內(nèi),最大形變值67 mm位于距離震中16 km 的蘆山站,距離震中約29 km 的天全站向南西位移21 mm.位移場(chǎng)運(yùn)動(dòng)方向顯示破裂以龍門(mén)山前山隱伏斷裂為分界,通過(guò)震中的斷裂縱剖面上的位移方向呈現(xiàn)垂直斷裂的對(duì)沖方式,上盤(pán)測(cè)站向南東運(yùn)動(dòng)指向震中,而下盤(pán)測(cè)站向北西運(yùn)動(dòng),顯示明顯的破裂逆沖方式.沿?cái)嗔炎呦驈恼鹬型鶅蓚?cè),下盤(pán)測(cè)站位移基本都垂直斷裂向北西運(yùn)動(dòng),但是上盤(pán)測(cè)站位移方向很快離開(kāi)垂直斷裂方向而向外側(cè)的南西或北東偏移,有可能這幾個(gè)站點(diǎn)已經(jīng)超出破裂斷層范圍,或者破裂方式發(fā)生了變化.(圖1)
圖1 蘆山地震構(gòu)造背景與同震形變場(chǎng)藍(lán)箭頭:陸態(tài)網(wǎng)絡(luò)和四川GPS連續(xù)站形變;紅箭頭:加密觀測(cè)站形變;黑箭頭:Jiang等(2014)測(cè)站形變;綠箭頭:強(qiáng)震儀測(cè)站形變;F01:大邑隱伏斷裂;F02:隱伏斷裂;F11:龍門(mén)山前山斷裂的大川—雙石斷裂;F12:灌縣—江油斷裂; F21:龍門(mén)山中央斷裂的鹽井—五龍斷裂;F22:映秀—北川斷裂;F31:龍門(mén)山后山斷裂的耿達(dá)—隴東斷裂;F32:汶川—茂汶斷裂;灰曲線:斷層.Fig.1 The tectonic setting and coseismic deformation of the Lushan earthquakeBlue arrow: deformations of Crustal Movement Observation Network of China and Sichuan GPS continusly operating reference stations;red arrow: deformations of repeated measure stations;black arrow:deformations of Jiang′s stations (2014); green arrow:deformations of seismic strong motion stations; F01:Dayi fault; F02:a blind fault; F11:Dachuan-Shuangshi fault; F12:Guanxian-Jiangyou fault; F21:Yanjing-Wulong fault; F22:Yingxiu-Beichuan fault; F31:Gengda-Longdong fault; F32:Wenchuan-Maowen fault;gray curve:fault.
一般假設(shè)地震區(qū)為彈性空間,則矩形斷層位錯(cuò)引起的地表變形可用彈性半空間(Okada,1985)或地殼分層(Wang et al.,2003)的彈性位錯(cuò)公式計(jì)算,主要與斷層幾何特征7參數(shù)(長(zhǎng)、寬、深、走向、傾向、水平投影坐標(biāo))和滑動(dòng)參數(shù)(走滑量、傾滑量、張性分量)有關(guān),以形變觀測(cè)值為約束可以反演斷層破裂參數(shù).如果假設(shè)破裂由少數(shù)幾個(gè)面積比較大的斷裂組成,則反演的目的主要是求得破裂的大體幾何形狀和平均意義的滑動(dòng)量.如果將破裂細(xì)分為更多的子斷層,則反演的目的是為了獲得精細(xì)的破裂滑動(dòng)分布.兩種反演方法都應(yīng)該滿足觀測(cè)數(shù)據(jù)擬合度和滑動(dòng)分布粗糙度最小(minimum),即
‖W(Gs-d)‖2+β2‖Ls‖2=min
(1)
這里d是形變觀測(cè)值,W是觀測(cè)值的權(quán)矩陣,是觀測(cè)值方差-協(xié)方差D的逆矩陣,即D=W-1W,G是格林函數(shù),s是子斷層滑動(dòng)矢量,β是平滑因子, 單位為m-1,L是拉普拉斯有限差分算子.
地表位移是斷層幾何形狀參數(shù)和滑動(dòng)參數(shù)的非線性函數(shù),可通過(guò)非線性優(yōu)化方法求解.如果假定斷層幾何形狀參數(shù)已知,則地表位移是滑動(dòng)參數(shù)的線性函數(shù),可用線性反演方法求解滑動(dòng)參數(shù),滑動(dòng)分布反演結(jié)果主要取決于破裂面先驗(yàn)幾何特征.蘆山地震破裂尚未出露地表,斷裂位置、走向、長(zhǎng)度、傾角、寬度等具有較大的不確定性.本文先用非線性反演的模擬退火算法,確定斷層幾何形狀參數(shù)及其置信區(qū)間,即將斷裂用數(shù)個(gè)斷裂段表示,用模擬退火算法確定破裂面的幾何形狀,獲得合理的、最優(yōu)的幾何形狀參數(shù),如斷裂傾角等,并估算其置信區(qū)間;然后,在最佳幾何形狀模型的基礎(chǔ)上,將破裂面離散化為更多的子斷層,用線性反演方法,如非負(fù)最小二乘法確定最終的破裂滑動(dòng)分布模型.
3.1 破裂面幾何形狀參數(shù)及置信區(qū)間
本節(jié)用幾個(gè)必要的矩形斷裂表征破裂面的幾何形狀,用模擬退火算法求解其幾何形狀參數(shù)和平均滑動(dòng)量,因而用彈性半空間矩形位錯(cuò)模型(Okada,1985)計(jì)算形變格林函數(shù)即可滿足精度要求.目標(biāo)函數(shù)采用(1)式的觀測(cè)數(shù)據(jù)擬合度和滑動(dòng)分布粗糙度之和,其中用平滑因子β來(lái)平衡兩者的權(quán)重.根據(jù)交叉驗(yàn)證的殘差平方和(CVSS)選取平滑因子(MatthewsandSegall,1993),其原理是:平滑因子在某個(gè)范圍內(nèi)依次取值;對(duì)于每個(gè)選取好的平滑因子,將每個(gè)測(cè)站依次從n個(gè)數(shù)據(jù)里排除,用剩下的n-1個(gè)測(cè)站反演最佳模型,并用該模型計(jì)算被漏掉的那個(gè)數(shù)據(jù)的模擬殘差,循環(huán)所有測(cè)站獲得全部數(shù)據(jù)的殘差平方和CVSS;CVSS最小值對(duì)應(yīng)的平滑因子就是合理的平滑因子.
首先嘗試假設(shè)蘆山地震破裂面是單個(gè)矩形斷層.由于地質(zhì)調(diào)查沒(méi)有發(fā)現(xiàn)破裂出露地表,所以需要在較大的范圍內(nèi)尋找破裂走向的最優(yōu)解.根據(jù)平滑因子與CVSS曲線選取β=3,反演獲得破裂最優(yōu)模型的走向約為212°,與地質(zhì)調(diào)查的龍門(mén)山南段走向接近.破裂長(zhǎng)度約16.8 km,寬度約24.7 km,以逆沖為主,兼有少量走滑.觀測(cè)值殘差均值為5.0 mm.
然后假設(shè)破裂面由上下相接的2個(gè)矩形斷裂段組成,選取β=3,最優(yōu)斷層模型走向?yàn)?11°,上層傾角大于40°,下層斷裂傾角小于30°,以逆沖為主,殘差均值為 4.7 mm,比單斷裂模型提高0.3 mm.
誤差影響的觀測(cè)值反演獲得的結(jié)果都有可能是破裂模型的最優(yōu)解,以此可以估算模型參數(shù)的置信區(qū)間.令某個(gè)觀測(cè)值的中誤差為mi,將所有觀測(cè)數(shù)據(jù)加上(-2mi,+2mi)之間的一個(gè)隨機(jī)數(shù),作為觀測(cè)值誤差干擾,分別重復(fù)反演單斷裂模型、雙斷裂模型1000次.然后將各參數(shù)結(jié)果按大小排列,去掉最小的2.5%和最大的2.5%,得到每個(gè)參數(shù)的95%置信區(qū)間(表1,表2).圖2是將單斷裂模型的10個(gè)參數(shù)配對(duì)繪制的參數(shù)分布,圖3是將雙斷裂模型的10個(gè)主要參數(shù)(第一斷裂段的長(zhǎng)、寬、深、傾角、走向、走滑量、第二斷裂段的寬、深、傾角、走滑量)配對(duì)繪制的參數(shù)分布.
從參數(shù)分布圖可以到斷裂走向、位置的變化.在圖2的單斷裂模型參數(shù)分布圖中,走向在207°~215°左右,傾角約39°~42°,下邊界深度約20~24 km,走滑量很小.雖然斷裂滑動(dòng)量變動(dòng)較大,但基本呈現(xiàn)逆沖的特征,逆沖量與破裂長(zhǎng)度成反相關(guān).長(zhǎng)度取15 km,則逆沖量有1 m;長(zhǎng)度取20 km,則逆沖量只有0.8 m.
圖2 單斷裂模型參數(shù)分布length:長(zhǎng);width:寬; depth:深;dip:傾角;strike:走向;east:X坐標(biāo);north:Y坐標(biāo);strike slip:走滑量;dip slip:傾滑量;open slip:張性分量.Fig.2 One-fault-model parameters distribution
圖3 雙斷裂模型參數(shù)分布1:雙斷裂上層;2:下層.Fig.3 Two-fault-model parameters distribution
表1 單斷裂模型參數(shù)的搜索上下界、最優(yōu)模型和95%置信區(qū)間Table 1 Searching ranges, optimum models and confidence intervals of one-fault-model parameters
注:斷層X(jué)、Y坐標(biāo)定位于斷層下邊界中點(diǎn).深度是下邊界埋深.坐標(biāo)原點(diǎn)的地理坐標(biāo)是103.18°E,30.21°N.
表2 雙斷裂模型參數(shù)的搜索上下界、最優(yōu)模型和95%置信區(qū)間Table 2 Searching ranges, optimum models and confidence intervals of two-fault-model parameters
注: N指該參數(shù)不獨(dú)立,沒(méi)有給出其搜索范圍進(jìn)行搜索,而是根據(jù)其他參數(shù)直接計(jì)算.
雙斷裂模型(圖3)與單斷裂模型具有相似的特征,斷裂走向和位置與單斷裂模型差異不大,走向約205°~214°.總深度與單斷層模型的基本一致,但分辨出兩個(gè)傾角,上層傾角大于43°,而下層傾角小于30°.長(zhǎng)度主要與上層逆沖量成反比相關(guān),似乎主要由上層逆沖量控制,而與下層逆沖量的相關(guān)性不是很明顯.上層破裂平均滑動(dòng)量大于0.5 m,破裂長(zhǎng)度為14~20 km,與單斷裂模型一致.可以分辨出下層具有大于0.3 m 的逆沖量.
基于雙斷裂模型上陡下緩特征、以及能量最小原理考慮,破裂面傾角也可能從地表的高傾角連續(xù)變化到近水平傾角狀態(tài),但是這種復(fù)雜模型未知參數(shù)多,需要更多的先驗(yàn)假設(shè),降低了形變觀測(cè)數(shù)據(jù)的分辨力,使反演難以收斂于穩(wěn)定解.由于雙斷裂模型幾何參數(shù)可以被形變觀測(cè)數(shù)據(jù)穩(wěn)健分辨,反演結(jié)果穩(wěn)定,并且與地質(zhì)調(diào)查和震源機(jī)制解比較吻合,我們認(rèn)為雙斷裂模型反映了破裂基本特征,可以做為后續(xù)研究的基礎(chǔ) .
3.2 破裂滑動(dòng)分布
本節(jié)假設(shè)前面的破裂幾何形狀參數(shù)為已知,將斷層分割為更多數(shù)目的子斷層,以獲得更精細(xì)的滑動(dòng)分布.以雙斷裂最優(yōu)模型為基礎(chǔ),固定斷裂的走向和傾角,考慮破裂可能超出雙斷裂模型范圍,新建模型的破裂面有意向兩側(cè)和深處延長(zhǎng)并超過(guò)余震分布范圍(長(zhǎng)54 km,寬56 km),下層斷裂長(zhǎng)度、走向與上層斷裂一致.最后新建破裂面都被離散化為約1.3×1.1 km的2142個(gè)矩形子斷層.
如此則地表位移是子斷層滑動(dòng)參數(shù)的線性函數(shù),遵循(1)式觀測(cè)數(shù)據(jù)擬合度和滑移分布粗糙度最小的原則,采用非負(fù)最小二乘法解算所有子斷層的走滑量和傾滑量.經(jīng)試算證實(shí):相同的平滑因子約束下,半空間的彈性位錯(cuò)模型與地殼分層的彈性位錯(cuò)模型反演的滑動(dòng)分布基本一致;如果平滑因子取值在7~14內(nèi),彈性半空間模型與地殼分層模型都可以明顯分辨出淺部和深部各有一個(gè)滑動(dòng)峰值區(qū),其中最大破裂都位于12~14 km深處.為簡(jiǎn)單起見(jiàn),采用彈性半空間位錯(cuò)模型計(jì)算地殼形變格林函數(shù).
Du等(1992)以分辨率矩陣R來(lái)評(píng)估模型的可靠性,其對(duì)角線元素rj(j=1,2,…,2n,n是子斷層數(shù))反映子斷層滑動(dòng)分量估算強(qiáng)度.如果該分量能夠完全分辨,其值為1;完全不能分辨,其值為0;通常對(duì)角元素 0 圖4 滑動(dòng)分布模型平滑因子與CVSS曲線Fig.4 Smooth factor VS. CVSS of slip distribution models 圖5 同震滑動(dòng)分布模型Fig.5 Coseismic slip distribution model (2) 這里α代表子斷層的線性尺度,在本模型為1.3 km.ρ代表了模型能夠可靠分辨的滑動(dòng)分布尺度大小.分析可知,模型上層十多公里深處分辨率達(dá)7 km左右,下層二十多公里深處分辨率也在10 km左右.因此可以認(rèn)為,模型顯示上層大于1 m的主滑動(dòng)區(qū),以及下層大于0.5 m的滑動(dòng)區(qū)是可靠的. 3.3 模型比較 以往基于遠(yuǎn)震地震波或大地測(cè)量資料反演蘆山地震破裂模型,揭示了破裂的基本特征.這些模型顯示在震源附近二十多公里范圍內(nèi),破裂以逆沖為主,尚未出露地表(表3). 表3所列以地震波資料約束的破裂模型假設(shè)震源機(jī)制的破裂幾何參數(shù)和震源位置為已知.由于各家采用的震源機(jī)制來(lái)源不同,預(yù)設(shè)的破裂走向、傾角、震源位置差異較大.這些破裂模型顯示了破裂滑動(dòng)的量級(jí),以及逆沖為主、滑動(dòng)區(qū)域較小等基本特征,滑動(dòng)分布的空間位置不具備可比性.張勇等(2013)假設(shè)USGS最早發(fā)布的破裂斷層走向219°、傾角33°、震中位置和深度為已知,反演結(jié)果顯示震級(jí)約MW6.8,兩次破裂子事件都發(fā)生在破裂起始點(diǎn)附近,集中在一個(gè)長(zhǎng)寬皆約為25 km的滑動(dòng)區(qū)域,最大滑動(dòng)量約1.3 m,破裂沒(méi)有大規(guī)模出露地表.趙翠萍等(2013)參考中國(guó)地震局地球物理研究所劉超等的矩張量反演結(jié)果(http:∥www.cea-gp.ac.cn/tpxw/266824.shtml), 假設(shè)斷層走向220°、傾角35°、震中位置和深度為已知,結(jié)果顯示震級(jí)MW6.8,破裂滑動(dòng)集中在起始破裂點(diǎn)之上40 km×30 km 的區(qū)域, 最大滑動(dòng)量達(dá)1.8 m.而王衛(wèi)民等(2013)綜合震源機(jī)制解以及發(fā)震斷層地質(zhì)構(gòu)造背景,假設(shè)走向205°、傾角38.5°為已知,結(jié)果顯示震級(jí)MW6.7,最大滑動(dòng)量1.6 m. 表3 蘆山地震破裂參數(shù)Table 3 Rupture parameters of the Lushan earthquake 相對(duì)而言,近場(chǎng)大地測(cè)量資料具有更好的空間分辨率.基于GPS和InSAR數(shù)據(jù)反演的破裂斷層幾何形狀參數(shù)(Jiang et al.,2014;劉云華等,2014)具有較好的一致性,破裂斷層走向、傾角、甚至破裂位置和深度的差異都較小,走向都在208°~212°間,傾角大于40°,最大破裂滑動(dòng)的深度約13 km,矩震級(jí)MW6.6與GCMT結(jié)果一致.本文不但分辨出上部破裂主要傾角46.8°,還分辨出下部16~25 km處破裂傾角低于26.5°,有可能是連接龍門(mén)山深部滑脫層的過(guò)渡斷層,也可能屬于滑脫層的一部分.上部16 km可能是彈性閉鎖層,釋放地震矩0.70×1019N·m, 占總地震矩的61.7%,相當(dāng)于矩震級(jí)MW6.50;下部低角度滑脫層的峰值滑動(dòng)深22 km,釋放地震矩0.44×1019N·m ,占總地震矩38.3%,相當(dāng)于震級(jí)MW6.37. 地震波反演的破裂過(guò)程(張勇等,2013;趙翠萍等,2013)也表明蘆山地震由2次破裂事件組成. 本文與地震波資料給出的最大滑動(dòng)量都超過(guò)1 m,而Jiang 等(2014)和劉云華等(2014)給出的結(jié)果都小于1 m.Jiang等(2014)給出的精細(xì)滑動(dòng)分布模型最大滑動(dòng)量為0.6 m.從物理意義上理解,均勻滑動(dòng)模型滑動(dòng)量應(yīng)該是所有滑動(dòng)量的大概平均,應(yīng)在0~0.6 m 間,但其給出的均勻滑動(dòng)量為0.7 m.觀測(cè)數(shù)據(jù)的空間分布、約束條件或平滑條件的影響,都可能影響均勻滑動(dòng)模型或滑動(dòng)分布模型的破裂面積大小,從而影響均勻滑動(dòng)量或者最大滑動(dòng)量的大?。畡⒃迫A等(2014)沒(méi)有分析平滑因子對(duì)滑動(dòng)分布的影響,但均勻滑動(dòng)模型滑動(dòng)量0.94 m 略大于滑動(dòng)分布模型最大滑動(dòng)量0.91 m,而且均勻滑動(dòng)模型矩震級(jí)與滑動(dòng)分布模型矩震級(jí)不一致,說(shuō)明至少有一個(gè)結(jié)果沒(méi)有反映約束數(shù)據(jù)所代表的地震能量. 破裂滑動(dòng)分布結(jié)果可與余震活動(dòng)分布相互驗(yàn)證.為此,根據(jù)本文最優(yōu)破裂模型計(jì)算蘆山地震靜態(tài)庫(kù)侖應(yīng)力(Toda et al.,2005)分布,將靜態(tài)庫(kù)侖應(yīng)力、破裂滑動(dòng)分布、余震分布分別投影到水平面(圖6a)和垂直剖面(圖6b)上,余震分布采用陳晨和胥頤(2013)重定位結(jié)果.本文有滑動(dòng)的最大深度是25~26 km,與余震重定位結(jié)果的最大深度一致(圖6b).大部分余震分布在5~24 km 深度內(nèi),與本文模型滑動(dòng)大于0.2 m 的區(qū)域一致,與房立華等(2013)、Jiang等(2014)的結(jié)果基本相同.但是,本文結(jié)果顯示,余震分布基本都落在最大滑動(dòng)量等值線外部、庫(kù)侖應(yīng)力增加的區(qū)域,與2008年汶川地震破裂類似,余震主要分布在高滑動(dòng)區(qū)外(Wang et al.,2011),而與Jiang 等(2014)展示的余震分布與高滑動(dòng)區(qū)相對(duì)應(yīng)的結(jié)果不同. 本文解算標(biāo)量地震矩與GCMT估算的標(biāo)量地震矩最接近,略大于Jiang等(2014)、劉云華等(2014)用大地測(cè)量資料反演得到的結(jié)果,比張勇等(2013)、王衛(wèi)民等(2013)、趙翠萍等(2013)利用遠(yuǎn)震地震波計(jì)算結(jié)果都?。惹肮嫉臄?shù)個(gè)模型,由于受到遠(yuǎn)場(chǎng)資料空間分辨率的制約,或者近場(chǎng)資料的稀少,都只能分辨出單一傾角斷層,以及一個(gè)比較集中的滑動(dòng)峰值區(qū).本文模型分辨出斷層由上部高傾角斷層下接低傾角斷層組成,與汶川地震破裂的鏟形幾何形狀類似. 圖7清晰地展示了汶川地震(譚凱等,2013)、蘆山地震破裂模型(本文結(jié)果)和地質(zhì)調(diào)查結(jié)果.汶川地震和蘆山地震均具有類似的鏟形破裂結(jié)構(gòu),深部低角度破裂深度都從16 km延伸到23~25 km,兩個(gè)模型低角度破裂頂部的水平投影大致平行,說(shuō)明龍門(mén)山斷裂南段與中北段具有類似的深部低角度滑脫結(jié)構(gòu).在汶川地震中,深部破裂幅度最大可達(dá)6~7 m,斷坡-滑脫斷層鏟形結(jié)構(gòu)的觀測(cè)值擬合度比高角度斷坡斷層的殘差平方和至少低16~17%,尤其是對(duì)斷層西南部觀測(cè)值的擬合.蘆山地震位于龍門(mén)山斷裂西南段,其深部低角度破裂幅度大于0.5 m,可能在中央斷裂的鹽井—五龍斷裂西北側(cè)(包括后山的耿達(dá)—隴東斷裂)的深部逐漸歸入近水平滑脫層.類似結(jié)構(gòu)在1999年臺(tái)灣集集MW7.6地震中也有發(fā)現(xiàn),其破裂模型也延伸至深部近水平的滑脫層(Johnson and Segall,2004),滑脫層上的滑動(dòng)量最大達(dá)2~4 m. 近年的川滇地區(qū)測(cè)震研究提供了關(guān)于深部物質(zhì)結(jié)構(gòu)更精細(xì)的資料.鄭勇等(2013)用2007—2009年區(qū)域?qū)掝l帶測(cè)震數(shù)據(jù),使用背景噪聲層析成像方法和遠(yuǎn)震接收函數(shù)法,獲得了汶川震區(qū)及周邊地區(qū)精細(xì)的S波速度結(jié)構(gòu).在淺部約5 km深度,四川盆地一側(cè)為低速區(qū),松潘甘孜地區(qū)為高速區(qū);在20 km 以下,四川盆地地區(qū)變成高速區(qū),而松潘甘孜地區(qū)變?yōu)榈退賲^(qū).深度13~17 km正好是波速變化劇烈的地區(qū),其下底邊界與本文模型深部滑脫層位置基本一致,破裂向上以高角度出露地表,與地殼厚度和莫霍面深度陡變位置對(duì)應(yīng).巖石成分變化和部分熔融可能導(dǎo)致波速比的變化,使得松潘甘孜上部具備與下部解偶滑脫的構(gòu)造條件;13~17 km深度下四川盆地的波速值逐漸高于松潘甘孜地塊,導(dǎo)致龍門(mén)山斷裂帶東側(cè)相對(duì)于西側(cè)表現(xiàn)的更加堅(jiān)固和穩(wěn)定.因此推測(cè),巴顏喀拉塊體向東沿滑脫層擠出,遇到堅(jiān)硬的四川盆地基底的阻擋積累了足夠的地震能量,爆發(fā)了汶川地震和蘆山地震. 圖7 蘆山地震(本文結(jié)果)、汶川地震(譚凱等,2013)破裂模型比對(duì)F01:大邑隱伏斷裂;F02:隱伏斷裂;F11:龍門(mén)山前山斷裂的大川—雙石斷裂;F12:灌縣—江油斷裂;F13:江油—廣元斷裂;F21:龍門(mén)山中央斷裂的鹽井—五龍斷裂;F22:映秀—北川斷裂;F23:茶壩—林庵寺斷裂;F31:龍門(mén)山后山斷裂的耿達(dá)—隴東斷裂;F32:汶川—茂汶斷裂;F33:青川—平武斷裂;虛線框:汶川和蘆山地震破裂模型邊界的水平投影;黑曲線:汶川地震地表破裂;灰曲線:斷層.Fig.7 The Lushan earthquake rupture model (in this paper) comparation with the one of the Wenchuan earthquake (Tan et al.,2013)F01:Dayi fault;F02:a blind fault;F11:Dachuan-Shuangshi fault;F12:Guanxian-Jiangyou fault; F13: Jiangyou-Guangyuan fault; F21:Yanjing-Wulong fault;F22:Yingxiu-Beichuan fault;F23:Chaba-Lin′ansi fault;F31:Gengda-Longdong fault;F32:Wenchuan-Maowen fault;F33:Qingchuan-Pingwu fault; dashed line: the fault model border of the Wenchuan and Lushan earthquake;black curve:the Wenchuan earthquake surface rupture;gray curve:fault. 跨青藏高原東緣及四川盆地的大地電磁測(cè)深剖面研究表明(張樂(lè)天等,2012),松潘—甘孜地塊淺部高阻、中下地殼低阻,四川盆地淺部低阻、中下地殼相對(duì)高阻.龍門(mén)山斷裂帶在電性結(jié)構(gòu)上表現(xiàn)為傾角較緩、北西傾向的逆沖高導(dǎo)低阻體,向青藏高原內(nèi)部延伸到20 km深度左右,與我們推測(cè)的滑脫層位置一致.其高導(dǎo)電性的成因可能是破碎帶內(nèi)的流體造成的.龍門(mén)山斷裂帶下方中下地殼范圍內(nèi)的高阻體很好地體現(xiàn)了揚(yáng)子地塊所應(yīng)具有的較為穩(wěn)定的高阻特征,說(shuō)明揚(yáng)子地塊基底在深部楔入龍門(mén)山下部. 李志偉等(2011)用四川地震臺(tái)網(wǎng)2000—2008年的地震波數(shù)據(jù)反演獲得的龍門(mén)山及鄰區(qū)地殼P波速度結(jié)構(gòu)與鄭勇等(2013)的結(jié)果類似,并且沿龍門(mén)山斷裂帶,地殼上部10 km 的蘆山段低速異常與映秀—北川段的低速異常被中間的高速異常所分隔,15 km 以下剛好相反,蘆山段和映秀—北川段的高速異常被中間的低速異常分開(kāi).南北差異的速度分布結(jié)構(gòu)說(shuō)明了介質(zhì)性質(zhì)的橫向變化,有利于地震能量分別積累于蘆山段和映秀—北川段15 km 深度. 蘆山地震和汶川地震破裂模式無(wú)疑給出深部滑脫層持續(xù)活動(dòng)的觀測(cè)證據(jù),進(jìn)一步支持地殼擠壓縮短增厚作為龍門(mén)山隆升的主要機(jī)制.龍門(mén)山推覆構(gòu)造前緣地帶代表了典型的中生代前陸盆地(許志琴等,1992),中地殼存在一個(gè)低角度的基底滑脫面,下部基底不變形.滑脫面在龍門(mén)山一側(cè)逐漸陡立,形成斷坡,呈現(xiàn)大致平行的疊瓦狀沖斷帶.沉積蓋層水平縮短、褶皺、增厚變形,山體抬升(許志琴等,1992),總體表現(xiàn)為薄皮構(gòu)造.與此對(duì)立的觀點(diǎn)是厚皮構(gòu)造,高角度深大沖斷帶的差異性垂直運(yùn)動(dòng)調(diào)節(jié)了青藏高原東緣相對(duì)四川盆地的抬升,但本文研究顯示,龍門(mén)山下部存在切割中下地殼高角度斷層的可能性不大,下地殼的層流不可能是龍門(mén)山隆升的唯一機(jī)制. 汶川地震撕裂了龍門(mén)山前山的灌縣—江油斷裂和中央的映秀—北川斷裂,說(shuō)明前山的灌縣—江油斷裂依然是青藏高原東緣與四川盆地的動(dòng)力學(xué)分界面.從地形一般可以推斷,同屬前山斷裂的大川—雙石斷裂可能曾經(jīng)也是高原東緣的動(dòng)力學(xué)分界面.然而,地質(zhì)調(diào)查并沒(méi)有發(fā)現(xiàn)蘆山地震在地表有明顯破裂,只是在邛崍高何鎮(zhèn)、雅安下里鎮(zhèn)有地面隆起(徐錫偉等,2013),與蘆山地震破裂模型的地表跡線一致.這暗示龍門(mén)山西南段與四川盆地的動(dòng)力學(xué)分界面已經(jīng)推移到高何—下里的隱伏斷裂,蘆山地震發(fā)震斷層與屬于前山斷裂的大川—雙石斷裂關(guān)系不大.蘆山地震模型低角度破裂頂部的水平投影幾乎與前山斷裂的大川—雙石斷裂重合,與中央斷裂的映秀—北川破裂平行,說(shuō)明大川—雙石斷裂深部與映秀—北川斷裂深部動(dòng)力構(gòu)造更為接近. 地質(zhì)調(diào)查也顯示龍門(mén)山斷裂帶中段發(fā)育的3條主要活動(dòng)斷裂和西南段發(fā)育的3條活動(dòng)斷裂不是完全一一對(duì)應(yīng)(張?jiān)罉蚝屠詈}垼?010),與本文推斷較為一致.雖然余震重定位結(jié)果都顯示余震震中主要沿雙石—大川斷裂地表出露位置的兩側(cè)分布(陳晨和胥頤,2013;房立華等,2013;張廣偉和雷建設(shè),2013),但是根據(jù)震源深度剖面推測(cè)蘆山地震的發(fā)震構(gòu)造可能為大川—雙石東側(cè)約10 km的一條盲逆沖斷層(陳晨和胥頤,2013;房立華等,2013),與大川—雙石斷裂直接相關(guān)(張廣偉和雷建設(shè),2013)的可能性不大. 因此,汶川地震和蘆山地震可能是在相同滑脫結(jié)構(gòu)的擠壓造山作用下,在相鄰區(qū)域引起的不同地震事件.汶川地震的發(fā)生對(duì)蘆山地震有觸發(fā)作用和解鎖作用,但不是主余震的關(guān)系.并且反演結(jié)果表明汶川地震由映秀往北東方向發(fā)展,而蘆山地震破裂由深部向淺部偏西南滑動(dòng).如果忽略微小的滑動(dòng)量,兩個(gè)模型在中央斷裂上有約40 km的破裂空區(qū),在前山斷裂上有約60 km的破裂空區(qū).破裂空區(qū)的形成,有可能是深度10 km以上的高速異常區(qū)阻礙了兩個(gè)地震破裂的傳播,其上部也許是個(gè)障礙體.破裂空區(qū)的形成也有可能與龍門(mén)山15 km下低速異常區(qū)沒(méi)有積累足夠的地震能量有關(guān),其深部可能是龍門(mén)山斷裂上比較“軟”的部位,使得龍門(mén)山蘆山段比較容易向高原東緣擠出.這與Pei等(2014)發(fā)現(xiàn)彭灌雜巖體與寶興雜巖體間較“軟”的低速異常區(qū)相對(duì)應(yīng).高何—下里隱伏斷裂、大川—雙石斷裂、鹽井—五龍斷裂、耿達(dá)—隴東斷裂的走向都與龍門(mén)山中北段相對(duì)應(yīng)的山前隱伏斷裂、灌縣—江油斷裂、映秀—北川斷裂、汶川—茂汶斷裂有微小的偏移,與對(duì)應(yīng)斷裂向東錯(cuò)位15~20 km以上,應(yīng)該是長(zhǎng)期東向擠出的結(jié)果. 我們綜合利用GPS靜態(tài)同震位移和強(qiáng)震臺(tái)加速度積分推算的同震位移約束蘆山地震幾何與運(yùn)動(dòng)學(xué)模型.近場(chǎng)同震資料的大地測(cè)量反演結(jié)果顯示,蘆山地震破裂發(fā)生在上陡下緩的鏟式斷層上, 其向下接入近水平的深部滑脫層.相比高角度(大于45°)深大斷裂為主要特征的厚皮構(gòu)造,以斷坡-滑脫斷層為特征的薄皮構(gòu)造作為破裂面擬合地表位移的殘差平方和至少低15%.?dāng)嗥聰嗔焉喜?傾角43°~50°)和下部(傾角<27°)各有一個(gè)滑動(dòng)幅度大于0.5 m 的峰值破裂區(qū).2013年蘆山地震揭示的龍門(mén)山斷裂西南段深部破裂,與2008年汶川地震展示的龍門(mén)山中北段深部破裂特征十分吻合,表明剛性的四川盆地基底沿滑脫層被動(dòng)插入龍門(mén)山之下, 龍門(mén)山?jīng)_斷帶內(nèi)疊瓦狀鏟式結(jié)構(gòu)的斷坡-滑脫斷層調(diào)節(jié)青藏高原東緣的擠壓變形, 導(dǎo)致龍門(mén)山隆升. 致謝 感謝中國(guó)地震局預(yù)測(cè)研究所的江在森研究員提供部分GPS數(shù)據(jù),國(guó)家強(qiáng)震動(dòng)臺(tái)網(wǎng)中心提供強(qiáng)震動(dòng)數(shù)據(jù).審稿專家的建議對(duì)本文質(zhì)量的提高有很大的幫助,在此一并致謝. Chen C, Xu Y. 2013. Relocation of the LushanMS7.0 earthquake sequence and its tectonic implication.ChineseJ.Geophys. (in Chinese), 56(12): 4028-4036, doi: 10. 6038/cjg20131208. Chen G G, Ji F J, Zhou R J, et al. 2007. Primary research of activity segmentation of Longmenshan fault zone since Late-Quaternary.SeismologyandGeology(in Chinese), 29(3): 657-673. Chen L C, Ran Y K, Chen J, et al. 2009. Geological and geomorphic evidence for repeatedly occurring severe earthquakes on fronter-range fault of Longmenshan, China.QuaternarySciences(in Chinese), 29(3): 472-483. Deng Q D, Chen S F, Zhao X L. 1994. Tectonics, seismicity and dynamics of Longmenshan Mountains and its adjacent regions.SeismologyandGeology(in Chinese), 16(4): 389-403. Du Y J, Aydin A, Segall P. 1992. Comparison of various inversion techniques as applied to the determination of a geophysical deformation model for the 1983 Borah Peak earthquake.Bull.Seism.Soc.Am., 82(4): 1840-1866. Fang L H, Wu J P, Wang W L, et al. 2013. Relocation of the mainshock and aftershock sequences ofMS7.0 Sichuan Lushan earthquake.ChineseScienceBulletin, 58(28-29): 3451-3459, doi: 10.1007/s11434-013-6000-2. Iwan W D, Moser M A, Peng C Y. 1985. Some observations on strong motion earthquake measurement using a digital accelerograph.Bull.Seism.Soc.Am., 75(5): 1225-1246. Jiang W L, Xie X S, Zhang J F, et al. 2009. Evidence of multistage late Quaternary strong earthquakes on typical segments of Longmenshan Active Fault Zone in Sichuan, China.Sci.ChinaSer.D-EarthSci., 52(9): 1412-1415, doi: 10.1007/s11430-009-0145-5. Jiang Z S, Wang M, Wang Y Z, et al. 2014. GPS constrained coseismic source and slip distribution of the 2013MW6.6 Lushan, China, earthquake and its tectonic implications.GeophysicalResearchLetters, 41(2): 407-413. Johnson K M, Segall P. 2004. Imaging the ramp-décollement geometry of the Chelungpu fault using coseismic GPS displacements from the 1999 Chi-Chi, Taiwan earthquake.Tectonophysics, 378(1-2): 123-139. Liao H, Xu R, Chen W F, et al. 2013. Property variation and statistical analysis of Sichuan GPS time series before and after Wenchuan Earthquake.ChinesJ.Geophys. (in Chinese), 56(4): 1237-1245, doi: 10. 6038/cjg20130418. Li Q, You X Z, Yang S M, et al. 2012. A precise velocity field of tectonic deformation in China as inferred from intensive GPS observations.ScienceChinaEarthSciences, 55(5): 695-698, doi: 10.1007/s11430-012-4412-5. Li Z W, Xu Y, Huang R Q, et al. 2011. Crustal P-wave velocity structure of the Longmenshan region and its tectonic implications for the 2008 Wenchuan earthquake.ScienceChinaEarthSciences, 54(9): 1386-1393, doi: 10.1007/s11430-011-4177-2. Liu Y H, Wang C S, Shan X J, et al. 2014. Result of SAR differential interferometry for the co-seismic deformation and source parameter of theMS7.0 Lushan Earthquake.ChineseJ.Geophys. (in Chinese), 57(8): 2495-2506, doi: 10.6038/cjg20140811. Matthews M V, Segall P. 1993. Statistical inversion of crustal deformation data and estimation of the depth distribution of slip in the 1906 earthquake.J.Geophys.Res., 98: 12153-12163. Okada Y. 1985. Surface deformation due to shear and tensile faults in a half-space.Bull.Seism.Soc.Am., 75(4): 1135-1154. Pei S P, Zhang H J, Su J R, et al. 2014. Ductile gap between the Wenchuan and Lushan earthquakes revealed from the two-dimensional Pg seismic tomography.ScientificReports, 4: 6489 doi: 10.1038/srep06489 Ran Y K, Chen L C, Chen G H, et al. 2008. Primary analyses of in-situ recurrence of large earthquake along seismogenic fault of theMS8.0 Wenchuan earthquake.SeismologyandGeology(in Chinese), 30(3): 630-643. Ran Y K, Chen L C, Chen J, et al. 2010. Paleoseismic evidence and repeat time of large earthquakes at three sites along the Longmenshan Fault zone.Tectonophysics, 491(1-4): 141-53. doi: 10.1016/j. tecto.2010.01.009. Ran Y K, Chen W S, Xu X W, et al. 2013. Paleoseismic events and recurrence interval along the Beichuan-Yingxiu Fault of Longmenshan Fault zone, Yingxiu, Sichuan, China.Tectonophysics, 584: 81-90. Shi X. 2009. Structural deformation style and study of the paleoearthquake at the Beichuan-Dengjia segment along the central fault, Longmen Mountain [Master′s thesis] (in Chinese). Beijing: Institute of Geology, CEA. Tan K, Yang S M, Qiao X J, et al. 2013. Rupture of ramp-décollement faults in the 2008 Wenchuan earthquake-Geodetic evidence for the uplift of the Longmen Shan driven by convergent deformation.ChineseJ.Geophys. (in Chinese), 56(5): 1506-1516, doi: 10.6038/cjg20130509. Toda S, Stein R S, Richards-Dinger K, et al. 2005. Forecasting the evolution of seismicity in southern California: Animations built on earthquake stress transfer.JournalofGeophysicalResearch, 110: B05S16, doi: 10.1029/2004JB003415. Wang H. 2011. The structure characteristics of the Wenchuan earthquake fault zone and its relationship with seismic activity [Master′s thesis] (in Chinese). Beijing: China University of Geosciences. Wang Q, Qiao X J, Lan Q G, et al. 2011. Rupture of deep faults in the 2008 Wenchuan earthquake and uplift of the Longmen Shan.NatureGeoscience, 4(9): 634-640. Wang R J, Lorenzo-Martín F, Roth F. 2003. Computation of deformation induced by earthquakes in a multi-layered elastic crust-FORTRAN programs EDGRN/EDCMP.ComputerandGeosciences, 29(2): 195-207. Wang W M, Hao J L, Yao Z X. 2013. Preliminary result for rupture process of Apr. 20, 2013, Lushan earthquake, Sichuan, China.ChineseJ.Geophys. (in Chinese), 56(4): 1412-1417, doi: 10.6038/cjg20130436. Xu X W, Wen X Z, Ye J Q, et al. 2008. TheMS8.0 Wenchuan earthquake surface ruptures and its seismogenic structure.SeismologyandGeology(in Chinese), 30(3): 597-629. Xu X W, Wen X Z, Han Z J, et al. 2013. LushanMS7.0 earthquake: A blind reserve-fault event.ChineseScienceBulletin, 58(28-29): 3437-3443, doi: 10.1007/s11434-013-5999-4. Xu Z Q, Hou L W, Wang Z X, et al. 1992. Orogenic Processes of the Songpan Ganze Orogenic Belt of China (in Chinese). Beijing: Geological Publishing House. Yang S M, Lan Q G, Nie Z S, et al. 2012. Coseismic displacement caused by the 2008 great Wenchuan earthquake derived from various types of geodetic data.ChineseJ.Geophys. (in Chinese), 55(8): 2575-2588, doi: 10.6038/j.issn.0001-5733.2012.08.009.Zhao B, Gao Y, Huang Z B, et al. 2013. Double difference relocation, focal mechanism and stress inversion of LushanMS7.0 earthquake sequence.ChineseJ.Geophys. (in Chinese), 56(10): 3385-3395, doi: 10.6038/cjg20131014.Zhao C P, Zhou L Q, Chen Z L. 2013. Source rupture process of LushanMS7.0 earthquake, Sichuan, China and its tectonic implication.ChineseScienceBulletin, 58(28-29): 3444-3450, doi: 10.1007/s11434-013-6017-6. Zhang G W, Lei J S. 2013. Relocations of Lushan, Sichuan Strong earthquake (MS7.0) and its aftershocks.ChineseJ.Geophys. (in Chinese), 56(5): 1764-1771, doi: 10.6038/cjg20130534. Zhang L T, Jin S, Wei W B, et al. 2012. Electrical structure of crust and upper mantle beneath the eastern margin of the Tibetan plateau and the Sichuan basin.ChineseJ.Geophys. (in Chinese), 55(12): 4126-4137, doi: 10.6038/j.issn.0001-5733.2012.12.025. Zhang P Z, Xu X W, Wen X Z, et al. 2008. Slip rates and recurrence intervals of the Longmen Shan active fault zone, and tectonic implications for the mechanism of the May 12 Wenchuan earthquake, 2008, Sichuan, China.ChineseJ.Geophys. (in Chinese), 51(4): 1066-1073, doi: 10.3321/j.issn:0001-5733.2008.04.015. Zhang Y, Feng W P, Xu L S, et al. 2009. Spatio-temporal rupture process of the 2008 great Wenchuan earthquake.ScienceinChinaSeriesD:EarthSciences, 52(2): 145-152. Zhang Y Q, Li H L. 2010. Late Quaternary active faulting along the SW segment of the Longmenshan Fault zone.QuaternarySciences(in Chinese), 30(4): 699-710. Zhang Y, Xu L S, Chen Y T. 2013. Rupture process of the Lushan 4.20 earthquake and preliminary analysis on the disaster-causing mechanism.ChineseJ.Geophys. (in Chinese), 56(4): 1408-1411, doi: 10.6038/cjg20130435. Zheng Y, Ge C, Xie Z J, et al. 2013. Crustal and upper mantle structure and the deep seismogenic environment in the source regions of the Lushan earthquake and the Wenchuan earthquake.ScienceChina:EarthSciences, 56(7): 1158-1168, doi: 10.1007/s11430-013-4641-2. 附中文參考文獻(xiàn) 陳晨, 胥頤. 2013. 蘆山MS7.0 級(jí)地震余震序列重新定位及構(gòu)造意義. 地球物理學(xué)報(bào), 56(12): 4028-4036, doi: 10.6038/cjg20131208.陳國(guó)光, 計(jì)鳳桔, 周榮軍等. 2007. 龍門(mén)山斷裂帶晚第四紀(jì)活動(dòng)性分段的初步研究. 地震地質(zhì), 29(3): 657-673. 陳立春, 冉勇康, 陳杰等. 2009. 龍門(mén)山前山斷裂大地震原地重復(fù)發(fā)生的地質(zhì)地貌證據(jù). 第四紀(jì)研究, 29(3): 472-483. 鄧起東, 陳社發(fā), 趙小麟. 1994. 龍門(mén)山及其鄰區(qū)的構(gòu)造和地震活動(dòng)及動(dòng)力學(xué). 地震地質(zhì), 16(4): 389-403. 房立華, 吳建平, 王未來(lái)等. 2013. 四川蘆山MS7.0 級(jí)地震及其余震序列重定位. 科學(xué)通報(bào), 58(20): 1901-1909. 江娃利, 謝新生, 張景發(fā)等. 2009. 四川龍門(mén)山活動(dòng)斷裂帶典型地段晚第四紀(jì)強(qiáng)震 多期活動(dòng)證據(jù). 中國(guó)科學(xué)D輯: 地球科學(xué), 39(12): 1688-1700. 廖華, 徐銳, 陳維鋒等. 2013. 汶川地震前后四川區(qū)域GPS時(shí)序特征演變及統(tǒng)計(jì)分析. 地球物理學(xué)報(bào), 56(4): 1237-1245, doi: 10.6038/cjg20130418. 李強(qiáng), 游新兆, 楊少敏等. 2012. 中國(guó)大陸構(gòu)造變形高精度大密度GPS 監(jiān)測(cè)—現(xiàn)今速度場(chǎng). 中國(guó)科學(xué): 地球科學(xué), 42(5): 629-632.李志偉, 胥頤, 黃潤(rùn)秋等. 2011. 龍門(mén)山地區(qū)的P波速度結(jié)構(gòu)與汶川地震的深部構(gòu)造特征. 中國(guó)科學(xué): 地球科學(xué), 41(3): 283-290. 劉云華, 汪馳升, 單新建等. 2014. 蘆山MS7.0 級(jí)地震InSAR 形變觀測(cè)及震源參數(shù)反演. 地球物理學(xué)報(bào), 57(8): 2495-2506, doi: 10.6038/cjg20140811. 冉勇康, 陳立春, 陳桂華等. 2008. 汶川MS8.0地震發(fā)震斷裂大地震原地重復(fù)現(xiàn)象初析. 地震地質(zhì), 30(3): 630-643. 史翔. 2009. 龍門(mén)山中央斷裂帶北川—鄧家地區(qū)構(gòu)造變形樣式及古地震研究[碩士論文]. 北京: 中國(guó)地震局地質(zhì)研究所. 譚凱, 楊少敏, 喬學(xué)軍等. 2013. 2008年汶川地震中斷坡-滑脫斷層破裂: 龍門(mén)山擠壓隆升的大地測(cè)量證據(jù). 地球物理學(xué)報(bào), 56(5): 1506-1516, doi: 10.6038/cjg20130509. 王煥. 2011. 汶川地震斷裂帶結(jié)構(gòu)特征及其與地震活動(dòng)性的關(guān)系[碩士論文]. 北京: 中國(guó)地質(zhì)大學(xué). 王衛(wèi)民, 郝金來(lái), 姚振興. 2013. 2013年4月20日四川蘆山地震震源破裂過(guò)程反演初步結(jié)果. 地球物理學(xué)報(bào), 56(4): 1412-1417, doi: 10.6038/cjg20130436. 徐錫偉, 聞學(xué)澤, 葉建青等. 2008. 汶川MS8.0地震地表破裂帶及其發(fā)震構(gòu)造. 地震地質(zhì), 30(3): 597-629. 徐錫偉, 聞學(xué)澤, 韓竹軍等. 2013. 四川蘆山7.0級(jí)強(qiáng)震: 一次典型的盲逆斷層型地震. 科學(xué)通報(bào), 58(20): 1887-1893. 許志琴, 侯立煒, 王宗秀等. 1992. 中國(guó)松潘一甘孜造山帶的造山過(guò)程. 北京: 地質(zhì)出版社. 楊少敏, 蘭啟貴, 聶兆生等. 2012. 用多種數(shù)據(jù)構(gòu)建2008年汶川特大地震同震位移場(chǎng). 地球物理學(xué)報(bào), 55(8): 2575-2588, doi: 10.6038/j.issn.0001-5733.2012.08.009. 趙博, 高原, 黃志斌等. 2013. 四川蘆山MS7.0地震余震序列雙差定位、震源機(jī)制及應(yīng)力場(chǎng)反演. 地球物理學(xué)報(bào), 56(10): 3385-3395, doi: 10.6038/cjg20131014. 趙翠萍, 周連慶, 陳章立. 2013. 2013年四川蘆山MS7.0 級(jí)地震震源破裂過(guò)程及其構(gòu)造意義. 科學(xué)通報(bào), 58(20): 1894-1900. 張廣偉, 雷建設(shè). 2013. 四川蘆山7.0級(jí)強(qiáng)震及其余震序列重定位. 地球物理學(xué)報(bào), 56(5): 1764-1771, doi: 10.6038/cjg20130534. 張樂(lè)天, 金勝, 魏文博等. 2012. 青藏高原東緣及四川盆地的殼幔導(dǎo)電性結(jié)構(gòu)研究. 地球物理學(xué)報(bào), 55(12): 4126-4137, doi: 10.6038/j.issn.0001-5733.2012.12.025. 張培震, 徐錫偉, 聞學(xué)澤等. 2008. 2008年汶川8.0級(jí)地震發(fā)震斷裂的滑動(dòng)速率、復(fù)發(fā)周期和構(gòu)造成因. 地球物理學(xué)報(bào), 51(4): 1066-1073, doi: 10.3321/j.issn:0001-5733.2008.04.015. 張勇, 馮萬(wàn)鵬, 許力生等. 2008. 2008年汶川大地震的時(shí)空破裂過(guò)程. 中國(guó)科學(xué)D輯: 地球科學(xué), 38(10): 1186-1194. 張?jiān)罉? 李海龍. 2010. 龍門(mén)山斷裂帶西南段晚第四紀(jì)活動(dòng)性調(diào)查與分析. 第四紀(jì)研究, 30(4): 699-710. 張勇, 許力生, 陳運(yùn)泰. 2013. 蘆山4.20地震破裂過(guò)程及其致災(zāi)特征初步分析. 地球物理學(xué)報(bào), 56(4): 1408-1411, doi: 10.6038/cjg20130435. 鄭勇, 葛粲, 謝祖軍等. 2013. 蘆山與汶川地震震區(qū)地殼上地幔結(jié)構(gòu)及深部孕震環(huán)境. 中國(guó)科學(xué): 地球科學(xué), 43(6): 1027-1037. (本文編輯 胡素芳) Rupture models of the 2013 Lushan earthquake constrained by near field displacements and its tectonic implications TAN Kai1, WANG Qi2*, DING Kai-Hua3, LI Heng1, ZOU Rong2, NIE Zhao-Sheng1,WANG Di-Jin1, YANG Shao-Min1, QIAO Xue-Jun1 1KeyLaboratoryofEarthquakeGeodesy,InstituteofSeismology,CEA,Wuhan430071,China2InstituteofGeophysics&Geomatics,ChinaUniversityofGeosciences,Wuhan430074,China3FacultyofInformationEngineering,ChinaUniversityofGeosciences,Wuhan430074,China Following the 2008 WenchuanMS8.0 earthquake in the Longmen Shan fault zone, the 2013 LushanMS7.0 earthquake struck its southwestern section. Analysis of seismic waves and coseismic surface displacements together with field investigations all show that this earthquake ruptured a compact areas at depths of 10~15 km, characterized by a predominating thrust motion that may not reach to the surface. Yet the results differ considerably in details about the geometry and distribution of coseismic slip, partly owing to the scarcity of near-field observation, it remains elusive whether the coseismic slip has extended downward in a fashion similar to what the Wenchuan earthquake did along a deep-seated décollement under the Longmen Shan. We used an improved data set of near-field coseismic deformation acquired by GPS receivers and accelerate-meter to constraint the rupture geometry and slip distribution. We show the existent of coseismic slip on a basal detachment fault, which suggests that the Longmen Shan is thickened at this part by crustal shortening between the Tibetan Plateau and the Sichuan Basin.We combined near field coseismic offsets acquired by continuously-operating and campaign GPS stations as well as seismic strong motion stations. We processed GPS data with GIPSY software, and integrate twice acceleration waveforms with zero correction method to obtain static offsets caused by the earthquake. We use these data to constrain the rupture geometry and slip distribution based on the half-space elastic dislocation model. Firstly, we used the simulated annealing algorithm to determine model parameters and corresponding confidence intervals assuming that coseismic slip is distributed on multiple rectangle planes. Then, we divide the optimal model planes further into numerous sub-fault patches, and use the nonnegative least squares to estimate slip values on these patches to outline collectively feature of the distribution of coseismic slip on this optimal rupture plane. In both inversions, we minimize misfits to surface displacements while maintaining smoothness of slip across the neighboring patches by the cross validation sum of squares.The updated data set of static offsets includes additional 30 stations, and are double in number the early one. All stations are located largely in a region in dimension of 80~100 km between Qionglai and Tianquan. The maximum offset of 67 mm is found at one GPS station 16 km away from the epicenter. The surface displacements are all directed to the epicenter, in consistent with thrusting mechanism. The uncertainty for the continuously-operating GPS stations is less than 3 mm, and less than 5 mm for the campaign stations whereas the uncertainty of horizontal components of a strong-motion station is estimated to be about 10 mm by comparing with the GPS data. Our modeling indicates that the strike of the rupture is about 212° parallel to the topographic front of the Longmen Shan. Our preferred model plane is approximated by a listric fault with a steeply-dipping ramp shallower than 16 km, which soles into agently-dipping detachment at 16~25 km. The dip angle is greater than 40° for the ramp and is reduced to less than 30° for the detachment fault. Our slip model consists of two main patches with local pick slip exceeding 0.5 m. One are located at the ramp with dip angles of 43°~50°. The maximum slip of 1.5 m is found at 13 km depth. Another is on the detachment faults that dips 27°. The slip on the detachment fault is extended downwards to the depths of 25~26 km where aftershocks ceased. The relocated aftershocks surround largely the main slip patches. The majority of the aftershocks were restricted to the areas that experienced an enhanced Coulomb stress due to the main shock. The geodetic moment is estimated to 1.14×1019N·m corresponding to the moment magnitudeMw6.6 in accordance with the GCMT solution. The seismic imaging of deep structures beneath the Longmen Shan reveals the P-wave velocities change abruptly at the depths of 13~17 km, coincident with the location of the décollement inferred from our model. The magnetotelluric observation shows that the Longmen Shan fault corresponds to as a gently-dipping conductor that extend downward to the depth of about 20 km, coincident with the location of the décollement.The deep slip caused by the 2013 Lushan earthquake is found along the southern section of the Longmen Shan, and is much consistent with that of central and northern Longmen Shan caused by the 2008 Wenchuan earthquake. The observation provides evidence for the existence of a ramp-décollement preferred by a thin-skined tectonics for the Longmen Shan. The opposite viewpoint is a thick-skinned tectonics that requires a steeper basal fault to uplift the Longmen Shan. Our modelling precludes such a deep fault extending into the mid- and lower crust. The ramp-décollement under the Longmen Shan highlights the crustal shortening across the eastern margin of the Tibetan Plateau to uplift the Longmen Shan. The rigid Sichuan Basin basement underthrusts under the Longmen Shan along the detachment layer. The ramp-décollement adjusts the compressional deformation on the eastern margin of the Tibetan Plateau, resulting in the uplift of the Longmen Shan. Lushan earthquake; Coseismic deformation; Fault geometry; Rupture slip distribution 譚凱, 王琪, 丁開(kāi)華等. 2015. 近場(chǎng)位移數(shù)據(jù)約束的2013年蘆山地震破裂模型及其構(gòu)造意義.地球物理學(xué)報(bào),58(9):3169-3182, 10.6038/cjg20150913. Tan K, Wang Q, Ding K H, et al. 2015. Rupture models of the 2013 Lushan earthquake constrained by near field displacements and its tectonic implications.ChineseJ.Geophys. (in Chinese),58(9):3169-3182,doi:10.6038/cjg20150913. 10.6038/cjg20150913 P228, P313 2015-01-09,2015-06-17收修定稿 中國(guó)地震局地震研究所所長(zhǎng)基金(IS201506220),國(guó)家自然科學(xué)基金(40974012,41304019),地震行業(yè)科研專項(xiàng)(201208006)資助. 譚凱,博士,研究員,主要從事大地測(cè)量與地球動(dòng)力學(xué)方面的研究.E-mail: whgpstan@163.com *通訊作者 王琪,博士,教授,主要從事大地測(cè)量與地球動(dòng)力學(xué)方面的研究. E-mail: wangqi@cug.edu.cn4 討論
5 結(jié)論