2010年智利馬烏萊地震滑移與安第斯俯沖帶震前閉鎖的相關(guān)性＊

Marcos Moreno，Matthias Rosenau，Onno Oncken

（Helmholtz Centre Potsdam，GFZ German Research Centre for Geosciences，Telegrafenberg，Potsdam 14473，Germany）

2010年智利馬烏萊地震滑移與安第斯俯沖帶震前閉鎖的相關(guān)性＊

Marcos Moreno，Matthias Rosenau，Onno Oncken

（Helmholtz Centre Potsdam，GFZ German Research Centre for Geosciences，Telegrafenberg，Potsdam 14473，Germany）

2010年2月27日馬烏萊（智利）M 8.8地震使安第斯大逆沖型俯沖帶的一部分區(qū)段破裂，該區(qū)段被認(rèn)為具有極高的潛在地震危險(xiǎn)性[1－6]。此次地震是使俯沖帶內(nèi)一個(gè)長期地震空區(qū)發(fā)生破裂的一次最大地震，在此之前，一個(gè)密集型空間大地測(cè)量臺(tái)網(wǎng)對(duì)該俯沖帶進(jìn)行了監(jiān)測(cè)。這就為評(píng)估震間閉鎖與同震滑移的空間相關(guān)性提供了板塊界面震前閉鎖狀態(tài)的前所未有的高分辨率圖像。鑒于閉鎖類似于滑動(dòng)的基本假設(shè)，震前閉鎖可被用于預(yù)測(cè)許多地震空區(qū)的未來破裂[6－12]。然而，沒有第一次“填補(bǔ)空白”的地震的發(fā)生，這一假說就無法得以驗(yàn)證。在此，我們提供的證據(jù)表明，2010年馬烏萊地震滑移分布狀況與根據(jù)之前10年間全球定位系統(tǒng)（GPS）觀測(cè)結(jié)果反演得到的震間閉鎖分布密切相關(guān)。此次地震成核于一個(gè)閉鎖程度很高的區(qū)域內(nèi)，并將1835年大地震以來該區(qū)聚集的應(yīng)力幾乎釋放一空。地震滑動(dòng)量很大的兩個(gè)區(qū)域（凸起體）在此次地震之前幾乎處于完全閉鎖狀態(tài)。在這兩個(gè)凸起體之間，此次破裂形成了一個(gè)處于震間蠕動(dòng)狀態(tài)的橋梁地帶，其同震滑動(dòng)量一貫較小。破裂中止于震前處于高度閉鎖狀態(tài)的區(qū)域，但其預(yù)應(yīng)力卻因20世紀(jì)重復(fù)發(fā)生的地震事件而大大降低。我們的結(jié)果顯示，單個(gè)凸起體尺度內(nèi)同震滑移的不均勻性應(yīng)該能夠說明未來大地震的潛在危險(xiǎn)，由此，可以根據(jù)大地測(cè)量結(jié)果對(duì)其進(jìn)行預(yù)測(cè)。

雖然我們無法預(yù)測(cè)下一次地震將在何時(shí)何地發(fā)生，但我們可以通過對(duì)今后數(shù)十年可能發(fā)生破裂的地帶進(jìn)行監(jiān)測(cè)，估算出未來地震事件的最大可信規(guī)模，關(guān)于這一點(diǎn)，根據(jù)板塊構(gòu)造理論以及從前發(fā)生的地震事件已經(jīng)得出過推論。這種易于發(fā)生破裂的地區(qū)被稱為“地震空區(qū)”[2－3，13－17]，因?yàn)樗鼈冊(cè)诘卣鸹顒?dòng)的時(shí)空分布上呈空白狀態(tài)。隨著空間大地測(cè)量技術(shù)的開展，現(xiàn)在，我們可以根據(jù)地震空區(qū)應(yīng)變（應(yīng)力）積累的近實(shí)時(shí)測(cè)量結(jié)果獲取有關(guān)地震危險(xiǎn)性的更多的約束條件[6－12]。然而，目前尚不清楚震間應(yīng)力積累究竟在多大程度上可以定量模擬未來地震的同震能量釋放。顯然，同震能量釋放在破裂區(qū)內(nèi)并不是均勻分布的，相反，破裂區(qū)是一個(gè)拼合而成的區(qū)域，其中某些地塊（凸起體）滑動(dòng)量大，某些地塊（障礙體）滑動(dòng)量小，甚至沒有滑動(dòng)[18]。的確，近期沿安第斯[10]和蘇門答臘[11－12]俯沖帶等地區(qū)的大地測(cè)量觀測(cè)結(jié)果表明，沿大逆沖型板塊界面可識(shí)別出一個(gè)由閉合地帶（潛在凸起體）拼合而成的類似區(qū)域，這些閉合地帶被蠕變區(qū)（潛在障礙體）所包圍。這一拼合區(qū)域在震間以滑移缺失的形式積累應(yīng)力，其空間模式是不均勻的，類似于未來的地震滑移分布，但這一假說還有待證實(shí)。

