青藏高原東西向差異形變與隆升機制

姜效典，李德勇，宮偉，秘叢永

1中國海洋大學海洋地球科學學院，青島 266100

2海底科學與探測技術教育部重點實驗室，青島 266100

1 引言

青藏高原內部構造形變本質上受形成于不同時期的三條邊界構造帶控制，由北至南分別為：柴達木地塊南緣的古生代昆侖斷裂縫合帶、羌塘地塊北緣的三疊紀金沙江斷裂縫合帶和拉薩地塊北緣的班公湖—怒江斷裂縫合帶（圖1），它們的地表形態、傾角及延伸長度、延展深度等基本要素是認識高原構造形變需要搞清的問題.其中之一是，與縫合帶有關的邊界斷層基底在哪兒？換句話說，它們的斷裂深度是多少？由于高昂的費用和二維地震成像精度的限制，它們的構造特征至今尚不清楚.然而，這些邊界斷裂的延展深度是理解青藏高原隆升過程和機制的重要基礎.如果邊界斷層切穿了整個巖石圈或者說上地幔巖石圈正處于被撕裂的準縫合狀態，如Airy地殼均衡理論所說的那樣，巖石圈為與印度和歐亞大陸之間的持續收斂運動相協調將會發生純粹增厚產生垂向補償，同時導致巖石圈強度逐漸減小.Dewey和Burke（1973）的上地殼增厚模型、Tapponnier等（2001）的巖石圈整體增厚模型以及England和Houseman（1989）的整體巖石圈垂向持續應變模型均強調巖石圈強度在青藏高原變形中的作用.與上述構造模式相對應，青藏高原90°E以西部分（圖1）可能與巖石圈整體增厚到一定程度導致的機械強度減小有關，定量估計其有效彈性厚度比較低，約5～20km（Braitenberg et al.，2003；Jiang et al.，2004；Jordan and Watts，2005）.

與西部不同，青藏高原90°E以東部分（圖1）巖石圈強度明顯增加，有效彈性厚度（Te）達30～50km（Braitenberg et al.，2003；Jiang et al.，2004；Jordan and Watts，2005）.即使應用最簡單的地熱模型，青藏高原地殼也達正常厚度的兩倍左右，因此僅用有效彈性厚度值（Te）難以有效描述衍生于地殼與上地幔之間流變熔融圈層的堅硬巖石圈特征（Jin et al.，2008）.由相干譜證明的兩種褶皺模式也均顯示完整的巖石圈并非由單層結構組成（Jin et al.，1994）.因此，關于青藏高原東部巖石圈的更合理解釋是：其機械強度為被黏塑性下地殼所分割的堅硬上地殼和堅硬上地幔兩部分巖石強度的疊加之和.青藏高原東部巖石圈的去耦模型與西部的整體增厚模型截然不同.Zhao和Morgan（1987）曾經提出一種液壓注入模型來解釋青藏高原的變形解耦過程.下地殼的流動性特征也已經被用來解釋上地殼與上地幔變形的不一致性（Royden et al.，1997；Clark and Royden，2000；Royden et al.，2008）.

柴達木盆地同樣位于90°E以東，其本是塔里木盆地原型的一部分，自漸新世開始逐漸由阿爾金斷裂裂離走滑所形成（Meng and Fang，2008）.該盆地并未因逆沖推覆作用發生明顯的巖石圈增厚，但已經達到了屈服應力的臨界點并且正在發生撓曲變形（Meng and Fang，2008）.受塔里木盆地的阻礙，柴達木盆地的壓縮量較小，說明其具有半解耦的巖石圈流變結構（Jin et al.，2008）.

