華南地區(qū)地幔過渡帶結(jié)構(gòu)及其動力學(xué)意義

張曉青, 徐濤, 陳立春, 李喆祥

1 桂林理工大學(xué)地球科學(xué)學(xué)院, 桂林 541004

2 中國科學(xué)院地質(zhì)與地球物理研究所, 中國科學(xué)院礦產(chǎn)資源研究重點(diǎn)實(shí)驗(yàn)室, 北京 100029

3 中國科學(xué)院地球科學(xué)研究院, 北京 100029

0 引言

華南大陸主要由揚(yáng)子克拉通和華夏古陸兩部分組成,位于亞歐板塊東南緣,同時受特提斯構(gòu)造域和太平洋構(gòu)造域的影響(圖1; Li et al., 2012, 2016; Faure et al., 2017; 吳福元等, 2020).自新元古代(900～800 Ma)揚(yáng)子與華夏塊體碰撞拼合以來,華南塊體持續(xù)受到周邊板塊(包括華北克拉通、太平洋板塊和印度板塊)的擠壓、巖漿侵入和殼幔相互作用等影響,構(gòu)造演化歷史復(fù)雜,至少經(jīng)歷了拼合后的快速裂解、兩期陸內(nèi)造山作用(Chu et al., 2012a,b)以及晚中生代的構(gòu)造體制轉(zhuǎn)換(舒良樹, 2012; 張國偉等, 2013)等.對華南地區(qū)地幔過渡帶開展系統(tǒng)研究,揭示不同區(qū)域板片俯沖與地幔過渡帶結(jié)構(gòu)變化的關(guān)系,對認(rèn)識太平洋構(gòu)造域和特提斯構(gòu)造域各自的深部影響范圍和探討華南地區(qū)構(gòu)造演化的深部動力學(xué)機(jī)制有重要意義.

20世紀(jì)60年代,地震學(xué)家通過體波走時數(shù)據(jù)發(fā)現(xiàn),在地幔內(nèi)部約410 km和660 km存在兩個速度間斷面(本文簡稱為410-km界面和660-km界面),這兩個間斷面之間的區(qū)域后來被定義為上、下地幔之間的過渡帶,即地幔過渡帶(Shearer, 2000; Nail, 2021).目前,地幔過渡帶已成為地球內(nèi)部結(jié)構(gòu)研究的熱點(diǎn),精細(xì)探測地幔過渡帶結(jié)構(gòu)對于揭示地球內(nèi)部溫度、組分以及動力學(xué)演化等問題具有重要作用(周曉亞等, 2015; 俞春泉等, 2023).由于橄欖石是地幔中的主要造巖礦物,基于橄欖石的高溫高壓物理實(shí)驗(yàn),多數(shù)學(xué)者認(rèn)為地幔過渡帶的強(qiáng)烈速度變化是由于橄欖石的相變引起的(Ringwood, 1975; Jackson, 1983; Ito and Takahashi, 1989).α相的橄欖石和γ相的林伍德石分別在410 km和660 km的溫壓條件下,會變成β相的瓦茲利石和鈣鈦礦及鎂質(zhì)方鐵礦(方鎂石、方鐵礦和超石英),兩者分別具有正的(dP/dT>0)和負(fù)的(dP/dT<0)相變克拉珀瓏斜率(Kreutzberger et al., 1986; Katsura and Ito, 1989; Irifune et al., 1998; Katsura et al., 2004; Zhu et al., 2021).目前,盡管由于缺乏基于深部樣品的礦物物理學(xué)和地球化學(xué)研究,地幔過渡帶的組分、結(jié)構(gòu)和形成機(jī)制尚懸而未決,但410-km界面、660-km界面以及地幔過渡帶的厚度對于這些問題的解決以及對區(qū)域地球動力學(xué)的研究仍然具有重要的參考價值.

