畢奔騰， 胡祥云*， 李麗清， 張恒磊， 劉雙， 蔡建超

1 中國(guó)地質(zhì)大學(xué)地球物理與空間信息學(xué)院， 武漢　430074 2 地球內(nèi)部多尺度成像湖北省重點(diǎn)實(shí)驗(yàn)室， 武漢　430074 3 中國(guó)地質(zhì)大學(xué)資源學(xué)院， 武漢　430074

?

青藏高原東北部多尺度重力場(chǎng)及其地球動(dòng)力學(xué)意義

畢奔騰1，2， 胡祥云1，2*， 李麗清3， 張恒磊1，2， 劉雙1，2， 蔡建超1，2

1 中國(guó)地質(zhì)大學(xué)地球物理與空間信息學(xué)院， 武漢430074 2 地球內(nèi)部多尺度成像湖北省重點(diǎn)實(shí)驗(yàn)室， 武漢430074 3 中國(guó)地質(zhì)大學(xué)資源學(xué)院， 武漢430074

摘要由于受到SN和EW雙向擠壓作用的綜合影響，處于年輕高原與古老地塊交接區(qū)的青藏高原東北部巖石圈強(qiáng)烈變形，構(gòu)造活動(dòng)十分活躍.為了整體性地了解青藏高原東北部重力場(chǎng)和深部動(dòng)力學(xué)機(jī)制，本文基于EGM2008全球重力場(chǎng)模型數(shù)據(jù)，利用小波多尺度分析獲得不同尺度的重力異常信息；同時(shí)反演了研究區(qū)的地殼厚度；通過(guò)構(gòu)建穿越龍門山造山帶和西秦嶺造山帶兩條剖面的巖石圈密度結(jié)構(gòu)模型，分析了地殼上地幔內(nèi)不同介質(zhì)的分布特征.研究結(jié)果表明，青藏高原東北部巖石圈顯示十分復(fù)雜的塑性體的特征，其重力異常走向多以EW或SSE為主，反映了高原巖石圈物質(zhì)向東運(yùn)移的趨勢(shì)；地殼厚度由西向東逐漸減薄，邊緣造山帶深部并沒(méi)有發(fā)現(xiàn)“山根”痕跡，結(jié)合該地區(qū)低重力的特征，推測(cè)西秦嶺—松潘構(gòu)造結(jié)巖石圈發(fā)生過(guò)大規(guī)模地幔流底侵作用；地幔流上涌的動(dòng)力可能來(lái)源于印度板塊向歐亞大陸板塊俯沖，激發(fā)了地幔流體側(cè)向移動(dòng)，在揚(yáng)子地臺(tái)和華北地臺(tái)附近受到堅(jiān)硬巖石圈的阻擋而被迫上移，并因此造成龍門山與西秦嶺的隆升.

關(guān)鍵詞青藏高原東北部； 重力場(chǎng)； 巖石圈結(jié)構(gòu)； 密度模型； EGM2008

1引言

自70—50 Ma以來(lái)，印度板塊向北與歐亞板塊持續(xù)強(qiáng)烈碰撞，至少有1360 km的南北向縮短的物質(zhì)被吸收(Yin and Harrison，2000)，物質(zhì)的積累導(dǎo)致了青藏高原的隆升(Molnar and Tapponnier，1975；England and Houseman，1986)，形成了全球海拔最高的獨(dú)特地域單元.青藏高原發(fā)育有厚達(dá)60～80 km的地殼(高銳等，2009)，由于其構(gòu)造的復(fù)雜性，碰撞過(guò)程的吸收機(jī)制與地殼增厚的方式也存在多種解釋.在高原腹地，一系列的俯沖、逆沖帶可以很好地對(duì)應(yīng)地殼縮短增厚模型.而在邊緣造山帶，一系列走滑斷裂的存在則與大陸塊體側(cè)向擠出模型一致(Tapponnier et al.，2001)；中下地殼顯示低速高導(dǎo)的性質(zhì)又為地殼流模型提供了依據(jù)(Royden et al.，1997；Clark and Royden，2000；Bai et al.，2010).作為印度歐亞板塊碰撞作用由南北向東西轉(zhuǎn)換的重要場(chǎng)所，青藏高原東北部(圖1)具有與一般線性造山帶不同的造山極性及構(gòu)造體制(許志琴等，1992).特別是西秦嶺—松潘大陸構(gòu)造結(jié)屬中國(guó)大陸完成其主體拼合過(guò)程中形成的地殼尺度上的巨型構(gòu)造，是青藏高原與周圍穩(wěn)定地體交接轉(zhuǎn)換的關(guān)鍵構(gòu)造域(張國(guó)偉等，2004).以往的地球物理深部探測(cè)在該區(qū)做了很多研究，包括人工反射地震(Liu et al.，2006；嘉世旭和張先康，2008；Zhang et al.，2011，2013；Gao et al.，2014)、天然地震(Vergne et al.，2003；薛光琦等，2003；Li et al.，2011, 2014；Pan and Niu，2011；Tian and Zhang，2013)、大地電磁測(cè)深(趙國(guó)澤等，2004；Xiao et al.，2013)、GPS形變監(jiān)測(cè)(Zhang et al.，2004；He et al.，2013)等，獲得了青藏高原東北部巖石圈的速度結(jié)構(gòu)、泊松比分布、SKS剪切波分裂、電性結(jié)構(gòu)以及水平和垂直位移分量等重要地球物理參數(shù)，并且對(duì)殼幔中的低速高導(dǎo)層及動(dòng)力學(xué)模式進(jìn)行了一系列的探討.由于受復(fù)雜的高原環(huán)境影響，上述研究多集中于沿特定地質(zhì)構(gòu)造方向的二維剖面測(cè)量或者是點(diǎn)距較大的固定臺(tái)站，這對(duì)于研究大尺度的區(qū)域地質(zhì)構(gòu)造稍顯不足.重力方面，前人分析了研究區(qū)的斷裂體系、變形特征等，取得了許多成果(Zhang et al.，2010；孟小紅等，2012；陳石等，2013).然而高原的地面重力觀測(cè)數(shù)據(jù)分布不均且十分有限，上述多集中于將龍門山造山帶、西秦嶺造山帶、青藏高原東北緣分隔開來(lái)單獨(dú)研究.事實(shí)上，受印度板塊向歐亞板塊俯沖遠(yuǎn)程效應(yīng)的影響，龍門山與西秦嶺的隆升時(shí)限及所處的構(gòu)造環(huán)境有諸多相似性，將其統(tǒng)一到青藏高原東北部向四川盆地、鄂爾多斯盆地周邊擴(kuò)展的動(dòng)力學(xué)背景下，可進(jìn)一步查明青藏高原東北部隆升變形機(jī)制以及川陜甘強(qiáng)震區(qū)的孕震構(gòu)造機(jī)理.

圖1　青藏高原東北部大地構(gòu)造輪廓圖左下圖紅色圖框?yàn)檠芯繀^(qū)位置.圖中黑色實(shí)線為文中涉及的斷裂，灰色細(xì)線為本文忽略斷裂.斷裂數(shù)據(jù)源自鄧起東等(2002)；SCB：四川盆地；SGB：松潘—甘孜地塊；WQL：西秦嶺造山帶；QLB：祁連山地塊；OB：鄂爾多斯盆地；JS：金沙江縫合帶；AKMS：阿尼瑪卿—昆侖—木孜塔格縫合帶；SQS：南祁連縫合帶；GY：甘孜—玉樹斷裂；XSH：鮮水河斷裂；LMS：龍門山斷裂；MJ：岷江斷裂；DBL：迭部—略陽(yáng)斷裂；NWQL：西秦嶺北緣斷裂；NKL：北昆侖斷裂；HY：海原斷裂.Fig.1　Outline map of the northeastern Tibetan plateau and the surrounding areaThe study region marked by the red frame in the below left map. The black solid lines represent faults involved, and the gray filaments represent faults omitted. Faults data depicted after Deng et al.(2002); SCB: Sichuan basin; SGB: Songpan-Garzê block; WQL: West Qinling Orogenic belt; QLB: Qilian block; OB: Ordos basin; JS: Jinsha River suture; AKMS: A′nyêmaqên-Kunlun-Muztag suture; SQS: South Qilian suture; GY: Garzê-Yushu fault; XSH: Xianshui River fault; LMS: Longmen Shan fault; MJ: Minjiang fault; DBL:Diebu-Lueyang fault; NWQL: Northern edge of the Western Qinling fault; NKL: North Kunlun fault; HY: Haiyuan fault.

