地球深內部探測及動力學研究*

傅容珊 黃建華 楊 亭

(蒙城地球物理國家野外科學觀測研究站，中國科學技術大學地球和空間科學學院，合肥 230026)

地球深內部探測及動力學研究*

傅容珊 黃建華 楊 亭

(蒙城地球物理國家野外科學觀測研究站，中國科學技術大學地球和空間科學學院，合肥 230026)

近20年來人類對于地球深內部的探測和認識已經上升到一個新的層次，滋生出許多新的研究熱點。一系列最新發現集中于地球深內部研究領域。簡述了這一領域研究的成果，探討其未來發展和研究的方向。

地球深內部;板塊構造;地幔動力學;地震層析成像;地球動力學

1 引言

自20世紀80年代開始，依賴于全球地震觀測臺網、數值技術和計算技術的飛速發展，地球科學已將人類對于地球深內部的探測和認識推進到了一個新的層次。地幔深部乃至地核的精細結構和橫向不均勻的探測和動力學研究滋生出許多新的熱點。可以毫不夸張地說，近20年來地球科學最新的發現都在地球深內部研究領域。例如，地球固體內核的超速轉動及復雜結構、地幔底部D”層復雜結構和動力學背景、CMB超低速區域以及深地幔過鈣鈦礦物質結構的發現等等。尤為重要的是，地球科學家有了一個強烈的共識:把地球看作一個完整的系統而不是一系列彼此孤立的、互不相關的部分。因此，研究固體地球子系統之間的相互耦合作用，探討地球整體系統演化和現代運轉過程，越來越受到地球科學的關注。與此同時，地球科學家已經將地球物理、地球化學、實際地球狀態下實驗礦物學以及地球動力學的相互交叉和互補作為地球深內部研究的主要途徑。本文將介紹近期地球深內部探測和動力學研究的重大發現和主要成果，同時探討該領域未來發展和研究的方向。

2 地球深內部探測及結構

地球內部結構的探測和研究始于上世紀初期，地震學家利用地震波震相的走時特征發現了地球的分層結構。近百年來，以地震學為主的地球內部結構探測的發展使得人類對于其內部構造的認識建立在兩個層次之上。第一層次的認識是發現地球具有層圈構造(圖1)。第二層次的認識則是探索和發現了地球內部的精細結構和橫向不均勻性。隨著地震觀測精度的提高，地殼、地幔甚至更深內部的精細結構逐漸被地震學家所揭示，例如上地幔中存在220 km、330 km、410 km、520 km以及下地幔中的710 km、900 km和1 050 km全球或者區域性層圈界面等等。特別對地幔底部核幔邊界(CMB)之上D”層(厚度為100～300 km)的發現以及對其精細結構的探測加深了人類對于地球的全新認識，進而全面開展了內部精細結構和橫向不均勻的探測和研究。

圖1 地球內部結構(根據CSEDI，1993編輯)Fig.1 Interior structure of the Earth(Edited from CSEDI，1993)

第二層次的研究起步較晚，始于1970年代。借助于現代計算機技術和數值地震觀測技術的高速發展，地震層析成像技術從全球和區域尺度揭示出地球內部的精細結構和橫向不均勻結構。以不同的方式和精度呈現出全新的地球模型或區域構造模型，極大地豐富和深化了人類對于地球的認識。

地震層析成像結果從三個方面展示出地球內部橫向不均勻結構。首先，全球和區域尺度，地殼、地幔、外核和內核都呈先出不同程度的橫向不均勻;在長波的范圍內P波和S波震相的結果均呈現出很好的相關性;在上地幔中呈現出以板塊構造為基本構架的巖石層和上地幔橫向不均勻，其格局以較為狹窄線性區域分布為基本特征。然而，在下地幔中，從1 600 km左右到核幔邊界橫向不均勻失去了狹窄線性區域分布，多表現為大尺度特點，不過觀測中還是可以尋求到上下地幔之間地震波異常連續性過渡的特征。第二，在固體內核不僅表現出地震波速呈軸對稱的各向異性，對稱軸和地球自轉軸之間的夾角為10°～11°，地震P波傳播速度沿對稱軸方向速度快而在垂直于該軸方向速度慢，同時表現出整個半球區域的差異。這一系列新的重大的發現不僅開拓了人類探測地球深部的嶄新領域，同時也為人類研究其深內部動力學過程提供了觀測基礎。第三，橫向不均勻的另一重要體現是地球內部層圈界面起伏，其包括內外核(ICB)、核幔邊界(CMB)、660km地震波速度不連續面、巖石層和地幔邊界(LMB)以及莫霍界面(MOHO)的起伏。地球內部邊界起伏是在動力學過程的驅動下其偏離流體靜力學平衡的體現，所以對于邊界形態的研究一直是地球動力學研究非常關注課題。早在1960年代Hide等就提出，核幔邊界(CMB)的起伏和由此而引起的核幔邊界的偶合可能是地球自轉的錢德勒晃動和十年尺度波動的原因，也可能是地球磁場反轉的激勵機。與此同時，地球科學家在地震學、重力學和天文等不同學科就采用不同的方式探測和確定CMB的形態。特別是地震層析數據給出了許多不同的模型。而其他邊界如MOHO、LMB、660 km和ICB的形態的探測也得到了一些初步的結果。以地震層析為主導的地球橫向不均勻研究展示了地幔結構的三維圖像和基本框架，而與全球構造相關的俯沖板塊和地幔熱柱在深地幔中的蹤跡，正是地幔熱動力系統演化過程的現代表觀。

核慢邊界區域精細結構的探測，特別是D”層復雜結構、地震波速度各向異性、超低速層和熱化學熱異常區域的發現都是現代地球科學的重大發現。在1950，1960年代D”層被認定為全球范圍內地幔底部地震波速度差異的圈層。但1980年代以來，地震學家們發現該區域的結構遠遠較過去的認識復雜。其不僅有大尺度的不均勻結構，而還有小尺度的橫向不均勻結構，存在地震波傳播速度各向異性。特別值得關注的是，這一層中存在地震S波超低速區域和熱化學橫向不均勻的結構區域。地震學發現，在核慢邊界上一些地域存在地震波超低速區域(ULVZ)。其主要特征是，在CMB上部5～50 km的范圍內，其地震P波速度迅速降低5～10%，而S波速度則降低10～30%。觀察圖2可以發現ULVZ的范圍非常大，其中厚度超過5 km的區域集中在中太平洋、非洲和冰島下方，與大型地幔熱柱密切相關。近來McNamara等［1］總結了地幔底部超低速異常區域的分布，勾畫出了更為詳細的圖像。這一圖像已經遠遠超出了過去地球科學的想象，因此其成因和動力學效應的研究很快成為深地幔研究的重要課題。與此同時，地震學在CMB上方一些區域例如南非下部存在地震S波的低速塊體其速度大約比正常速度低3%左右。區域延伸和橫跨度在1 200 km左右，其形態傾斜，在高度方向伸展1 500 km左右(圖3)。地震波速異常還顯示，其塊體邊緣尖銳，異常邊界區域的寬度約為50 km。基于地幔熱動力學分析，這一塊體被認定為熱-化學異常體。