本文應用歐亞大陸中部、中國、蒙古及南亞的高分辨率重力數據（圖2）、DEM高程數據以及深源地震震源深度數據（圖3）對青藏高原的巖石圈結構進行研究，提出一種高原巖石圈分區差異形變模式，其東西向分區界限在90°E附近.在模型中，構造邊界斷裂在青藏高原90°E以西部分往往切穿整個巖石圈，而90°E以東部分在柔軟、黏塑性的下地殼既已終止.差異的邊界斷裂深度說明青藏高原東西兩部分分別具有不同的隆升機制.

2 布格重力數據

本研究中所用的布格重力數據分別來自中國、蒙古、前蘇聯及印度（圖2）.中國境內的原始重力數據由地質礦產部地球物理地球化學勘查研究所1989年主編的中華人民共和國1∶2500000布格重力異常圖數字化而來，其數據累計誤差在中國東部（102°E以東）約為±1mGal，在西部約±2mGal.原始重力測點存在幾個大的未覆蓋空白帶，特別是在西南部的青藏高原和北部邊緣.為使青藏高原地區達到全覆蓋，研究中補充了部分新的重力測點.新增加的重力測點有四個來源，分別為：中石油1994—1999年所做的測量（高瑞祺和趙政璋，2001）；王謙身和安玉林于1998年在青藏高原東部測量所得的349個數據；中國海洋大學在1997年和1998年橫穿青藏高原北部測量所得的468個數據（Jiang et al.，2004）以及在2009年橫穿龍門山褶皺帶測量所得的231個數據.

圖1 青藏高原及其鄰區巖石圈有效彈性厚度（Te）和主要構造邊界圖中Te等值線間距為5km，MFT：前緣逆沖帶；MCT：中央俯沖帶；YZS：雅魯藏布江縫合帶；KRF：喀喇昆侖斷裂；LS：拉薩地塊；BNS：班公湖—怒江縫合帶；QT：羌塘地塊；JRS：金沙江縫合帶；SG：松潘—甘孜地塊；KF：昆侖山斷裂；KS：昆侖縫合帶；WKF：西昆侖斷裂；ATF：阿爾金斷裂；ATT：阿爾金逆沖帶.Fig.1 Effective elastic thickness（Te）of the Tibetan plateau and its vicinity shown by the colored background and scaled in color bars The interval of Te contours is 5km.The abbreviations from south to north are Main Frontal Thrust（MFT），Main Central Thrust（MCT），Yarlung Zangbo Suture（YZS），Karakorum Fault（KRF），Lhasa Block（LS），Bangong—Nujiang Suture（BNS），Qiangtang Block（QT），Jinsha River Suture（JRS），Songpan GarzêBlock（SG），Kunlun Fault（KF），Kunlun Suture（KS），West Kunlun Fault（WKF），Altun Fault（ATF），and Altun Thrust（ATT）.

圖2 青藏高原及鄰區觀測布格重力異常Fig.2 Observed Bouguer gravity anomalies（mGal）in the Tibetan Plateau and its vicinity

圖3 青藏高原及其鄰區地形圖和深源天然地震分布地形高程數據來自DEM，等值線是沉積盆地的沉積物厚度，單位km，間距1km.圓點表示地震震源中心位置，地震記錄來自NEIC（NationalEarthquakeInformationCenter）和USGS（UnitedStatesGeologicalSurvey）1973—2008年觀測到的深度大于80km的深源地震，其中黑色圓點代表震源中心深度在80～100km之間，藍色代表100～200km，紅色代表200～300km.Fig.3 Topographic image of the Tibetau Plateau and its vicinity from DEM Contours are the basement depths in km in the sedimentary basins over Tibet and its vicinity.Solid circles show deep earthquakes acquired by NEIC and USGS from 1973to 2008with the focal depths larger than 80km.Black solid circles have a depth range from 80～100km，blue for100～200km and red for 200～300km

印度、尼泊爾、緬甸及不丹境內的布格重力數據由Balakrishnan（1997，2003）編譯的印度次大陸及毗連區布格重力異常圖數字化得到.前蘇聯及蒙古境內的的重力數據由俄羅斯武裝部隊地形服務處于20世紀90年代測量，該重力數據分辨率達到了10km，精度高于0.5mGal（Kogan and McNutt，1993）.

