4地球物理大地測(cè)量聯(lián)合反演理論與應(yīng)用.ppt

1、地球物理反演理論地球物理反演理論課程組Balancing the plate motion budget in the South Island/New Zealand using GPS, geological and seismological data基于GPS、地質(zhì)學(xué)和地震數(shù)據(jù)對(duì)新西蘭南島板塊運(yùn)動(dòng)的平衡預(yù)算Laura M. Wallace， Geophys. J. Int. (2006) doi: 10.1111/j.1365-246X.2006.03183.x地球物理大地測(cè)量(聯(lián)合)反演理論與應(yīng)用之五主要內(nèi)容1，背景知識(shí)介紹2，反演模型與方法3，反演結(jié)果與分析4，討論與結(jié)論GPS網(wǎng)速

2、度的作用與不足：GPS數(shù)據(jù)能夠定量給出累積的應(yīng)變量及滑移的分解狀況 Wallace et al. 2004a證實(shí)，有很多的構(gòu)造解釋模型都能很好的擬合GPS數(shù)據(jù)；這種不唯一的原因是： 板塊邊界區(qū)域內(nèi)的不痛類型速度場(chǎng)存在較強(qiáng)的折； 俯沖區(qū)的復(fù)雜構(gòu)造及其活動(dòng)斷層；1、背景知識(shí)介紹改進(jìn)途徑和模型方法：利用這3種數(shù)據(jù)，采用擬合估計(jì)法：彈性旋轉(zhuǎn)塊體模型擬合法 (McCaffrey 1995, 2002). 斷層滑移速率和方位：塊體間的相對(duì)運(yùn)動(dòng) (Morgan 1968); GPS測(cè)量得到的速度場(chǎng)：方位，參考框架，塊體間斷層的 間震期耦合。對(duì)NEW Zealand的已有研究： 活動(dòng)斷層研究Active fa

3、ulting studies: (Wellman 1953; Lensen 1968; Berryman et al. 1992; Van Dissen & Berryman 1996)； 地震學(xué)家們: 估計(jì)震源參數(shù)；本研究的重點(diǎn)： 太平洋澳大利亞板塊間的相對(duì)運(yùn)動(dòng) (Beavan et al. 2002)； 斷層滑移速率，間震期的耦合程度；同時(shí)對(duì)其中的三個(gè)方面做了探討。俯沖/純走滑會(huì)聚擠壓右旋走滑純走滑旋轉(zhuǎn)俯沖backarc rifting活動(dòng)構(gòu)造背景數(shù)據(jù)分析1994至2004年間，15個(gè)區(qū)域GPS測(cè)量站；有地震的區(qū)域，則做相關(guān)處理；獲取站速度及其不確定性，將其轉(zhuǎn)換至澳大利亞塊體參考框架下，處

4、理方法同Wallace et al. (2004a);2個(gè)取舍標(biāo)準(zhǔn)：觀測(cè)時(shí)間與周期、方差協(xié)方差判斷其他數(shù)據(jù)： 太平洋板塊上和澳大利亞板塊上的GPS觀測(cè)數(shù)據(jù)： 建立澳大利亞板塊參考框架； 對(duì)運(yùn)動(dòng)學(xué)模型提供邊界條件； 斷層上的地震滑動(dòng)矢量方位角：對(duì)形變的運(yùn)動(dòng)學(xué)提供約束 地質(zhì)觀測(cè)資料：對(duì)運(yùn)動(dòng)學(xué)模型提供更強(qiáng)的約束二、反演模型方法基于構(gòu)造板塊長(zhǎng)期旋轉(zhuǎn)而被觀測(cè)到的GPS速度和板塊邊界斷層耦合而產(chǎn)生的彈性應(yīng)力而建立的模型：An elastic、rotating block approach (McCaffrey 1995, 2002)。在最優(yōu)擬合GPS速度、地震滑動(dòng)矢量和地質(zhì)斷層滑動(dòng)速率和方位角的條件上，獲