圖2 全球超低速區ULVZ的分布(藍色區域為缺乏超低速結構，深紅色為超低速區［1］)Fig.2 Global distribution of ULVZ based on seismic studies(Blue areas in the foreground indicate probed areas lacking evidence for ULVZ structure，while red patches in the foreground mark regions with detected ULVZs［1］)

圖3 南部非洲核幔邊界上方熱化學異常體［2］Fig.3 Thermo-chemical anomaly up on the CMB of the south Africa［2］

地球流體外核被地球科學家認定為地球磁場的源區。對于人類，地磁起源一直神秘莫測，而由于數學和物理學的困難，地磁起源的理論研究一直徘徊不前，大量研究僅僅停留在簡單的地磁發電機模型之上的定性討論。然而，這一領域由于其具有挑戰性，吸引了許多優秀的理論地球物理學家。1990年代中期，基于大規模計算機數值模擬的地磁起源的研究取得了突破性的進展。文獻［3］利用當時的超級計算機以磁流體基本方程為基礎，模擬了地球磁場的成長和演化，甚至觀察到了磁極的反轉。無疑，他們的成果使地磁起源的研究看到了未來，在人類揭示地磁場的起源的漫長旅途上邁出了一大步。

地球固體內核的探測和動力學最令人興奮和吃驚。在地球固體內核地震波速各向異性的基礎上，文獻［4，5］利用地震學分析發現地球內核相對于地幔和地殼的快速差異旋轉，其差異轉動的速率在1°/年左右。盡管一些研究對此仍持異議，特別是差異轉動速率的確定。近幾年來，一些學者用相似度極高的地震事件分析了地球內核差異轉動，他們的結果支持地球內核相對于地幔和地殼的快速差異旋轉。固體地球內核超速轉動的發現和確認，無論是對于地球系統的認識，還是探討地球深內部的動力學過程都將揭開新的一頁。基于地震學的探測，有關地球固體內核的結構在過去的幾年不斷有新的發現。例如:固體內核地震波衰減也存在各向異性，其吸收強(弱)的方向和地震波速度高(低)方向一致。利用這一結果和地震波速各向異性對于理解地球固體內核的礦物學以及地球動力學有重要的價值;固體內核的表層，其地震波速度呈現出東西半球的差異;Koper等［6］提供了內核存在小尺度橫向不均勻的地震學證據，其不均勻體的尺度在數十千米范圍;針對固體內核和流體外核之間密度差異存在0.25～1.0 g/cm3低跳躍和0.6～1.8 g/cm3高跳躍模型，Krasnoshchekov等［7］提出了內核表面的鑲嵌構造模型，即從固體內核向流體外核過渡包含了可部分流動的薄層(密度差異低跳躍)，但其鑲嵌著一些尖銳躍變(密度)的區域;Wen［8］利用1993年12月1日和2003年9月6日地震事件的波形分析，發現了該時間段內在中非下部內核表面的半徑有0.98到1.75 km的增大。當然這一變化可用內核差異轉動或該內核表面區域地形的快速增高來解釋。

3 深地幔物質及性質

地球深內部物質成分和巖石礦物組成，地球深內部物質物理特性，一直是地球科學家急切需要認識的。但由于人類還沒有能力獲取原位(真實深度上壓力和溫度)條件下巖石的樣本，更無法測定該條件下的巖石礦物的物理-力學性質等等。所以，在許多情況下只能依賴科學的方法推斷深內部的物質狀態。現代有三種途徑可以獲取地球深內部物質成分和性質的信息:第一，通過地球物理場觀測，特別是地震觀測。因為地震波能夠穿透地球深內部，其不僅能攜帶地球內部地震波速度結構的信息，而且其波形的變化能攜帶大量地球內部物質物理-力學性質的信息;第二，通過高溫高壓條件下地幔巖石的礦物學研究，推斷地球深內部條件下的巖石礦物的組分，相變條件及其物理-力學性質;第三，通過大規模計算機數值模擬的方式，在基本的物理-力學定律的基礎上，探討地球深內部的巖石礦物狀態等等。

在過去一個世紀中，地球物理學家們發展了以PREM為代表的地球結構和物理-力學性質模型，地球化學家們也提供了地球內部巖石礦物組成的基本模型。就地球內部巖石礦物組成而言，廣泛地認識到地球流體外核和固體內核的物質組成是以鐵鎳合金為主體，內核僅僅是在該深度上溫度和壓力下鐵鎳的相變結果。地幔的巖石礦物組成則比較復雜，特別是660 km地震波速度不連續面的性質(化學界面或相變面)長期以來倍受關注。不過大量觀測數據支持該地震波速度不連續面為相變面。就巖石礦物而言，地殼巖石礦物由于容易直接獲取，所以其研究程度相對較高。但地幔就不同了，其研究困難程度遠遠大于地殼物質。不過經過長期探索，地球科學給出了地幔巖石礦物的基本模型(表1)。如表1所顯示，660 km地震波不連續面為一相變面，通過這一界面，下地幔巖石礦物的主要成分為鈣鈦礦結構，大量地球深內部研究都以此為據而開展。不過，在探討地幔底部特別是D”層復雜結構核動力學問題時，其巖石礦物用簡單的鈣鈦礦結構解釋面臨許多困難，有的甚至是難以逾越的障礙。

表1 地幔礦物相模式(根據CSEDI SCIENCE PLAN，1993)Tab.1 Mineral phase model of the mantle(CSEDI SCIENCE PLAN，1993)

