陳愛華， 許鶴華

（1. 國土資源部海底礦產(chǎn)資源重點實驗室 廣州海洋地質(zhì)調(diào)查局， 廣州 510760； 2. 中國科學院邊緣海地質(zhì)重點實驗室 中國科學院南海海洋研究所， 廣州 510301）

0 引言

通過分析俯沖板塊年齡和熱流值， 馬尼拉俯沖帶屬于熱俯沖[5-11]。 陳愛華等[12]嘗試利用熱俯沖探討馬尼拉俯沖帶的地震分布規(guī)律。 但馬尼拉俯沖帶的俯沖速度和俯沖角度比其它熱俯沖大得多，會對俯沖洋殼的熱結構分布產(chǎn)生很大影響[13]。 馬尼拉俯沖帶經(jīng)歷了復雜的地質(zhì)演化過程， 俯沖板塊深度較深， 同時存在陸殼、 洋陸過渡帶和洋殼，且不同區(qū)域的地質(zhì)條件相差較大[14]。 高翔等[15]選取馬尼拉俯沖帶南、 北兩條測線， 根據(jù)兩處俯沖板塊的角度、年齡和運動特征的差異， 分析了摩擦熱和剪切熱對俯沖帶熱結構的影響， 討論俯沖板塊發(fā)生脫水和部分熔融的位置。 由于研究的問題不同， 建立的模型和考慮的因素也不同， 必然會導致模擬結果產(chǎn)生差異。 近年來， 通過對馬尼拉俯沖帶以及呂宋島弧地球物理測量， 獲得了馬尼拉俯沖帶各構造單元莫霍面的形態(tài)、 走向、 梯度和地殼性質(zhì)， 為本文研究馬尼拉俯沖帶的熱結構提供了最新的地球物理資料[16]。 本文綜合馬尼拉俯沖帶的地震、 古地磁、 地熱與層析成像等多種最新地球物理資料， 選取3 條典型剖面（位置見圖1）， 構建地質(zhì)模型和二維數(shù)值模型， 根據(jù)剖面板塊俯沖角度、 速度等， 應用有限元方法求解非線性瞬態(tài)熱方程， 得到馬尼拉俯沖帶的熱結構特征，為進一步探討馬尼拉俯沖帶地震成因機制提供熱力學依據(jù)。

1 區(qū)域地質(zhì)背景

馬尼拉俯沖帶是菲律賓海板塊向南海板塊仰沖所形成的匯聚邊界， 南起民都洛島西南陸架的海底峽谷， 北部延伸至北緯21°15′， 并且逐漸與臺灣造山帶西部的變形相連接， 空間上呈南北向弧形展布， 斷面呈不對稱的“V”形， 整體表現(xiàn)為向西“凸”形， 形態(tài)為向東傾的貝尼奧夫帶（圖1）[17-18]。

圖1 馬尼拉俯沖帶地形圖Fig.1 Topographies maps of Manila subduction zone

馬尼拉俯沖帶主要受近東西向的擠壓， 大多數(shù)地殼內(nèi)地震的主壓應力軸接近E-W 方向； 菲律賓斷層以北以NW 向擠壓應力為主， 顯示北部以擠壓逆沖為特征； 菲律賓海槽及海溝地區(qū)， 主壓應力方向為NWW 到NEE 之間； 南部較為復雜，存在NW、 NE 和近NS 向擠壓應力， NW 向的擠壓應力主要來自菲律賓板塊繼續(xù)向西北的運動， 以順時針旋轉(zhuǎn)為特征， 且存在左旋走滑斷層[19-20]。