研究區(qū)域尺度上地幔速度間斷面和速度結(jié)構(gòu)的常用地震學(xué)方法主要有遠(yuǎn)震層析成像、P波三重震相波形擬合、接收函數(shù),以及基于背景噪聲互相關(guān)函數(shù)的共中心點(diǎn)反射波疊加方法等.前人在華南及鄰區(qū)已經(jīng)開展了大量基于上述方法的相關(guān)研究,其中全球或區(qū)域?qū)游龀上瘾@得了俯沖、滯留板片及可能的海南地幔柱的深部形態(tài)和范圍(Fukao et al., 2001; Zhao, 2001, 2004; Lebedev and Nolet, 2003; Montelli et al., 2004, 2006; Huang and Zhao, 2006; Lei et al., 2009; Li and van der Hilst, 2010; Zhao et al., 2011; Fukao and Obayashi, 2013; Chen et al., 2015; Huang et al., 2015; Lei and Zhao, 2016; 張昌榕等, 2018; 曲平等, 2020; 徐峣等, 2022; Feng et al., 2022).P波三重震相波形擬合則獲得了準(zhǔn)確的區(qū)域一維P波速度結(jié)構(gòu),不同方位的一維速度模型揭示了華南地區(qū)上地幔結(jié)構(gòu)的橫向差異以及410-km界面頂部低速層的存在及其含水量(Revenaugh and Sipkin, 1994;李國輝等, 2014; 周曉亞等, 2015; 呂苗苗等, 2017).相比于上述兩種方法,遠(yuǎn)震P波接收函數(shù)具有更高的橫向和垂向分辨率,該方法在華南地區(qū)的應(yīng)用也多有報(bào)道,但是研究區(qū)域均比較分散,一般僅聚焦于特提斯構(gòu)造域或太平洋構(gòu)造域的局部地區(qū),缺乏對華南地區(qū)地幔過渡帶結(jié)構(gòu)的整體研究和對其淺部響應(yīng)的聯(lián)合分析(Ai et al., 2007; Wang and Niu, 2011; 葉卓等, 2013; Huang et al., 2014;Zhang et al., 2017; Lin et al., 2021).

圖1 東亞地區(qū)構(gòu)造簡圖(修改自Li et al., 2016)華南塊體周緣分別是華北克拉通、青藏高原東部、思茅—印支塊體東南部和西太平洋板塊;黑色箭頭表示太平洋板塊和印度板塊漂移的方向.

圖2 華南寬頻帶地震臺站分布圖圖中黃色正方形表示寬頻帶地震臺站,藍(lán)色叉號和紅色叉號分別表示遠(yuǎn)震P波在410 km和660 km深度處的透射點(diǎn)位置,透射點(diǎn)依據(jù)IASP91模型進(jìn)行計(jì)算,黑色實(shí)線表示截取的CCP疊加剖面的位置.

基于上述原因,我們收集了華南地區(qū)的寬頻帶固定臺站數(shù)據(jù),利用接收函數(shù)方法開展了覆蓋整個華南地區(qū)的地幔過渡帶結(jié)構(gòu)研究,并討論不同構(gòu)造域的影響范圍和地幔過渡帶區(qū)域差異的原因.此外,為進(jìn)一步分析深部構(gòu)造的淺部響應(yīng),我們還需要參考華南地區(qū)的地殼厚度和平均縱橫波速度比.考慮到這部分工作前人已經(jīng)利用接收函數(shù)H-κ掃描方法做了大量研究和總結(jié)(Chen et al., 2010; 葉卓等, 2013; He et al., 2014; Li et al., 2014; Zhang et al., 2021; 趙延娜等, 2017; 韓如冰等, 2019; 張永謙等, 2019),因此本文不再重復(fù),而將直接予以引用.

圖3 遠(yuǎn)震事件分布紅色三角形坐標(biāo)為(110°E,25°N),大致對應(yīng)華南地區(qū)中部;黃色圓圈表示本研究所用的地震事件,其震中距按臺站和地震事件的實(shí)際經(jīng)緯度進(jìn)行計(jì)算和挑選(30°～90°).

1 數(shù)據(jù)與方法