圖4 過鈣鈦礦晶體結構在［001］，［100］，和［010］方向上的投影以及層狀結構的三維視像(八面體為圍繞Si原子的O原子，圓球為Mg2+離子。黑線表示一個晶包單元［9］Fig.4 Projections of post-perovskite crystal structure［9］in the direections［001］，［100］and［010］and the 3D view of sandwich

然而，經過了長期的疑惑，地球科學探索戲劇性地發現，在核幔邊界附近的溫度和高壓條件下，鈣鈦礦MgSiO3結構可以轉化為新的高壓相結構，即過鈣鈦礦結構。研究表明，在125 GPa和2 500 K溫度下(相應于核慢邊界附近2 700 km深度)鈣鈦礦Mg-SiO3轉化為SiO6八面體的二維片狀結構，其密度有1.0 to 1.2%的跳躍(圖4)。這一結構與鈣鈦礦Mg-SiO3的三維結構明顯不同。他們還推斷，D”地震波不連續面可能就是這一轉換的結果。而由于這一結構的晶體具有明顯的各向異性，可以導致D”層中地震波強烈的各向異性。與此同時一系列理論和實驗研究都證實了鈣鈦礦結構的存在。

由于這一新的發現，地球深內部特別D”層所表現出來的一些怪異特征或許可以得到合理的解釋。Wookey等［10］利用過鈣鈦礦模型解釋了D”層的形態和位置，解釋了P波和S波呈現出來的差異特征及其在深地幔中不連續性的觀測結果。與此同時，他們還通過計算展示了下地幔中S波各向異性區域性的變化和MgSiO3鈣鈦礦相變之間的一致性等等。而在D”層中的觀測到的特殊結構，一些區域呈現出成對的地震波不連續特征，Hernlund等［11］則用鈣鈦礦的相變邊界來解釋，并賦予其地幔動力學內涵(圖5)。當冷的俯沖板塊進入下地幔，在該深度的溫、壓狀態下由MgSiO3鈣鈦礦結構(Pv)相變為過鈣鈦礦結構(pPv)，并堆積在CMB之上。它們與源于D”層的地幔熱柱的上升運動的物質(鈣鈦礦結構)之間形成了冷-溫暖-熱的橫向熱結構(圖5 (c))，也正是這一下地幔的溫度結構結構形成了成對的地震波不連續特征。對于地幔底部一些地震波速度的超低速區域(UVLZ)也有了新的解釋。Mao等利用實驗室觀測到，在過鈣鈦礦相中包含了40 mol% 的FeSiO3的富鐵過鈣鈦礦具有很低P波和S波速度及高泊松比。富鐵過鈣鈦礦在CMB上的地幔底部的壓力-溫度及化學條件下是穩定的，這一礦物相可以由地幔和地核的反應形成。而地幔動力學則驅動了這些物質在一些區域的積累成為觀測到的地震波的超低速區。顯然這一解釋和過去將ULVZ理解為該區物質的部分熔融明顯不同。

過鈣鈦礦的發現為推進人類理解和認識核幔邊界上部D”區域的實質邁出了很大的一步，但是并不意味所有問題都解決了，因為D”區域復雜的結構和動力學背景還仍然模糊不清。源于該區域的地幔熱柱如何形成?俯沖板塊在該區域如何堆積、演化和重循環到上地幔，甚至到我們這個星球的表層? D”區域中可能存在的大規模化學異常體的來源? D”層是如何和地球外核相互作用的?等等重大的科學問題有待我們去深入探討。

4 地幔對流和地幔熱柱

自1930年代提出地幔對流假想以來，地球動力學家們力圖從地球內部追尋地球表層大規模構造運動的驅動力源。現在地幔對流不再是少數動力學家的假想，包括地幔熱柱在內的地幔對流的深入研究不僅成為研究地幔熱動力系統演化的主線，也成為研究大陸形成和演化驅動機理的主線。可以理解，在深度2 900 km以上的固體地幔的對流是驅動板塊運動和地球表層地質構造、活動，包括大陸漂移、地震、火山和造山運動的主要力源。也可以理解，地幔對流如同熱機一樣在幾十億年的演化過程中不斷將其內部熱量輸送到表層，散發到空間。使得我們這個行星得以緩慢地冷卻而不至于在放射性元素的加熱下過熱熔融。與此同時，地幔對流如同一個巨型攪拌機使原始的不均一的地幔混合、均勻化。可以說，地幔對流(包括地幔熱柱)、板塊構造和仍然神秘的D”層猶如基本框架將固體地球科學的許多分支學科聯系起來。

圖5 深地幔等溫線和鈣鈦礦(Pv)相變的關系((a)三條示意等溫線和過鈣鈦礦(pPv)邊界，(b)與前面相應的Vs速度剖面，(c)可能的下地幔結構，圖中淺色部分為pPv區域，箭頭表示物質流動方向，點劃線表示溫暖，冷和熱地幔物質剖面［11］Fig.5 Relations of the deep mantle isotherms［11］with perovskite phase change

包括地幔熱柱在內的地幔對流是在一個復雜系統中運轉的，我們面對的是比想象要復雜得多的地球系統。由于數學分析能力的局限，大規模數值模擬已成為模擬地幔熱動力系統的主要手段之一。而數值模擬的時空域也從二維發展向三維，從直角坐標面向球域，從穩定對流的描述轉向對流系統演化的探討，特別是從假設模型的研究轉向面對實際地幔。例如，研究地幔分層結構對地幔對流格局的影響;復雜系統對流演化;地幔對流格局和尺度隨熱動力學參數變化的關系;板塊運動和地幔對流的相互作用以及地幔對流和俯沖板片的互動關系等等。

地幔復雜性的表征之一是，其結構的復雜性，物質成分、相變(礦物組成)和流變學特性(黏滯性)等等。通常，物質成分的差異將阻止流體穿越界面，此時對流將分層，形成兩個不相混合的對流系統而形成分層對流的格局。660 km地震波不連續面性質的認定長期以來倍受地球科學家的關注，不過近20年來地震層析成像揭示該界面不能阻止俯沖板片進入深地幔，從而推定了該界面不是化學分層。盡管有的研究提出，在地幔底部CMB以上1 000 km范圍內可能存在化學分層，并以此為據提出了一些對流模型。地幔條件下巖石礦物表明，地幔巖石的流變性強烈依賴溫度的變化及其礦物成分，同時相變過程吸熱和放熱或地幔各分層物質成分變化使得其黏滯性產生差異。一些研究推測，地幔中可能存在兩個低黏滯區(LVZ，或軟流層)，其中第一個LVZ區在100～400 km之間，具有全球分布的特征。第二個LVZ區在660～1 000 km之間，在三個大洋下部都能追蹤到此低速層的存在。與之相應的對流3D數值模擬給出，當地幔Rayleigh數達到2×107時，控制對流模式的基本因素仍然是660km界面，表現為局部或全球性的分層對流，第二LVZ的存在只是促進分層對流而已。地幔熱柱可能在轉換區域產生，在第二LVZ區域出現小尺度的對流格局。值得注意的是，一些研究還通過設定不同的條件來研究俯沖板片和地幔熱柱在660 km界面附近的行為，用以解釋在該界面附近的地震波速不連續面的圖像。

地幔復雜性的表征之二是，由于其物理-力學狀態和物理環境的復雜性，與之相應對流的格局呈現出的復雜和多變特征。觀測數據反映了地幔熱柱在內的地幔熱對流系統的多尺度性的特征。Anderson［12］分析了地幔對流的尺度后指出，1階球諧函數對應的對流可能反映了超級大陸及其分裂，2階呈現的上升和下降流動則對應了地球表層構造的基本格局(冷的地幔部分、大陸和俯沖帶)，6階對流對應了熱點分布。而上地幔小尺度對流的波長在400～1 000 km之間，它受到多種因素的控制，例如，上地幔吸熱相變面深度，以及上地幔低速層厚度等等。3D直角坐標、分層黏滯性的熱動力系統中多種尺度對流的數值模型顯示，當系統Rayleish數從105增加到1.75×107時，上部邊界層呈現出不穩定狀態，小尺度對流的格局也隨Rayleigh數的增加而變化，同時被大尺度對流水平地輸送。分層流變性可能是解釋地幔中同時存在多尺度對流的最好理由，而流變特性隨深度的變化可能在激發巖石層下部小尺度對流中起關鍵作用。當然這一實驗可以定性地說明同一熱動力系統中多尺度對流同時存在的可能性，以及存在的可能條件，不過這離實際地幔的討論還相差甚遠。可以預期，諸如球型幾何特征，如地幔可壓縮性、固體相變、地幔流變性隨溫度和壓力變化、地幔熱動力參數隨深度變化以及部分熔融等都可能是影響對流格局和演化的因素。

地幔復雜性的表征之三是，現代地幔對流系統和以板塊構造為主導的巖石層和在CMB之上的D”層之間既統一又相互獨立。地球是我們觀測到的太陽系中行星表層構造中唯一顯示為板塊構造的星球。地幔對流、巖石層板塊運動和D”層內的物質運移可以被看為統一的動力學系統，因為海洋巖石層板塊被證明可能是對流系統上部冷卻的熱邊界層，而D”也常常被想象為對流系統的下部邊界層。然而，由于他們之間流變性質的巨大差異又使得其各自表現出相對的獨立性。所以可以理解它們為相對獨立、相互偶合的子系統。地球動力學家在模擬這一既統一又相對獨立的復雜系統時，面臨許多困惑，采用了一些不同的流變模型。實驗室研究表明，巖石黏滯性和溫度之間存在如下的關系，

式中，E是與物質原子結構相關的常數，T為溫度，而T0和η0則為選取的參考溫度和相應的黏滯性。例如，將地幔和巖石層分成流變性差異很大的兩層，或者將巖石層視為伴隨弱化帶的高黏滯性區域。一些研究則從更為復雜的非牛頓流變性出發，整體地模擬巖石層板塊的形成。早期模型顯示了一個能夠自適應生成類巖石層板塊結構的地幔對流系統，不過其板塊是固定不動的，通常稱此類模型為自適應板塊-地幔模型。近年來，發展了一種考慮形變和軟弱性的機制，廣泛地應用于自適應板塊-地幔模型，其依賴于時間的三維對流數值模擬獲得了類板塊構造特征。例如，在模擬中僅僅考慮巖石的屈服強度而不顧及其應變的弱化。文獻［13］給出的三維自適應板塊-地幔模型的數值模擬，其模型采用了三維直角坐標系，對應于11 600×11 600×2 900 (km)的實際地幔幾何框架。其模型從底部加熱，采用溫度-應變率依賴關系的流變性，同時結合應力屈服點的概念。模擬結果顯示，對流系統形成了一個高強度的頂層(類似板塊)，當應力足夠高時該層將破裂。該模型還產生一些類板塊構造特征的現象，如高強度層的俯沖和在流體表面類似板塊運動等等。然而，這一模擬也出現了與地球表層構造運動不相似的圖像。例如，伴隨新的高強度頂層形成而偶爾產生的消減流，這一特征更像在金星上對流的格局。Tackley［14］的數值模型采用了更徹底的參數化空間，包括統一的黏滯度-溫度依賴關系(變化可達5個數量級)，統一的屈服強度以及溫度達到固化點時其十倍變化的黏滯性等等。模擬顯示出對流系統的行為很大地依賴于屈服強度，低屈服強度時可以得到類板塊構造特征，但是板塊很弱;中等屈服強度時可以得到強的板塊和弱的邊界;而高屈服強度時則會形成一個剛性的蓋層。數值模擬得到了更加連續和平滑演化的類板塊構造特征，同時在板塊下部存在一個低黏滯性區域，呈現出可滑動的擴張中心，其分布更加區域化，減少了隨機性等。然而，模擬未能給出一個純粹的轉換邊界，而每個邊界都是發散和轉換運動的混合等等。對此，正如 Tackley［14］指出的“盡管取得了很大的進展，但是這一類板塊特征并不像地球的板塊構造”。所以，真正地實現模擬地球板塊構造還有很多工作要做，這也是今后全球構造-地幔對流數值模擬的一個重要的方向性課題。可喜的是，有關此類研究的大規模數值模擬已經開展，有的研究探討了深地幔過程和超級大陸之間的耦合關系，有的研究探討了地幔底部熱化學異常體與熱點和大規模火山區域之間的關系。不過，特別需要思考的是，以黏滯性為代表的巖石層地幔和D”層的力學性質不僅依賴作用力的時間和空間物理環境，還依賴其巖石礦物組成。與此同時，在同樣的溫度壓力下巖石礦物部分熔融將改變其力學特性，影響包括巖石層在內的地幔物質運移的方式和格局，特別是在探討巖石層和下伏地幔之間的耦合時這一因素可能起到關鍵的控制性作用。

地幔復雜性的表征之四是，在地幔礦物分層結構的基本成分的構架基礎上，地球化學家發現了地幔物質橫向不均勻分布的復雜特征。數據顯示，源于地幔巖漿巖中不相容示蹤元素分成兩大類:一類來自大洋中脊玄武巖(MORB)，呈現了部分元素的虧損，被解釋為上地幔巖石部分熔融的結果，稱之為虧損或耗散了的地幔(DMM)。而另一類則是洋島玄武巖(OIB)，如夏威夷熱點火山巖，其被解釋為來自更深地幔的地幔熱柱，具有原始地幔的地球化學表征。早期地球化學研究將上地幔認定為耗散地幔源，而下地幔則為原始地幔源。然而，隨著地球化學數據的不斷增加和積累，這一簡單的劃分無法解釋新的數據，一些地球化學家更趨向于MORB和OIB來自不同地幔化學儲體的推測。OIB除了可能來自深地幔D”層的原始地幔外，有的還可能來自別的不同類型地球化學儲體的物質和MORB對應的耗散地幔(DMM)混合的結果。例如，一些年齡大約為10～20億年的地幔巖樣品，它們不是原始的地幔，可能是大洋和大陸地殼物質再循環進入地幔形成的化學儲體。有的地幔巖樣品則認為其源于高3He/4He值的地幔儲體。由于3He被認為與地球同時形成的，而4He是U、Th放射性衰變的結果，所以高3He/4He可以理解為原始地球的物質儲體。如此，地幔對流的問題就顯得更加復雜，更加撲朔迷離。Takley［14］列舉了一些依賴于地球化學數據的框架模型:

1)典型的上下地幔(660 km)化學分層模型;

2)典型全地幔對流模型，除了地幔底部存在富集重循環地殼(ERC)和原始地幔外，其余為均勻地幔;