根據(jù)南海海底磁異常對比分析， 多數(shù)學者認為， 南海中央次海盆具有東西向?qū)ΨQ分布的磁條帶(11～5 d)[6]。 大陸裂谷及早期海底擴張發(fā)生于32～26 Ma， 期間伸展應力場由SE 向轉(zhuǎn)為SSE 向；24 Ma 后， 轉(zhuǎn)為NW-SE 向擴張， 最后終止于15 Ma[8,21]。 在中新世， 菲律賓海板塊沿馬尼拉海溝以7.0 cm/a 的速度、 NW55°方向斜向仰沖[22]。 Bautista等[23]的研究也表明， 菲律賓板塊在呂宋島北端以約7 cm/a， 向南逐漸增大， 在棉蘭老島的東南部最大， 達到9 cm/a， 歐亞板塊則以1 cm/a 的速度沿相近方向做同向運動， 由于菲律賓板塊速度比歐亞板塊運動速度快， 二者表現(xiàn)出強烈會聚。 隨著俯沖作用繼續(xù)和菲律賓海微板塊向西北方向移動， 呂宋島弧和歐亞板塊發(fā)生碰撞[22]。 菲律賓海板塊持續(xù)地向西北運動， 碰撞作用加劇并向南擴展， 馬尼拉海溝的北端逐漸進入碰撞造山階段， 形成了非火山?。ㄔ錾ǎ?弧前盆地（北呂宋海槽和西呂宋海槽）-火山?。▍嗡位鹕交。嬙旖M合[24，25]。

2 計算模型

2.1 剖面特征與初始條件

馬尼拉俯沖帶是南海東部的俯沖邊界， 北部受臺灣島碰撞帶的影響， 南部由于巴拉望地塊和菲律賓板塊的作用， 使得地質(zhì)構造復雜。

AA′剖面位于臺灣島南部附近， 沿20°N 的縱波速度結構剖面顯示， 東沙隆起（116°～118°E）處地殼厚度約32 km， 巖石圈厚度約95 km； 往東至120°10′E 處， 地殼減薄至20 km， 此段西部為陸坡區(qū)， 東部南海東北部海盆是洋殼[26]。 由此可以推斷， 剖面AA′俯沖板塊為陸殼， 地殼厚度取為25 km， 巖石圈厚度取95 km。

BB′剖面處于洋陸過渡帶， 往北逐漸為陸殼，往南則為南海中央海盆。 沿18°N 的P 波速度結構剖面表明， 海南島之下的巖石圈厚度在80～85 km左右， 從這里開始至南海中央海盆東部邊緣（118°E）， 巖石圈厚90～95 km。 在南海西北海盆至中央海盆之下（114°～119°E）， 地殼厚度減薄到5～6 km， 莫霍面深度10 km， 在馬尼拉海溝（119°30′E）之下， 莫霍面深度從10 km 向40 km 下降， 但巖石圈厚度未見變化[26]。

尾礦資源化和減量化是指將尾礦直接作為資源進行回收利用或者對尾礦進行再生利用，其方法包括尾礦再選和尾礦制備材料等。尾礦的資源化和減量化利用不僅可提高我國礦業(yè)開發(fā)綜合利用水平，同時還可減少尾礦和廢水等礦業(yè)固體廢棄物的排放，減少礦業(yè)開發(fā)過程對環(huán)境的影響。

CC′剖面總體屬于中央海盆區(qū)域， 也處于古俯沖帶的洋中脊附近， 俯沖板塊為洋殼。 根據(jù)上地幔三維S 波速度結構結果：

（1）沿北緯16°的速度剖面上， 在東經(jīng)114°附近， 軟流層頂界面深度最淺， 在60 km 左右， 向東逐步變深， 在東經(jīng)120°最深， 達到70 km 左右，然后又變淺， 達到60 km 左右。

（2）東經(jīng)115°的速度剖面上， 在北緯16°附近， 出現(xiàn)一個隆起， 此處軟流層頂界面深度最淺，在60 km 左右， 北緯16°以北部分， 向北逐漸變深， 最深達到100 km 左右， 北緯16°以南部分，向南先是變深， 在65 km 左右， 然后又逐步變淺，最淺達到60 km 左右[27]。