2010年2月27日襲擊智利中南部馬烏萊和比奧比奧地區(qū)的MW8.8地震（時(shí)間：06：34：14 UTC；震中：35.909°S，72.733°W；深度：35 km；http：∥earthquake.usgs.gov／earthquakes／eqinthenews／2010／us2010tfan／；參考文獻(xiàn)[1]）是現(xiàn)代空間大地測(cè)量監(jiān)測(cè)網(wǎng)捕捉到的第一次大地震，而且有可能是填補(bǔ)地震空區(qū)的一次大地震。由此，我們能夠?qū)φ痖g應(yīng)力積累模式與潛在災(zāi)難性事件期間同震釋放之間的相似性進(jìn)行空前詳盡的檢測(cè)和評(píng)估。沿安第斯板塊邊緣，納斯卡（Nazca）海洋板塊向南美大陸板塊下方俯沖（圖1a），此處大型逆沖區(qū)地震發(fā)生的頻率為每100～200年一次，地點(diǎn)可位于該邊緣的任一特定區(qū)段[2－4，19－20]。2010年馬烏萊地震使一個(gè)長期地震空區(qū)——康塞普西翁（Concepción）空區(qū)——破裂，根據(jù)歷史事件推測(cè)，這一地區(qū)被疑為地震的高危險(xiǎn)區(qū)[2－5]。其高地震危險(xiǎn)性最近得到了空間大地測(cè)量結(jié)果的支持[6]，該結(jié)果顯示出震前數(shù)十年內(nèi)積累了大量應(yīng)力的一個(gè)閉鎖板塊界面。

圖1 研究區(qū)的構(gòu)造背景、數(shù)據(jù)、觀測(cè)及結(jié)果。（a）智利中南部安第斯俯沖帶的暈渲地形圖。沿邊緣的地震分段用橢圓表示，這些橢圓顯示出歷史地震的大概破裂區(qū)范圍（更新自文獻(xiàn)[4－6]）。小插圖（矩形）示出了圖a相對(duì)于南美大陸的位置。（b）2010年馬烏萊地震前GPS觀測(cè)的穩(wěn)定的南美地表速度匯總（1996—2008）（參考數(shù)據(jù)見在線補(bǔ)充資料“方法詳述”）。帶箭頭的橢圓代表95%的置信界限。（c）根據(jù)GPS＋FEM模擬獲得的2010年馬烏萊地震以前的10年中沿安第斯大逆沖型俯沖帶的界面閉鎖（板塊會(huì)聚部分）分布。2010年馬烏萊地震的震中（白色五角星，USGS／NEIC）和震源機(jī)制標(biāo)示于圖a和c中

然而，就地震能量釋放分布而言，現(xiàn)有的模型未能揭示出預(yù)測(cè)未來地震危險(xiǎn)性所需的單一凸起體和障礙體閉鎖分布的詳細(xì)狀況?？等瘴魑炭諈^(qū)的上次大地震是由達(dá)爾文報(bào)道的[21]，發(fā)生在175年前的1835年，估算震級(jí)為M≈8.5（見文獻(xiàn)[19]，圖1a）。自那以來，該空區(qū)以北（1906年M 8.4，文獻(xiàn)[4]；1985年M7.8，文獻(xiàn)[22]）和以南（1960年M 9.5，文獻(xiàn)[23－24]）地區(qū)共發(fā)生了3次“大型逆沖區(qū)”地震。康塞普西翁空區(qū)內(nèi)的北部曾于1928年發(fā)生過一次“大型逆沖區(qū)”地震（M≈8；參考文獻(xiàn)[4]），據(jù)推測(cè)只釋放了自1835年以來積累的小部分應(yīng)力。