隨著衛(wèi)星技術(shù)的發(fā)展，快速高精度測(cè)量全球重力場(chǎng)成為可能，繼而建立了各種高分辨率的重力場(chǎng)模型.EGM2008(Earth Gravitational Model 2008)是近年來(lái)應(yīng)用較為廣泛的全球超高階地球重力場(chǎng)模型(Pavlis et al.，2008)，它的球諧展開階次為2160，空間分辨率約為5′(9 km).其構(gòu)建過(guò)程是用GRACE衛(wèi)星數(shù)據(jù)建立前60階球諧系數(shù)，加上EGM96重力模型的61～360階球諧系數(shù)，再加上剩余地形重力模型的361～2160階球諧系數(shù)計(jì)算得到的地球重力異常.與實(shí)測(cè)地面重力數(shù)據(jù)相比，EGM2008模型數(shù)據(jù)在青藏高原東北部的誤差約為11.6 mGal(章傳銀等，2009)，其異常的變化趨勢(shì)上大體一致(付廣裕等，2013).值得說(shuō)明的是，青藏高原地面重力點(diǎn)位分布不均，其反映的異常信息遠(yuǎn)不如模型，地球重力模型網(wǎng)格數(shù)據(jù)的精度被低估了(楊金玉等，2012).

為了能從宏觀上整體反映研究區(qū)深部構(gòu)造特征，本文基于EGM2008模型的自由空氣重力異常數(shù)據(jù)，研究青藏高原東北部不同尺度的重力場(chǎng)分布，推測(cè)橫縱向地殼上地幔的物質(zhì)成分及構(gòu)造特征，進(jìn)而探討其動(dòng)力學(xué)涵義.

2青藏高原東北部及周緣地區(qū)重力異常

圖2顯示了青藏高原東北部及周緣地區(qū)自由空氣重力異常，可以看出，其形態(tài)與地形圖相似.在山脈分布的高地勢(shì)地區(qū)顯示正異常，尤其是在劇烈切割的高山峽谷地區(qū)，由于山體形態(tài)復(fù)雜，受密度不均勻地質(zhì)體的影響，自由空氣重力異常特征變得十分復(fù)雜；而盆地中因在近地表分布有中、新生界的低密度地層，因而呈現(xiàn)負(fù)值，負(fù)異常的幅值和形態(tài)與盆地內(nèi)沉積層的厚薄和分布呈正相關(guān)關(guān)系.四川盆地和鄂爾多斯盆地顯示弱的負(fù)異常，變化較為平緩，表明地質(zhì)構(gòu)造比較穩(wěn)定，負(fù)異常也說(shuō)明盆地在均衡補(bǔ)償過(guò)程中質(zhì)量虧損.高原腹地的正異常走向由EW向向東轉(zhuǎn)為SE向，幅值在0～100 mGal內(nèi)變化，而在邊緣造山帶異常值在100～200 mGal之間，龍門山斷裂帶附近達(dá)到最大300 mGal.印度大陸持續(xù)向北推擠造成了青藏高原的抬升，形成高的地勢(shì)，正是由于高原內(nèi)部軟流圈流體受到了地形對(duì)巖石圈的垂向重力均衡作用，側(cè)向流動(dòng)，遇到揚(yáng)子剛性地塊的阻擋，在龍門山斷裂帶附近積累大量物質(zhì)，因而形成了強(qiáng)烈的正異常帶.

圖2　研究區(qū)自由空氣重力異常圖Fig.2　Free air gravity anomaly map of the study region

在自由空氣重力數(shù)據(jù)的基礎(chǔ)上結(jié)合ETOPO1數(shù)據(jù)地形改正后，獲得了整個(gè)研究區(qū)的布格重力異常分布(見圖3).由圖可見，研究區(qū)內(nèi)布格重力異常以負(fù)值為主，范圍在50～-500 mGal之間，自東向西逐漸衰減.青藏高原布格重力異常在-300 mGal以下，而四川盆地、鄂爾多斯盆地在-100 mGal以上.在西寧—汶川一線有一條明顯的重力梯度帶，特別是梯度帶南部龍門山斷裂帶附近，由西向東重力值劇烈遞增，在不到100 km范圍內(nèi)重力值增加了100 mGal；而梯度帶北段則相對(duì)平緩得多，西寧—蘭州—海原一線地區(qū)呈現(xiàn)出過(guò)渡帶的性質(zhì).

3小波多尺度分析

重力異常反映的是地表下不同深度、不同規(guī)模、不同密度的地質(zhì)體對(duì)應(yīng)的重力異常疊加效應(yīng).小波多尺度分析方法是將信號(hào)分解成各種不同頻率成分，并且聚焦到任意細(xì)節(jié)上加以分析(Mallat，1989；侯遵澤和楊文采，1997).為了得到目的地質(zhì)體對(duì)應(yīng)的重力特征，可以利用小波多尺度分析的方法分離出反映其橫向、縱向上不同尺度地質(zhì)體產(chǎn)生的重力異常場(chǎng)(Jiang et al.，2012；孟小紅等，2012).本文通過(guò)對(duì)EGM2008模型的布格重力異常應(yīng)用Daubechies(db3)小波提取各個(gè)尺度的重力場(chǎng)細(xì)節(jié)信息；對(duì)于分解后的異常，應(yīng)用功率譜分析的原理計(jì)算其徑向?qū)?shù)功率譜，估計(jì)每一階異常的場(chǎng)源深度，所得結(jié)果見表1.1階和2階細(xì)節(jié)異常反映的是近地表的密度變化，在此不予討論.

3.1小尺度重力異常

3階和4階細(xì)節(jié)異常圖(圖4a、4b)反映了研究區(qū)上地殼的布格重力異常特征.四川盆地和鄂爾多斯盆地分別處于揚(yáng)子地臺(tái)和華北地臺(tái)，構(gòu)造穩(wěn)定，其基底之上覆蓋巨厚的中新生代沉積物，因此顯示出均勻的弱異常特征.松潘—甘孜地塊中部巴顏喀拉、可可西里地區(qū)廣泛分布三疊紀(jì)復(fù)理石沉積，厚度多在10 km以上，由于后期構(gòu)造活動(dòng)強(qiáng)烈，基底褶皺變形使沉積物重新改造，其重力異常形態(tài)表現(xiàn)為略微雜亂的弱異常的特征.與高原腹地相比，邊緣地區(qū)的異常形態(tài)更加復(fù)雜.阿尼瑪卿—昆侖—木孜塔格縫合帶(AKMS)作為東古特提斯洋在晚古生代向北俯沖閉合的場(chǎng)所(Zhu and Helmberger，1998)，東西綿延近1500 km，在異常圖上顯示其東段為負(fù)的線性異常，一直延伸到瑪沁斷裂以東.松潘—甘孜地塊東緣及東南緣則表現(xiàn)出一系列正負(fù)相伴的條帶異常.東緣是青藏高原向東運(yùn)動(dòng)與揚(yáng)子地臺(tái)碰撞的場(chǎng)所，其雜亂無(wú)序的異常方向體現(xiàn)了物質(zhì)受擠壓褶皺變形或破碎的特征，而四川盆地受到強(qiáng)烈的擠壓其重力異常形態(tài)依然穩(wěn)定，顯示出剛性塊體的特征.松潘甘孜東南緣接近于喜馬拉雅東構(gòu)造結(jié)，是青藏高原物質(zhì)流向東南方向轉(zhuǎn)折的場(chǎng)所，其不僅在地形上表現(xiàn)為強(qiáng)烈的褶皺變形、地殼縮短以及一系列走滑斷層，在重力異常圖上亦顯現(xiàn)SSE向和SN向的串珠狀異常，負(fù)異?？赡苁菐r石擠壓破碎的反映，而正異常則可能是由于巖漿巖沿破碎帶侵入形成的.鮮水河斷裂帶一直被認(rèn)為是現(xiàn)今東亞大陸最活躍的斷裂帶之一，其包括甘孜—玉樹(GY)、鮮水河(XSH)、小江、則木河和安寧河等次級(jí)斷裂.在3階和4階圖上，鮮水河斷裂和甘孜—玉樹斷裂表現(xiàn)為連續(xù)的SE和SSE向的負(fù)異常條帶，并且兩側(cè)伴有正異常.