3)原始地幔泡狀塊體伴隨富集重循環地殼(ERC)模型;

4)完全重循環模型;

5)原始殘留塊體模型;

6)深原始地幔層模型。

其中 Kelloge等［15］的模型假設，從深度大約1 600 km可能存在一個極不規則的化學界面(圖6)，其頂部深度可能為1 600 km，由下層中上升的流動將其推動，并在上面層產生向上的可能是熱柱的流動。而上層中的下降流(俯沖板片)迫使該界面向下移動至核幔邊界(CMB)附近，兩層之間不存在物質交換。他們的數值模擬圖形也顯示了存在此種對流模型的可能性。不過，現今的地震觀測的精度還不能夠探測在地幔深部是否存在橫向如此大起伏的全球性的物質(化學)界面。Lava燈模型是在實驗室模擬和地球物理及地球化學數據的基礎上，以地幔熱柱形成機為中心建立新的地幔熱動力學框架。基于Davaille的實驗認識到，在地幔中即使很細小的密度分層(1%左右)都將產生很重要的動力學結果。實驗模擬了底部加熱頂部冷卻的密度(化學)分層的黏滯流體的對流特征。結果表明，當流體中層間的密度(化學)差異和由于熱膨脹產生的密度異常之間的比率在0.3到0.5之間時，其密度分界面變形，成為很寬闊的圓頂以至上升到層的表面，而冷的物質又會下沉，這一過程猶如Lava燈一樣。Davaille［16］假設，地球開始可能是分層對流，那時其浮力的比率大于1，后來有上升的熱柱和片狀下降的流動使兩層的地幔逐漸混合，其之間浮力的比率下降到0.5以下。該模型可能有助于我們調和諸多觀測之間的矛盾，特別是地球物理(地震層析成像)和地球化學數據的矛盾。不過，這一模型不僅要三維球域的實驗或數值模擬研究，而且還需要更加接近于地幔實際狀態的力學框架和熱動力學參數，更需要地震層析等直接觀測數據的證明。看起來問題又回到了原點，地幔的自然屬性是化學雙層地幔還是化學統一的單層地幔?不同的只是化學分層從660 km地震不連續面轉向1 600 km深度，而后者在目前技術水平上，難以用地震學方法探測。