本文使用的區域布格重力異常數據與基于EGM2008模型的衛星重力數據（BGI，http：／／bgi.omp.obs－mip.fr／data－products／Toolbox／Anomaly－mapsvisualization）存在很大的差異.地表區域重力數據較衛星重力數據更為平滑，衛星重力明顯蘊含著與地形地貌有關的短波信息，而由地表重力數據計算的長波布格重力異常則可以更好地反映了巖石圈深部結構和莫霍面起伏特征.

3 布格重力約束的三維空間域巖石圈撓曲形變模擬

巖石圈撓曲模型可以表征地球巖石圈機械強度的橫向非均勻性變化，進而達到揭示板塊構造活動性質的目的.目前，撓曲模擬的方法可以歸納為三類：協函數技術、小波表示及數值求解撓曲方程.Forsyth（1985）提出的協函數技術對巖石圈上部負載和下部負載具有良好的表征效果（McNutt，1983；Simons，2000），其中有效彈性厚度（Te）的計算依賴于譜窗口的最優估計（Bechtel et al.，1990；Jin et al.，1994；Ojeda and Whitman，2002；Audet and Mareschal，2004）.小波表示方法則假設地球巖石圈撓曲譜可以用長橢球波動方程的本征函數表示（Thomson，1982；Swain and Kirby，2006；Kirby and Swain，2008）.多窗口技術在芬諾斯堪地盾、北歐、南美以及安第斯山邊緣巖石圈強度橫向非均勻性研究中都取得了很好的應用效果（Pérez－Gussinyé et al.，2004，2007，2008）.

本研究應用3D有限差分法在空間域對強度橫向非均勻性的地球巖石圈進行撓曲模擬.這種方法可以有效消除前述的頻率域算法中每一個單元Te需要保持定值的局限，使結果與構造動力過程更加貼近.但由于數值求解難度大，研究和應用實例極少.Van Wees和Cloetingh（1994），Stewart和 Watts（1997）以及Torne等（2000）簡單介紹了該方法，金煜和姜效典（2002）給出了其差分方程的完整表達，Jordan和Watts（2005）在青藏做了實驗.考慮巖石圈橫向非均勻性及上部地形負載的撓曲基本方程表達式為：

其中w（x，y）為3D空間域的巖石圈撓曲形變，N為單位長度上的背景構造應力，g為重力加速度，ht為高程，v為泊松比，ρc為地殼密度，Δρ地幔與地殼的密度差，D為撓曲強度，可由下式計算：

式中E為楊氏模量，h為所模擬巖石圈的有效彈性厚度（Te）.邊界條件由2D撓曲模擬結果確定.模擬中所用到的參數取值為：地殼密度ρc＝2670kg·m－3，地幔密度ρm＝3260kg·m－3，楊氏模量E＝1011Pa，泊松比ν＝0.25，水平負載＝1013N·m－1.由地形引起的上部負載層密度和由盆地沉積及地幔熱物質膨脹上涌導致的下部負載層密度分別賦值為2670kg·m－3、2450kg·m－3和3350kg·m－3（Jiang et al.，2004；Jiang and Jin，2005）.理論布格重力異常由Parker公式（Parker，1973）計算.

地形和盆地的沉積巖被作為上部負載和內部負載加載到撓曲主控方程中.依據不同深度界面的布格異常譜函數不同，我們求解Parker公式的功率譜分布，以確定巖石圈的內部和下部負載的界面，具體模擬工作流程則如圖4所示.