圖3　研究區(qū)布格重力異常圖A-A′和B-B′兩條實(shí)線代表二維密度模型剖面切割的位置.Fig.3　Bouguer gravity anomaly map of the study region Solid lines (A-A′ and B-B′) represent the location of density modeling profiles.

階次場(chǎng)源似深度(km)1122354165346125

3.2中尺度重力異常

5階細(xì)節(jié)異常圖(圖4c)反映了研究區(qū)中下地殼的布格重力異常特征.由圖可知，青藏高原東北部大都顯示負(fù)異常，而周圍塊體，如四川盆地、鄂爾多斯盆地、祁連地塊以及柴達(dá)木盆地等都表現(xiàn)為突出的正異常，依舊為剛性塊體的特征.西秦嶺造山帶北緣、松潘—甘孜地塊東緣則呈現(xiàn)明顯的負(fù)異常，推測(cè)是巖石受壓應(yīng)力破碎或者存在地殼熔融物質(zhì).四川盆地正異常區(qū)向西延伸，越過(guò)龍門山斷裂帶，反映了松潘—甘孜地塊巖石圈向四川盆地仰沖的構(gòu)造形態(tài).

圖4　布格重力異常3～6階小波變換細(xì)節(jié)圖Fig.4　The 3rd to 6th order detailed images of Bouguer gravity anomaly from multi-scale wavelet analysis method

3.3大尺度重力異常

對(duì)數(shù)功率譜估計(jì)的6階細(xì)節(jié)異常所反映的平均場(chǎng)源深度為125 km，體現(xiàn)了研究區(qū)上地幔布格重力異常特征(圖4d).鮮水河斷裂帶和金沙江斷裂帶所圍限的正異常區(qū)，表現(xiàn)為剛性塊體的性質(zhì)，與周圍的負(fù)異常帶形成對(duì)比，推測(cè)該區(qū)曾經(jīng)為揚(yáng)子地臺(tái)的一部分，其被順時(shí)針旋轉(zhuǎn)的地幔流沖斷.研究區(qū)最為明顯的重力低值區(qū)出現(xiàn)在松潘甘孜東緣和西秦嶺造山帶，范圍較為連續(xù)，其東北緣幅值最低，向兩翼逐漸增大.鄂爾多斯盆地西南緣分布弧形的重力梯度帶，對(duì)應(yīng)于地理上的海原斷裂帶，分割了低重力的青藏高原東北緣和高重力的鄂爾多斯盆地，形成深達(dá)巖石圈尺度的超殼斷裂，而青藏高原殼幔物質(zhì)擠入到阿拉善地臺(tái)和鄂爾多斯盆地之間，在這樣一個(gè)狹小的空間內(nèi)勢(shì)必會(huì)積累巨大的應(yīng)變能，因此該區(qū)成為第四系以來(lái)構(gòu)造變動(dòng)十分強(qiáng)烈的地震活動(dòng)區(qū).

4地殼厚度模型

地球外部重力場(chǎng)是由地球內(nèi)部物質(zhì)分布所決定的，由于地殼和地幔存在著較大的密度差異，因此可以利用重力資料來(lái)確定地殼厚度(Feng et al.，2014).基于上述結(jié)論，本文對(duì)布格重力異常利用Parker-Oldenburg界面反演算法(Parker，1973；Oldenburg，1974)計(jì)算了青藏高原東北部及周緣地區(qū)的地殼厚度，其結(jié)果與GRACE衛(wèi)星提供的GX-OG-2-GCM模型所計(jì)算的結(jié)果基本一致(段虎榮等，2010).

由圖5可見，青藏高原東北部地殼厚度在整體上由西向東逐漸減薄，其變化趨勢(shì)以SE向?yàn)橹?，顯示了青藏高原物質(zhì)向東滑移引起地殼縮短變形的趨勢(shì).青藏高原巖石圈物質(zhì)向東運(yùn)動(dòng)，受到揚(yáng)子地臺(tái)和華北地臺(tái)的阻擋，引起物質(zhì)的堆積；然而，在邊緣造山帶并沒(méi)有發(fā)現(xiàn)“山根”的痕跡，重力異常多尺度分解結(jié)果顯示該地殼減薄區(qū)存在明顯的負(fù)異常，體現(xiàn)了深部物質(zhì)上升并與巖石圈底部相互作用導(dǎo)致部分熔融的可能性(圖6).研究區(qū)莫霍面最深處在松潘甘孜腹地，約59 km，而最淺處在四川盆地，約37 km.以102°E—104°E為界，以西地殼厚度在50 km以上，以東在45 km以下，尤其在龍門山附近，存在一條NNE向的地殼陡變帶，其兩側(cè)莫霍面有5 km以上的錯(cuò)斷，龍門山兩側(cè)地殼厚度的顯著差異反映了松潘—甘孜地塊巖石圈與揚(yáng)子地臺(tái)巖石圈結(jié)構(gòu)的差異，這可能是龍門山形成的深部原因.總之，青藏高原東北部的地殼厚度東西向有顯著差異，呈現(xiàn)出地殼內(nèi)部受到了東西方向的強(qiáng)烈擠壓的趨勢(shì)；在印度板塊的擠壓環(huán)境下，南北向的莫霍面起伏卻較為平緩，推測(cè)下地殼和上地幔經(jīng)歷過(guò)強(qiáng)烈的伸展運(yùn)動(dòng).

圖5　研究區(qū)地殼厚度分布圖Fig.5　Variation of the crustal thickness of the study region

圖6　松潘—甘孜地塊東緣深部動(dòng)力學(xué)示意圖據(jù)Wang et al.(2010)修改. 剖面位置見圖5黑色虛線.Fig.6　The cartoon showing the dynamic process under the eastern margin of the Songpan-Garzê block The result was derived from Wang et al.(2010). The location was labeled with black dash line in Fig.5.

圖7　A-A′剖面巖石圈密度模型不同塊體和層位的平均密度已標(biāo)注在模型中，單位為g·cm-3.Fig.7　Lithospheric density structure along the A-A′ profile Densities of different positions are labeled in the model, and the unit of density is g·cm-3.

圖8　B-B′ 剖面巖石圈密度模型不同塊體和層位的平均密度已標(biāo)注在模型中，單位為g·cm-3.Fig.8　Lithospheric density structure along the B-B′ profile Densities of different positions are labeled in the model, and the unit of density is g·cm-3.

5地殼上地幔二維密度結(jié)構(gòu)