圖6 下地幔固有化學分層動力學模型［15］Fig.6 Inherent chemical stratification dynamics model of the lower mantle［15］

正如本節前面談及，地幔對流如同一個巨型攪拌機使原始的不均一的地幔混合、均勻化。但是，一些模型，特別是化學分層地幔模型顯示，地幔對流不斷地將深地幔化學差異的物質帶到地球表層，如此在對流循環過程中會增加地幔不均勻性。那么地幔均勻嗎?早期的研究表明，對于黏滯分層地幔而言，如果下地幔黏滯性100倍于上地幔時，全地幔對流仍然能使得下地幔保持地幔原始狀態。直角坐標系中地幔混合的3D數值模擬表明，在環形運動存在的條件下，即使是穩定對流也能夠產生有效、無序的混合。而近期計算更表明，在地球存在如此活躍的對流時，盡管存在一系列不利因素，如地幔具有黏滯性-溫度依賴關系，在660 km界面的下地幔有很高的黏滯性，而且660 km界面為吸熱相變面等等，在幾十億年的時間內下地幔仍然能很好地混合和去氣(outgassed)。文獻［17］首次研究了3D球域中現代地幔對流條件下地幔混合的效率問題。他們的數值模擬實驗中的示蹤元素顯示，在現代地幔對流(包括板塊運動)的驅動下，存在多尺度的地幔混合。只有一些確定的地域呈現出單一對流細胞的片狀混合，而多數區域表現出螺旋狀混合。如此，對流可以將示蹤元素從一個地點輸送到很遠的地方而呈現出無序混合狀態。他們得出，主導現代地幔混合效率的關鍵因素是，現代俯沖板片(特別是環太平洋區域)的驅動力和現代巖石層板塊運動環型速度場;相對而言，現代地球已經有效地混合了，而經過了長期演化的地幔中不可能保留大的孤立塊體。傅等［18］定義了地幔快體置換度和對流混合濃度，他們的研究結果同樣顯示，即使在穩態對流的狀態下，經過40億年演化后的地幔基本均勻，其地幔塊體置換度超過80%。然而，一些研究仍然指出，地幔混合的效率還會由于地幔黏滯性橫向變化或者原始地幔中高黏滯性(高于正常地幔10到100倍)殘余泡狀塊體所降低，這些泡狀塊體將保留幾十億年而保留下來不被混合。可以理解，地幔中的化學異常體殘存時間不僅取決于異常體的黏滯性，還取決于其和周圍地幔的浮力比等等，所以這一問題仍然需要更多的關注和研究。

圖7 地幔系統示意和地幔最底部邊界層［19］(a)地幔剖面及關結構單元:包括板塊構造近地表的熱化學邊界層和由地震學提供的下地幔大尺度圖像，(b)熱邊界層模型(TBL)，(c)化學邊界層模型(CBL)，(d)集中在上升流下部的化學殘渣模型，(e)俯沖板塊堆積模型，(f)均勻下地幔相變模型(PC)，(g)b-f的混雜模型)Fig.7 Sketch of the mantle system and its bottom boundary layer model［19］

5 俯沖板塊和地幔熱柱的深地幔追蹤

地球科學已經確認，板塊構造主導了地球表層構造運動。來自地幔的物質在大洋中脊上升，形成新的海洋巖石層，巨大的巖石層塊體(板塊)在不斷背向洋中脊運動的過程中不斷地冷卻、增厚，而在海溝區域俯沖，消減進入地幔。俯沖板塊去向何方? Walson循環的關鍵環節迷失了。然而，當我們為雄偉壯觀的火山噴發而贊嘆大自然的巨大力量時，自然會問到這些物質來自何處?它們屬于Walson循環環節嗎?自Morgan［20］提出地幔熱柱的概念以及它在板塊動力學中的作用以來，地幔熱柱就作為地幔演化動力學一個重要的子系統而得到廣泛的關注。正如本文前面已經說明，卓有成效的地球深內部探測揭示，地幔底部存在一個D”層，其結構和成分的復雜性并不亞于巖石層。或許它是俯沖板塊的歸宿?或許它同時又是地幔熱柱的發源地?是俯沖板塊和地幔熱柱將相距2 900 km的巖石層和D”層連接起來和地幔構成了統一的熱動力系統?對此，Lay等給出了一個簡單的框架(圖8)和推測D”層的基本模式，并討論了這些模式。對于D”層而言，他們詳細分析和推薦了熱-化學邊界層(TCBL)模型。不過，這一框架中既沒有顧及相變的因素(Mg-SiO3過鈣鈦礦結構的發現而改變了地球科學對D”層的基本認識)，同時也缺乏將活躍的巖石層和復雜的D”層聯系起來的元素。事實上近十來年在深地幔探測的研究中已經確立起這一聯系巖石層和D”層的主要元素，那就是來自地球表層的巖石層俯沖板塊和源于D”層的地幔熱柱。

地球是現代仍然活動的演化的行星。現代所獲取的地球觀測數據都是這個星球演化過程的現代表觀。地震層析成像數據在許多地區明顯地呈現出現代俯沖板片可能進入下地幔之中，有的可達1 300～2 000 km的深度，甚至堆砌在核慢邊界之上，解決了長期以來困擾地學界的有關660 km地震波不連續面的性質問題。Bijwaard等［21］列出了全球各地區俯沖帶和包括了印度板塊和歐亞大陸碰撞帶在內的19個地震層析成像剖面。他們發現，在愛情海、阿留申弧、堪察加、馬尼亞納海溝等一系列的剖面上，地震波高速異常都穿過了660 km間斷面，有的延續到下地幔深部，甚至2 000 km的深度上。只有少數剖面顯示高速異常或中斷于660 km間斷面或沿該界面延伸，如Izu-Bonin島弧和Scotta Sea俯沖帶等地區。不過Zhao等對全球特別是東亞-西太平洋地區俯沖帶的地震層析成像的研究顯示出一些不同的特征，即較強的、擴散開的地震波高速異常停留在俯沖帶的相變區域(660 km)，而由于相變帶的強烈的重力不穩定性，其最終沉入下地幔。這一差異可能改變以對流為中心的地幔動力過程的一些故事情節，但并沒有改變下地幔乃至D”層可能是俯沖板塊的最后歸宿的結論。