本文從國家測震臺網(wǎng)數(shù)據(jù)備份中心搜集了華南地區(qū)(18°N—30°N,98°E—122°E)共計(jì)302個寬頻帶固定臺站(圖2)的連續(xù)波形數(shù)據(jù),數(shù)據(jù)采集區(qū)間為2015年9月1日至2016年12月31日.經(jīng)過挑選,震中距位于30°～90°之間且震級≥5.5級的有效地震事件共計(jì)355個(圖3).數(shù)據(jù)處理主要包括地震事件的三分量波形截取(理論直達(dá)P波到時前20 s至后100 s)、去線性趨勢、去均值、尖滅、帶通濾波(0.01～0.24 Hz)、坐標(biāo)旋轉(zhuǎn)(由ENZ坐標(biāo)系旋轉(zhuǎn)至RTZ坐標(biāo)系),以及時間域迭代反褶積計(jì)算(Ligorría and Ammon, 1999).最后,經(jīng)挑選共獲得13032條高質(zhì)量的P波徑向接收函數(shù)(高斯系數(shù)為0.5).對于410-km界面和660-km界面,由于界面深度和P波入射角均較大,無法利用基于多次波的H-κ掃描方法(Zhu and Kanamori, 2000)來估計(jì)界面深度.為此,本文直接采用接收函數(shù)一次波共轉(zhuǎn)換點(diǎn)(CCP)疊加方法(Zhu, 2000)對地幔過渡帶進(jìn)行成像,并基于CCP疊加圖像拾取410-km和660-km界面的實(shí)際深度,然后進(jìn)一步計(jì)算地幔過渡帶厚度.CCP疊加采用基于全波形反演的三維速度模型FWEA18(Tao et al., 2018).

上述帶通濾波參數(shù)、高斯系數(shù)和速度模型的選擇是經(jīng)過實(shí)際參數(shù)測試決定的.首先,我們以一條北西—南東向的長剖面(圖2中MM′剖面)為例,在速度模型相同的情況下(IASP91模型,Kennett and Engdahl, 1991),比較了基于三組不同濾波參數(shù)(帶通濾波0.01～1.24 Hz,高斯系數(shù)2.5;帶通濾波0.01～0.50 Hz,高斯系數(shù)1.0;帶通濾波0.01～0.24 Hz,高斯系數(shù)0.5)的CCP成像結(jié)果,分別如圖4a—c所示.結(jié)果表明,低頻段數(shù)據(jù)(圖4c)的成像結(jié)果信噪比最高,410-km和660-km間斷面的Ps震相連續(xù)可追蹤.因此,本文預(yù)處理階段帶通濾波的頻帶范圍為0.01～0.24 Hz,接收函數(shù)計(jì)算時高斯系數(shù)為0.5.

為了進(jìn)一步測試速度模型對地幔過渡帶結(jié)構(gòu)特征的影響,我們同樣以MM′剖面為例,在濾波參數(shù)均相同的情況下,對比了基于三種不同速度模型的CCP成像結(jié)果(圖4c—e),分別為全球一維IASP91模型、基于伴隨成像的全球三維模型GLAD-M25(Lei et al., 2020)以及基于全波形反演的區(qū)域三維模型FWEA18.結(jié)果表明,410-km和660-km界面的絕對深度、界面本身的相對深度變化都較大程度的依賴于速度模型.為此,我們采用分辨率最高的三維全波形反演模型FWEA18進(jìn)行CCP疊加成像.

基于三維的CCP疊加數(shù)據(jù)體,由100°E至122°E以2°為采樣間隔,截取了12條剖面(圖2黑色實(shí)線),然后分別從12條剖面中提取410-km和660-km界面的實(shí)際深度(圖5).由于遠(yuǎn)震P波在410 km和660 km的透射點(diǎn)分布范圍均遠(yuǎn)大于臺站分布范圍(圖2),對華南北部、東南沿海、海南島等地區(qū)也實(shí)現(xiàn)了較好的數(shù)據(jù)覆蓋,因此接收函數(shù)對地幔過渡帶的成像范圍要略大于地震臺站的分布范圍.

需要指出的是,410-km和660-km界面的成像深度與CCP疊加時所采用的速度模型有關(guān).當(dāng)平均速度偏高時成像界面偏深,當(dāng)平均速度偏低時則成像界面偏淺.盡管我們采用了分辨率較高的三維全波形反演速度模型,但有限的數(shù)據(jù)覆蓋并不能保證所有區(qū)域的絕對速度均很準(zhǔn)確,因此仍然可能存在由于速度模型不確定而引起的界面絕對深度或相對深度的局部變化.為此,本文討論時將重點(diǎn)關(guān)注地幔過渡帶的厚度及其頂、底界面的主要特征.

2 華南大陸地幔過渡帶結(jié)構(gòu)特征