從布格重力異常平面等值線上切割兩條剖面(位置見圖3中實(shí)線A-A′、B-B′)，由于地形質(zhì)量產(chǎn)生的短波高頻信息對(duì)深部構(gòu)造研究會(huì)造成嚴(yán)重干擾(楊金玉等，2012)，因此對(duì)原始的布格重力數(shù)據(jù)采用巴特沃斯低通濾波過(guò)濾掉高頻成分，獲得了平穩(wěn)的區(qū)域異常曲線.基于前人的研究結(jié)果作為分層結(jié)構(gòu)(Liu et al.，2006；Wang et al.，2010，2013)，采用二度半體重力異常人機(jī)交互算法正演模擬了兩條剖面的密度結(jié)構(gòu)(圖7，圖8).在構(gòu)建過(guò)程中，對(duì)密度模型的形態(tài)與密度等參數(shù)進(jìn)行多次擬合，力圖達(dá)到最好的效果.

5.1A-A′剖面密度結(jié)構(gòu)

A-A′剖面由99.65°E，33.75°N至106.114°E，30.83°N，全長(zhǎng)697 km，方位為112.7°，該剖面跨越松潘—甘孜地塊中部，以SE向穿越龍門山造山帶，最終抵達(dá)四川盆地腹地.從橫向上看，龍門山兩側(cè)巖石圈物質(zhì)和結(jié)構(gòu)差異明顯(見圖7).上地殼中，松潘—甘孜地塊密度要大于四川盆地.由于早二疊紀(jì)松潘甘孜洋向昆侖柴達(dá)木地塊俯沖，使陸緣弧地帶抬升成山，在三疊紀(jì)時(shí)期松潘甘孜洋上沉積了巨厚海相復(fù)理石堆積(Yin and Harrison，2000)，這種高致密的沉積物甚至被學(xué)者稱為“中生代基底”(嘉世旭等，2009).

在中下地殼，松潘—甘孜地塊平均密度要明顯低于四川盆地，表明處于青藏高原東北部地域的松潘—甘孜地塊中下地殼物質(zhì)相對(duì)柔軟.大地電磁資料顯示青藏高原東部中下地殼范圍內(nèi)存在大規(guī)模的高導(dǎo)體(張樂(lè)天等，2012；Wang et al.，2014)；地震資料也揭示了龍門山造山帶以西中下地殼速度遠(yuǎn)小于平均值(Wang et al.，2010)；地?zé)崮M的松潘—甘孜地塊東部地下溫度，中地殼接近700～800 ℃，下地殼達(dá)到了1000 ℃(Jiménez-Munt et al.，2008)，如此高的溫度足以使地殼物質(zhì)熔融；GPS測(cè)量顯示青藏高原中部和北部地殼物質(zhì)存在著向東移動(dòng)的現(xiàn)象，并且自西向東運(yùn)動(dòng)的速度越來(lái)越小(Zhang et al.，2004).在松潘—甘孜地塊東緣龍門山造山帶，中下地殼的密度達(dá)到橫向上的最低值，且界面起伏較大，顯示該區(qū)受到過(guò)強(qiáng)烈的擠壓致使基底褶皺變形.下地殼流模型似乎可以解釋這些地球物理現(xiàn)象，然而地殼流意味著殼幔解耦，并且有足夠大范圍內(nèi)的下地殼存在流體.布格重力異常多尺度分解結(jié)果(圖4c)顯示在松潘—甘孜地塊內(nèi)部并沒(méi)有出現(xiàn)明顯的低重力異常，巖石圈密度模型(圖7)中松潘—甘孜地塊內(nèi)部下地殼密度為2.80 g·cm-3，密度略低，但并未達(dá)到熔融的程度.

在上地幔結(jié)構(gòu)中，松潘—甘孜地區(qū)的物質(zhì)密度仍小于四川盆地，龍門山造山帶最低.由于青藏高原內(nèi)部存在豐富的地幔熱流活動(dòng)，在龍門山造山帶附近，地幔熱流已充斥到下地殼，導(dǎo)致下地殼物質(zhì)部分熔融，形成塑性流體，因而密度達(dá)到橫向上的最小值.松潘—甘孜地塊向四川盆地仰沖的過(guò)程中，脆性的上地殼相互擠壓碰撞，在接觸部位形成破碎帶，而塑性的巖石圈物質(zhì)在東流過(guò)程中，遇到剛硬的四川盆地巖石圈阻擋，一部分轉(zhuǎn)向東南方，一部分被迫上移，侵入到破碎巖石中并且抬升地殼，因而造就了龍門山的隆升.

5.2B-B′剖面密度結(jié)構(gòu)

B-B′剖面由99.65°E，33.75°N至107.64°E，35.71°N，全長(zhǎng)772 km，方位為67.3°，該剖面跨越松潘—甘孜地塊中部以NE向穿越西秦嶺造山帶，抵達(dá)鄂爾多斯盆地.B-B′剖面密度結(jié)構(gòu)模型呈現(xiàn)出松潘甘孜—西秦嶺—鄂爾多斯地區(qū)的巖石圈密度結(jié)構(gòu)具有橫向分塊，縱向分層的特點(diǎn)(見圖8).上地殼各地塊的平均密度由西南到東北逐漸減?。缓穸壬衔髑貛X造山帶最薄，往兩邊逐漸增厚.

在中下地殼，西秦嶺造山帶與松潘—甘孜地塊密度值相同，反映了中下地殼物質(zhì)成分相同，二者可能統(tǒng)屬一個(gè)基底，這與前人的深地震反射結(jié)果有良好的一致性(高銳等，2006).中下地殼較低的密度值顯示其柔軟的性質(zhì)，也與該區(qū)其他測(cè)量結(jié)果相對(duì)應(yīng)：大地電磁測(cè)深剖面揭示該區(qū)地表下13～20 km到40～50 km之間存在低阻異常體，其電阻率在幾歐姆米到幾十歐姆米(金勝等，2012)；深地震反射剖面結(jié)果顯示該區(qū)域在地下20～45 km深度范圍內(nèi)具有高低速相間、低速為主的結(jié)構(gòu)，與周邊穩(wěn)定地塊(阿拉善地臺(tái)、鄂爾多斯盆地、四川盆地)的中下地殼速度相比，平均低0.3～0.6 km·s-1(嘉世旭和張先康，2008).相比于龍門山，青藏高原東北緣大部分地區(qū)的地殼物質(zhì)并不具備高泊松比的性質(zhì)，單純的地殼流體模型解釋該區(qū)的低密度體有一定的局限，對(duì)此有學(xué)者認(rèn)為青藏高原東北緣中下地殼含有較多的長(zhǎng)英質(zhì)組分(Pan and Niu，2011).長(zhǎng)英質(zhì)組分提升了地殼物質(zhì)的熔點(diǎn)，致使西秦嶺造山帶下地殼物質(zhì)的熔融程度弱于龍門山造山帶.

上地幔中，西秦嶺造山帶密度最小，且莫霍面上隆.結(jié)合前人研究成果推斷西秦嶺地區(qū)上地幔低密度體可能是來(lái)自于地幔流物質(zhì)上涌加熱弱化了巖石圈底部，關(guān)于地幔流物質(zhì)的來(lái)源和成因，下一節(jié)會(huì)詳細(xì)討論.

6地球動(dòng)力學(xué)意義

巖石圈厚度分布揭示青藏高原地區(qū)有兩處明顯減薄區(qū)：羌塘地體以北到昆侖造山帶之間的藏北高原地區(qū)和青藏高原東北緣西秦嶺松潘構(gòu)造結(jié)附近(An and Shi，2006).在藏北高原，存在著大面積低Pn波速和Sn波缺失的事實(shí)(馬開義等，1996)，同時(shí)廣泛發(fā)育新生代堿性火山巖.印度板塊向北俯沖到班公湖—怒江縫合帶附近(潘桂棠等，2004；Kumar et al.，2006)，其巖石圈地幔前緣進(jìn)入亞洲軟流圈深部，造成深部軟流圈地幔的熱物質(zhì)上涌(Owens and Zandt，1997；Tilmann et al.，2003).SKS分裂波顯示該區(qū)具有青藏高原目前所知最強(qiáng)的各向異性區(qū)，其快波極化方向揭示了軟流圈地幔物質(zhì)向東流動(dòng)(楊曉松等，2002).