地幔動力學的主要任務之一是將研究的視野由近向遠地推向更遙遠的地質時代，從現代觀測的地震層析成像數據中去尋找殘存在地幔中的古老的俯沖板片的痕跡，推求古老的地質塊體漂移、拼合的歷史，以追尋地球演化的軌跡。Van der Voo等［22，23］研究了西太平洋俯沖削減帶，包括 Kurile島，日本及南部地域以及古亞洲特提斯(Tethys)邊緣區域的地震層析成像結果，聯系晚侏蘿紀135～155 Ma及早白堊紀95～135 Ma以來西伯利亞，蒙古等地質塊體的運動拼合以及Mongol-Okhotsk和Kula-Nero洋的閉合史，發現了可能殘存于地幔中的中生代俯沖板片的蹤跡。他們標示出與現代活動的西太平洋俯沖板片相聯系的地震波高速異常帶以及與其近鄰的深入到約1 400 km以下的Mongol-Okhotsk洋俯沖巖石層板片殘留體。他們推斷，當時介于西伯利亞(Siberia)和連接蒙古-華北(Mongolia-North China)塊體及東北西伯利亞(Omolon)塊體之間的Mongol-Okhotsk和Kula-Nero洋閉合了。確定該殘片位于貝加爾湖西部約深度在2 500 km的地幔深處，延至核幔邊界，估計其年齡至少為150 Ma或更長時間。以此為據，他們推估、對比了古西伯利亞活動邊緣位置和在地幔中與其相應的地震波高速異常帶，勾劃出該區域自晚侏羅紀以來板塊俯沖以及地幔演化圖像。與此同時，他們還分析了東亞和中亞以及印度洋區域從地幔地震層析成像數據中推演出該區域古老俯沖板塊在現代地幔中的殘存蹤跡(板片化石fossil slab)。探討了特提斯海(Tethys)俯沖板塊的消減、古特提斯海閉合以及印度板塊、拉薩(Lhasa)塊體的碰撞過程。他們推定，拉薩地塊和古亞洲大陸之間的中生代特提斯海(Meso-Tethys)大陸俯沖后碰撞斷離的巖石層板片(殘片I)，位于塔吉克斯坦(Tajeikistan)、北阿富汗(northern Afghanistan)、西藏和喜馬拉雅下方，深度為1 000～1 900 km之間，年齡為130 Ma。推斷出新生代特提斯海(Neo-Tethys)的兩部分，北部向拉薩地塊俯沖以及南部則為在該海中部的洋-洋俯沖帶。其俯沖終結于印度板塊和歐亞板塊碰撞后在地幔中留下的殘片(II和III)。其南部洋-洋俯沖殘留板片(III)較早脫離巖石層，時代大約為中生代白堊記后期85～70 Ma之間，在地幔中的深度為1 600～2 300 km。印度大陸和歐亞板塊碰撞后的殘留板片(II)的年齡推測為66～65 Ma，在地幔中的深度為1 000～1 900 km。由于印度大陸的向北推進，其殘片II和III分別位于印度大陸的下部。最后，他們還發現在印度板塊和歐亞板塊碰撞帶下方現代俯沖板塊的蹤跡(IV)，其深度為600 km左右。Hall和Spakman詳細分析了澳大利亞北部和東北部地幔地震層析所顯示的地幔中一系列高速異常體，這些高速異常體與俯沖板塊相關，在 Tonga-Kermadec地區可以深達1 500 km。他們利用現代板塊構造重建方法，分析解釋這些高速異常和現代及過去板塊俯沖系統之間的關系得到，許多很強的異常區和現代俯沖板塊有關。推斷的俯沖板塊長度和位置則與在 New Britain和Halmahera島弧下面晚第三紀的俯沖相一致，也與自10 Ma以來曾經急速后退的 Tonga和 New Hebride海溝位置一致。同時，他們還解釋了一系列更深的異常，并將這些異常與25 Ma，45～25 Ma過去該區域的俯沖帶聯系起來，推測了該地區的地幔動力學過程。值得提出的是，一種俯沖帶4D時空演化模式的研究在東亞這一構造復雜的地帶得到應用。他們以板塊構造重建印度板塊-歐亞板塊，澳大利亞-歐亞板塊碰撞歷史為依據(圖8)，選擇了7個剖面，通過匹配板塊俯沖-碰撞邊界位置的方式，追蹤50 Ma以來俯沖板塊在地幔中的殘片或區域(圖 9)。他們發現板片的輪廓與亞洲邊緣形變相匹配。同時，他們還估計了俯沖板片在地幔中沉降速度在700 km以下為2 cm/a，而在700 km以上為5 cm/a，估計由于亞洲巖石層變形產生的吸收匯聚總量至少有1 500 km。

為了更好地認識和理解地球表層(巖石層)物質通過俯沖消減進入地幔深部的動力學過程，一些研究用數值模擬的方式追蹤了全球俯沖帶物質運移。Steinberger利用俯沖板片密度異常驅動地幔對流模型，顧及到120 Ma以來全球板塊邊界演化運動的歷史;Gordon和Jurdy研究了俯沖板片在地幔中的運動演化。模型顯示，俯沖板片在100 Ma時間內插入地幔的平均深度為1 700 km，而在水平方向上的運動為600～700 km。這一結果和觀測到的俯沖板片的位置和下地幔橫向不均勻結構有很好的吻合。Fu等［24］利用板塊運動和地幔密度異常驅動地幔對流模型探討了俯沖板片在地幔中運動的軌跡，發現多數俯沖板片都能俯沖到地幔的底部，而且由于板塊運動環型場的作用，其俯沖板片的軌跡遠遠比單一極型場驅動復雜得多。利用這一方式，我們可以對一些特定的俯沖板片詳細地追蹤和研究。俯沖板片動力學模型一直受到關注，一些研究利用數值模擬的方式，從地幔熱動力學模型出發探討俯沖板片所面臨的特殊問題。如，俯沖板片的非對稱性，影響俯沖板片消減強度的因素及俯沖板片與相變區域、黏滯分層或地幔中可能的化學壁壘之間的相互作用等。文獻［25，26］用2D和3D數值模擬探討了俯沖板片在地幔中殘存時間，發現板片殘存時間由兩個基本因素確定:板片的初始長度和地幔黏滯分層。殘存時間隨板片長度增加而減小，隨分層黏滯度的增加而加大。當通過660 km地震不連續面黏滯度增加100倍時，地幔深部板片在150 Ma期間拖動其上部向下地幔運動，殘留在地幔1 500 km以上的淺的板片殘存時間超過150 Ma。

圖8 印度板塊和歐亞板塊相對于穩定的西伯利亞在50、40、30、15、10、5 Ma位置的重建Fig.8 Reconstruceion of the India plate and Eurasian plate relative to the stable Siberia in 50，40，30，15 and 5 Ma positiongs

圖9 印度和東南亞地區100～1 600 km地幔地震波速度結構［27］(6個試驗的相關性顯示重建亞洲板塊邊界的位置和形狀(圖9，點線)，以及相應深度上高速異常板片:5 Ma，200 km;10 Ma，400 km;15 Ma，600 km;30 Ma，900 km;40 Ma，1100 km;and 50 Ma，1300 km)Fig.9 Seismic velocity structure［27］from 100 to 1 600 km in India and Southeast Asia of the mantle