圖5顯示的是研究區布格重力異常功率譜曲線，各密度（負載）界面對應的深度被模擬.8km界面可以解釋為沉積巖底面，19km是康氏面的平均深度，47km是莫霍面的平均深度，220km為巖石圈底界面.與大地電磁測深和廣角地震顯示的在青藏高原200km以下有異常物質面相一致（Kind et al.，2002；Zhang et al.，2010）.反演四個密度界面產生的重力異常（圖6）顯示平均深度47km和220km的密度界面起伏的重力異常構成了實測重力數據（圖2）的90%.因此，我們把導致青藏高原和鄰區巖石圈形變的內部和下部的異常熱物質限定在大于47km的深部.按照上下負載準相關理論（Jin et al.，1994），我們得到了兩個熱物質界面的起伏展布（圖7和圖8），它們被作為內部負載加載到撓曲主控方程參加模擬，Airy均衡做方程的邊界條件.

求解板塊彈性形變模型的主控方程，實測布格重力異常用于約束，得到莫霍面深度的理論計算值（圖9）.當巖石圈處于彈性形變時，莫霍面深度的變化反映了巖石圈撓曲形變的強度和形態.假定印度大陸和歐亞巖石圈的初始彈性厚度分別是80km和40km（圖10），用Parker公式計算理論重力異常，由實測布格重力（圖2）做約束，根據剩余重力異常調整巖石圈彈性強度、進行疊代.剩余布格重力異常最小化會導致上部負載函數的收斂速度減小，用其反演替代下部負載會使模擬迭代過程明顯收斂.當布格異常觀測值和理論模擬值之間的均方根差最小時即為有效彈性厚度（Te）的最佳擬合值.據此，我們得到了青藏高原與鄰區的巖石圈有效彈性厚度分布（圖1）.

明顯的巖石圈機械強度的差異展現在圖中：在經度90°以西的區域，Te大多小于30km，在西昆侖斷裂、喀喇昆侖斷裂和雅魯藏布江縫合帶數值減小到大約5km，部分接近于0，預示斷裂切穿莫霍面進入上地幔，地殼達到局部均衡.而在經度90°以東的區域Te多在40～50km，反映了中等剛度的巖石圈強度，盡管青藏高原的地殼厚度在60～75km，是平均地殼的兩倍，但其有效彈性厚度僅40～50km，只能解釋為塑性柔軟的下地殼使剛度大的上地殼和上地幔解耦，減弱了整個巖石圈的機械強度.反射地震和大地電磁測深得到的地層波速（圖11）支持我們的模擬結果：經度小于90°的高原西南部，縱波速度明顯小于其他地區，當地層深度大于50km后這個特征更加顯著.

圖4 三維空間域巖石圈撓曲形變模擬流程Fig.4 Flow chart of 3Dspatial flexural modeling of lithosphere

圖5 青藏高原及鄰區布格異常功率譜隨深度的變化Fig.5 Power spectra of Bouguer anomalies versus depth over the Tibetan Plateau and its vicinity

4 剩余布格重力異常討論

觀測數據與理論模型數值之間的差異可以幫助我們評估本文模擬結果的有效性.圖12展示了布格異常觀測值（圖2）與理論模擬數值之間的剩余重力異常特征，中國境內的剩余異常約為±15mGal，與實測重力異常值的絕對誤差約為3%，該剩余異常精度印證了我們模型的可靠性，完全滿足本次構造動力學研究的需求.剩余異常的極大和極小值主要分布在東印度板塊和中亞中國大陸西部，這些地區往往缺少重力測點數據的控制，只有對以上地區做一些實地測量以填補其數據空白才可以提高模擬精度.

5 來自撓曲模擬的地球模型及深源地震震源解釋

應用實測布格重力異常數值約束的巖石圈撓曲形變模型顯示了青藏高原及鄰區莫霍面的起伏變化（圖9）.NEIC（National Earthquake Information Center）和USGS（United States Geological Survey）觀測到的位于研究區內的深部地幔地震的震源深度均大于80km，表明地震主要發生在上地幔巖石圈內.為說明青藏高原東西區域隆升機制的差異，本文建立了一個新的覆蓋青藏高原及其鄰區的3D地球巖石圈結構模型.

