東北地區至華北北緣地殼結構的區域差異：地殼厚度與波速比的聯合約束

危自根，陳 凌

中國科學院地質與地球物理研究所 巖石圈演化國家重點實驗室，北京 100029

1 引 言

東北地區位于西伯利亞和華北克拉通之間，屬于中亞造山帶東部，主要由松遼盆地、吉黑褶皺帶和興蒙造山帶三個構造單元組成[1］.隨著晚古生代古亞洲洋和晚中生代蒙古洋的相繼閉合，東北地區于白堊紀由多個小陸塊拼合成為一個完整陸塊[2-3］.地質地球化學資料表明東北地區地殼在新生代受到拉伸減薄，其中生代花崗巖和火山巖時代存在著明顯的東西差異[4-5］.與東北地區緊鄰的華北克拉通北緣，其巖石圈在中新生代的克拉通構造活化中受到了顯著的改造與破壞[6-12］.地震學與其它學科資料觀測結果表明，這一地區巖石圈被改造的程度與過程存在著橫向差異[13-17］.

作為華北克拉通的邊界帶和相鄰區域，遼東臺隆、燕山帶和東北地區的巖石圈結構和演化必然與華北克拉通在顯生屆的構造活化和破壞緊密相關.因此，對這些區域的地殼結構研究能為研究華北克拉通破壞的過程與機制提供重要約束.在這方面前人已經做了大量的地球物理觀測研究[1，18-29］，包括：兩條地學斷面[1，22］的地球物理綜合研究、人工源地震勘探[26］、天然地震流動臺陣[15，23］和固定臺站[28-29］資料的研究等.然而，目前天然地震流動臺陣觀測和人工源地震勘探只能約束線性區域的地殼結構，而基于天然地震固定臺站的研究主要集中在華北北緣部分，從而不能從整體上反映包括東北地區和華北北緣在內的地殼結構區域變化.

利用兩條密集線性流動臺陣地震波形數據，危自根等[15］采用接收函數 H-κ疊加方法[30］獲得了華北克拉通北部邊界區域燕山帶、遼東臺隆和興蒙造山帶臺站下方的地殼厚度（H）與平均波速比（Vp／Vs，κ），發現這三個構造區域的地殼結構存在差異.相比于遼東臺隆，燕山帶的地殼在中新生代的華北克拉通活化中可能經歷了更為復雜的后期改造[16］.然而，由于數據空間分布的限制（臺陣對研究區線性采樣），由它們的結果尚不足以形成對這三個構造區域地殼結構的全面認識.

本研究將綜合流動地震臺陣與固定臺網的觀測資料，對比研究東北和華北北緣的地殼結構和性質.通過增加分布在興蒙造山帶、吉黑褶皺帶、燕山帶和遼東臺隆的127個中國地震局固定臺站記錄的波形數據，并結合危自根等[15］對97個流動臺站資料的分析結果，我們共獲得了224個地震臺站（圖1）下方的地殼厚度與平均波速比信息.結合區域地質構造，我們對比分析了上述四個構造區的地殼結構特征及其差異性，并探討了造成這些結構差異的原因.

2 數據與方法

本文研究數據包括兩部分.新增加的127個固定臺站分布在東北和華北北緣的廣大區域，對興蒙造山帶、吉黑褶皺帶、燕山帶及遼東臺隆形成較好的空間覆蓋（圖1）.其數據來自中國地震局地球物理研究所“國家數字測震臺網數據備份中心”，包括從2007年7月到2010年7月的遠震資料[31］.其余97個流動地震臺站組成的兩條線性臺陣為中國科學院地質與地球物理研究所地震臺陣實驗室執行的“華北內部結構計劃”第三期（NCISP III，2003—2004年）和第六期（NCISP VI，2007—2008年）在華北東北邊界及其附近區域的觀測實驗，分別沿N—S和NW—SE向穿越興蒙造山帶、燕山帶和遼東臺隆（圖1）.危自根等[15］采用接收函數 H-κ疊加方法已經獲得了臺陣下方的地殼厚度與波速比分布.本文所選用的遠震事件震中距為30°到90°（圖1右上角插圖），震級都大于5.5級.大多數地震事件分布在背方位角110°到230°之間.

接收函數是由遠震徑向分量和垂直分量反褶積得到的時間序列，由于扣除了震源效應和傳播路徑的影響，因此非常適合用來研究地殼和上地幔間斷面的結構信息.接收函數波形記錄中，緊隨直達P波之后的來源于Moho面的轉換震相主要為Ps轉換波和PpPs及PpSs+PsPs地表多次波.利用這些轉換震相和地表多次波與直達P波的到時關系可以約束Moho間斷面的深度與地殼平均波速比[30］.

構造一個H-κ平面內的疊加函數s（H，κ），

其中，H為地表到Moho面之間的厚度，κ為地殼平均P波和S波速度之比，r（t）為徑向接收函數，tPs，tPpPs，tPpSs+PsPs為不同震相的理論到時，WPs，WPpPs，WPpSs+PsPs為加權系數，且WPs+WPpPs+WPpSs+PsPs=1.最佳的地殼厚度與波速比對應H-κ平面上最大能量值，詳情可參考[16，30］.

電子商務專業是一個大步前進并且有別于傳統的專業，屬于新興學科。在院校教學從事電子商務的教師，很多是非本專業畢業的。很多教師是由類似專業過來的。他們中有的電腦知識儲備不足，有的業務水平和實習培訓經驗不足，給教學帶來了很多問題，學生的專業技能的訓練也不能很好地展開。

上述方法即為Zhu等[30］提出的接收函數的H-κ方法，該方法的主要優點是不用人工挑選震相走時，故而可以避免挑選震相時人為因素的影響.該方法產生的誤差可以用下列公式估計[30］：

其中，σs為疊加函數s（H，κ）的均方差，σH和σk分別為H 和κ的均方差.

本文采用最大熵譜反褶積方法提取接收函數[32-33］，然后在0.03～0.5Hz范圍內濾波（該頻段充分平衡了高頻噪聲干擾和低頻較低分辨率的影響）.在人工挑選接收函數過程中，由于不同臺站數據的質量存在一定差別，我們根據Moho面Ps轉換波和多次反射波的波形特征和識別程度，把臺站數據分為三類（圖2）.如果Ps和PpPs震相清晰并且具有很好的一致性，這些臺站列入“良好”、如臺站LNWFD；如果Ps震相容易識別，但是PpPs波卻受到某種程度干擾，則為“中等”，例如臺站HLYIL；假如Ps和PpPs都受到干擾，則為“較差”，例如臺站LNLHT.通過分類，良好、中等和較差的臺站數目分別為112，87，25（圖1和表1），其中絕大多數數據較差的臺站位于盆地邊緣和構造帶交界處.

經過挑選，接收函數都超過50個的地震臺站被用來做進一步研究.參考前人在研究區域人工源地震勘探[1，22，26］和天然地震 成像結 果[23］，根據區 域地質情況，選取6.35km／s和6.15km／s分別作為海拔高于100m和低于100m的地殼平均P波速度來估計最終的H-κ疊加結果.基于三震相信噪比順序減小的特點，經過大量測試，選取0.6，0.3，0.1分別作為Ps，PpPs和PpSs+PsPs的加權值（不同加權值組合對結果影響非常小）.

3 結 果

本文通過對分布相對均勻的固定臺站下方接收函數的H-κ疊加分析，并結合前期對兩條線性密集流動臺陣的研究結果[15］，獲得了東北地區和華北克拉通北緣地殼厚度與平均波速比的區域分布（圖3，表1）.

如圖3所示，研究區域地殼厚度整體上從松遼盆地（～30km）向四周造山帶逐漸增加.具體表現為：遼東臺隆地殼中間厚四周薄，主要在29～35km變化，平均值～32km；燕山帶地殼厚度由東往西逐漸增加，主要在28～40km變化，平均值～34km；興蒙造山帶地殼東薄西厚，主要在31～39km變化，平均值為～35km；吉黑褶皺帶則表現為地殼南東厚北西薄，主要在28～37km變化，平均值為～33km.相比于東側吉黑褶皺帶和遼東臺隆＜37km厚的地殼，西側興蒙造山帶和燕山帶的最大地殼厚度接近40km.地殼平均波速比在研究區變化顯著，尤其是在松遼盆地和渤海灣盆地周邊以及構造邊界處和郯廬斷裂帶附近變化比較復雜，并且出現局部增高（＞1.87）.在興蒙造山帶和遼東臺隆，波速比變化相對較小，都主要在1.71～1.83范圍內變化，平均值都約為1.77；而在燕山帶和吉黑褶皺帶變化則相對較大，分別在1.70～1.91和1.72～1.89范圍內變化，平均值也相對稍高（～1.79）.

表1 研究區域臺站下方地殼厚度（已去掉地形）和Vp／Vs比（κ）的分布Table 1 Crustal thickness（altitude removed）and Vp／Vsratio（κ）under the stations

續表1

續表1

圖5 研究區域不同構造帶的地殼厚度（a，b，c，d）與地表地形以及Vp／Vs（e，f，g，h）的相關性Fig.5 Altitude versus crustal thickness（a，b，c，d）and Vp／Vsratio versus crustal thickness（e，f，g，h）in different tectonic units of the study region

其中，H代表地殼厚度（去掉地形），A表示海拔，單位為km，相關系數R為0.63.研究區地殼厚度與波速比沒有明顯的相關性（圖4c）.圖5展示了四個構造帶各自的地殼厚度與海拔以及波速比的分布.從圖中可以看出，興蒙造山帶、遼東臺隆和燕山帶三個區域的地殼大致隨著海拔增高而逐漸增厚，其相關系數R分別為0.57，0.60和0.31，而吉黑褶皺帶地殼厚度與海拔則顯示出不相關性質（R=0.13），主要表現為地殼厚度從～28km到～37km的顯著變化，而海拔變化較小（～100-600m）.興蒙造山帶、遼東臺隆和燕山帶的地殼厚度與波速比表現出弱的負相關關系，而吉黑褶皺帶則沒有這種關系.（盡管本文采用的吉黑褶皺帶臺站數量相對較少，但分布比較均勻，且在海拔相對較高和較低區域都有覆蓋（圖1），因此本文的研究結果基本反映了吉黑褶皺帶的地殼結構特征）.

4 討 論

4.1 與前人結果對比

我們通過對比基于不同方法獲得的地殼厚度分布（圖6），對研究區的地殼結構做了進一步分析.對比結果顯示，采用接收函數H-κ疊加方法得到的東北地區地殼厚度變化圖像與前人基于人工源資料獲得的結果具有較好的一致性.Li等[26］采用深地震測深資料獲得了中國大陸的地殼厚度分布（圖6b），發現在本文研究區，燕山帶最西端和興蒙造山帶內部（～45°N）的地殼較厚，可達～40km.我們在這兩個區域同樣得到了約40km的地殼厚度（圖6a）.沿滿州里—綏芬河地學斷面，廣角地震測深資料顯示，該區域普遍存在縱波波速＞7.2km／s的殼幔過渡帶[34］（圖6c中虛線和莫霍面之間）.本文采用接收函數H-κ疊加方法得到的沿斷面的地殼厚度均落在這一殼幔過渡帶深度范圍，其中在興蒙造山帶（滿洲里—大慶）內位于殼幔過渡帶底邊界附近，而在吉黑褶皺帶內（哈爾濱—穆棱）則接近于殼幔過渡帶頂部（圖6c中紅點）.這一特征可能反映了殼幔過渡帶波速隨深度的變化及其相關的殼幔相互作用存在區域差異.此外，本文沿NCISP III和VI兩條流動地震臺陣獲得的地殼厚度與接收函數波形反演[35］和偏移成像[36］得到的結果也基本一致（圖6d，6e）.

采用不同方法得到的區域地殼厚度分布的一致性（圖6）進一步證明了本文以及前人研究結果的可靠性.對于某些臺站，不同方法得到的地殼厚度存在較大的差別，可能來自于以下幾方面的原因：（1）采用不同資料或方法揭示的是地殼結構不同方面的特征.比如，高頻人工源反射波對尖銳的波速和密度跳變（間斷面）結構敏感，而遠震接收函數轉換波則反映相對寬緩的速度隨深度的變化；（2）小尺度的復雜結構變化會增加研究結果的不確定性，增大不同結果之間的差異.例如當存在厚的殼幔轉換帶結構時，由不同方法、不同頻率資料獲得的地殼厚度估計往往存在明顯差別[37-38］.

圖6 采用不同地震學方法獲得的研究區地殼厚度對比（a）接收函數H-κ疊加得到的地殼厚度分布（本文）；（b）由人工地震資料得到的地殼厚度分布[26］；（c）中國滿州里—綏芬河地學斷面下方的地殼P波速度[34］，紅點表示接收函數H-κ疊加結果；（d）NCISP III臺陣下方接收函數H-κ疊加結果[15］（圓圈）和波形反演結果[35］（紅點）的對比.（e）NCISP VI臺陣下方接收函數 H-κ 疊加結果[15］（圓圈）和偏移成像結果[36］（紅點）的對比.Fig.6 Comparison of crustal thickness by different seismological methods（a）Distribution of Hby H-κstacking of RFs in the study region.（b）Distribution of Hby active seismic exploration[26］.（c）P wave velocity under the Manzhouli-Suifenhe Geoscience Transect[34］，red dots represent Hby H-κstacking of RFs.（d）Hon the profile NCISP III from H-κstacking of RFs[15］（circles）and waveform inversion of RFs[35］（red dots）.（e）Hon the profile NCISP VI from H-κstacking of RFs[15］（circles）and receiver function migration[36］（red dots）.

4.2 地殼結構特征與橫向密度變化

本文接收函數H-κ疊加結果表明，研究區域地殼結構特征與區域地質構造密切相關.地殼厚度從松遼盆地向四周造山帶逐漸增厚（圖3），與地表地形大致對應.地殼總體較薄，尤其是在松遼盆地周邊區域薄于30km，具有典型大陸伸展裂陷區的地殼厚度.這與這些區域在中新生代經歷了廣泛的構造伸展并且伸展中心位于松遼盆地的特征是一致的.總體來看，地殼在重力梯度帶西側比東側明顯增厚，與地表地形相對應，這種現象在東北地區尤其明顯，可能與兩側區域顯生屆以來截然不同的淺表地質構造與深部地幔動力學過程有關[21，24，34］.遼東臺隆、燕山帶、吉黑褶皺帶和興蒙造山帶四個構造帶的地殼厚度與波速比呈現出不同的變化特征（圖5），可能反映了不同構造單元在區域動力學背景下具有不同的構造響應.

艾利均衡模型認為山區地殼存在山根，地表起伏引起的重力異常由山根厚度變化來補償，因此地殼隨著地表海拔的增加而逐漸增厚.上文所述（圖4b），研究區地殼大致隨著海拔增高逐漸增厚，相關系數R為0.63，決定系數R2=0.4，暗示地殼厚度中只有40%的信息可以被海拔高程解釋（假設地殼達到均衡狀態），表明地殼厚度與海拔的相關性相對較弱.研究區地殼厚度與波速比也沒有明顯的相關性（圖4c）.這些現象揭示了該區域地殼結構的橫向不均勻性，暗示地殼存在著橫向密度變化.地殼密度的不均勻性得到了重力[4，21］和磁場資料[4］的證明.研究區域大范圍分布的不同時期的花崗巖與火成巖[5-6］，由于各自密度差異可能會產生地殼密度橫向變化.中—新生代的熱構造事件可能不同程度地改變了局部區域地殼厚度與成分，從而造成地殼厚度與海拔之間相對較弱的相關關系和波速比值的復雜變化.

4.3 地殼結構區域差異

本文通過接收函數H-κ疊加方法得到的地殼厚度與波速比的分布特征表明東北地區至華北北緣地殼結構存在明顯的東西差異（圖3，5）.這種差異同樣被其 它 地 球 物 理[4，26，39］和 地 球 化 學[5-6］資 料 所證明.我們研究發現，地殼結構的東、西差異在東北與華北北緣兩個區域表現出截然相反的特征.東北地區東側吉黑褶皺帶相比于西側興蒙造山帶具有較薄的地殼和復雜的波速比，而華北北緣地殼東側遼東臺隆相比于西側燕山帶具有小范圍變化的地殼厚度和相對穩定的波速比.

4.3.1 東北地區東西差異

東北地區主要包括興蒙造山帶、松遼盆地和吉黑褶皺帶（圖1）.由于盆地厚的沉積層多次波對深部結構震相的干擾，本文無法有效約束松遼盆地內部的地殼結構.然而，我們獲得了松遼盆地邊界區域（包括位于華北北緣部分）相對可靠的地殼厚度與波速比.研究結果顯示（圖3，5），這些盆地邊界區域的地殼都比較薄，最薄處只有～28km，這與兩條近SE—NE向穿過東北地區的地學斷面得到的結果一致[1，22］.波速比在邊界處變化復雜，高低值共存，暗示了地殼結構強烈的橫向不均勻性.

研究結果顯示東北地區東側吉黑褶皺帶和西側興蒙造山帶地殼結構存在明顯差異（圖3，5）.整體上興蒙造山帶地殼相對較厚，波速比相對穩定，地殼厚度與海拔相關系數為0.57；而吉黑褶皺帶地殼相對較薄，波速比變化范圍大，地殼厚度與海拔相關系數僅為0.13.東北地區地殼結構的東西差異已被其它地球物理學觀測資料所證明.郝天珧等[18］通過對重力場、磁場、電性結構與深地震測深資料的研究認為中國東部及其鄰域具有東西分帶、南北分塊的特征.地球化學研究發現東北地區花崗巖的形成年代存在東西差異，西側興蒙造山帶的花崗巖主要形成于早白堊紀，而東側吉黑褶皺帶的花崗巖則主要形成于侏羅紀[5-6］.這反映兩個區域中生代以來的地殼構造演化過程有所不同.

地質資料表明[40］，興蒙造山帶主要由相對簡單的額爾古納和興安古陸塊組成；而吉黑褶皺帶則是由多個小古陸塊拼合的塊體（例如松嫩、張廣才嶺、佳木斯、老爺嶺、興凱、布列亞等古陸塊），因而自形成時地殼結構可能更加復雜.侏羅紀以來，東北地區主要受到太平洋板塊向歐亞大陸俯沖的影響[5-6，41］.侏羅紀時期，東北地區遭到了強烈的陸內擠壓造山過程；白堊紀蒙古洋閉合后，該地區作為一個構造整體遭到了強烈的陸內伸展作用，造成巖石圈和地殼的減薄以及廣泛的巖漿作用[2-3，5-6］.與興蒙造山帶相比，吉黑褶皺帶由于更靠近太平洋板塊前緣從而可能受到其更大的影響.這與本文獲得的興蒙造山帶和吉黑褶皺帶下方差異顯著的地殼結構是一致的.興蒙造山帶較厚的地殼和相對穩定的波速比，東部吉黑褶皺帶較薄的地殼和復雜的波速比以及地形與地殼厚度的非鏡像關系，可能代表了這兩部分對太平洋板塊俯沖不同程度的響應.相比于西側興蒙造山帶，東側吉黑褶皺帶可能由于其自身的地殼結構復雜性和更靠近太平洋板塊前緣從而在中新生代遭受到了更為嚴重的減薄與改造.

4.3.2 華北北緣東西差異

危自根等[15］通過對兩條線性密集流動臺陣下方地殼結構的研究認為華北北緣西側燕山帶比東側遼東臺隆可能經歷了更為強烈的后期改造.通過在該區域增加近均勻分布的固定地震臺站（臺站數目接近危自根等[15］三倍），我們得到了更全面的地殼結構信息.研究表明燕山帶和遼東臺隆的地殼厚度與波速比呈現出弱的負相關性并且比危自根等的結果[15］更加明顯（圖5）.Ji等[42］認為這種關系可能是由相對富集長英質（較低波速比）的中上地殼減薄和下地殼鐵鎂質物質（較高波速比）底侵共同引起.研究區域中新生代火山巖和變質核雜巖的出露[14-15，43］證實了這兩種作用的存在.研究結果表明燕山帶和遼東臺隆地殼結構差異非常明顯（圖3，5）.遼東臺隆地殼相對較薄且厚度變化范圍小，波速比相對穩定，地殼厚度與海拔相關系數為0.60；而燕山帶地殼相對較厚且由東向西厚度顯著增加，波速比分布范圍也較大，地殼厚度與海拔相關系數僅為0.31.這種地殼厚度的東西差異和變化趨勢與盧造勛等[1］得到的內蒙古東烏珠穆沁旗至遼寧東溝地學斷面結果一致.

地殼結構差異特征表明，相比于遼東臺隆，燕山帶地殼結構更加復雜，經歷的后期改造可能更加強烈和不均勻.復雜的地殼改造過程可能不同程度地改變了局部區域的地殼厚度與成分，從而造成燕山帶地殼厚度和波速比值的復雜變化.與東北地區單一的NNE構造走向不同（圖1，3），重力梯度帶在華北北緣呈現出E—W和NNE走向共存，表現為該區域重力梯度帶附近地殼結構的復雜變化（圖3）.前人研究表明[14-15］，華北北緣在中生代經歷了與中亞造山帶增生（E—W 走向）和太平洋板塊俯沖（NNE走向）相關的復雜的擠壓和伸展變形并伴隨著強烈的巖漿活動.這兩種構造過程可能造成了華北北緣地殼結構的不均勻改造，相對于燕山帶，東側的遼東臺隆可能主要受太平洋板塊俯沖影響，地殼改造相對較為均勻.

5 結 論

通過對分布相對均勻的127個固定臺站下方接收函數的H-κ疊加分析，并結合危自根等[15］對兩條線性密集流動臺陣的研究結果，我們獲得了東北地區和華北克拉通北緣地殼厚度與平均波速比信息.

（1）研究區域地殼總體較薄，波速比變化復雜，地殼厚度與海拔的相關關系存在區域差異.這些結果表明了地殼結構的強烈橫向不均勻性，暗示著地殼在中—新生代可能經歷了強烈且不均勻的減薄與改造.

（2）東北地區地殼結構存在明顯的東西差異.西側興蒙造山帶地殼較厚且大致隨著海拔增高逐漸增厚，波速比變化相對平穩.東側吉黑褶皺帶地殼較薄且地殼厚度與海拔不成鏡像關系，波速比變化明顯.這種差異可能來自于這兩部分對太平洋板塊俯沖不同程度的響應.相比于西側興蒙造山帶，東側吉黑褶皺帶由于其自身的地殼結構復雜性和更接近太平洋板塊前緣從而在中新生代遭受到了更為嚴重的減薄與改造.

（3）華北北緣地殼結構呈現出明顯的東西向差異.東側遼東臺隆地殼相對較薄且厚度變化范圍小，波速比相對穩定；西側燕山帶地殼相對較厚且由東向西厚度顯著增加，波速比分布范圍較大.這一特征表明了這兩部分對中亞造山帶增生（E—W走向）和太平洋板塊俯沖（NNE走向）構造過程的不同響應.與遼東臺隆相比，燕山帶地殼經歷了更復雜的改造.致 謝 中國科學院地質與地球物理研究所地震臺陣實驗室提供流動地震數據資料；中國地震局地球物理研究所“國家數字測震臺網數據備份中心”提供固定臺站地震波形數據；鄭天愉和王炳瑜分別提供地殼結構的接收函數波形反演和偏移成像結果；評審專家給出了建設性意見和建議，在此一并致謝！

（References）

[1]盧造勛，夏懷寬.內蒙古東烏珠穆沁旗至遼寧東溝地學斷面.北京：地震出版社，1992.Lu Z X，Xia H K.Geoscience Transect From Dong Ujimqinqi，Inner Mongolia（in Chinese）.Beijing：Seismological Press，1992.

[2]Li S T，Mo X X，Yang S G.Evolution of circum Pacific basins and volcanic belts in East China and their geodynamic background.Journal of China University of Geosciences，1995，6（1）：48-58.

[3]張岳橋，趙越，董樹文等.中國東部及鄰區早白堊世裂陷盆地構造演化階段.地學前緣，2004，11（3）：123-133.Zhang Y Q，Zhao Y，Dong S W，et al.Tectonic evolution stages of the Early Cretaceous rift basins in Eastern China and adjacent areas and their geodynamic background.Earth Science Frontiers （in Chinese），2004，11（3）：123-133.

[4]Wu F Y，Sun D Y，Ge W C，et al.Geochronology of the Phanerozoic granitoids in northeastern China.Journal of Asian Earth Sciences，2011，41（1）：1-30.

[5]Zhang J H，Gao S，Ge W C，et al.Geochronology of the Mesozoic volcanic rocks in the Great Xing′an Range，northeastern China：Implications for subduction-induced delamination.Chemical Geology，2010，276（3／4）：144-165.

[6]Menzies M A，Fan W M，Zhang M.Palaeo-zoic and Cenozoic lithoprobes and the loss of＞120km of Archean lithosphere，Sino-Korean craton，China.／／Prichard H M，Alabaster T，Harris N B M eds.Magmatic Processes and Plate Tectonics：Geological Society.London：Special Publication，1993，76：71-81.

[7]Griffin W L，Zhang A，O′Reilly S Y，et al.Phanerozoic evolution of the lithosphere beneath the Sino-Korean Craton，in Flower.／／Martin F J，Chung S L，Lo C H eds.Mantle Dynamics and Plate Interactions in East Asia.American Geophysical Union Geodynamics Series，1998，27：107-126.

[8]Chen L.Concordant structural variations from the surface to the base of the upper mantle in the North China Craton and its tectonic implications.Lithos，2010，120（1-2）：96-115.

[9]Yang J H，Wu F Y，Chung S L，et al.Rapid exhumation and cooling of the Liaonan metamorphic core complex：Inferences from40Ar／39Ar thermochronology and implications for Late Mesozoic extension in the eastern North China Craton.Geol.Soc.Am.Bull.，2007，119（11-12）：1405-1414.

[10]Liu J L，Davis G A，Lin Z Y，et al.The Liaonan metamorphic core complex，Southeastern Liaoning Province，North China：A likely contributor to Cretaceous rotation of Eastern Liaoning，Korea and contiguous areas.Tectonophysics，2005，407（1／2）：65-80.

[11]Chen L，Wang T，Zhao L，et al.Distinct lateral variation of lithospheric thickness in the Northeastern North China Craton.Earth and Planetary Science Letters，2008，267（1／2）：56-68.

[12]陳凌，程騁，危自根.華北克拉通邊界帶區域深部結構的特征差異性及其構造意義.地球科學進展，2010，25（6）：571-581.Chen L，Cheng C，Wei Z G.Contrasting structural features at different boundary areas of the North China Craton and its tectonic implications.Advances in Earth Science（in Chinese），2010，25（6）：571-581.

[13]Davis G A，Zhang Y D，Wang C，et al.Mesozoic tectonic evolution of the Yanshan fold and thrust belt，with emphasis on Hebei and Liaoning Provinces，northern China.GSA Memoir，2001，194：171-197.

[14]鄭亞東，Davis G A，王琮等.燕山帶中生代主要構造事件與板塊構造背景問題.地質學報，2000，74（4）：289-302.Zheng Y D，Davis G A，Wang C，et al.Major Mesozoic tectonic events in the Yanshan belt and the plate tectonic setting.Acta Geologica Sinica （in Chinese），2000，74（4）：289-302.

[15]危自根，陳凌，楊小林.遼東臺隆、燕山帶和興蒙造山帶臺站下方地殼厚度和平均波速比（Vp／Vs）的橫向變化及其構造意義.地球物理學報，2011，54（11）：2799-2808.Wei Z G，Chen L，Yang X L.Transverse variations of crustal thickness and Vp／Vsratio under the stations in the Liaodong anteclise-Yanshan belt-Xingmeng orogenic belt and their tectonic implications.Chinese J.Geophys.（in Chinese），2011，54（11）：2799-2808.

[16]何正勤，葉太蘭，丁志峰.華北東北部的面波相速度層析成像研究.地球物理學報，2009，52（5）：1233-1242.He Z Q，Ye D L，Ding Z F.Surface wave tomography for the phase velocity in the northeastern part of North China.Chinese J.Geophys.（in Chinese），2009，52（5）：1233-1242.

[17]陳石，王謙身，徐偉民等.華北地區熱均衡、重力均衡與深部構造.地球物理學報，2011，54（11）：2864-2875.Chen S，Wang Q S，Xu W M，et al.Thermal isostasy of North China and its gravity isostasy and deep structure.Chinese J.Geophys.（in Chinese），2011，54（11）：2864-2875.

[18]郝天珧，劉伊克，段昶.中國東部及其鄰域地球物理場特征與大地構造意義.地球物理學報，1997，40（5）：677-690.Hao T Y，Liu Y K，Duan C.The characteristic of geophysical field in the east China and adjacent regions and its tectonics significance.Chinese J.Geophys.（in Chinese），1997，40（5）：677-690.

[19]李明明，何玉梅.利用瑞雷面波反演華北克拉通東北部邊界的巖石圈結構.地震學報，2011，33（2）：143-155.Li M M，He Y M.Lithospheric structure beneath northeastern boundary region of the North China Craton from Rayleigh wave dispersion inversion.Acta Seismologica Sinica（in Chinese），2011，33（2）：143-155.

[20]盧造勛，姜德祿，白云等.東北地區地殼上地幔結構的探測與研究.東北地震研究，2005，21（1）：1-8.Lu Z X，Jiang D L，Bai Y，et al.Exploration and research on the structure of the crust and upper mantle in northeast China.Seismological Research of Northeast China（in Chinese），2005，21（1）：1-8.

[21]楊寶俊，劉萬崧，王喜臣等.中國東部大興安嶺重力梯級帶域地球物理場特征及其成因.地球物理學報，2005，48（1）：86-97.Yang B J，Liu W S，Wang X C，et al.Geophysical characteristics of Daxing’anling gravitational gradient zone in the East China and its geodynamic mechanism.Chinese J.Geophys.（in Chinese），2005，48（1）：86-97.

[22]金旭，楊寶俊.中國滿州里—綏芬河地學斷面地球物理場及深部構造特征研究.北京：地震出版社，1994.Jin X，Yang B J.Manzhouli-Suifenhe Geoscience Transect，China（in Chinese）.Beijing：Seismological Press，1994.

[23]Zheng T Y，Chen L，Zhao L，et al.Crustal structure across the Yanshan belt at the northern margin of the North China Craton.Physics of the Earth and Planetary Interiors，2007，161（1-2）：36-49.

[24]Chen L.Lithospheric structure variations between the eastern and central North China Craton from S-and P-receiver function migration.Earth Planet Interiors，2009，173（3-4）：216-227.

[25]Tang Q S，Chen L.Structure of the crust and uppermost mantle of the Yanshan Belt and adjacent regions at the northeastern boundary of the North China Craton from Rayleigh wave dispersion analysis.Tectonophysics，2008，455（1-4）：43-52.

[26]Li S L，Mooney W D，Fan J C.Crustal structure of mainland China from deep seismic sounding data.Tectonophysics，2006，420（1-2）：239-252.

[27]Zhao L，Zheng T Y.Complex upper-mantle deformation beneath the North China Craton：implications for lithospheric thinning.Geophys.J.Int.，2007，170（3）：1095-1099.

[28]許衛衛，鄭天愉.渤海灣盆地北西盆山邊界地區泊松比分布.地球物理學報，2005，48（5）：1077-1084.Xu W W，Zheng T Y.Distribution of Poisson′s ratios in the northwestern basin-mountain boundary of the Bohai Bay Basin.Chinese J.Geophys.（in Chinese），2005，48（5）：1077-1084.

[29]葛粲，鄭勇，熊熊.華北地區地殼厚度與泊松比研究.地球物理學報，2011，54（10）：2538-2548.Ge C，Zheng Y，Xiong X.Study of crustal thickness and Poisson ratio of the North China Craton.Chinese J.Geophys.（in Chinese），2011，54（10）：2538-2548.

[30]Zhu L P，Kanamori H.Moho depth variation in southern California from teleseismic receiver functions.J.Geophys.Res.，2000，105（B2）：2969-2980.

[31]鄭秀芬，歐陽飚，張東寧等."國家數字測震臺網數據備份中心"技術系統建設及其對汶川大地震研究的數據支撐.地球物理學報，2009，52（5）：1412-1417.Zheng X F，Ouyang B，Zhang D N，et al.Technical system construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake.Chinese J.Geophys.（in Chinese），2009，52（5）：1412-1417.

[32]吳慶舉，曾融生.用寬頻帶遠震接收函數研究青藏高原的地殼結構.地球物理學報，1998，41（5）：669-679.Wu Q J，Zeng R S.The crustal structure of Qinghai-Xizang plateau inferred from broadband teleseismic waveform.Chinese J.Geophys.（in Chinese），1998，41（5）：669-679.

[33]Ai Y S，Zheng T Y，Xu W W，et al.A complex 660km discontinuity beneath northeast China.Earth Planet.Sci.Lett.，2003，212（1-2）：63-71.

[34]劉殿秘，韓立國，翁愛華等.松遼盆地西北邊界部分地球物理特征.地球物理學進展，2007，22（6）：1722-1727.Liu D M，Han L G，Weng A H，et al.Partial geophysical features of the northwest boundary in Songliao basin.Progress in Geophysics （in Chinese），2007，22（6）：1722-1727.

[35]Zheng T Y，Chen L，Zhao L，et al.Crust-mantle structure difference across the gravity gradient zone in North China Craton：Seismic image of the thinned continental crust.Physics of the Earth and Planetary Interiors，2006，159（1-2）：43-58.

[36]王炳瑜，陳凌，艾印雙等.華北克拉通東北部及鄰區地殼和地幔轉換帶厚度研究.地球物理學報，2013，56：待刊.Wang B Y，Chen L，Ai Y S，et al.Crustal structure and mantle transition zone thickness beneath the northeastern area of the North China Craton and adjacent region.Chinese J.Geophys.（in Chinese），2013，56：in Press.

[37]Rumpfhuber E M，Keller G R，Sandvol E，et al.Rocky Mountain evolution：Tying Continental Dynamics of the Rocky Mountains and Deep Probe seismic experiments with receiver functions.J.Geophys.Res.，2009，114：B08301，doi：10.1029／2008JB005726.

[38]Wei Z G，Chen L，Xu W W.Crustal thickness and Vp／Vs ratio of the central and western North China Craton and its tectonic implications.Geophys.J.Int.，2011，186（2）：385-389.

[39]劉光鼎.中國海區及鄰域地質-地球物理系列圖 （1∶500萬）.北京：地質出版社，1989.Liu G D.China Seas and Adjacent Regions，Geophysical Series Maps（1∶500 000）（in Chinese）.Beijing：Geological Publishing House，1989

[40]張興洲，楊寶俊，吳福元等.中國興蒙—吉黑地區巖石圈結構基本特征.中國地質，2006，33（4）：816-823.Zhang X Z，Yang B J，Wu F Y，et al.The lithosphere structure in the Hingmong-Jihei（Hinggan-Mongolia-Jilin-Heilongjiang）region，northeastern China.Geology in China（in Chinese），2006，33（4）：816-823.

[41]Ren J Y，Tamaki K，Li S T，et al.Late Mesozoic and Cenozoic rifting and its dynamic setting in Eastern China and adjacent areas.Tectonophysics，2002，344（3-4）：175-205.

[42]Ji S C，Wang Q，Salisbury M H.Composition and tectonic evolution of the Chinese continental crust constrained by Poisson′s ratio.Tectonophysics，2009，463（1-4）：15-30.

[43]Meng Q R.What drove late Mesozoic extension of the northern China-Mongolia tract？Tectonophysics，2003，369（3-4）：155-174.

[44]Zhao G C，Sun M，Wilde S A，et al.Late Archean to Paleoproterozoic evolution of the North China Craton：key issues revisited.Precambrian Research，2005，136（2）：177-202.