近十余年來，熱柱的深地幔追蹤一直是地球動力學的重要課題。因為自Mogan提出地幔熱柱的概念和理論以來，它們在地幔中的現代表觀就成為該理論能否成立的核心。有關熱柱的討論已經很多，其主要集中在:什么是熱柱?熱柱的形態是什么?熱柱源于何處?熱柱形成的動力學機制是什么?熱柱和板塊構造之間的關系等等。現代地震層析成像以及由地球化學提供的有關海島玄武巖(OIB)的數據說明，來自下地幔深處乃至D”層的地幔熱柱可能駕馭著地幔深部熱物質上升流動，其部分熱柱可能從核幔邊界直接穿越地幔而到達地球表面，參與、驅動著地球表層構造和動力學演化。這一理念來源于大量地震、地熱和地球化學證據。地球科學將熱柱，熱點和海島玄武巖(OIB)聯系起來并用地震層析成像的數據去發現和尋找地幔中熱柱的蹤跡。然而，只有在高精度的地幔地震成像技術的支持下，熱柱假說的實踐支柱才有可能得到確認。

Ritsema等［27］從地震層析成像數據分析追蹤非洲和冰島下部地幔熱物質上涌流動(圖10)。他們從深度為2 500 km到700 km的地幔地震層析切片和剖面，追蹤成長于核幔邊界(CMB)的地震S波低速異常。他們發現，在該區域地幔中，一地震剪切波低速異常(通常被認定為高溫物質)從南大西洋核幔邊界(CMB)一直延伸到東非上地幔，到達Afar區域和活動的東非裂谷下部，成為新生代玄武巖火山區域。體積大約為4 000×200×200 km3的一低速高溫異常物質源于西南大西洋下方的核幔邊界，其向上伸展到中非洲、西北非洲和南印度洋下部。其在水平方向上被彎曲跨越了達45°之巨。Goes等［28］研究中歐地區地震層析成像發現，中歐火山活動和地幔中深度為600～2000 km的低速地震波速度異常相關。據此推定，中歐下部下地幔深處熱的物質上升支持該區域上地幔中小的幔熱柱，以及相應的火山活動。他們推測，歐洲新生代的裂谷和火山和地幔中的主動和被動的上升流相關。而對于夏威夷下面的地幔熱柱更成為其深地幔追蹤的熱點，許多研究顯示夏威夷熱點下面的熱柱源于地幔底部。與此同時，Zhao［29］在討論全球范圍內地震層析數據和地幔熱柱分布特征時指出，在熱點下部的地幔中能夠看到低速異常，但熱柱在地幔中可能被地幔流所扭曲;CMB處的低速異常和地球表面的熱點之間存在很好的相關性(圖11)，說明那些主要熱點對應的強的熱柱可能源于CMB;巨大的低速異常在非洲和中-南太平洋的下部地幔之中(圖10)，相應有兩個源于CMB的超級熱柱，但夏威夷不屬于這兩個超級熱柱。Lei和Zhao［30］給出了該區域地幔層析數據顯示了夏威夷熱柱的確源于地幔底部(圖12)。

盡管熱柱的研究還存在許多疑問，很多細節需要更精密的觀測數據，但有理由認為，駕馭著地幔深部熱物質上升流動的地幔熱柱可能是深部地幔物質快速進入上地幔以及出露地球表面的主要渠道。

6 地幔熱動力學框架

圖10 地震層析成像顯示的超級熱柱源于地幔底部(右:中-南太平洋超級地幔熱柱［27］，A:非洲超級地幔熱柱;B:中-南太平洋超級地幔熱柱)［29］Fig.10 Seismic tomographic images showing super plumes originating in the lowermost mantle［27，29］

在過去20年中，地球深內部領域研究取得了重大的進展，為人類認識和研究地球深內部的結構和動力學提供了一系列重要的成果。值得提出的是，地球科學家統一認識到地球系統及系統研究方法的重要性，1993年美國固體地球科學與社會一書提出:“研究地球一系列過程的一種新方法已經形成，這種方法就是把地球看作一個完整的動力系統，而不是各部分彼此孤立的一種集合體”。Roman和Smith［31］概括到:“地球受一個動力的整體系統控制，這些動力作用在所有組成她的物態之上，包括氣體、流體和固體。這些動力和物態組合成一個復雜的包含許多未知數的圖式”。所以現代地球深內部的探索和研究不再是孤軍奮戰，而形成了一支由地球動力學、地球物理學、地質學、地球化學、實驗物理學和計算機科學等學科交叉研究的領域。地球科學家力圖將核-幔邊界以D”層為代表的深內部結構和動力學過程與以板塊構造為代表，以大陸演化動力學為新的方向為代表的地表構造運動聯系起來，這就是地幔熱動力系統。

圖11 地幔底部P波層析成像和地表熱點分布［29］(紅色代表低速擾動，藍色代表高速擾動百分比，實三角表示地表熱點)Fig.11 P wave tomography of the lower mantle and the hot spots distribution at the surface of the Earth［29］

圖12 夏威夷熱點地區在不同深度上P波地震層析成像［30］(紅色代表低速擾動，藍色代表高速擾動百分比，速度參考模型為ISPA91。三角形為夏威夷熱點)Fig.12 P wave seismic tomography in different depth of Hawaii hot spots area［30］

Hilst基于地震學、地球動力學和礦物物理學觀測、實驗的基本數據，以運動的視覺構建了運動的現代地幔結構的框圖(圖13)。綜合分析這些模型有著共同的特征:在板塊運動的會聚區域，俯沖板塊及其周圍物質的下降流動是地球表層物質消減(重循環)進入地幔以至到地幔深層的主要渠道;在板塊運動發散區域以大洋中脊為中心，其下部深至200～300km呈現出地幔熱物質上涌的格局，也是地幔物質通向地表的主要通道;自D”層的熱柱駕馭著地幔深部熱物質上升流動，這一流動可能被地幔中的流動所扭曲或湮沒，但是至少有部分物質(直接或重循環)仍然可以作為火山熔巖噴發通向地球表面，凝固成新的地殼，成為熱點。

圖13 運動中的現代地幔圖像(上面部分:地球表面到核幔邊界(～2 890 km)從中美洲到非洲P速度變化。左邊藍色結構描述了Farallon板塊俯沖到地幔深部，右邊紅色結構描述了西非深部地幔的低地震波速分布。下面部分:不完整數據族的最佳契合并不唯一，通常受到先驗假設和調整的強烈影響，也不一定要很精確)Fig.13 Mantle view in motion(Top:P-wave velocity variations in the mantle between Central America and Africa，from Earth’s surface to the core-mantle boundary(～2 890 km depth).The blue structure on the left depicts the slab of subducted lithosphere of the former Farallon tectonic plate，which dips deeply into the mantle，whereas the red structures on the right mark the low wave velocity deep beneath West Africa.Bottom:Optimal fits to incomplete data sets(stars) are not unique，often strongly influenced by prior assumptions and regularization，and not necessarily accurate)