圖5是沿經(jīng)度方向且以2°為采樣間隔提取的12條CCP疊加剖面,圖中藍(lán)色和紅色長虛線分別表示410 km和660 km深度,藍(lán)色和紅色短虛線則分別表示實(shí)際的410-km和660-km界面.結(jié)果顯示,除川滇地塊和滇南地塊下方數(shù)據(jù)質(zhì)量略低以外,華南地區(qū)的410-km和660-km界面整體上清晰連續(xù),且410-km界面平均深度為423 km(圖6a),660-km界面平均深度為684 km(圖6b),均較IASP91全球平均模型偏深;平均地幔過渡帶厚度為261 km(圖6c),較IASP91全球平均模型偏厚.其中,華南大陸東南緣區(qū)域的地幔過渡帶結(jié)構(gòu)特征與葉卓等(2013)的結(jié)果基本一致.

華南地區(qū)內(nèi)部,地幔過渡帶結(jié)構(gòu)表現(xiàn)出顯著的橫向差異,且410-km和660-km界面的特征也各不相同.410-km界面主要呈現(xiàn)出南北差異,北部揚(yáng)子準(zhǔn)克拉通下方(圖6a藍(lán)色虛線區(qū)域),410-km界面深度接近416 km,與IASP91模型相當(dāng);而南部的海南島地區(qū)(圖6a紅色虛線區(qū)域),410-km界面平均深度約為435 km,較IASP91模型明顯偏深;東南沿海的華夏地塊和青藏高原東南緣下方410-km界面深度居中,平均深度約為426 km,較IASP91模型亦偏深.

660-km界面不僅具有南北差異,還呈現(xiàn)出東西差異,大致可以分為圖6b所示的A、B、C、D、E五個區(qū)域,其中E區(qū)660-km界面深度最淺,平均深度約為673 km,且這一深度較全球平均模型亦略偏深;B區(qū)660-km界面深度居中,平均深度約為685 km;以B區(qū)為界,其兩側(cè)的A區(qū)、C區(qū)和D區(qū)660-km界面明顯偏深,A區(qū)大致與青藏高原東南緣對應(yīng),C區(qū)大致位于南北重力梯度帶以東;值得注意的是D區(qū)660-km界面向700 km深度逼近,結(jié)合遠(yuǎn)震層析成像結(jié)果(Fukao et al., 2001; Huang and Zhao, 2006; Li and van der Hilst, 2010; Zhao et al., 2011; Chen et al., 2015)可知,該區(qū)域大致與西太平洋滯留板片的位置相對應(yīng).

基于410-km和660-km界面的深度,我們進(jìn)一步計(jì)算得到了華南地區(qū)地幔過渡帶的厚度分布圖(圖6c),其整體特征與Wang和Niu(2011)獲得的中國東部地幔過渡帶厚度特征基本一致,但本文成像結(jié)果的分辨率更高,這得益于在臺站覆蓋密度、遠(yuǎn)震事件數(shù)量以及所采用的速度模型的分辨率等方面均有所改善.結(jié)果顯示,華南地區(qū)平均地幔過渡帶厚度(～261 km)較IASP91模型偏厚,僅海南島及其北部地幔過渡帶厚度偏薄,與圖6a和圖6b暗紅色虛線所示區(qū)域具有相關(guān)性,但并不完全重合.

3 討論

本文對華南地區(qū)地幔過渡帶結(jié)構(gòu)的整體成像,為研究特提斯構(gòu)造域和太平洋構(gòu)造域的深部影響范圍及其動力學(xué)意義提供了全面的視角.結(jié)果表明,660-km界面(圖6b)的深度主要呈現(xiàn)出東西差異(除海南島地區(qū)以外),與上述兩大構(gòu)造域的作用范圍有很好的對應(yīng)關(guān)系;410-km界面的深度(圖6a)和地幔過渡帶的厚度(圖6c)分布則體現(xiàn)了揚(yáng)子克拉通和海南地幔柱分別作為穩(wěn)定陸塊和地?；钴S地區(qū)的地幔過渡帶結(jié)構(gòu)特征.下面將分區(qū)域予以討論：

3.1 特提斯構(gòu)造域的影響范圍及其深部構(gòu)造演化