在此，我們利用2010年馬烏萊地震以前十幾年間（1996—2008）的表面速度GPS觀測(cè)結(jié)果，得到了震前閉鎖的安第斯大逆沖型俯沖帶的詳細(xì)圖像。在利用匯總的研究區(qū)GPS觀測(cè)資料（圖1b、補(bǔ)充圖1、方法詳述和補(bǔ)充信息）進(jìn)行反演時(shí)，我們使用了安第斯俯沖帶的一個(gè)球形的、層狀有限元模型（FEM），該模型與根據(jù)近期地球物理大斷面匯編[25]得到的模型一樣，也包含了地形學(xué)、測(cè)深學(xué)和真實(shí)的板塊形態(tài)等信息（見方法詳述和補(bǔ)充圖2）。通過除去由其他一些源——如1960年M9.5地震之后區(qū)域粘彈性震后地幔弛豫、地殼斷裂運(yùn)動(dòng)和弧前裂塊運(yùn)動(dòng)（圖1b和補(bǔ)充圖1b）——產(chǎn)生的地表形變信號(hào)，我們對(duì)GPS速度進(jìn)行了校正，然后，對(duì)殘余模式進(jìn)行解釋（補(bǔ)充圖1c），以單獨(dú)反映存在非均勻閉鎖板塊界面的情況下由納斯卡－南美板塊會(huì)聚產(chǎn)生的震間彈性應(yīng)變積累。在此，沿大逆沖型俯沖帶的震前閉鎖在運(yùn)動(dòng)學(xué)上被定義為地震事件前作為滑移缺失而積累的板塊會(huì)聚的部分。

將GPS獲得的震前閉鎖分布與2010年馬烏萊地震的兩個(gè)初始同震滑移模型進(jìn)行比較，這兩個(gè)模型一個(gè)源自遠(yuǎn)震體波反演（文獻(xiàn)[26－27]；圖2a、2b），另一個(gè)作為滑移模型約束，是利用圣地亞哥和康塞普西翁附近IGS－GPS臺(tái)站的同震GPS速度進(jìn)行遠(yuǎn)震－GPS聯(lián)合反演獲得的（文獻(xiàn)[28]；圖2c）。據(jù)此，我們檢驗(yàn)了震間閉鎖與同震滑移之間的相似性假說。最后，我們利用GEOFON擴(kuò)展虛擬網(wǎng)絡(luò)（GEOFON Extended Virtual Network——GEVN）提供的余震分布[29]對(duì)震后行為進(jìn)行詮釋（圖3a、3b）。

雖然這里分析的滑移分布只是初步結(jié)果，而且其復(fù)雜程度也不盡相同，但它們顯示出震中以南和震中以北兩個(gè)滑動(dòng)量很大的地塊（凸起體）的一級(jí)模式，兩個(gè)地塊被一個(gè)50～100 km寬的滑動(dòng)量很小的地帶分開（圖2a—c）。敏感性分析[1]表明這一主體模式是此次地震遠(yuǎn)震反演的一個(gè)穩(wěn)定特征。因此，滑移成核于從前高度閉鎖（板塊會(huì)聚＞0.75的部分）區(qū)域外圍破裂區(qū)段的中心部位。此處的震前閉鎖程度朝板塊會(huì)聚方向呈迅速下傾式降低（圖1c），很可能強(qiáng)加了一個(gè)局部剪切應(yīng)力梯度，使破裂成核。雙邊破裂[1]向南傳播，南邊的凸起體破裂，由此連通了從前的北邊的一個(gè)蠕滑帶，而后使北邊的又一個(gè)凸起體破裂。該破裂雖然連通了數(shù)十公里寬的一個(gè)蠕變區(qū)段（即潛在的松弛障礙體[10－12，18]），但它卻在馬烏萊地震前（圖2a—c）反直覺地中止于震前閉鎖程度很高的區(qū)域（即潛在凸起體），這似乎自相矛盾，關(guān)于這一點(diǎn)我們以后再討論。