作者在1990年代初對地幔對流總結了五種不同的對流型態，提出地幔是多尺度、多形態、復雜的熱動力系統。同時還提出統一熱動力系統中不同源的概念，這一概念是建立在原始地幔非均勻分異的假設上。化學成分的非均勻分異、相變和對流的結果既使得地幔呈現層狀結構，又產生橫向不均勻性。漂浮在地幔之上的大陸就是非均勻分異的最直接證據。既然地幔能分異出非均勻的大陸，那么也沒有理由反對地幔非均勻分異同樣會在地幔的底部D”層中沉淀、聚合一個或幾個冷或熱的沉淀的“大陸”。這幾個沉淀的“大陸”猶如現在漂浮的大陸一樣，在地球幾十億年的演化中沒被對流的攪動所混合而均勻化，它們漂浮在核慢邊界之上，在 D”層中，它們保留了其固有的特征。也許可以推想一些超級地幔熱柱和這些沉沒的“大陸”密切相關，就像俯沖板塊總是和漂浮的大陸相關一樣。圖14是作者對2005年提出的經過修改后的地幔熱動力系統模型，這一模型具備了上面闡述的地幔物質運移的共同特征。在以地幔對流和地幔熱柱為主導的地幔熱動力學系統中，其主要的載體——地幔在不斷地被對流攪動而均勻化，同時也被地幔熱柱攜帶的原始或重循環的物質，被俯沖板塊消減的來自巖石層(地殼)的物質破壞其均勻性。所以，一個成分基本均勻的地幔和基本不均勻的巖石層(特別是地殼)和D”層的同時存在，正是我們這顆行星現代仍然保持活力，仍然充滿大規模的突發事件和令人神秘莫測的魅力所在。

圖14 地幔熱動力學模型［32］Fig.14 Thermo-dynamical model of the mantle［32］

致歉 文中引用的部分參考文獻，由于文章過長由編輯刪除，在此表示歉意!有興趣的讀者可以電郵索取未刪減稿。

1 McNamara A K，Garnero E J and Rost S.Tracking deep mantle reservoirs with ultra-low velocity zones［J］.Earth and Planetary Science Letters，2010，299:1-9.

2 Ni Sidao，et al.Sharp sides to the African superplume［J］.Science，2002，(296):1 850-1 852.

3 Glatzmaier G A and Roberts PH.A 3-dimensional convective dynamo solution with rotating and finitely conducting inner-core and mantle［J］.Phys Earth Planet Inter.，1995，91(1-3):63-75.

4 Song X D and Richards P G.Seismological evidence for differential rotation of the Earth’s inner core［J］.Nature，1996，(382):221-224.

5 Su W，Dziewonski A M and Jeanloz R.Planet within a planet:rotation of the inner core of the Earth［J］.Science，1996，(274):1 883-1 887.

6 Koper Keith D，et al.Evidence for small-scale heterogeneity in Earth’s inner core from a global study of PKiKP coda waves［J］.Earth and Planetary Science Letters，2004，(228):227-241.

7 Krasnoshchekov，et al.Seismological evidence for mosaic structure of the surface of the Earth’s inner core［J］.Nature，2005，(435):483-487.

8 Wen Lianxing.Localized temporal change of the Earth’s inner core boundary［J］.Science，2006，(314):967-970.

9 Murakami，et al.Post-perovskite phase transition in MgSiO3［J］.Science，2004，(304):855-858.

10 Wookey J，et al.Efficacy of the post-perovskite phase as an explanation for lowermost-mantle seismic properties［J］.Nature，2005，(438):1 004-1 007.

11 Hernlund John W，et al.A doubling of the post-perovskite phase boundary and structure of the Earth’s lowermost mantle［J］.Nature，2005，(434):882-886.

12 Anderson and Don L.The scales of mantle convection［J］.Tectonophysics，1998，(284):1-17.

13 Trompert R and Hansen U.Mantle convection simulations with rheologies that generate plate-like behaviour［J］.Nature，1998，(395):686-689.

14 Tackley P J.Mantle convection and plate tectonics:Toward an integrated physical and chemical theory［J］.Science，2000，(288):2 002-2 007.

15 Kelloge B，et al.Stratification in the deep mantle［J］.Science，1999，(283):1 881-1 884.

16 Davaille A.Simultaneous generation of hotspots and superswells by convection in a heterogeneous planetary mantle［J］.Nature，1999，(402):756-760.

17 Van Keken P and Shijie Zhong.Mixing in a 3D spherical model of present-day mantle convection［J］.Earth and Planetary Science Letters，1999，(171):533-547.

18 傅容珊，等.現代板塊驅動地幔塊體置換度及地幔混合的研究［J］.地球物理學報，2007，50(5):1 409-1 417.(Fu Rongshan，et al.Study of the mantle mixing driving by plate motions［J］.Chinese Journal of Geophysics，2007，5-(5):1 409-1 417)

19 Lay T，et al.Partial melting in a thermo-chemical boundary layer at the base of the mantle［J］.Physics of the Earth and Planetary Interiors，2004，(146):441-467.

20 Morgan W J.Convection plumes in the lower mantle［J］.Nature，1971，(230):42-43.

21 Bijwaard H and Spakman W.Closing the gap between regional and global travel time Tomography［J］.J Geophy Res.，1998，(103):30 055-30 078.

22 Van der Voo R，Spakman W and Bijiwaard H.Tethyan subducted slab under India［J］.Earth.Plant Sci Lett.，1999a，(171):7-20.

23 Van der Voo R.Mesozonic subducted slab under Siberia［J］.Nature，1999b，(21):246-249.

24 Fu R S，et al.Effect of rheology structure to flow patterns driving by density anomaly in mantle［J］.Chinese J Geophys.，2005，48(4):900-910.

25 Jarvis G T and Lowman J P.Sinking slabs below fossil subduction zones［J］.Phys Earth Planet Inter.，2005，(152): 103-115.

26 Jarvis Gary T，Julian P and Lowman.Survival times of subducted slab remnants in numerical models of mantle flow［J］.Earth and Planetary Science Letters，2007，(26):23-36.

27 Ristema J，et al.Complex shear wave velocity structure imaged beneath Africa and Iceland［J］.Science，1999，(286):1 925-1 928.

28 Goes S，Spakman W and Bijwaard H.A low source for central European volcanism［J］.Science，1999，(286):1 928 -1 932.

29 Zhao D.Global tomographic images of mantle plumes and subducting slabs:insight into deep Earth dynamics［J］.Physics oftheEarth and PlanetaryInteriors，2004，(146):3-34.

30 Lei Jianshe and Dapeng Zhao.A new insight into the Hawaiian plume［J］.Earth and Planetary Science Letters，2006，(241):438-453.

31 Roman L and Smith H J.The terrestrial web［J］.Science，2000，(288):1 983

32 傅容珊，等.地幔對流與深部物質運移研究的新進展［J］.地球物理學進展，2005，20(1):170-179.(Fu Rongshan，et al.Advancements in the study of mantle convection and the material movements in the deep earth interior［J］.Progress in Geophysics，20-05，20(1):170-179)

EXPLORATIONS OF THE EARTH’S DEEP INTERIORS AND ITS DYNAMICS

Fu Rongshan，Huang Jianhua and Yang Ting
(National Geophysical Observatory at Mengcheng，School of Earth and Space Sciences，University of Science and Technology of China，Hefei 230026)

The geoscience has focused on the explorations of the Earth’s deep interiors and its dynamics for a long time.In the past twenty years the study in this area reached a new stage and a series of new discoveries had been published.In this article，we review the main results obtained from the studies so far and explore its developments in the future.

earth’s deep interiors;plate tectonics;mantle dynamics;seismic tomography;geodynamics

1671-5942(2012)03-0139-14

2012-04-02

國家自然科學基金(41074071)

傅容珊，男，1942年生，教授，博士生導師，1965年畢業于中國科學技術大學地震專業，主要從事地球重力場、地幔對流、巖石層應力場以及大陸演化動力學等方面的研究.E-mail:frs@ustc.edu.cn

P542

A