如圖6b所示,以B區(qū)為界,其兩側(cè)的A、C、D區(qū)域660-km界面都不同程度下沉,且A、D區(qū)域660-km界面最深.其中,A區(qū)位于華南以西的青藏高原東南緣,410-km和660-km界面均顯著下沉,地幔過渡帶厚度略大于全球平均水平,整體特征區(qū)別于B區(qū)及其以東地區(qū).此外,A區(qū)所對應(yīng)的410-km界面深度也較其東部的揚(yáng)子克拉通地區(qū)有明顯差異.據(jù)此,我們認(rèn)為A區(qū)屬于特提斯構(gòu)造域的深部影響范圍.A區(qū)410-km和660-km界面的絕對深度分別表明,該區(qū)上地幔溫度可能偏高,而地幔過渡帶底部溫度可能偏低.

圖4 采用不同濾波參數(shù)和不同速度模型的CCP疊加結(jié)果(a)—(c) 計(jì)算接收函數(shù)的濾波頻段分別為0.01～1.24 Hz、0.01～0.50 Hz和0.01～0.24 Hz,高斯系數(shù)分別為2.5、1.0和0.5,CCP疊加所采用的速度模型均為IASP91模型(Kennett and Engdahl, 1991);(d)和(e)計(jì)算接收函數(shù)的濾波頻段均為0.01～0.24 Hz,高斯系數(shù)均為0.5,CCP疊加所采用的速度模型分別為GLAD-M25(Lei et al., 2020)和FWEA18(Tao et al., 2018). 藍(lán)色和紅色長虛線分別表示410 km和660 km深度線,藍(lán)色和紅色短虛線分別表示基于實(shí)際CCP疊加剖面解譯得到的410-km界面和660-km界面深度. 剖面位置對應(yīng)圖2中的MM′.

圖5 沿不同經(jīng)度方向截取的CCP疊加剖面藍(lán)色和紅色長虛線分別表示410 km和660 km深度線,藍(lán)色和紅色短虛線分別表示基于實(shí)際CCP疊加剖面解譯得到的410-km界面和660-km界面深度,各剖面對應(yīng)的經(jīng)度見左下角.

圖6 華南地區(qū)地幔過渡帶結(jié)構(gòu)(a) 410-km界面深度分布圖; (b) 660-km界面深度分布圖; (c) 地幔過渡帶厚度分布圖. 灰色粗實(shí)線表示華南地塊邊界(據(jù)張培震等, 2013),灰色細(xì)實(shí)線為各次級塊體邊界(據(jù)朱介壽等, 2005;張培震等, 2013).

當(dāng)然,溫度并不是影響410-km和660-km界面深度的唯一因素,古老的俯沖滯留板片或拆沉的巖石圈物質(zhì)也可能造成地幔過渡帶內(nèi)的化學(xué)成分差異,進(jìn)而影響地幔過渡帶頂、底界面的深度.我們注意到,A區(qū)410-km和660-km界面整體向西傾.這一特征與以往的研究結(jié)果是一致的,包括遠(yuǎn)震層析成像結(jié)果、接收函數(shù)圖像和背景噪聲互相關(guān)函數(shù)共反射點(diǎn)疊加波形均顯示(Hall and Spakman, 2015; Zhang et al., 2017; Feng et al., 2022),青藏高原東南緣地幔過渡帶中存在一個西傾的高速異常體.對于該地幔過渡帶內(nèi)高速異常的成因至今尚未達(dá)成共識,早期有學(xué)者將其解釋為斷離的印度板塊(Huang et al., 2015; Lei and Zhao, 2016),但考慮到該高速異常的傾向(西傾)與印度板塊的俯沖方向相反,近年來有越來越多的學(xué)者提出該異?？赡苁切律鸪恋拇箨憥r石圈,這一觀點(diǎn)可以得到地幔過渡帶上方減薄的地幔巖石圈和廣泛分布的鉀質(zhì)巖漿巖的支持(Zhang et al., 2017; 黃周傳等, 2021; Feng et al., 2022).根據(jù)本文的成像結(jié)果,我們更傾向于后者的觀點(diǎn).

3.2 太平洋構(gòu)造域的影響范圍及其橫向差異