南海海盆的巖石圈厚度為65 km 左右， 最淺可達到60 km， 向東到菲律賓海盆平均地殼厚度25.7 km， 平均巖石圈厚度80 km[28]。 BB′和CC′兩條剖面處的磁條帶異常對比大概為12（約35 Ma）和6（約20 Ma）， 根據(jù)半空間冷卻模型進行計算，兩處剖面的巖石圈厚度約為77 km 和58 km [7,29]。 海溝兩側(cè)的臺灣-呂宋島地區(qū)的地殼厚度約為34 km， 巖石圈厚度為70～80 km[15，30]。 根據(jù)層析成像結果可以推測， 三條剖面處的板塊俯沖深度均超過300 km[31-32]。

根據(jù)南海及鄰區(qū)莫霍面深度反演計算， 南海板塊在BB′剖面處地殼中的俯沖角度很平緩， 大約3°， 中央海盆莫霍面深度較厚， 約為19～22 km，馬尼拉俯沖帶處地殼厚度為15～20 km； 由于CC′剖面處靠近古俯沖洋中脊， 中央海盆莫霍面深度減薄至10 km 左右， 馬尼拉俯沖帶地殼厚度約為10～15 km， 地殼中的俯沖角度不大， 約4°左右[16]。 B autista 等[23]根據(jù)地震分布和震源機制解， 推斷南海俯沖板片表面形態(tài)， 馬尼拉俯沖帶南海俯沖板片的形態(tài)復雜， 底部的俯沖角度變化非常大， AA′剖面俯沖板塊在淺部的俯沖角度約5°， 深部約60°，至18°逐漸變緩， 向南俯沖角度又逐漸變陡， 至最南端的民都洛島處， 俯沖角度近乎垂直。

模型中選取的尺度較大， 洋殼沉積層厚度較薄， 且本文的模擬結果主要針對中深源地震， 為方便建模和計算， 忽略洋殼沉積層。 在BB′剖面菲律賓以下的南海板塊俯沖角度達到45°； CC′剖面南部剖面深部約65°， 綜合層析成像和地震剖面資料， 本文中CC′剖面取為60°。

綜合上述地質(zhì)資料及分析， 地質(zhì)模型參數(shù)選取如表1 所示。

2.2 模型參數(shù)

根據(jù)馬尼拉俯沖帶的地質(zhì)資料， 可以建立理論數(shù)值模型（圖2）。 為了研究俯沖板塊的熱結構及減少邊界影響， 本文中設定計算區(qū)域馬尼拉海溝橫向?qū)挾葹? 000 km， 模擬底邊界深度為600 km，板片俯沖深度為300 km。 陸殼模型的初始溫度由穩(wěn)態(tài)熱傳導方程給出， 洋殼俯沖模型的初始溫度根據(jù)菲律賓板塊的年齡及厚度， 通過半空間冷卻模型的計算方法給出。 上邊界溫度取為0℃， 底邊界處溫度取為1 600 ℃， 海洋板塊巖石圈邊界和上覆板塊巖石圈邊界均采用半空間冷卻模型計算出的溫度分布作為邊界條件， 軟流圈地幔左右邊界均為自由邊界。 模型各參數(shù)如下表2。

表1 地質(zhì)模型參數(shù)Table 1 Parameters of geology model

圖2 數(shù)值模型示意圖Fig.2 Schematic diagrams of numerical model

表2 數(shù)值模型參數(shù)Table 2 Parameters of numerical model

根據(jù)臧紹先等[33]的研究結果， 俯沖洋殼的剪切生熱以及地幔中的放射性物質(zhì)的產(chǎn)熱量A 相對其他熱源項要小很多， 忽略它們作為一種近似是合理的。 由于AA′剖面俯沖板塊為陸殼， 考慮到在同一圈層差異很小， 也為方便計算， 放射性生熱率取為1*10-6W/m3。

3 模擬結果