5、得彈性板塊的角速度和斷層的耦合系數(shù)，應(yīng)用的方法是模擬退火法的極小化。最小化數(shù)據(jù)殘差：由凱方的減小來確定。該方法也能使我們把多種GPS速度數(shù)據(jù)統(tǒng)一到同一參考框架 (澳大利亞板塊框架)內(nèi)。依據(jù)已知的活動(dòng)斷層劃分成很多構(gòu)造板塊:對(duì)于那些由于觀測(cè)的問題，而無法區(qū)分的多重?cái)鄬訁^(qū)域，采用近似的地殼形變分布劃分;有證據(jù)表明,有些構(gòu)造板塊內(nèi)部的斷層(無法確定的)有很小的地殼形變反應(yīng)到斷層上，但由于太小(1mm/yr)而用GPS技術(shù)無法檢測(cè)到。南島的劃分區(qū)域劃分的結(jié)果是一個(gè)純運(yùn)動(dòng)量，用來描述摩擦導(dǎo)致的斷層長(zhǎng)期的不動(dòng)的過程。 表達(dá)式：其中：V 長(zhǎng)期滑動(dòng)速率(經(jīng)歷了很多地震周期)，Vc為短期蠕動(dòng)。=0，表示只有蠕動(dòng)

6、;=1，表示在震間周期內(nèi)沒有蠕動(dòng)。 注：估計(jì)的只是網(wǎng)點(diǎn)上的值，其他的可以通過雙線性內(nèi)插獲得。我們不考慮它隨時(shí)間的變化。耦合系數(shù)對(duì)值的限制新西蘭島的地殼斷層尚未發(fā)掘震間蠕動(dòng)的證據(jù)相對(duì)移動(dòng)矢量矩陣；滑動(dòng)速率秩虧矩陣；幾個(gè)計(jì)算的量后滑模型 (Savage 1983),；彈性半空間的位錯(cuò)模型 (Okada 1985) ；二者均需要定義耦合可能發(fā)生出的位錯(cuò)幾何參數(shù)計(jì)算彈性形變對(duì)速度場(chǎng)的貢獻(xiàn)就最優(yōu)模型而言，2= 1.17，這表明數(shù)據(jù)擬合程度在不確定性水平內(nèi)。使用了818個(gè)GPS水平速度、30個(gè)地震滑動(dòng)矢量或斷層滑動(dòng)方位角以及35個(gè)地質(zhì)學(xué)滑動(dòng)速率，對(duì)217個(gè)參數(shù)(包括北島的138個(gè)角速度的自由參數(shù)和斷層耦合

7、系數(shù)，Wallace et al. 2004a) 做了估計(jì)。反演的自由參數(shù)包括：旋轉(zhuǎn)參數(shù)、應(yīng)變率張量、。四、反演結(jié)果與分析斷層滑動(dòng)速率虧損的估計(jì)地震學(xué)估計(jì)結(jié)果：12 km (Leitner et al. 2001)；大地測(cè)量學(xué)估計(jì)結(jié)果：58 km (Beavan et al. 1999)； 原因：可能是沒有考慮滑動(dòng)速率虧損沿走勢(shì)是變化的， 同時(shí)也受到臨近斷層的影響。 考慮以上兩點(diǎn)，對(duì)于大多數(shù)的the Alpine Fault，估計(jì)獲得了一個(gè)18KM的震間耦合深度(=0.70.85)。 可能=1的深度小于18KM而但仍然大于地震研究估計(jì)的12KM。 南島構(gòu)造板塊間可能的內(nèi)部應(yīng)力為了估計(jì)CANT