C、D區(qū)域的660-km界面均較B區(qū)有所下沉,且與B區(qū)大致以南北重力梯度帶為界,據(jù)此我們認(rèn)為C、D區(qū)域主要受控于太平洋構(gòu)造域.但值得注意的是,C區(qū)下沉較少,D區(qū)下沉強(qiáng)烈,兩者大致以28°N為界.參考已有層析成像結(jié)果(Fukao et al., 2001; Huang and Zhao, 2006; Chen et al., 2015;Tao et al., 2018)可知,地幔過渡帶中高速的太平洋滯留板片主要位于南北重力梯度帶以東、28°N以北,而華夏地塊下方地幔過渡帶中未發(fā)現(xiàn)顯著的高速異常體,且P波三重震相波形擬合獲得的速度結(jié)構(gòu)特征(Revenaugh and Sipkin, 1994;李國輝等, 2014; 周曉亞等, 2015; 呂苗苗等, 2017)與層析成像結(jié)果一致.據(jù)此我們認(rèn)為,盡管C區(qū)和D區(qū)均受太平洋板塊俯沖影響,但俯沖之后的動力學(xué)過程并不一樣.

D區(qū)與整個華北克拉通東部相似,西太平洋俯沖板片并未穿透660-km界面,而是停滯平躺在地幔過渡帶中,且偏冷的滯留板片可能導(dǎo)致了660-km界面下沉(Chen and Ai, 2009; 呂苗苗等, 2017).但C區(qū)有所不同,可能只是早期經(jīng)歷了太平洋板片俯沖,隨后板片發(fā)生快速后撤,大洋板片并未在地幔過渡帶中滯留,已有地殼結(jié)構(gòu)研究也支持這一推斷.比如：從整體上看,位于東南地區(qū)的華夏地塊地殼厚度明顯薄于揚(yáng)子地塊(He et al., 2014; Li et al., 2014),這與660-km界面所指示的太平洋板塊俯沖的影響范圍是相一致的;其次,在東南沿海地區(qū),VP/VS顯著偏高,局部地區(qū)>1.8(Chen et al., 2010; 葉卓等, 2013; He et al., 2014; Li et al., 2014; 趙延娜等, 2017; 韓如冰等, 2019; 張永謙等, 2019),顯示該局部地區(qū)經(jīng)歷了強(qiáng)烈的鎂鐵質(zhì)巖漿活動.據(jù)此我們推測,華南東南部與華北克拉通東部類似,晚中生代經(jīng)歷了古西太平洋板塊的俯沖作用,軟流圈物質(zhì)上涌并形成強(qiáng)烈的鎂鐵質(zhì)巖漿活動,使地殼減薄、平均波速比顯著升高,地表大規(guī)模出露的中-新生代巖漿巖(Faure et al., 2017)就是最直接的證據(jù),但這一強(qiáng)烈的巖漿活動并未進(jìn)入華南大陸腹地,表明俯沖板片隨后發(fā)生了后撤.此外,巖石學(xué)和地球化學(xué)研究也支持這一時期俯沖作用類型的轉(zhuǎn)變(從平俯沖到陡俯沖;Li et al., 2012).

無論具體的動力學(xué)過程如何,但至少可以表明太平洋板塊的西向俯沖及其在東亞地幔過渡帶的停滯存在區(qū)域性差異,且均對中國東部地殼和上地幔結(jié)構(gòu)產(chǎn)生了重要影響.俯沖板片或滯留板片脫水,可能導(dǎo)致上地幔發(fā)生部分熔融以及軟流圈物質(zhì)大量上涌,并最終導(dǎo)致中國東部巖石圈不同程度減薄.

3.3 海南地幔柱的證據(jù)和起源深度

不同于660-km界面,華南地區(qū)410-km界面(圖6a)和地幔過渡帶厚度(圖6c)主要呈現(xiàn)出南北差異.其中,410-km界面的最淺區(qū)域(圖6a藍(lán)色虛線區(qū)域)與地表揚(yáng)子準(zhǔn)克拉通地塊對應(yīng)較好,并與660-km界面的B區(qū)域(圖6b)大致對應(yīng),其地幔過渡帶厚度接近IASP91模型.以上特征表明揚(yáng)子地塊上地幔結(jié)構(gòu)與全球平均水平相當(dāng),整體較穩(wěn)定,太平洋板塊俯沖和印度板塊俯沖均未對該地區(qū)的地幔結(jié)構(gòu)造成顯著影響.而在揚(yáng)子克拉通周緣,410-km界面明顯變深,表明其周緣上地幔溫度較克拉通下方可能偏高,當(dāng)然這也不能排除周邊地區(qū)由于板片滯留等引起化學(xué)成分差異進(jìn)而導(dǎo)致410-km界面的深度發(fā)生變化.