圖6 不同深度界面的理論布格重力異常（a）研究區平均深度8km密度界面產生的重力異常，界面上下密度差180kg·m－3；（b）平均深度19km密度界面產生的重力異常，界面上下密度差170kg·m－3；（c）平均深度47km密度界面產生的重力異常，界面上下密度差400kg·m－3；（d）平均深度220km密度界面產生的重力異常，界面上下密度差250kg·m－3.Fig.6 Theoretical Bouguer anomalies at different depths（a）Gravity anomalies at 8km depth with density difference of 180kg·m－3 between the upper and lower layers；（b）Gravity anomalies at 19km depth with density difference of 170kg·m－3；（c）Gravity anomalies at 47km depth with density difference of 400kg·m－3；（d）Gravity anomalies at 220km deep depth with density difference of 250kg·m－3.

圖13展示的即為青藏高原及其鄰區的3D地球模型，模型垂向幅度較水平方向放大了15倍.模型頂層為地形；中層為莫霍面，彩色圓點表示了1973年至2008年間發生的、深度介于81～288km的1125次地幔地震的震源中心位置；底部的高傾角層面系震源中心的擬合面，在這個由震源中心展布約束的區帶以下巖層發生了脆裂形變.

我們的模擬顯示，青藏高原內部的昆侖縫合帶、金沙江縫合帶和班公湖—怒江縫合帶處的斷裂活動僅僅局限于地殼內部，沒有深切至地幔.同時，深部地震的分布（圖3）清楚地顯示：絕大多數深源地震分布在青藏高原的西部，少數深度介于80～100km的地震事件發生于雅魯藏布江縫合帶（圖3），在經度大于90°的高原東部沒有地幔地震發生.

6 結論

據高精度重力約束的巖石圈撓曲模擬和深源天然地震觀測結果，得出如下結論：連接青藏高原西構造結的主要邊界斷裂切穿莫霍面進入地幔，嚴重削弱了該區域的巖石圈機械強度，并發生純剪切構造形變；這里的地殼接近局部均衡，厚皮逆沖是地形隆升的主要因素.而與西部不同，青藏高原中東部縫合帶邊界斷層的構造形變僅僅局限于地殼內部，這里的巖石圈由于塑性柔軟下地殼的作用發生了構造解耦；地殼處于區域均衡狀態，下地殼熱物質的流動膨脹是地殼隆升的主控要素.

圖7 加載至巖石圈撓曲形變主控方程模擬的內部負載分布，來自下地殼的融熔膨脹的熱物質Fig.7 Surface loading distribution stemmed from molten and expanding hot materials of the lower crust，which are loaded to the lithosphere flexural modeling

圖8 加載至巖石圈撓曲形變主控方程的內部負載分布，來自上地幔熱物質上涌Fig.8 Surface loading distribution associated with upwelling hot materials of upper mantle，which are loaded to the lithosphere flexural modeling

圖9 實測布格重力約束的青藏高原及鄰區理論莫霍面撓曲形變Fig.9 The oretical flexure of Moho constrained by observed Bouguer gravity anomalies over the Tibetan Plateau and its vicinity

圖10 初始巖石圈有效彈性厚度模型Fig.10 Initial model of lithosphere effective elastic thickness

圖11 青藏高原地殼縱波速度（李松林等，2002；李秋生等，2003；Jiang et al.，2006；劉明軍等，2008；Zhang et al.，2011a，2011b，2013；Teng et al.，2013；Wang et al.，2013）（a）40km深度地殼縱波速度；（b）50km深度地殼縱波速度；（c）60km深度地殼縱波速度；（d）70km深度地殼縱波速度.Fig.11 Compressional velocity of crust over the Tibetan Plateau（Li et al.，2002；Li et al.，2003；Jiang et al.，2006；Liu et al.，2008；Zhang et al.，2011a，2011b，2013；Teng et al.，2013；Wang et al.，2013）（a）Compressional velocity at 40km depth of crust；（b）Compressional velocity at 50km depth of crust；（c）Compressional velocity at 60km depth of crust；（d）Compressional velocity at 70km depth of crust.