8、和 SALP塊體間Alpine斷層上可能的永久應(yīng)力分布，我們?cè)谧罴褦M合模型下，反演了這些板塊的水平應(yīng)變率張量；估計(jì)的有些板塊間的永久應(yīng)力是彈性的，由于構(gòu)造板塊內(nèi)的耦合斷層。然而，考慮GPS數(shù)據(jù)的不確定性，很難區(qū)分無震滑動(dòng)板塊內(nèi)斷層的永久地殼形變與震間耦合。估計(jì)應(yīng)力速率后，對(duì)擬合結(jié)果的改進(jìn)作用不明顯。 就CANT塊體而言，塊體增厚導(dǎo)致約0.1 to 0.2 mm/yr 的抬升；在有些構(gòu)造板塊，利用速度殘差估計(jì)殘余應(yīng)變率，估計(jì)的結(jié)果表明: 假設(shè)塊體內(nèi)部未發(fā)生形變是合理的，符合模型假設(shè)要求。但也有少數(shù)例外。由于Marlborough臨近斷層的滑動(dòng)速率秩虧引起的彈性應(yīng)力，使得只依靠GPS數(shù)據(jù)很難確定它

9、的滑動(dòng)速率，這一區(qū)域的旋轉(zhuǎn)具有很大的不確定性。地質(zhì)資料作為硬性條件被應(yīng)用到反演中，這對(duì)于解決有些參數(shù)是極其關(guān)鍵的，因?yàn)樗鼈儾皇軓椥詰?yīng)力成分的影響，不像GPS數(shù)據(jù)。 GPS速度的研究表明大約一半的總滑動(dòng)速率可以由地質(zhì)研究解釋，余下的滑動(dòng)，肯定發(fā)生在難以描述和發(fā)現(xiàn)的斷層中。 PPAFZ的速率是之前用地質(zhì)資料估計(jì)的速率的2倍，這是因?yàn)樗麄冎肮烙?jì)的時(shí)候，只考慮了大斷層，然而PPAFZ是一個(gè)復(fù)雜的、廣闊的斷層區(qū)域，還有很多包含右旋的斷層(e.g. Pettinga et al. 2001)。 Marlborough地區(qū)的板塊運(yùn)動(dòng)和斷層滑動(dòng)速率古地磁和地質(zhì)數(shù)據(jù)表明，800W年來，Marlborough塊

10、體的北部相對(duì)于太平洋和澳大利亞板塊邊界有35的順時(shí)針方向的旋轉(zhuǎn)。Marlborough 塊體表現(xiàn)出與北島塊體很大的不同的旋轉(zhuǎn)點(diǎn)，表明Marlborough 塊體中部和 北島塊體之間是不連續(xù)貫通的；Alpine Fault的滑動(dòng)速率總體來說，與地質(zhì)資料估計(jì)的值比較一致。從北端往南，滑動(dòng)速率在不斷的增加。 到Haast南部，運(yùn)動(dòng)又變?yōu)榧冏呋俾氏陆禐椋?24 mm/yr. 這與地質(zhì)資料預(yù)測(cè)的差不多。這種減少，表明有些滑動(dòng)從Alpine Fault 中部轉(zhuǎn)化到了 Fiordland 俯沖區(qū)。 我們預(yù)測(cè)Alpine Fault中部的收斂運(yùn)動(dòng)是5 mm /yr，比地質(zhì)資料估計(jì)的要少12 mm/yr。從

11、我們的模型中，可以獲知， Alpine Fault中部的走滑運(yùn)動(dòng)成分可以解釋80%的太平洋澳大利亞板塊間的相對(duì)運(yùn)動(dòng)。Southern Alps的地殼形變分布討論了一個(gè)對(duì)Central Southern Alps region可能更適合的模型.我們估計(jì)Alpine Fault的運(yùn)動(dòng)的走滑處于，活動(dòng)斷層研究估計(jì)的平均值的最大值。對(duì)不活動(dòng)的主分離斷層區(qū)的構(gòu)造觀測(cè)，可以獲得額外的證據(jù)，表明Southern Alps的地殼形變分布發(fā)生在過去。就Alpine 斷層?xùn)|側(cè)100千米的斷層，包含有巨大的剪切形變，可能不能正確的模擬出這種彈性剪切形變，這是因?yàn)閠he Alpine Fault的震間滑動(dòng)速率秩虧的存