在主要破裂區(qū)（～34°S—37°S），從視覺上觀察震前閉鎖與同震滑移分布之間的相似性（圖2a—c），可以發(fā)現(xiàn)高度的一致性?；瑒?dòng)量＞5～10 m的兩個(gè)凸起體與高度閉鎖（＞0.75）互相關(guān)聯(lián)，其間滑動(dòng)量小的地帶（＜5～10 m）與相對(duì)較低程度（＜0.75）的震前閉鎖區(qū)域吻合。我們得出了震后滑移缺失殘余圖（圖2d—f），力求對(duì)震間與同震滑移模式的相似性進(jìn)行定量分析。為此，我們對(duì)每一子斷層段中的同震滑移（180個(gè)來自USGS[26]模型的子斷層段，220個(gè)來自UCSB[27]模型的子斷層段，276個(gè)來自Caltech／JPL[28]模型的子斷層段）和1835年最后一次大地震以來理論上積累的滑移缺失進(jìn)行抗衡分析，并假定此次地震事件以前的10年間觀測(cè)到的閉鎖模式代表著先前的整個(gè)震間時(shí)期。結(jié)果（圖2d—f）表明，2010年馬烏萊地震的主要破裂區(qū)的殘差很低（USGS[26]和UCSB[27]滑移模型中＜2 m，Caltech／JPL[28]滑移模型中～5 m）。這一定量關(guān)系說明，在平行于海溝的＞200～300 km距離范圍內(nèi)（圖3c），震前滑移缺失幾乎完全被釋放，同時(shí)，馬烏萊地震可能意味著康塞普西翁地震空區(qū)已完全破裂。兩個(gè)閉鎖地塊的破裂導(dǎo)致了完全破裂，而且還有一個(gè)蠕變區(qū)被連通，這一觀點(diǎn)也得到了較大（M≥5）余震模式的支持（圖3a、3b）。較早的（最初48小時(shí)）和晚些時(shí)候的（最初3個(gè)月）余震呈雙峰分布模式，其特點(diǎn)是在從前的閉鎖區(qū)密度非常高，它們?cè)诘卣鹌陂g起到了凸起體的作用。相比之下，從前的蠕變區(qū)余震活動(dòng)性不強(qiáng)，說明那里震后滑移較小并（或）呈現(xiàn)出無震狀態(tài)。

圖2 同震滑移與震間閉鎖的相似性。（a—c）疊加于2010年馬烏萊地震初步同震滑移分布上的研究區(qū)震前閉鎖分布（等值線）。（a）USGS遠(yuǎn)震模型（文獻(xiàn)[26]）；（b）UCSB遠(yuǎn)震模型（文獻(xiàn)[27]）；（c）Caltech／JPL遠(yuǎn)震＋GPS模型（文獻(xiàn)[28]）；（d—f）震后滑移缺失圖，顯示出2010年馬烏萊地震后存在于主要破裂區(qū)的殘差很低。通過計(jì)算震后滑移缺失，得出了2010年同震滑移與滑移缺失（即自1835年最后一次大地震以來同震滑移分布的每一子斷裂段假定積累的滑移缺失）之間的殘差。圖中的白色五角星和海濱球分別代表震中（據(jù)USGS／NEIC資料）和震源機(jī)制（據(jù)GCMT資料），紅色三角形代表火山。d—f中的殘差分別源自a—c中的滑移模型

圖3 震前、同震和震后形變圖像的關(guān)系。（a）疊加于震前閉鎖區(qū)域圖像（紅色陰影部分≥0.75）之上的2010年（藍(lán)色等值線；USGS滑移模型[26]）和1960年（綠色等值線；據(jù)文獻(xiàn)[30]）地震期間的同震滑移分布，以及早期的（震后最初48小時(shí)）M≥5余震地點(diǎn)（震級(jí)用大小不等的灰色圓圈表示；GEOFON數(shù)據(jù)[29]）。（b）沿北—南剖面早期的（最初48小時(shí)；共80個(gè)地震事件）和晚期的（最初3個(gè)月；共168個(gè)地震事件）余震密度柱狀圖（GEOFON數(shù)據(jù)[29]，M≥5）。（c）2010年馬烏萊地震后沿北—南剖面觀測(cè)的1835年以來的剩余滑移缺失（左欄，據(jù)USGS滑移模型[26]）。中欄和右欄示出重疊的20世紀(jì)地震對(duì)滑移缺失的影響（黑線是由多項(xiàng)式擬合數(shù)據(jù)得到的）。彩色數(shù)據(jù)點(diǎn)及日期示出按年表示的地震