對建立的兩個數(shù)值模型賦予相同的邊界條件，為驗證本模擬結果， 采用熱流作為檢驗標準。 AA′剖面屬陸殼， 地殼厚度較厚， 海溝西側(cè)熱流值為66 mW/m2， 最低值約為45 mW/m2， 海溝東側(cè)約為72 mW/m2， 陸殼本身熱流值較低， 加上角落流、巖漿活動和火山活動等因素的影響， 上覆板塊熱流值會相對上升， 使得海溝東側(cè)熱流值高于西側(cè)；BB′剖面中， 在俯沖較淺位置， 熱流值為68 mW/m2， 海溝最低處熱流值為45 mW/m2， 在菲律賓板塊， 熱流值約為65.8 mW/m2； CC′剖面由于巖石圈厚度較薄， 淺部達到70 mW/m2， 而海溝最低熱流值約為20 mW/m2， 與BB′剖面處的俯沖角度和俯沖速度有關， 菲律賓板塊熱流值與BB′剖面基本一致。 模型中的熱流值和熱流趨勢與南海熱流特征一致[9,10,34]， 因此， 此模擬結果可以較好地反應馬尼拉俯沖帶的熱結構。

圖3 AA′、BB′和CC′剖面計算熱結構圖Fig.3 Calculated thermal structure for AA′、BB′ and CC′cross-section

根據(jù)圖3 可知， 在俯沖的淺部， 由于俯沖角度平緩， 俯沖板塊的溫度主要受熱傳導作用控制，基本與俯沖板塊形狀保持一致， 隨俯沖深度加深而增加。 當深度加深， 俯沖角度逐漸增大， 縱向的分速度也逐漸增大， 較冷的俯沖板塊運動至較深的深度。 至軟流圈邊界時， 開始受軟流圈對流影響， 俯沖板塊溫度變化明顯， AA′剖面俯沖板塊頂部約為570 ℃， BB′剖面洋殼頂部溫度比AA′高，約為680 ℃， 可能是由于BB′剖面厚度較AA′剖面薄， 相對更 “熱”。 CC′剖面的俯沖角度和俯沖速度均較BB′剖面大， 溫度也較低， 洋殼頂部僅有約450 ℃。 當洋殼至軟流圈后， 洋殼溫度快速升高，在90 km 處， AA′剖面板塊頂部溫度達到730 ℃，底部增加較緩， 約為650 ℃； BB′剖面表面溫度約為805 ℃， 底部溫度約為673 ℃； CC′剖面表面溫度約為885 ℃， 底部溫度比BB′剖面低的多， 約為210 ℃。 由于CC′剖面的俯沖速度和角度均比AA′和BB′剖面大， 使得軟流圈對洋殼內(nèi)部“加熱”作用比AA′和BB′剖面弱的多， 保持較冷的狀態(tài)到更深的深度。 經(jīng)過軟流圈的加熱后， 洋殼表面和底部溫度均有較大升高， 在160 km 時， AA′剖面溫度約875 ℃， 底部溫度約為671 ℃； BB′剖面洋殼頂部溫度升至950 ℃， 洋殼底部溫度也較高， 達到840℃； CC′剖面洋殼頂部溫度更高， 達到1 050 ℃， 洋殼底部溫度也升高到500 ℃。 三條剖面的熱結構對比也表明， 俯沖帶熱結構主要受俯沖角度、 俯沖速度和俯沖板塊本身地質(zhì)條件如軟流圈深度、厚度、年齡等因素影響。

4 地震意義

Peacock 等[35]和Kirby 等[36]將熱流值＞75 mW/m2，且年齡＜15～25 Ma 的洋殼俯沖稱為熱俯沖， 其它熱流較低、年齡較老的俯沖為冷俯沖。 BB′和CC′剖面的地溫梯度值比日本西南部熱俯沖的地溫梯度還高， 因此， 本模擬結果也進一步證明BB′和CC′剖面屬于熱俯沖[35]。 根據(jù)熱俯沖的研究， 地震成因機制主要有以下兩種觀點：