圖12 觀測布格異常和3D空間域巖石圈撓曲模擬的理論布格異常間的剩余重力異常Fig.12 Residual gravity anomalies between the observed Bouguer anomalies and the theoretical Bouguer anomalies derived from 3Dspatial flexural modeling of the lithosphere

圖13 青藏高原及其鄰區的3D地球模型（垂向比例幅度較水平方向放大15倍）模型頂層為地形高程，中部為莫霍面，由南至北的高角度傾斜面代表印度巖石圈俯沖，俯沖帶分布受莫霍面以下深源地震震源中心深度（彩色圓點）的控制.1125個地震來自NEIC和USGS記錄的1973年—2008年間震源深度在81～288km的深地震.Fig.13 3DEarth model of the Tibetan plateau and its vicinity（The vertical exaggeration is 15times the horizontal scale）The top horizon is topography and the middle horizon is Moho.The high angle diving horizon from south to north represents the subducting Indian mantle lithosphere，and the diving horizon is constrained by the earthquake focal depths below Moho（balls in color）.There are 1125deep earthquakes from 1973to 2008with focal depths from 81to 288km.The earthquake data are from NEIC and USGS.

Audet P，Mareschal J C.2004.Variations in elastic thickness in the Canadian Shield.EarthandPlanet.Sci.Lett.，226（1－2）：17－31.

Balakrishnan T S.1997.Major tectonic elements of the Indian subcontinent and contiguous areas：ageophysical view.GeologicalSocietyofIndia，38：18.

Balakrishnan T S.2003.Impact of gravity and other geophysical data on the geology of Indian subcontinent.Journalofthe VirtualExplorer，12：83－92.

Bechtel T D，Forsyth D W，Sharpton V L，et al.1990.Variations in effective elastic thickness of the North American lithosphere.Nature，343（6259）：636－638.

Braitenberg C，Wang Y，Fang J，et al.2003.Spatial variations of flexure parameters over the Tibet－Quinghai plateau.Earth Planet.Sci.Lett.，205（3－4）：211－224.

Clark M K，Royden L H.2000.Topographic ooze：Building the eastern margin of Tibet by lower crustal flow.Geology，28（8）：703－706.

Dewey J F，Burke K C.1973.Tibetan，Variscan，and Precambrian basement reactivation：Products of continental collision.TheJournalofGeology，81（6）：683－692.

England P，Houseman G.1989.Extension during continental convergence，with application to the Tibetan Plateau.J.Geophys.Res.：SolidEarth（1978—2012），94（B12）：17561－17579.

Forsyth D W.1985.Subsurface loading and estimates of the flexural rigidity of continental lithosphere.J.Geophys.Res.：Solid Earth（1978—2012），90（B14）：12623－12632.

Gao R Q，Zhao Z Z.2001.The Frontier Petroleum Exploration in China－Petroleum Geology of Qinghai－Tibetan Plateau （in Chinese）.Beijing：Petroleum Industry Press，28－30.

Jiang M，GalvéA，Hirn A，et al.2006.Crustal thickening and variations in architecture from the Qaidam basin to the Qang Tang （North－Central Tibetan Plateau）from wide－angle reflection seismology.Tectonophysics，412（3－4）：121－140.

Jiang X D，Jin Y，McNutt M K.2004.Lithospheric deformation beneath the Altyn Tagh and West Kunlun faults from recent gravity surveys.J.Geophys.Res.：SolidEarth（1978—2012），109（B5）：B05406，doi：10.1029／2003JB002444.

Jiang X D，Jin Y.2005.Mapping the deep lithospheric structure beneath the eastern margin of the Tibetan Plateau from gravity anomalies.J.Geophys.Res.：SolidEarth（1978—2012），110（B7）：B07407，doi：10.1029／2004JB003394.