然而，破裂為什么會(huì)側(cè)向中止于閉鎖程度很高而且沒有明顯的幾何障礙（如北邊的Juan Fernández洋脊）可以阻止其傳播的區(qū)域（圖1a、3a）？兩種似乎合理的障礙機(jī)制是[18]：破裂中止是由于——至少在此次事件期間——局部強(qiáng)界面的存在（強(qiáng)度障礙體），也可能由于從前的地震降低了2010年破裂外圍的預(yù)應(yīng)力（低應(yīng)力障礙體）。第一種假說無從驗(yàn)證，但對(duì)低應(yīng)力障礙體假說可通過下述方式做進(jìn)一步研究，即在計(jì)算1835年以來積累的震后滑移缺失時(shí)將20世紀(jì)中重疊的大型逆沖區(qū)地震記錄考慮在內(nèi)。

1960年M9.5地震使南部鄰近的大逆沖型安第斯區(qū)段破裂，在此期間，阿勞科半島（Arauco Peninsula）下方的滑移逐漸停止[30]（37°S—38°S；圖3a）。相應(yīng)地，它也抵消了2010年破裂區(qū)南部邊緣多達(dá)10 m的視滑移缺失，從而實(shí)際上降低了此高度閉鎖區(qū)積累的應(yīng)力（圖3c，中欄）。因此，2010年地震似乎向南傳播到了一個(gè)預(yù)應(yīng)力很低的環(huán)境，這種環(huán)境足以有效地阻止破裂。北部可能也存在一個(gè)類似的低應(yīng)力障礙體。那里，1985年M＝7.8地震期間的平均滑移約為1.6 m（文獻(xiàn)[22]以及其中的參考文獻(xiàn)），由震源尺度規(guī)律關(guān)系估算[31]的1906年M＝8.4和1928年M≈8地震事件的滑移分別約為3 m和1.5 m，這些滑移使2010年馬烏萊地震北端的視滑移缺失有效地降低至0±5 m（圖3c，右欄）。現(xiàn)在，沿安第斯俯沖帶整個(gè)研究區(qū)段的震后滑移缺失似乎很低，在37°S緯度附近的最大值約為5 m（圖3c）。后者可能歸因于遠(yuǎn)震反演未發(fā)現(xiàn)的同震滑動(dòng)分量（或由于網(wǎng)絡(luò)配置[1]，或由于滑動(dòng)速度低），也可能引發(fā)沿阿勞科半島下方大逆沖型斷層的震后余滑。

我們得出的結(jié)論是，2010年馬烏萊地震定量性地彌合了康塞普西翁地震空白，這與震前觀測(cè)是一致的——因此，通過震前觀測(cè)，它在某種程度上是可以預(yù)測(cè)的。然而，地震往往很復(fù)雜，它們（如2010年馬烏萊地震）連通的蠕變區(qū)在其他時(shí)候[32]或其他地點(diǎn)[10－12]可能起到障礙體作用，地震也可能以小震序列的形式僅釋放出閉鎖地塊聚集的部分滑移缺失[9，11－12]。將來，通過對(duì)震間閉鎖和大震期間的同震滑移進(jìn)行進(jìn)一步觀測(cè)，我們可以更加詳細(xì)地監(jiān)測(cè)沿大逆沖型斷層的應(yīng)力聚集與釋放，從而對(duì)俯沖帶的地震危險(xiǎn)性做出更加可靠的評(píng)估。

方法概述

本文使用的GPS數(shù)據(jù)集包含1996—2008年間觀測(cè)的232個(gè)已發(fā)表的速度數(shù)據(jù)（見補(bǔ)充圖1a和方法詳述）。利用文獻(xiàn)[30]中詳細(xì)描述的FEM技術(shù)對(duì)GPS速度進(jìn)行反演。該FEM是三維的、球形的、層狀的（彈性巖石圈、粘彈性軟流圈），它包括地形學(xué)、測(cè)深學(xué)以及真實(shí)的板片資料，也包括莫霍界面（Moho）幾何學(xué)資料（見文獻(xiàn)[25]；補(bǔ)充圖2）?，F(xiàn)時(shí)觀測(cè)到的智利中南部地表形變場(chǎng)模擬結(jié)果源自不同的地震構(gòu)造過程[30，33]（圖1b，補(bǔ)充圖1）：（1）1960年瓦爾迪維亞（Valdivia）地震后，本研究區(qū)以南持續(xù)很久的震后地幔弛豫，（2）阿勞科半島地區(qū)和沿Liqui?e－Ofqui斷裂帶的局部地殼形變，（3）由于震間閉鎖的大逆沖型俯沖帶的存在而產(chǎn)生的震間應(yīng)變積累。我們?cè)诖怂榻B的模型框架中，有閉鎖界面存在情況下的會(huì)聚（每年66 mm，取向N077°E；文獻(xiàn)[34]）是通過一個(gè)后滑模型進(jìn)行運(yùn)動(dòng)學(xué)描述的?！胺椒ㄔ斒觥焙脱a(bǔ)充圖1－6中對(duì)于GPS反演方法有更詳盡的描述，包括分辨率測(cè)試和敏感性分析等。研究區(qū)內(nèi)（32°～40°S，71°～75°W，0～100 km深）2010年馬烏萊地震的余震資料是通過GEVN服務(wù)器下載的（文獻(xiàn)[29]）。經(jīng)統(tǒng)計(jì)檢驗(yàn)后，我們?yōu)樵摂?shù)據(jù)集指定了一個(gè)M≈5的完整震級(jí)，并僅利用這一閾值以上的數(shù)據(jù)進(jìn)行進(jìn)一步分析。