12、在。這種無法辨別的壓力可能恰好解釋了，為什么我們的最佳擬合模型獲得了更深的震間耦合深度，與震源層相比。 為了確定Southern Alps的地殼形變分布，把SALP分成了：Southern Alps and Southland 塊體。添加2個(gè)限制： Alpine斷層中部的走滑速率為：22 and 27 mm/yr； Southern Alps and Canterbury/Otago 塊體邊界的總相對(duì)運(yùn)動(dòng)少于7 mm/yr。這2個(gè)限制條件的使得可以很好的擬合the central South Island GPS數(shù)據(jù)。在這種條件下的模型的 = 1.63，比最佳模型的 =1.17要大?？紤]一個(gè)更

13、復(fù)雜的參數(shù)應(yīng)力模式，可能對(duì)數(shù)據(jù)的擬合會(huì)和前面的最佳擬合模型一致。Southern Alps東邊山麓斷層的地殼形變運(yùn)動(dòng)學(xué)模型預(yù)測(cè)SALP/CANT 邊界有36 mm/yr的右旋收斂速率。前人做的研究，太多是地質(zhì)學(xué)研究，都致力于估計(jì)聚斂速率。我們認(rèn)為，在這個(gè)區(qū)域的斷層中，存在右旋走滑和聚斂成分，使得總的速率比地質(zhì)研究估計(jì)的值要大23倍。 地質(zhì)學(xué)研究獲得小的值的原因可能是：沒有估計(jì)右旋速率，或者是考慮的活動(dòng)斷層數(shù)量要少一些。本文的運(yùn)動(dòng)學(xué)模型也預(yù)測(cè)，在CANT 塊體 和太平洋板塊間有12 mm/yr的收縮(有點(diǎn)右旋); 而有人預(yù)測(cè)是3mm/yr. 這種差別，可以從 Alpine 和Marlboroug