青藏高原東北緣亦表現(xiàn)出與藏北高原相似的地質(zhì)地球物理性質(zhì).地質(zhì)調(diào)查表明在高原東北部廣泛發(fā)育新生代火山巖(喻學(xué)惠等，2009；湯慶艷等，2012；賴紹聰?shù)龋?014)，有些地區(qū)巖漿活動(dòng)仍在進(jìn)行.地殼下巖漿上升至地表噴發(fā)的前提之一是存在地殼破裂作為其運(yùn)移通道，而破裂源于地殼底部的張應(yīng)力作用(熊熊等，2007).深地震反射結(jié)果揭示了西秦嶺造山帶巖石圈經(jīng)歷了強(qiáng)烈的伸展運(yùn)動(dòng)(高銳等，2006)，說(shuō)明該地區(qū)巖石圈深部確實(shí)存在張應(yīng)力作用.西秦嶺地區(qū)新生代超鉀質(zhì)火山巖的地球化學(xué)特征揭示了其源區(qū)為地?；蜍浟魅?，并且具有地幔柱源的特征(喻學(xué)惠等，2009).火山巖出露位置恰位于上地幔重力低值區(qū)內(nèi)，亦是巖石圈減薄區(qū)，殼幔界面不明顯，同時(shí)周邊有熱泉分布，天水地區(qū)莫霍面附近存在高導(dǎo)低速層預(yù)示著該區(qū)莫霍面以下局部熔融的可能性(林長(zhǎng)佑等，1995).一系列證據(jù)說(shuō)明了該區(qū)巖石圈底部軟流圈熱物質(zhì)通過(guò)底侵作用上涌，強(qiáng)烈的上拱力造成巖石圈的整體形變，地幔熱物質(zhì)進(jìn)入巖石圈致使巖石圈底部的物質(zhì)被加熱弱化，甚至擠入地殼部分，造成下地殼物質(zhì)局部減壓熔融.

地幔物質(zhì)上涌的動(dòng)力學(xué)機(jī)制來(lái)源于青藏高原中心軟流圈地幔流的東向移動(dòng).GPS精準(zhǔn)測(cè)量結(jié)果顯示高原東部物質(zhì)向東北和東南方向運(yùn)動(dòng)(見圖9)；莫宣學(xué)等(2007)根據(jù)青藏高原新生代火山巖的時(shí)空分布關(guān)系并綜合其他地質(zhì)地球物理資料，推斷三江地區(qū)和西秦嶺地區(qū)深部分別存在一條物質(zhì)流通道，而高原腹地下軟流圈地幔物質(zhì)正是沿“軟流圈通道”向東流動(dòng).在西秦嶺松潘構(gòu)造結(jié)，巖石圈上地幔分布著明顯的低密度異常，其重力異常區(qū)呈現(xiàn)側(cè)臥的馬鞍狀形態(tài)，前緣朝NE向凸出，而火山巖出露的位置恰位于前緣位置，體現(xiàn)了地幔熱物質(zhì)并不是垂直上涌，而是存在NE方向的位移趨勢(shì).由于西秦嶺松潘構(gòu)造結(jié)周緣為穩(wěn)定的大陸克拉通地區(qū)，東移的軟流圈地幔流受到揚(yáng)子地臺(tái)和華北地臺(tái)堅(jiān)硬巖石圈的阻擋，在陸續(xù)的擠壓之下，被迫上涌即產(chǎn)生底侵作用，物質(zhì)和能量積累，影響到龍門山和西秦嶺的隆升和近年來(lái)川陜甘地區(qū)地震的誘發(fā).

7結(jié)論

本文將EGM2008模型的重力數(shù)據(jù)應(yīng)用于青藏高原東北部殼幔結(jié)構(gòu)的研究之中.通過(guò)對(duì)多尺度重力場(chǎng)、地殼厚度、巖石圈密度結(jié)構(gòu)的分析以及其他地球物理、地質(zhì)等研究的比較，得到以下幾點(diǎn)認(rèn)識(shí)：

(1) 青藏高原東北部的殼幔物質(zhì)結(jié)構(gòu)與周邊地區(qū)相比有很大的差異：四川盆地、鄂爾多斯盆地以及祁連山地塊整體上表現(xiàn)出均勻的剛性塊體性質(zhì)；而青藏高原東北部則呈現(xiàn)出復(fù)雜的塑性介質(zhì)的特征，其異常的走向多以 EW 或 SSE 為主，反映了巖石圈物質(zhì)向東運(yùn)移的趨勢(shì).

圖9　青藏高原東部地球動(dòng)力學(xué)模式圖GPS數(shù)據(jù)資料取自Liang等(2013)，青藏高原新生代火山巖的時(shí)空分布及地幔流向引自莫宣學(xué)等(2007).Fig.9　The geodynamic model image of the eastern Tibetan plateauGPS data is compiled from Liang et al. (2013), the distribution of the Cenozoic volcanism with time is quoted from Mo et al. (2007).

(2) 重力場(chǎng)多尺度分解結(jié)果并不支持青藏高原東北部下地殼流動(dòng)模型.根據(jù)青藏高原東北部的地殼厚度由西向東逐漸減薄，邊緣造山帶深部無(wú)“山根”痕跡，推測(cè)地幔熱流上涌，抬升了莫霍面，致使巖石圈底部物質(zhì)部分熔融，因而高原東緣、東北緣的巖石圈底部顯示出低密度的性質(zhì).

(3) 龍門山和西秦嶺的隆升可能同時(shí)受大規(guī)模地幔流底侵作用的影響.深部物質(zhì)上涌的動(dòng)力受控于青藏高原向東北擴(kuò)張的過(guò)程中，高原內(nèi)部的軟流圈物質(zhì)側(cè)向流動(dòng)，受到揚(yáng)子地臺(tái)和華北地臺(tái)堅(jiān)硬巖石圈的阻擋而被迫上移.這一動(dòng)力學(xué)過(guò)程或許也是川陜甘地區(qū)地震孕育的重要條件.

本文的研究結(jié)果對(duì)深化認(rèn)識(shí)青藏高原周緣造山帶隆升機(jī)制及深層動(dòng)力學(xué)過(guò)程有著重要意義.如果今后能結(jié)合地面實(shí)測(cè)重力數(shù)據(jù)對(duì)模型進(jìn)行有效校正，則有望為青藏高原的巖石圈物質(zhì)分布、地震孕育等提供更豐富細(xì)致的信息.

致謝感謝兩位審稿專家提出的寶貴意見.文中部分圖件使用了GMT軟件(Wessel and Smith, 1991)進(jìn)行繪制.

References

An M J, Shi Y L. 2006. Lithospheric thickness of the Chinese continent.Phys.EarthPlanet.Int., 159(3-4): 257-266.

Bai D H, Unsworth M J, Meju M A, et al. 2010. Crustal deformation of the eastern Tibetan plateau revealed by magnetotelluric imaging.NatureGeoscience, 3(5): 358-362. Chen S, Xu W M, Shi L, et al. 2013. Gravity field and lithospheric mechanical properties of Longmenshan fault zone and its surrounding areas.ActaSeismologicaSinica(in Chinese), 35(5): 692-703.

Clark M K, Royden L H. 2000. Topographic ooze: Building the eastern margin of Tibet by lower crustal flow.Geology, 28(8): 703-706.

Deng Q D, Zhang P Z, Ran Y K, et al. 2002. Basic feature of China active structures.ScienceinChina(SeriesD) (in Chinese), 32(12): 1020-1030.