Jin Y，Jiang X D.2002.Lithosphere Dynamics （in Chinese）.Beijing：Sciences Press，113－119.

Jin Y，McNutt M K，Zhu Y S.1994.Evidence from gravity and topography data for folding of Tibet.Nature，371：669－674.

Jin Y，Wang E，Jiang X D.2008.The dynamic support and decoupling process of the Tibetan lithosphere based on the integration of flexural modeling with other geological and geophysical studies.／／Burchfiel B C，Wang E.Investigations into the Tectonics of the Tibetan Plateau，GSA Special Paper，444：89－104.

Jordan T A，Watts A B.2005.Gravity anomalies，flexure and the elastic thickness structure of the India－Eurasia collisional system.EarthandPlanetaryScienceLetters，236（3－4）：732－750，doi：10.1016／j.epsl.2005.05.036.

Kind R，Yuan X，Saul J，et al.2002.Seismic images of crust and upper mantle beneath Tibet：Evidence for Eurasian plate subduction.Science，298（5596）：1219－1221.

Kirby J F，Swain C J.2008.An accuracy assessment of the fan wavelet coherence method for elastic thickness estimation.Geochemistry，Geophysics，Geosystems，9（3）：Q03022，doi：10.1029／2007GC001773.

Kogan M G，McNutt M K.1993.Gravity field over northern Eurasia and variations in the strength of the upper mantle.Science，259（5094）：473－479.

Li Q S，Peng S P，Gao R，et al.2003.Seismic evidence of the basement uplift in the Bayan Har tectonic belt，Qinghai，and its tectonic significance.GeologicalBulletinofChina（in Chinese），22（10）：782－788.

Li S L，Zhang X K，Zhang C K，et al.2002.A preliminary study on the crustal velocity structure of Maqin－Lanzhou－Jingbian by means of deep seismic sounding profile.ChineseJ.Geophys.（in Chinese），45（2）：210－217.

Liu M J，Li S L，Fang S M，et al.2008.Study on crustal composition and geodynamics using seismic velocities in the northeastern margin of the Tibetan Plateau.ChineseJ.Geophys.（in Chinese），51（2）：412－430.

McNutt M K.1983.Influence of plate subduction on isostatic compensation in northern California.Tectonics，2（4）：399－415.

Meng Q R，Fang X.2008.Cenozoic tectonic development of the Qaidam basin in the northeastern Tibetan Plateau.GSASpecial Paper，444：1－24.

Ojeda G Y，Whitman D.2002.Effect of windowing on lithosphere elastic thickness estimates obtained via the coherence method：results from northern South America.J.Geophys.Res.：Solid Earth（1978—2012），107（B11）：ETG 3－1－ETG 3－12，doi：10.1029／200JB000114.

Parker R L.1973.The rapid calculation of potential anomalies.Geophys，J.R.Astr.Soc.，31（4）：447－455.

Pérez－GussinyéM，Lowry A R，Phipps Morgan J，et al.2008.Effective Elastic thickness variations along the Andean margin and their relationship to subduction geometry.Geochemistry，Geophysics，Geosystems，9（2）：Q02003，doi：10.1029／2007GC001786.

Pérez－GussinyéM，Lowry A R，Watts A B，et al.2004.On the recovery of effective elastic thickness using spectral methods：examples from synthetic data and from the Fennoscandian Shield.J.Geophys.Res.：SolidEarth（1978—2012），109（B10），doi：10.1029／2003JB002788.

Pérez－GussinyéM，Lowry A R，Watts A B.2007.Effective elastic thickness of South America and its implications for intracontinental deformation.Geochemistry，Geophysics，Geosystems，8（5）：Q05009，doi：10.1029／2006GC001511.

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.