（1）俯沖深度比較淺[37]。

（2）含水相脫水和榴輝巖的形成， 在俯沖板塊與軟流圈接觸之前停止[36]。

一般認為， 淺源地震是由脆性斷裂導致的，大洋俯沖帶深部地震基本由巖石相變引起[38-39]。 當俯沖洋殼開始俯沖時， 上覆巖石壓力和溫度增大，大多數(shù)的孔隙被擠壓消失， 排出大量孔隙水。 在俯沖的早期階段， 隨著變質(zhì)程度從沸石相的遞增，相應的主要含水礦物為沸石、 綠纖石等， 巖石中的含水量可達到8%～9%[40]。 隨著進一步俯沖， 洋殼在約20 km 的深度進入綠簾石藍片巖相區(qū)域，在50 km 的深度進入角閃石榴輝巖相[41]。 根據(jù)本文模擬結果， 當俯沖至軟流圈邊界時， 開始受軟流圈的影響， 俯沖板塊溫度迅速升高， 自由水降低200～300 ℃反應溫度， 使俯沖洋殼大量脫水， 孔隙壓力繼續(xù)增加， 上覆巖石壓力減小， 巖石脫水脆化發(fā)生破裂， 生成的水促進周圍無水輝長巖和玄武巖相變， 形成大量榴輝巖， 伴隨著密集的地震發(fā)生[36,41-43]。 榴輝巖的形成也使洋殼密度增大， 板塊應力改變， 負浮力增加。 由于BB′剖面的俯沖角度和俯沖速度比CC′剖面小， 使得BB′剖面發(fā)生地震的深度更淺。 俯沖洋殼頂部溫度較高， 在南海洋殼到達軟流圈之前， 含水相脫水和榴輝巖的形成基本停止， 使軟流圈之上的洋殼發(fā)生密集的地震活動。 洋殼底部溫度更低， 需要在更深的深度才能形成榴輝巖， 地震活動也持續(xù)到更大的深度，當洋殼全部進入干榴輝巖相時， 地震活動也基本停止， 這種推測與馬尼拉俯沖帶地震震源深度分布結果吻合。 但進入榴輝巖相的具體深度， 則需要結合實驗巖石學結果， 詳細分析俯沖洋殼的相變。

5 結論

通過本次研究， 得出以下幾點認識和結論：

（1）模擬計算結果表明， 俯沖帶熱結構主要受俯沖角度、 俯沖速度和俯沖板塊本身地質(zhì)條件等因素影響。

（2）根據(jù)模擬結果的地溫梯度數(shù)據(jù)， 進一步證明BB′剖面和CC′剖面屬于熱俯沖。

（3）當洋殼俯沖至軟流圈邊界時， 受軟流圈的影響， 俯沖板塊溫度迅速升高， 容易形成大量榴輝巖， 產(chǎn)生地震活動。

致謝：感謝中國科學院南海海洋研究所張佳政博士給予的幫助。

[1] 張健， 汪集旸. 南海北部大陸邊緣深部地熱特征[J]. 科學通報， 2000， 45 (10)： 1095-1100.

[2] 宋茂雙， 謝洪森， 徐有生. 俯沖帶的熱結構研究進展[J].地質(zhì)地球化學， 1996（2）： 18～22.

[3] 臧紹先， 寧杰遠， 景志成. 俯沖帶流變性質(zhì)的研究[J].中國科學：D 輯， 2001， 31（09）： 705-711.

[4] 周蒂， 王萬銀， 龐雄， 等. 地球物理資料所揭示的南海東北部中生代俯沖增生帶[J]. 中國科學： D 輯, 2006,36 (03)： 209-218.

[5] Taylor B, Hayes D E. The tectonic evolution of the South China Basin[J]. The tectonic and geologic evolution of Southeast Asian seas and islands， 1980， 1： 89-104.