Duan H R, Zhang Y Z, Liu F, et al. 2010. Crustal thickness computed by using satellite gravity data in China and its adjacent regions.ProgressinGeophys. (in Chinese), 25(2): 494-499, doi: 10.3969/j.issn.1004-2903.2010.02.016.

England P, Houseman G. 1986. Finite strain calculations of continental deformation: 2. Comparison with the India-Asia collision zone.J.Geophys.Res., 91(B3): 3664-3676.

Feng R H, Dong P, Wang L S, et al. 2014. Inversion of Moho interface in northeastern China with prior information.JournalofEarthScience, 25(1): 146-151, doi: 10.1007/s12583-014-0407-9.

Fu G Y, Zhu Y Q, Gao S H, et al. 2013. Discrepancies between free air gravity anomalies from EGM2008 and the ones from dense gravity/GPS observations at west Sichuan Basin.ChineseJ.Geophys. (in Chinese), 56(11): 3761-3769, doi: 10.6038/cjg20131117.

Gao R, Ma Y S, Li Q S, et al. 2006. Structure of the lower crust beneath the Songpan block and West Qinling orogen and their relation as revealed by deep seismic reflection profiling.GeologicalBulletinofChina(in Chinese), 25(12): 1361-1367.

Gao R, Xiong X S, Li Q S, et al. 2009. The Moho depth of Qinghai-Tibet Plateau revealed by seismic detection.ActaGeoscientiaSinica(in Chinese), 30(6): 761-773.

Gao R, Wang H Y, Zeng L S, et al. 2014. The crust structures and the connection of the Songpan block and West Qinling orogen revealed by the Hezuo-Tangke deep seismic reflection profiling.Tectonophysics, 634: 227-236.

He J K, Lu S J, Wang W M. 2013. Three-dimensional mechanical modeling of the GPS velocity field around the northeastern Tibetan plateau and surrounding regions.Tectonophysics, 584: 257-266.

Hou Z Z, Yang W C. 1997. Wavelet transform and multi-scale analysis on gravity anomalies of China.ChineseJ.Geophys. (ActaGeophysicaSinica) (in Chinese), 40(1): 85-95.

Jia S X, Zhang X K. 2008. Study on the crust phases of deep seismic sounding experiments and fine crust structures in the northeast margin of Tibetan plateau.ChineseJ.Geophys. (in Chinese), 51(5): 1431-1443.

Jia S X, Zhang X K, Zhao J R, et al. 2010. Deep seismic sounding data reveal the crustal structures beneath Zoige basin and its surrounding folded orogenic belts.ScienceChinaEarthSciences, 53(2): 203-212, doi: 10.1007/s11430-009-0166-0.

Jiang W L, Zhang J F, Tian T, et al. 2012. Crustal structure of Chuan-Dian region derived from gravity data and its tectonic implications.Phys.EarthPlanet.Int., 212-213: 76-87.

Jiménez-Munt I, Fernàndez M, Vergés J, et al. 2008. Lithosphere structure underneath the Tibetan Plateau inferred from elevation, gravity and geoid anomalies.EarthPlanet.Sci.Lett., 267(1-2): 276-289.Jin S, Zhang L T, Jin Y J, et al. 2012. Crustal electrical structure along the Hezuo-Dajing profile across the Northeastern Margin of the Tibetan Plateau.ChineseJ.Geophys. (in Chinese), 55(12): 3979-3990, doi: 10.6038/j.issn.0001-5733.2012.12.010.Kumar P, Yuan X H, Kind R, et al. 2006. Imaging the colliding Indian and Asian lithospheric plates beneath Tibet.J.Geophys.Res.-SolidEarth, 111: B06308, doi:10. 1029/2005JB003930.

Lai S C, Qin J F, Zhao S W, et al. 2014. Geochemistry and its continental dynamic implication of the Cenozoic sodic picritic basalt from the Liuping area, northeastern margin of the Tibetan Plateau.ActaPetrologicaSinica(in Chinese), 30(2): 361-370.

Li H Y, Shen Y, Huang Z X, et al. 2014. The distribution of the mid-to-lower crustal low-velocity zone beneath the northeastern Tibetan Plateau revealed from ambient noise tomography.J.Geophys.Res., 119(3): 1954-1970.

Li Y H, Wu Q J, Zhang F X, et al. 2011. Seismic anisotropy of the Northeastern Tibetan Plateau from shear wave splitting analysis.EarthPlanet.Sci.Lett., 304(1-2): 147-157.

Liang S M, Gan W J, Shen C Z, et al. 2013. Three-dimensional velocity field of present-day crustal motion of the Tibetan Plateau derived from GPS measurements.J.Geophys.Res., 118(10): 5722-5732.

Lin C Y, Wu Y X, Yang C F. 1995. On the high conductive layer near the Moho and the transitional zone between the crust and mantle in Tianshui region.ActaSeismologicaSinica(in Chinese), 17(2): 230-236.Liu M J, Mooney W D, Li S L, et al. 2006. Crustal structure of the northeastern margin of the Tibetan plateau from the Songpan-Ganzi terrane to the Ordos basin.Tectonophysics, 420(1-2): 253-266.

Ma K Y, Xue G Q, Dong Y J, et al. 1996. Uplift of Qinghai-Xizang (Tibet) Plateau and mantle deformation due to partial melting.ActaGeoscientiaSinica(in Chinese), 17(2): 214-223. Mallat S G. 1989. A theory for multiresolution signal decomposition: the wavelet representation.IEEETransactionsonPatternAnalysisandMachineIntelligence, 11(7): 674-693.Meng X H, Shi L, Guo L H, et al. 2012. Multi-scale analyses of transverse structures based on gravity anomalies in the northeastern margin of the Tibetan Plateau.ChineseJ.Geophys. (in Chinese), 55(12): 3933-3941, doi: 10.6038/j.issn.0001-5733.2012.12.006.

Mo X X, Zhao Z D, Deng J F, et al. 2007. Migration of the Tibetan Cenozoic potassic volcanism and its transition to eastern basaltic province: implications for crustal and mantle flow.Geoscience(in Chinese), 21(2): 255-264.

Molnar P, Tapponnier P. 1975. Cenozoic tectonics of Asia: effects of a continental collision.Science, 189(4201): 419-426.

Oldenburg D W. 1974. The inversion and interpretation of gravity anomalies.Geophysics, 39(4): 526-536.

Owens T J, Zandt G. 1997. Implications of crustal property variations for models of Tibetan plateau evolution.Nature, 387(6628): 37-43.

Pan G T, Zhu D C, Wang L Q, et al. 2004. Bangong Lake-Nu River suture zone-the northern boundary of Gondwanaland: Evidence from geology and geophysics.EarthScienceFrontiers(in Chinese), 11(4): 371-382.

Pan S Z, Niu F L. 2011. Large contrasts in crustal structure and composition between the Ordos plateau and the NE Tibetan plateau from receiver function analysis.EarthPlanet.Sci.Lett., 303(3-4): 291-298.

Parker R L. 1973. The rapid calculation of potential anomalies.Geophys.J.Int., 31(4): 447-455.

Pavlis N, Kenyon S, Factor J, et al. 2008. Earth gravitational model 2008.∥ SEG Technical Program Expanded Abstracts.LasVegas,Nevada:SEG, 761-763.

Royden L H, Burchfiel B C, King R W, et al. 1997. Surface deformation and lower crustal flow in eastern Tibet.Science, 276(5313): 788-790.

Tang Q Y, Zhang M J, Li X Y, et al. 2012. The chemical and carbon isotopic compositions of volatiles in Cenozoic high-potassic basalts in western Qinling, China and their mantle geodynamic implications.ActaPetrologicaSinica(in Chinese), 28(4): 1251-1260.

Tapponnier P, Xu Z Q, Roger F, et al. 2001. Oblique stepwise rise and growth of the Tibet plateau.Science, 294(5547): 1671-1677.