14、h 斷層系統(tǒng)存在的震間滑動(dòng)秩虧得到解釋。Fiordland 區(qū)域運(yùn)動(dòng)學(xué)模型對(duì)Fiordland 區(qū)域的估計(jì)： 受到Fiordland 俯沖區(qū)的GPS速度、地震滑動(dòng)矢量的限制。南島構(gòu)造板塊獨(dú)立性的統(tǒng)計(jì)檢驗(yàn)如果假定SALP, CANT, 及 FIOR塊體獨(dú)立旋轉(zhuǎn)，則觀測(cè)數(shù)據(jù)的擬合結(jié)果不如另一種旋轉(zhuǎn)模式好。SALP/FIOR: 基于GPS速度場(chǎng)，尚不能確定FLOR塊體是否與SALP塊體相關(guān)。如果FLOR和SALP塊體間發(fā)生形變，則其形變必然很小。基于這些檢驗(yàn)，Alpine斷層的東南島并未表明是剛性太平洋板塊的一部分。與最優(yōu)擬合模型比較，假定(FIOR,SALP及CANT塊體作為太平洋板塊的一部分進(jìn)行

15、檢測(cè)，對(duì)數(shù)據(jù)的擬合上出現(xiàn)了很大的偏差；SALP 和 CANT 塊體作為太平洋板塊的一部分的檢測(cè)，同樣出現(xiàn)了很差的擬合效果。如果CANT塊體作為太平洋板塊的一部分，F(xiàn)檢驗(yàn)表明可以對(duì)觀測(cè)數(shù)據(jù)集擬合的較好，但其擬合效果仍差于本文給出的最優(yōu)擬合效果。這些似乎表明CANT塊體較其他兩個(gè)塊體的運(yùn)動(dòng)更加近似于太平洋板塊的運(yùn)動(dòng)，但仍然與太平洋板塊的運(yùn)動(dòng)失相關(guān)。討論與結(jié)論南島滑動(dòng)分割樣式的變化 板塊邊界類型的變化，從Hikurangi 邊緣的俯沖轉(zhuǎn)換為Marlborough邊緣的純走滑，南Alps的斜滑使得新西蘭板塊邊界區(qū)域滑動(dòng)分解模式發(fā)生十分巨大的變化。北島的情況：北島主要發(fā)生構(gòu)造旋轉(zhuǎn)，而走滑斷層效應(yīng)為次級(jí)作

16、用；北島區(qū)域發(fā)生滑移分解現(xiàn)象，從垂直于Hikurangi俯沖邊緣的俯沖地震(Webb & Anderson 1998)和上覆板塊區(qū)域內(nèi)的走滑斷層的地質(zhì)學(xué)證據(jù)可以得到證實(shí) (Cashman et al. 1992)；本研究認(rèn)為北島南端部分較大的走滑速率是由用于調(diào)整平行于板塊邊界的旋轉(zhuǎn)運(yùn)動(dòng)速率降低引起的。南島的情況：Marlborough區(qū)域: 相對(duì)板塊運(yùn)動(dòng)的平行于和垂直于邊緣分量在很大程度上是體現(xiàn)為分離結(jié)構(gòu)發(fā)生的。Marlborough 中央和南部區(qū)域: 板塊旋轉(zhuǎn)幾乎可以不考慮進(jìn)滑動(dòng)劃分；Marlborough北部的斷層:由于受到來自北島的走勢(shì)的作用，可能包含了很大的收縮成分。(e.g. Jor

17、dan Thrust, Van Dissen & Yeats 1991; northern Awatere Fault, Little et al. 1998; London Hills Fault,Townsend&Little 1998; Clarence Fault, Nicol & Van Dissen 2002).PPAFZ區(qū)域的一部分?jǐn)鄬?，包含右旋扭壓。小?00W年的古地磁獲得的northeastern Marlborough的順時(shí)針旋轉(zhuǎn)證據(jù)表明: 旋轉(zhuǎn)在這個(gè)區(qū)域的滑動(dòng)劃分上起到了作用；雖然利用GPS技術(shù)，目前還無法證實(shí)是否這些旋轉(zhuǎn)還在繼續(xù)（這可能是由于這個(gè)區(qū)域比較小的原因）。假

18、如古地磁獲得的這種旋轉(zhuǎn)還在繼續(xù)，Marlborough 北東段可能是 北島構(gòu)造邊界的過渡區(qū)域；同時(shí)在Marlborough中部也可以觀測(cè)到接近純走滑的運(yùn)動(dòng)。 Alpine斷層：包含了大約75%的相對(duì)板塊運(yùn)動(dòng)，也有收聚和走滑成分；其余的量值分化在Alpine斷層?xùn)|部的走滑、傾斜滑動(dòng)和收縮的斷層中。大約63%的相對(duì)板塊運(yùn)動(dòng)的平行于邊緣的分量，發(fā)生在Alpine斷層；總體來說，北島的滑動(dòng)分割過程，是以板塊旋轉(zhuǎn)為特征的；南島這樣的旋轉(zhuǎn)由于是以遠(yuǎn)點(diǎn)為中心的而導(dǎo)致了低速度梯度。因此，南島的相對(duì)板塊運(yùn)動(dòng)更多的是，直接體現(xiàn)到斷層上，有一些是比較間接的。New Zealand的板塊旋轉(zhuǎn)的驅(qū)動(dòng)機(jī)制南島構(gòu)造板塊的垂

19、直軸旋轉(zhuǎn)速率，在數(shù)量和形式(相對(duì)于澳大利亞板塊是逆時(shí)針的)上和與太平洋板塊是相似的。這意味著南島的運(yùn)動(dòng)，受到太平洋板塊的影響更大。 有趣的是，SALP 塊體相對(duì)澳大利亞的垂直速率 (0.35/Myr, anti-clockwise) 比 CANT 塊體相對(duì)澳大利亞塊體的旋轉(zhuǎn)速率(0.99/Myr, anti-clockwise)要小。這表明：越接近澳大利亞板塊，由于邊壓的增加，板塊的運(yùn)動(dòng)受到澳大利亞板塊的影響越大。驅(qū)動(dòng)原因，可能是該區(qū)域底下的地幔運(yùn)動(dòng)所致.Northeastern Marlborough: 古地磁和地質(zhì)研究表明, 在過去的800W年里, 累積發(fā)生了巨大的順時(shí)針旋轉(zhuǎn)運(yùn)動(dòng), 達(dá)到了