[6] Taylor B, Hayes D. Origin and history of the South China Sea basin [J]. The tectonic and geologic evolution of Southeast Asian seas and islands， Part， 1983， 2： 23-56.

[7] Briais A， Patriat P， Tapponnier P. Updated interpretation of magnetic anomalies and seafloor spreading stages in the South China Sea： implications for the Tertiary tectonics of Southeast Asia [J]. Journal of Geophysical Research，1993， 98 (B4)： 6299-6328.

[8] 孫珍， 孫龍濤， 周蒂， 等. 南海巖石圈破裂方式與擴張過程的三維物理模擬[J]. 地球科學(中國地質(zhì)大學學報)， 2009 (03)： 435-447.

[9] He L， Wang K， Xiong L， et al. Heat flow and thermal history of the South China Sea[J]. Physics of the Earth and Planetary Interiors， 2001， 126 (3-4)： 211-220.

[10] Shi X， Qiu X， Xia K， et al. Characteristics of surface heat flow in the South China Sea[J]. Journal of Asian Earth Sciences, 2003， 22（3）： 265-277.

[11] 劉紹文， 施小斌， 王良書， 等. 南海成因機制及北部巖石圈熱-流變結構研究進展[J]. 海洋地質(zhì)與第四紀地質(zhì)， 2006， (04)： 117-124.

[12] 陳愛華， 許鶴華， 馬輝， 等. 馬尼拉俯沖帶缺失中深源地震成因初探[J]. 華南地震， 2011， 31（04）： 98-107.

[13] Yang T F， Lee T， Chen C H， et al. A double island arc between Taiwan and Luzon: consequence of ridge subduction[J]. Tectonophysics， 1996， 258（1-4）： 85-101.

[14] Chemenda A, Yang R, Hsieh C, et al. Evolutionary model for the Taiwan collision based on physical modelling[J].Tectonophysics， 1997, 274 (1-3)： 253-274.

[15] 高翔， 張健， 孫玉軍， 等. 馬尼拉海溝俯沖帶熱結構的模擬研究[J]. 地球物理學報， 2012， 55（1）： 117-125.

[16] 秦靜欣， 郝天珧， 徐亞， 等. 南海及鄰區(qū)莫霍面深度分布特征及其與各構造單元的關系[J]. 地球物理學報，2011， 54 (12)： 3171-3183.

[17] Hsu S K， Yeh Y， Doo W B， et al. New bathymetry and magnetic lineations identifications in the northernmost South China Sea and their tectonic implications[J]. Marine Geophysical Research, 2004, 25 (1)： 29-44.

[18] Pautot G， Rangin C. Subduction of the south China sea axial ridge below Luzon ( Philippines) [J]. Earth and Planetary Science Letters， 1989， 92 (1)： 57-69.

[19] 臧紹先， 陳奇志， 黃金水. 臺灣南部-菲律賓地區(qū)的地震分布、 應力狀態(tài)及板塊的相互作用[J]. 地震地質(zhì)，1994， 16（01）： 29-37.

[20] 朱俊江， 丘學林， 詹文歡， 等. 南海東部海溝的震源機制解及其構造意義[J]. 地震學報， 2005 (03)： 260-268.

[21] 孫珍， 鐘志洪， 周蒂， 等. 南海的發(fā)育機制研究： 相似模擬證據(jù)[J]. 中國科學： D 輯， 2006（09）： 797-810.

[22] 丁巍偉， 程曉敢， 陳漢林， 等. 臺灣增生楔的構造單元劃分及其變形特征[J]. 熱帶海洋學報， 2005， 24（05）： 53-59.

[23] Bautista B C， Bautista M L P， Oike K， et al. A new insight on the geometry of subducting slabs in northern Luzon, Philippines[J]. Tectonophysics， 2001， 339 (3-4)：279-310.