Tian X B, Zhang Z J. 2013. Bulk crustal properties in NE Tibet and their implications for deformation model.GondwanaRes., 24(2): 548-559.

Tilmann F, Ni J, INDEPTH III Seismic Team. 2003. Seismic imaging of the downwelling Indian lithosphere beneath central Tibet.Science, 300(5624): 1424-1427.

Vergne J, Wittlinger G, Farra V, et al. 2003. Evidence for upper crustal anisotropy in the Songpan-Ganze (northeastern Tibet) terrane.Geophys.Res.Lett., 30(11): 1552, doi:10.1029/2002GL016847.Wang X B, Zhang G, Fang H, et al. 2014. Crust and upper mantle resistivity structure at middle section of Longmenshan, eastern Tibetan plateau.Tectonophysics, 619-620: 143-148.

Wang Y X, Mooney W D, Yuan X C, et al. 2013. Crustal structure of the northeastern Tibetan Plateau from the Southern Tarim Basin to the Sichuan Basin, China.Tectonophysics, 584: 191-208.

Wang Z, Zhao D P, Wang J A. 2010. Deep structure and seismogenesis of the north-south seismic zone in southwest China.J.Geophys.Res., 115(B12): B12334.

Wessel P, Smith W H F. 1991. Free software helps map and display data.Eos.Trans.AGU, 72(41): 441-446.

Xiao Q B, Zhang J, Zhao G Z, et al. 2013. Electrical resistivity structures northeast of the Eastern Kunlun Fault in the Northeastern Tibet: Tectonic implications.Tectonophysics, 601: 125-138.Xiong X, Wang J Y, Teng J W. 2007. Deep mechanical background for the Cenozoic volcanism in the Tibetan plateau.EarthScience(JournalofChinaUniversityofGeosciences) (in Chinese), 32(1): 1-6.

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.

Xue G Q, Qian H, Jiang M, et al. 2003. Studies on the velocity structure of crust-upper mantle beneath northeast Qinghai-Tibet Plateau, China.ActaGeoscientiaSinica(in Chinese), 24(1): 19-26.

Yang J Y, Zhang X H, Zhang F F, et al. 2012. On the accuracy of EGM2008 earth gravitational model in Chinese Mainland.ProgressinGeophys. (in Chinese), 27(4): 1298-1306, doi: 10.6038/j.issn.1004-2903.2012.04.003.

Yang X S, Jin Z M, Ma J, et al. 2002. Genesis of SKS splitting in the north-central Qinghai-Xizang plateau: melt alignment enhanced lithosphere anisotropy.ChineseJ.Geophys. (in Chinese), 45(6): 821-831.

Yin A, Harrison T M. 2000. Geologic evolution of the Himalayan-Tibetan orogen.Annu.Rev.EarthPlant.Sci., 28: 211-280.

Yu X H, Mo X X, Zhao Z D, et al. 2009. Two types of Cenozoic potassic volcanic rocks in West Qinling, Gansu Province: Their petrology, geochemistry and petrogenesis.EarthScienceFrontiers(in Chinese), 16(2): 79-89.

Zhang C Y, Guo C X, Chen J Y, et al. 2009. EGM 2008 and its application analysis in Chinese mainland.ActaGeodaeticaetCartographicaSinica(in Chinese), 38(4): 283-289.

Zhang G W, Guo A L, Yao A P. 2004. Western Qinling-Songpan continental tectonic node in China′s continental tectonics.EarthScienceFrontiers(in Chinese), 11(3): 23-32.

Zhang J S, Gao R, Zeng L S, et al. 2010. Relationship between characteristics of gravity and magnetic anomalies and the earthquakes in the Longmenshan range and adjacent areas.Tectonophysics, 491(1-4): 218-229.

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, Shen Z K, Wang M, et al. 2004. Continuous deformation of the Tibetan Plateau from global positioning system data.Geology, 32(9): 809-812.

Zhang Z J, Klemperer S, Bai Z M, et al. 2011. Crustal structure of the Paleozoic Kunlun orogeny from an active-source seismic profile between Moba and Guide in East Tibet, China.GondwanaRes., 19(4): 994-1007.

Zhang Z J, Bai Z M, Klemperer S L, et al. 2013. Crustal structure across northeastern Tibet from wide-angle seismic profiling: Constraints on the Caledonian Qilian orogeny and its reactivation.Tectonophysics, 606: 140-159.

Zhao G Z, Tang J, Zhan Y, et al. 2005. Relation between electricity structure of the crust and deformation of crustal blocks on the northeastern margin of Qinghai-Tibet plateau.ScienceChinaEarthSciences, 48(10): 1613-1626.

Zhu L P, Helmberger D V. 1998. Moho offset across the northern margin of the Tibetan Plateau.Science, 281(5380): 1170-1172.

附中文參考文獻(xiàn)

陳石, 徐偉民, 石磊等. 2013. 龍門山斷裂帶及其周邊地區(qū)重力場(chǎng)和巖石層力學(xué)特性研究. 地震學(xué)報(bào), 35(5): 692-703.

鄧起東, 張培震, 冉勇康等. 2002. 中國(guó)活動(dòng)構(gòu)造基本特征. 中國(guó)科學(xué)(D輯: 地球科學(xué)), 32(12): 1020-1030.

段虎榮, 張永志, 劉鋒等. 2010. 利用衛(wèi)星重力數(shù)據(jù)研究中國(guó)及鄰域地殼厚度. 地球物理學(xué)進(jìn)展, 25(2): 494-499, doi: 10.3969/j.issn.1004-2903.2010.02.16.

付廣裕, 祝意青, 高尚華等. 2013. 川西地區(qū)實(shí)測(cè)自由空氣重力異常與EGM2008模型結(jié)果的差異. 地球物理學(xué)報(bào), 56(11): 3761-3769, doi: 10.6038/cjg20131117.

高銳, 馬永生, 李秋生等. 2006. 松潘地塊與西秦嶺造山帶下地殼的性質(zhì)和關(guān)系——深地震反射剖面的揭露. 地質(zhì)通報(bào), 25(12): 1361-1367.

高銳, 熊小松, 李秋生等. 2009. 由地震探測(cè)揭示的青藏高原莫霍面深度. 地球?qū)W報(bào), 30(6): 761-773.

侯遵澤, 楊文采. 1997. 中國(guó)重力異常的小波變換與多尺度分析. 地球物理學(xué)報(bào), 40(1): 85-95.

嘉世旭, 張先康. 2008. 青藏高原東北緣深地震測(cè)深震相研究與地殼細(xì)結(jié)構(gòu). 地球物理學(xué)報(bào), 51(5): 1431-1443.

嘉世旭, 張先康, 趙金仁等. 2009. 若爾蓋盆地及周緣褶皺造山帶地殼結(jié)構(gòu)——深地震測(cè)深結(jié)果. 中國(guó)科學(xué)(D輯: 地球科學(xué)), 39(9): 1200-1208.

金勝, 張樂(lè)天, 金永吉等. 2012. 青藏高原東北緣合作—大井剖面地殼電性結(jié)構(gòu)研究. 地球物理學(xué)報(bào), 55(12): 3979-3990, doi: 10.6038/j.issn.0001-5733.2012.12.010.

賴紹聰, 秦江鋒, 趙少偉等. 2014. 青藏高原東北緣柳坪新生代苦橄玄武巖地球化學(xué)及其大陸動(dòng)力學(xué)意義. 巖石學(xué)報(bào), 30(2): 361-370.

林長(zhǎng)佑, 武玉霞, 楊長(zhǎng)福. 1995. 天水地區(qū)莫霍面附近的高導(dǎo)層及殼幔過(guò)渡帶. 地震學(xué)報(bào), 17(2): 230-236.