Royden L H，Burchfiel B C，van der Hilst R D.2008.The geological evolution of the Tibetan Plateau.Science，321（5892）：1054－1058.

Simons F J，van der Hilst R D，Zuber M T.2000.Spatiospectral localization of isostatic coherence anisotropy in Australia and its relation to seismic anisotropy：Implications for lithospheric deformation.J.Geophys.Res.：SolidEarth（1978—2012），108（B5）：19163－19184.

Stewart J， Watts A B.1997.Gravity anomalies and spatial variations of flexural rigidity at mountain ranges.J.Geophys.Res.：SolidEarth（1978—2012），102（B3）：5327－5352.

Swain C J，Kirby J F.2006.An effective elastic thickness map of Australia from wavelet transforms of gravity and topography using Forsyth′s method.Geophys.Res.Lett.，33（2），doi：10.1029／2005GL025090.

Tapponnier P，Xu Z Q，Roger F，et al.2001.Oblique stepwise rise and growth of the Tibet Plateau.Science，294（5547）：1671－1677.

Teng J W，Zhang Z J，Zhang X K，et al.2013.Investigation of the Moho discontinuity beneath the Chinese mainland using deep seismic sounding profiles.Tectonophysics，608：202－216，doi：10.1016／j.tecto.2012.11.024.

Thomson D J.1982.Spectrum estimation and harmonic analysis.ProceedingsoftheIEEE，70（9）：1055－1096.

Torne M，Fernàndez M，Comas M C，et al.2000.Lithospheric structure beneath the Alboran Basin：results from 3Dgravity modeling and tectonic relevance.J.Geophys.Res.：Solid Earth（1978—2012），105（B2）：3209－3228.

Van Wees J D，Cloetingh S.1994.A finite－difference technique to incorporate spatial variations in rigidity and planar faults into 3－D models for lithospheric flexure.Geophys.J.Int.，117（1）：179－195.

Wang Q S，An Y L.2001.Gravity field and deep structure of Madoi－Shama region in eastern Qinghai－Xizang （Tibetan）plateau.ProgressinGeophysics（in Chinese），16（4）：4－10.

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.

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，doi：10.1016／j.tecto.2013.02.040.

Zhang Z J，Deng Y F，Teng J W，et al.2011a.An overview of the crustal structure of the Tibetan plateau after 35years of deep seismic soundings.JournalofAsianEarthSciences，40（4）：977－989.

Zhang Z J，Klemperer S，Bai Z M，et al.2011b.Crustal structure of the Paleozoic Kunlun orogeny from an active－source seismic profile between Moba and Guide in East Tibet，China.GondwanaResearch，19（4）：994－1007.

Zhang Z J，Yuan X H，Chen Y，et al.2010.Seismic signature of the collision between the east Tibetan escape flow and the Sichuan Basin.EarthandPlanetaryScienceLetters，292（3－4）：254－264.

Zhao W L，Morgan W J.1987.Injection of Indian crust into Tibetan lower crust：a two－dimensional finite element model study.Tectonics，6（4）：489－504.

附中文參考文獻

高瑞祺，趙政璋.2001.中國油氣新區勘探－第六卷－青藏高原石油地質.北京：石油工業出版社，28－30.

金煜，姜效典.2002.巖石圈動力學.北京：科學出版社，113－119.

李秋生，彭蘇萍，高銳等.2003.青藏高原北部巴顏喀拉構造帶基底隆起的地震學證據.地質通報，22（10）：782－788.

李松林，張先康，張成科等.2002.瑪沁—蘭州—靖邊地震測深剖面地殼速度結構的初步研究.地球物理學報，45（2）：210－217.

劉明軍，李松林，方盛明等.2008.利用地震波速研究青藏高原東北緣地殼組成及其動力學.地球物理學報，51（2）：412－430.

王謙身，安玉林.2001.青藏高原東部瑪多—沙馬地區的重力場與深部構造.地球物理學進展，16（4）：4－10.