[24] 陳志豪， 李家彪， 吳自銀， 等. 馬尼拉海溝集合形態(tài)特征的構造演化意義[J]. 海洋地質(zhì)與第四紀地質(zhì)，2009， 29（2）： 59-65.

[25] Clift P, Vannucchi P. Controls on tectonic accretion versus erosion in subduction zones: Implications for the origin and recycling of the continental crust[J]. Reviews of Geophysics, 2004， 42， RG2001.

[26] 姚伯初， 萬玲. 南海巖石圈厚度變化特征及其構造意義[J]. 中國地質(zhì)， 2010（04）： 888-899.

[27] 曹小林， 朱介壽， 趙連鋒， 等. 南海及鄰區(qū)地殼上地幔三維S 波速度結構的面波波形反演[J]. 地震學報,2001（02）： 113-124+223.

[28] 蔡學林， 朱介壽， 曹家敏， 等. 東亞西太平洋巖石圈三維結構及其地幔動力學[J]. 地學前緣， 2007（03）：21-38.

[29] Dl T，G S，Eds.. Geodynamics[M]. Lambridge： Cambridge University Press， 2002.

[30] 李志偉， 胥頤， 郝天珧， 等. 南海東北部及其鄰近地區(qū)地殼上地幔P 波速度結構[J]. 地學前緣， 2009， 16（04）： 252-260.

[31] Lallemand S， Font Y， Bijwaard H， et al. New insights on 3-D plates interaction near Taiwan from tomography and tectonic implications[J]. Tectonophysics， 2001， 335(3-4)： 229-253.

[32] 瞿辰， 周蕙蘭， 趙大鵬. 使用縱波和橫波走時層析成像研究菲律賓海板塊西邊緣帶和南海地區(qū)的深部結構[J]. 地球物理學報， 2007， 50（6）： 1757-1768.

[33] 臧紹先， 宋惠珍， 寧杰遠.日本海俯沖帶的熱結構及熱源的影響[J]. 地球物理學報， 1993， 36（2）： 164-173.

[34] 施小斌， 丘學林， 夏戡原， 等. 南海熱流特征及其構造意義[J]. 熱帶海洋學報， 2003（2）： 63-73.

[35] Peacock S M， Wang K. Seismic consequences of warm versus cool subduction metamorphism: Examples from southwest and northeast Japan[J]. Science， 1999， 286(29)： 937-939.

[36] Kirby S, Engdahl E, Denlinger R. Intermediate-depth intraslab earthquakes and arc volcanism as physical expressions of crustal and uppermost mantle metamorphism in subducting slabs [J]. Geophysical monograph， 1996（96）： 195-214.

[37] Sugi N, Uyeda S. Subduction of young oceanic plates without deep focus earthquakes [J]. Bull. Soc. Geol.France， 1984， 26： 245-254.

[38] Liu L. Phase transformations, earthquakes and the descending lithosphere [J]. Physics of the Earth and Planetary Interiors， 1983， 32（3）： 226-240.

[39] Dobson D， Meredith P， Boon S. Simulation of subduction zone seismicity by dehydration of serpentine[J]. Science， 2002， 298 (5597)： 1407-1410.

[40] Peacock S M. The importance of blueschist to eclogite dehydration reactions in subducting oceanic crust [J].Geol Soc Am Bull， 1993， 105： 684-694.

[41] 魏春景， 張穎慧. 俯沖大洋巖石圈的相轉(zhuǎn)變與俯沖帶巖漿作用[J]. 科學通報， 2008， 53（20）： 2449-2459.

[42] Wang K， Hyndman R D， Yamano M. Thermal regime of the Southwest Japan subduction zone: effects of age history of the subducting plate [J]. Tectonophysics，1995， 248 (1-2）： 53-69.

[43] Hyndman R， Wang K. Thermal constraints on the zone of major thrust earthquake failure: The Cascadia subduction zone [J]. Journal of Geophysical Research，1993， 98（B2）： 2039-2060.