（注：①原圖均為彩圖；②方法詳述、相關(guān)參考文獻(xiàn)及補(bǔ)充圖可通過此文的在線版獲取，網(wǎng)址：http：∥www.nature.com／nature／journal／v467／n7312／full／nature09349.html。）

譯自：Nature，9 September 2010，Vol.467，198－202

原題：2010 Maule earthquake slip correlates with pre－seismic locking of Andean subduction zone

（四川機(jī)電職業(yè)技術(shù)學(xué)院 徐寶學(xué) 譯；鄭需要 校）

（譯者電子信箱，徐寶學(xué)：xbx0505＠163.com）

[1] Lay T，Ammon C J，Kanamori H，et al.Teleseismic inversion for rupture process of the 27 February 2010 Chile（MW8.8）earthquake.Geophys.Res.Lett.，2010，37，L13301.doi：10.1029／2010GL043379

[2] Nishenko R.Seismic potential for large and great intraplate earthquakes along the Chilean and Southern Peruvian margins of South America：a quantitative reappraisal.J.Geophys.Res.，1985，90：3 589－3 615

[3] Nishenko S P.Circum－Pacific seismic potential：1989－1999.Pure Appl.Geophys.，1991，135（2）：169－259

[4] Beck S，Barrientos S，Kausel E，et al.Source characteristics of historic earthquakes along the central Chile subduction zone.J.South Am.Earth Sci.，1998，11（2）：115－129

[5]Campos J，Hatzfeld D，Madariaga R，et al.A seismological study of the 1835 seismic gap in southcentral Chile.Phys.Earth Planet.Inter.，2002，132（1－3）：177－195

[6] Ruegg J C，Rudloff A，Vigny C，et al.Interseismic strain accumulation measured by GPS in the seismic gap between Constitución and Concepción in Chile.Phys.Earth Planet.Inter.，2009，175（1－2）：78－85

[7] Savage J C，Lisowski M，Prescott W H.Strain accumulation in the Shumagin and Yakataga seismic gaps，Alaska.Science，1986，231（4 738）：585－587

[8] Bürgmann R，Kogan M G，Steblov G M，et al.Interseismic coupling and asperity distribution along the Kamchatka subduction zone.J.Geophys.Res.，2005，110，B07405.doi：10.1029／2005JB003648

[9] Murray J，Langbein J.Slip on the San Andreas Fault at Parkfield，California，over two earthquake cycles，and the implications for seismic hazard.Bull.Seism.Soc.Am.，2006，96：282－303

[10] Perfettini H，Avouac J－P，Tavera H，et al.Seismic and aseismic slip on the Central Peru megathrust.Nature，2010，465：78－81

[11] Chlieh M，Avouac J P，Sieh K，et al.Heterogeneous coupling of the Sumatran megathrust constrained by geodetic and paleogeodetic measurements.J.Geophys.Res.，2008，113，B05305.doi：10.1029／2007JB004981

[12] Konca O，Avouac J－P，Sladen A，et al.Partial rupture of a locked patch of the Sumatra megathrust during the 2007 earthquake sequence.Nature，2008，456：631－635

[13] Fedotov S A.Regularities of the distribution of strong earthquakes in Kamchatka，the Kurile Islands，and northeastern Japan.Trudy Inst.Phys.Earth.Acad.Sci.USSR，1965，36：66－93