馬開義, 薛光琦, 董英君等. 1996. 地幔內(nèi)異常熱熔變與青藏高原的隆升. 地球?qū)W報(bào), 17(2): 214-223.

孟小紅, 石磊, 郭良輝等. 2012. 青藏高原東北緣重力異常多尺度橫向構(gòu)造分析. 地球物理學(xué)報(bào), 55(12): 3933-3941, doi: 10.6038/j.issn.0001-5733.2012.12.006.

莫宣學(xué), 趙志丹, 鄧晉福等. 2007. 青藏新生代鉀質(zhì)火山活動(dòng)的時(shí)空遷移及向東部玄武巖省的過(guò)渡: 殼幔深部物質(zhì)流的暗示. 現(xiàn)代地質(zhì), 21(2): 255-264.

潘桂棠, 朱弟成, 王立全等. 2004. 班公湖—怒江縫合帶作為岡瓦納大陸北界的地質(zhì)地球物理證據(jù). 地學(xué)前緣, 11(4): 371-382.

湯慶艷, 張銘杰, 李曉亞等. 2012. 西秦嶺新生代高鉀質(zhì)玄武巖流體組成及其地幔動(dòng)力學(xué)意義. 巖石學(xué)報(bào), 28(4): 1251-1260.

熊熊, 王繼業(yè), 滕吉文. 2007. 青藏高原新生代火山活動(dòng)的深部力學(xué)背景. 地球科學(xué)(中國(guó)地質(zhì)大學(xué)學(xué)報(bào)), 32(1): 1-6.

許志琴, 侯立瑋, 王宗秀等. 1992. 中國(guó)松潘—甘孜造山帶的造山過(guò)程. 北京: 地質(zhì)出版社.

薛光琦, 錢輝, 姜枚等. 2003. 青藏高原東北部天然地震探測(cè)與巖石圈深部特征. 地球?qū)W報(bào), 24(1): 19-26.

楊金玉, 張訓(xùn)華, 張菲菲等. 2012. EGM2008地球重力模型數(shù)據(jù)在中國(guó)大陸地區(qū)的精度分析. 地球物理學(xué)進(jìn)展, (4): 1298-1306, doi: 10.6038/j.issn.1004-2903.2012.04.003

楊曉松, 金振民, 馬瑾等. 2002. 青藏高原北部異常SKS分裂成因的初步探討——被熔體強(qiáng)化的巖石圈各向異性. 地球物理學(xué)報(bào), 45(6): 821-831.

喻學(xué)惠, 莫宣學(xué), 趙志丹等. 2009. 甘肅西秦嶺兩類新生代鉀質(zhì)火山巖: 巖石地球化學(xué)與成因. 地學(xué)前緣, 16(2): 79-89.

章傳銀, 郭春喜, 陳俊勇等. 2009. EGM 2008地球重力場(chǎng)模型在中國(guó)大陸適用性分析. 測(cè)繪學(xué)報(bào), 38(4): 283-289.

張國(guó)偉, 郭安林, 姚安平. 2004. 中國(guó)大陸構(gòu)造中的西秦嶺—松潘大陸構(gòu)造結(jié). 地學(xué)前緣, 11(3): 23-32.

張樂(lè)天, 金勝, 魏文博等. 2012. 青藏高原東緣及四川盆地的殼幔導(dǎo)電性結(jié)構(gòu)研究. 地球物理學(xué)報(bào),55(12): 4126-4137, doi: 10.6038/j.issn.0001-5733.2012.12.025.

趙國(guó)澤, 湯吉, 詹艷等. 2004. 青藏高原東北緣地殼電性結(jié)構(gòu)和地塊變形關(guān)系的研究. 中國(guó)科學(xué)(D輯: 地球科學(xué)), 48(10): 908-918.

(本文編輯何燕)

基金項(xiàng)目國(guó)家自然科學(xué)基金(41274077，41474055)資助.

作者簡(jiǎn)介畢奔騰，男，1990年生，研究生在讀，研究方向?yàn)榈厍蛭锢頂?shù)據(jù)處理與解釋.E-mail:bibenteng@163.com *通訊作者胡祥云，男，1966年生，教授，博士生導(dǎo)師，主要從事地球物理方法理論與應(yīng)用研究.E-mail:xyhu@cug.edu.cn

doi:10.6038/cjg20160213 中圖分類號(hào)P312，P541

收稿日期2015-04-28，2015-07-30收修定稿

Multi-scale analysis to the gravity field of the northeastern Tibetan plateau and its geodynamic implications

BI Ben-Teng1,2，HU Xiang-Yun1,2*，LI Li-Qing3，ZHANG Heng-Lei1,2，LIU Shuang1,2, CAI Jian-Chao1,2

1InstituteofGeophysicsandGeomatics，ChinaUniversityofGeosciences，Wuhan430074，China2HubeiSubsurfaceMulti-scaleImagingKeyLaboratory，Wuhan430074，China3FacultyofEarthResources，ChinaUniversityofGeosciences，Wuhan430074，China

AbstractAffected by both N—S and E—W directed compression，the lithosphere of the northeastern Tibetan plateau, located in the junction region of the youngest plateau and ancient blocks，has experienced active tectonic motions and intense deformation. Such dynamic processes should have left signatures in the gravity field of this region. To understand the gravity field and its relation to deep structure of the northeastern Tibetan plateau，this work made a multi-scale wavelet to the data from the Earth Gravitation Model 2008 (EGM2008) of this region. Our analysis separated the Bouguer gravity anomalies into different parts of shallow to deep and then estimated the corresponding source depths by logarithmic power spectrum technique. Moreover，we calculated the crust thickness and described the Moho topography. The two lithosphere density models respectively across the Longmen Shan orogenic belt and the western Qinling orogenic belt have been constructed, which permit to reveal the distribution features of the different materials in the crust and uppermost mantle beneath the study area.

The results show that the crust-mantle structure in northeastern Tibetan plateau is different from that of the peripheral areas. The Sichuan basin，the Ordos basin and the Qilian block have homogeneous and rigid crust-mantle structure on the whole. While the crust-mantle materials in the Tibetan plateau are heterogeneous and plastic，where gravity anomalies strike dominantly in EW and SSE. It is inferred that the lithosphere materials extrude towards east in the interior of the Tibetan plateau，and then turn to south due to resistance of the Yangtze and the North China platforms in the margin of the plateau. The image of the crust thickness indicates that the northeastern Tibetan plateau gradually thins from west to east. However，there is no trace of “mountain roots” under the peripheral orogenic belts. Instead，the Moho below these areas uplifts. Combined with the low gravity anomalies，this study suggests that the upwelling mantle heat-flow interacts with the substances in upmost mantle and lower crust，uplifting the Moho or forming the new shallow Moho.

We infer that under the influence of the northward motion of the Indian plate，thermal activity may have occurred in the Qiangtang block lithosphere，and then induced thermal mantle flow. This mantle flow shifts toward east along the asthenospheric channel. Due to the obstruction of the rigid lithosphere in the Yangtze platform and the North China platform, the mantle flow is forced to rise up in the western Qinling-Songpan tectonic node. The uplift of the Longmen Shan and the western Qinling as well as the generation mechanism of earthquakes in the western Sichuan and Gansu province may be all related to this dynamic process.

KeywordsNortheastern Tibetan plateau； Gravity field； Lithosphere structure； Density model； EGM2008

畢奔騰， 胡祥云， 李麗清等. 2016. 青藏高原東北部多尺度重力場(chǎng)及其地球動(dòng)力學(xué)意義.地球物理學(xué)報(bào)，59(2):543-555,doi:10.6038/cjg20160213.

Bi B T, Hu X Y, Li L Q, et al. 2016. Multi-scale analysis to the gravity field of the northeastern Tibetan plateau and its geodynamic implications.ChineseJ.Geophys. (in Chinese),59(2):543-555,doi:10.6038/cjg20160213.