20、35度。在北島東部, 古地磁和大地測(cè)量資料表明, 有快速的順時(shí)針旋轉(zhuǎn)。Wallace et al. (2004a)認(rèn)為，這是由于太平洋俯沖浮力向南增加，這一俯沖過程的影響。目前, 一般認(rèn)為L(zhǎng)amb(1988)的浮塊模型解釋了 Marlborough北東區(qū)域的順時(shí)針旋轉(zhuǎn)。但是, 我們更相信，這是由于相對(duì)俯沖過程的作用, 與導(dǎo)致北島構(gòu)造板塊旋轉(zhuǎn)的類似(Wallace et al. 2004a)。假如正如所假設(shè)一樣, 那么北島的構(gòu)造板塊通過Cook Strait與Marlborough北東部應(yīng)該是相鏈接的。然而, 當(dāng)前在 Marlborough北東部的GPS速度測(cè)量的分辨率和精度, 都不能讓我們確定

21、這些順時(shí)針旋轉(zhuǎn)是否還在繼續(xù)。板塊邊界區(qū)域演化的主要影響因素Chatham Rise沿著新西蘭板塊邊線以2837 mm/yr的速率向西南遷移。也使得Hikurangi 俯沖區(qū)相對(duì)Australia向南蔓延。這種蔓延, 可能證實(shí)了導(dǎo)致MFS活動(dòng)性向南遷移, 這也導(dǎo)致了新的斷層出現(xiàn), 比如PPAFZ (Berryman et al. 1992; Little & Roberts 1997)。同時(shí), 俯沖帶以及隨后的MFS區(qū)域的斷層向南遷移, 從MFS地區(qū)南部發(fā)現(xiàn)的新斷層這一情況得到證實(shí), 比如PPAFZ (Cowan et al. 1996)。此外, 也出現(xiàn)了一個(gè)伴隨的滑動(dòng)轉(zhuǎn)移, 從Alpine斷層

22、向南東向傳播到; 這也有效的減短了Alpine 斷層中間部分的快速滑動(dòng)的長(zhǎng)度(Berryman et al. 1992)。 根據(jù) Little & Roberts (1997)， 我們額外認(rèn)為: Alpine 斷層/MFS 區(qū)域以相似速率向南西遷移；因此，Alpine斷層的最南端部分相對(duì)Marlborough斷層而言, 與200W年前相比, 遷移了大約5674KM；假如以上的確實(shí)發(fā)生了, 那么Marlborough 斷層系統(tǒng)在大約1.41.9百萬前開始形成有新的斷層。這與該區(qū)域的一些觀測(cè)研究結(jié)果相一致。Chatham Rise 和 Hikurangi 高原的南向遷移, 同樣也影響了北島的構(gòu)造過程。Wallace et al. (2004a)認(rèn)為: 北島東部的快速順時(shí)針旋轉(zhuǎn)的發(fā)生, 是由于Chatham Rise (a continental fragment)俯沖到Hikurangi 高原 (buoyant oceanic plateau)這一變化導(dǎo)致的。這也導(dǎo)致相對(duì) Chatham Rise這一“固定點(diǎn)”(相對(duì)于太平洋板塊的)發(fā)生旋轉(zhuǎn), 以及TVZ地區(qū)的后弧分裂。而由于相對(duì)于澳大利亞板塊的旋轉(zhuǎn)點(diǎn)位置的存在, 也使得TVZ地區(qū)發(fā)生分裂運(yùn)動(dòng)。