[14] Mogi K.Some features of recent seismic activity in and near Japan（2）.Activity before and after great earthquakes.Bull.Earthq.Res.Inst.，Univ.Tokyo，1969，47：395－417

[15] Sykes L R.Aftershock zones of great earthquakes，seismicity gaps，and earthquake prediction for Alaska and the Aleutians.J.Geophys.Res.，1971，76：8 021－8 041

[16] Kelleher J，Sykes L，Oliver J.Possible criteria for predicting earthquake locations and their application to major plate boundaries of the Pacific and the Caribbean.J.Geophys.Res.，1973，78：2 547－2 585

[17] Mc Nally K C.Seismic gaps in space and time.Ann.Rev.Earth Planet.Sci.，1983，11：359－369

[18] Kanamori H.Mechanics of Earthquakes.Ann.Rev.Earth Planet.Sci.，1994，22：207－237

[19] Lomnitz C.Grandes terremotos y tsunamis en Chile durante el periodo 1535－1955.Geofis.Panam.，1971，1：151－178

[20] Kelleher J.Rupture zones of large South American earthquakes and some predictions.J.Geophys.Res.，1972，77：2 087－2 103

[21] Darwin C.Journal and Remarks 1832－1836：Narrative of the Surveying Voyages of His Majesty’s Ships Adventure and Beagle between the Years 1826 and 1836，Describing their Examination of the Southern Shores of South America，and the Beagle’s Circumnavigation.Vol.3，370－381（Henry Colburn，1839）

[22] Barrientos S E.Dual seismogenic behavior：the 1985 Central Chile earthquake.Geophys.Res.Lett.，1995，22：3 541－3 544

[23] Plafker G，Savage J C.Mechanism of the Chilean earthquake of May 21 and 22，1960.Geol.Soc.Am.Bull.，1970，81：1 001－1 030

[24] Cifuentes I L.The 1960 Chilean earthquake.J.Geophys.Res.，1989，94（B1）：665－680

[25] Tassara A，Goetze H－J，Schmidt S，et al.Three－dimensional density model of the Nazca plate and the Andean continental margin.J.Geophys.Res.，2006，111，B09404.doi：10.1029／2005JB003976

[26]Hayes G.Finite Fault Model.Updated Result of the Feb 27，2010 MW8.8 Maule，Chile Earthquake.（http：∥earthquake.usgs.gov／earthquakes／eqinthenews／2010／us2010tfan／finite＿fault.php）（US Geological Survey／NEIC，2010）

[27] Shao G，Li X，Liu Q，et al.Preliminary slip model of the Feb 27，2010 MW8.9 Maule，Chile Earthquake.（http：∥www.geol.ucsb.edu／faculty／ji／big＿earthquakes／2010／02／27／chile＿2＿27.Html）（University of California Santa Barbara，2010）

[28] Sladen A，Owen S.Preliminary Model Combining Teleseismic and GPS Data 02／27／2010（MW8.8），Chile.（http：∥tectonics.caltech.edu／slip＿h(yuǎn)istory／2010＿chile／prelim－gps.html）（Caltech／JPL，2010）

[29] Hanka W.GEOFON Extended Virtual Network.（http：∥geofon.gfz－potsdam.de／geofon／）（GFZ，German Research Centre for Geosciences，2010）

[30] Moreno M S，Bolte J，Klotz J，et al.Impact of megathrust geometry on inversion of coseismic slip from geodetic data：Application to the 1960 Chile earthquake.Geophys.Res.Lett.，2009，36，L16310.doi：10.1029／2009GL039276

[31] Cloos M.Thrust－type subduction－zone earthquakes and seamount asperities：A physical model for seismic rupture.Geology，1992，20：601－604

[32] Kaneko Y，Avouac J－P，Lapusta N.Towards inferring earthquake patterns from geodetic observations of interseismic coupling.Nature Geosci.，2010，3：363－396

[33] Moreno M S，Klotz J，Melnick D，et al.Active faulting and heterogeneous deformation across a megathrust segment boundary from GPS data，south central Chile（36°S－39°S）.Geochem.Geophys.Geosyst.，2008，9，Q12024.doi：10.1029／2008GC002198

[34] Angermann D，Klotz J，Reigber C.Space－geodetic estimation of the Nazca－South America Euler vector.Earth Planet.Sci.Lett.，1999，171（3）：329－334

P315；

A；

10.3969／j.issn.0235－4975.2011.12.005

2011－09－22。