與揚(yáng)子克拉通不同的是,在華南南部的海南島地區(qū),410-km界面最深(圖6a紅色虛線區(qū)域),660-km界面最淺(圖6b,E區(qū)域),且地幔過渡帶厚度明顯偏薄(圖6c黑色虛線區(qū)域).已有遠(yuǎn)震層析成像研究表明(Huang and Zhao, 2006;Lebedev and Nolet, 2003; Fukao and Obayashi, 2013),海南島地區(qū)地震波速度明顯偏低,且低速異常的深部可達(dá)下地幔.對比可知,該區(qū)偏薄的地幔過渡帶現(xiàn)象與遠(yuǎn)震層析成像結(jié)果所反映的低速異常相對應(yīng).考慮到該偏薄的地幔過渡帶異常與太平洋構(gòu)造域和特提斯構(gòu)造域均無明顯的空間耦合關(guān)系,結(jié)合對該區(qū)新生代玄武質(zhì)巖漿(Ho et al., 2000)的認(rèn)識,我們推測該異?？赡芘c近年來提出的海南地幔柱相關(guān)(Lei et al., 2009; 曲平等, 2020).

目前,對于該地幔柱的起源和結(jié)構(gòu)形態(tài)尚有爭議.大尺度層析成像研究認(rèn)為,海南地幔柱起源于1900 km的中地幔深度或更深至地幔底部,其直徑約為80 km(Zhao, 2001, 2004; Montelli et al., 2004, 2006; Lei and Zhao, 2006; Lei et al., 2009).本文接收函數(shù)成像雖然只能約束地幔過渡帶結(jié)構(gòu),但通過分析410-km、660-km界面深度和地幔過渡帶厚度的橫向差異及其影響因素,依然可以對海南地幔柱可能的起源深度加以約束.首先,海南島地區(qū)410-km和660-km界面分別較鄰區(qū)偏深和偏淺,符合熱的地幔柱穿過地幔過渡帶以后應(yīng)有的現(xiàn)象.其次,該地區(qū)地幔過渡帶厚度最薄,最接近IASP91模型,而其他地區(qū)均略大于全球平均水平.雖然該地區(qū)的地幔過渡帶厚度并沒有明顯薄于IASP91模型,但已有研究表明,正常大陸地區(qū)的地幔過渡帶厚度確實(shí)略大于全球平均值(Gu and Dziewonski, 2002).據(jù)此我們認(rèn)為,海南島地區(qū)以外的區(qū)域,其地幔過渡帶的結(jié)構(gòu)僅受上地幔動力學(xué)過程的影響,因此地幔過渡帶厚度整體上變化不大,僅存在板片滯留的地區(qū)地幔過渡帶厚度明顯偏厚;但海南島地區(qū)偏薄的地幔過渡帶,還受到上涌的下地幔熱物質(zhì)的影響,即海南地幔柱至少起源于中地幔.

4 結(jié)論

華南地區(qū)的12條沿經(jīng)度方向的CCP疊加剖面顯示,華南大陸410-km界面和660-km界面連續(xù)可追蹤,兩個界面的成像深度均略大于全球一維平均模型.410-km界面和地幔過渡帶厚度主要呈現(xiàn)為南北差異,660-km界面既有南北差異又有東西差異.綜合多學(xué)科資料和本文成像結(jié)果,將華南及鄰區(qū)的上地幔大致劃分為5個構(gòu)造域,分別為揚(yáng)子克拉通地區(qū)、華南以西地區(qū)(青藏高原東南緣)、華南東北部、華南東南部和海南島地區(qū),其中揚(yáng)子克拉通地區(qū)相對穩(wěn)定,其余4個區(qū)域則可能分別受到印度板塊俯沖與拆沉巖石圈滯留、太平洋板塊俯沖與滯留、太平洋板塊俯沖與快速回撤,以及海南地幔柱的影響.受深部地幔過渡帶結(jié)構(gòu)的橫向變化影響,淺部地殼結(jié)構(gòu)和地表地質(zhì)也有顯著響應(yīng).

致謝謹(jǐn)以此文祝賀滕吉文先生90華誕暨從事地球物理工作70年.中國地震局地球物理研究所國家數(shù)字測震臺網(wǎng)數(shù)據(jù)備份中心(SEISDMC, doi： 10.11998/SeisDmc/SN, http：∥www.seisdmc.ac.cn; Zheng et al., 2010)為本研究提供了地震波形數(shù)據(jù).