川滇地區(qū)巖石圈溫度、化學(xué)組分結(jié)構(gòu)研究及對巖石圈減薄過程的認(rèn)識(shí)

周萬里, 單斌, 熊熊, 肖陽

中國地質(zhì)大學(xué)(武漢)地球物理與空間信息學(xué)院, 地球內(nèi)部多尺度成像湖北省重點(diǎn)實(shí)驗(yàn)室, 武漢 430074

0 引言

川滇地區(qū)地處青藏高原東南緣,是特提斯—喜馬拉雅造山系的東部轉(zhuǎn)折點(diǎn)之一(Deng et al., 2003),在經(jīng)歷了古-新特提斯漫長的歷史演化過程后,其強(qiáng)烈的構(gòu)造變形特征與青藏高原其他邊緣地帶形成鮮明對比,在解釋和建立大陸巖石圈物質(zhì)變形、青藏高原物質(zhì)東流模型中扮演著重要角色(Clark and Royden, 2000; Bai et al., 2010; Bao et al., 2015),而對該地區(qū)巖石圈結(jié)構(gòu)的深入認(rèn)識(shí)是了解該地區(qū)構(gòu)造變形機(jī)制的關(guān)鍵.目前的研究表明(Yang et al., 2017; Feng et al., 2022),印度板塊與亞歐板塊間的碰撞和擠壓導(dǎo)致川滇地區(qū)南部存在巖石圈減薄的現(xiàn)象,但關(guān)于川滇地區(qū)南部巖石圈減薄機(jī)制問題仍存在較大的爭議.例如:Huang等(2015)在川滇地區(qū)開展的地震體波成像結(jié)果顯示川滇地區(qū)南部存在大范圍的上地幔低速異常,被解釋為印度板塊沿緬甸弧東向俯沖導(dǎo)致川滇地區(qū)南部上地幔熱物質(zhì)上涌,巖石圈發(fā)生減薄,巖石圈減薄的機(jī)制和動(dòng)力來源分別為熱侵蝕作用和印度板塊的東向平俯沖作用(Huang et al., 2015; Lei and Zhao, 2016).Chen等(2017)在川滇地區(qū)開展的巖石學(xué)、地球化學(xué)研究發(fā)現(xiàn)在川滇南部地區(qū)出露部分富鉀鎂鐵質(zhì)巖石,被解釋為印度板塊與亞歐板塊南北向碰撞導(dǎo)致川滇地區(qū)南部巖石圈增厚拆沉,深部的富鉀鎂鐵質(zhì)巖漿也噴出地表,而巖石圈減薄的機(jī)制和動(dòng)力學(xué)來源分別為巖石圈增厚拆沉作用和印度板塊與亞歐板塊的南北向碰撞作用(Chung et al., 1997; Chen et al., 2017; Huang et al., 2019; Feng et al., 2022).這些爭議性問題的存在無疑制約著我們對川滇地區(qū)巖石圈減薄過程及其深部動(dòng)力學(xué)背景等研究的深入認(rèn)識(shí).

巖石圈溫度、化學(xué)組分結(jié)構(gòu)是理解大陸巖石圈形成、變形和破壞的重要屬性(Carlson et al., 2005),因此對川滇地區(qū)巖石圈溫度、化學(xué)組分結(jié)構(gòu)的認(rèn)識(shí)是理解該地區(qū)巖石圈減薄過程及深部動(dòng)力學(xué)背景的基礎(chǔ).巖石圈溫度結(jié)構(gòu)主要是指巖石圈熱厚度及溫度分布,是控制巖石圈巖石學(xué)、流變學(xué)和動(dòng)力學(xué)的重要參量之一(臧紹先等, 2002).目前對巖石圈溫度結(jié)構(gòu)的研究主要來源于地?zé)釋W(xué)、地震學(xué)和重力學(xué)的研究.地?zé)釋W(xué)研究通過求解穩(wěn)態(tài)熱傳導(dǎo)方程得到巖石圈溫度結(jié)構(gòu),地表熱流值作為方程的邊界條件(Hu et al., 2000).但是一方面,大陸地區(qū)地表熱流測量十分昂貴,因此川滇地區(qū)所擁有的地表熱流十分稀少且分布不均勻;另一方面,求解穩(wěn)態(tài)熱傳導(dǎo)方程所需的巖石熱參數(shù)(生熱率和熱傳導(dǎo)系數(shù))取值缺乏約束,求解方程后得到的巖石圈溫度結(jié)構(gòu)不確定度較高.而地震、重力學(xué)研究(Priestley and McKenzie, 2006, 2013; An and Shi, 2006; Fullea et al., 2007; 單斌等, 2008)基于實(shí)驗(yàn)室模擬的經(jīng)驗(yàn)轉(zhuǎn)換關(guān)系得到與實(shí)際地震/重力資料一致的溫度結(jié)構(gòu)(Zeyen and Fernàndez, 1994; Priestley and McKenzie, 2006, 2013; Jiménez-Munt et al., 2008; Ritsema et al., 2009; Cammarano et al., 2011),但是由于巖石圈波速/密度結(jié)構(gòu)受到溫度和化學(xué)組分結(jié)構(gòu)的共同作用,導(dǎo)致波速/密度結(jié)構(gòu)與溫度結(jié)構(gòu)間的轉(zhuǎn)換關(guān)系不準(zhǔn)確,加之相關(guān)參數(shù)難以約束,其所造成的誤差使得不同研究結(jié)果間存在較大爭議(Priestley and McKenzie, 2006, 2013; Fishwick, 2010; Becker, 2012).上述這些因素共同制約了地?zé)釋W(xué)、地震學(xué)、重力學(xué)方法對巖石圈溫度結(jié)構(gòu)精確度和空間分辨率的確定.此外,先前的研究(Priestley and McKenzie, 2006, 2013; Yang et al., 2017; Feng et al., 2022)基于不同研究方法得到川滇地區(qū)南部巖石圈較薄,表明該地區(qū)經(jīng)歷了巖石圈減薄的過程,但是這些研究難以回答川滇地區(qū)南部巖石圈減薄機(jī)制及其動(dòng)力來源等問題,而化學(xué)組分結(jié)構(gòu)能夠提供巖石圈減薄過程的更多信息.

巖石圈化學(xué)組分結(jié)構(gòu)是指構(gòu)成地殼、地幔巖石中各種礦物和元素的含量,是控制巖石圈密度、波速等物理性質(zhì)的重要參數(shù),也是巖石圈地幔演化研究和礦產(chǎn)資源開發(fā)的地球化學(xué)依據(jù)(Rudnick et al., 1998).目前,對川滇地區(qū)巖石圈化學(xué)組分的認(rèn)識(shí)主要來自地球化學(xué)或者巖石學(xué)研究.通過對川滇地區(qū)南部出露的富鉀鎂鐵質(zhì)巖石的分析,Chen等(2017)認(rèn)為由于印度板塊與亞歐板塊的南北向碰撞導(dǎo)致青藏高原物質(zhì)東向流出,而青藏高原物質(zhì)的東向擠壓作用導(dǎo)致川滇地區(qū)南部巖石圈增厚拆沉,富鉀鎂鐵質(zhì)巖漿也噴出地表,因此巖石圈增厚拆沉作用可能是控制川滇地區(qū)南部巖石圈減薄的主要機(jī)制.此外,Liu等(2013)通過分析馬關(guān)地區(qū)(圖1)幔源捕虜體的化學(xué)組分,認(rèn)為右江盆地巖石圈地幔以新生代地幔物質(zhì)為主,說明馬關(guān)地區(qū)巖石圈在減薄的同時(shí),巖石圈地幔物質(zhì)也發(fā)生了替換,熱化學(xué)侵蝕作用可能是主導(dǎo)川滇地區(qū)南部巖石圈減薄的動(dòng)力學(xué)機(jī)制.然而,地球化學(xué)和巖石學(xué)方法所使用的巖石樣品時(shí)空分布十分有限,并且這些巖石樣品在出露到地表的過程中其組分結(jié)構(gòu)不可避免地會(huì)發(fā)生改變.因此,巖石學(xué)、地球化學(xué)觀測難以得到川滇地區(qū)巖石圈全面、可靠的化學(xué)組分結(jié)構(gòu)認(rèn)識(shí),而要解答川滇地區(qū)南部巖石圈減薄過程的問題只能聯(lián)合巖石學(xué)、地球化學(xué)、地球物理的綜合成像研究.

隨著地球物理和地球化學(xué)數(shù)據(jù)的不斷豐富,計(jì)算能力、實(shí)驗(yàn)室測量技術(shù)和多學(xué)科交叉融合研究的迅速發(fā)展,為聯(lián)合多學(xué)科觀測資料研究巖石圈溫度、化學(xué)組分結(jié)構(gòu)提供了數(shù)據(jù)和技術(shù)支持.多觀測概率密度反演能夠充分發(fā)揮不同地球物理觀測資料(地形、地表熱流、大地水準(zhǔn)面高、地震波、大地電磁等)和地球化學(xué)觀測對巖石圈溫度、化學(xué)組分結(jié)構(gòu)不同的敏感度和解析度,進(jìn)而獲取更為全面、可靠的巖石圈溫度、化學(xué)組分結(jié)構(gòu).本文綜合利用川滇地區(qū)地形、大地水準(zhǔn)面高、地表熱流、瑞雷面波相速度頻散曲線等地球物理觀測資料,通過熱動(dòng)力學(xué)正演模擬和概率密度反演的方法,得到川滇地區(qū)巖石圈溫度、化學(xué)組分結(jié)構(gòu),并以此為基礎(chǔ)分析川滇地區(qū)南部巖石圈的減薄過程.

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

本文聯(lián)合多種地球物理觀測反演巖石圈溫度、化學(xué)組分結(jié)構(gòu),反演采用的是Afonso等(2008, 2013a,b)開發(fā)的Litmod_4INV軟件(https:∥www.juanafonso.com/software),該軟件基于熱動(dòng)力學(xué)平衡方程,利用概率密度反演的方法,將不同的地球物理觀測資料融合在一個(gè)反演體系中.本節(jié)將簡單介紹反演所用到的地球物理數(shù)據(jù)和Litmod_4INV軟件的反演方法.

1.1 觀測數(shù)據(jù)

本文聯(lián)合反演使用的地球物理觀測數(shù)據(jù)包括:地表大地?zé)崃?Jiang et al., 2019)、大地水準(zhǔn)面高(Pavlis et al., 2012)、地形(http:/lta.cr.usgs.gov/GTOPO30)和基階瑞雷面波相速度頻散數(shù)據(jù)(Shen et al., 2016;Zhang et al., 2020)(如圖2所示).Shen等(2016)使用背景噪聲成像給出了8～70 s、分辨率為1°×1°的相速度頻散數(shù)據(jù),Zhang等(2020)使用遠(yuǎn)震面波層析成像給出了15～150 s、分辨率為1°×1°的相速度頻散數(shù)據(jù),本文采用兩者的觀測數(shù)據(jù),并在15～70 s的重疊周期內(nèi)取兩者的加權(quán)平均值,獲得每個(gè)地理坐標(biāo)下的8～150 s的相速度頻散曲線,不同周期相速度的標(biāo)準(zhǔn)差根據(jù)蒙特卡洛誤差傳播計(jì)算得到.所采用的地表大地?zé)崃?SHF)數(shù)據(jù)來源于Jiang等(2019),并去除研究區(qū)域<20 mW·m-2和>150 mW·m-2熱流觀測的異常值(Shan et al., 2014; Guo et al., 2016; Qashqai et al., 2016).考慮到研究區(qū)域SHF較為稀少且分布極不均勻(測點(diǎn)大多集中于川滇菱形塊體和四川盆地),本文收集了川滇地區(qū)及其鄰區(qū)所有SHF并通過曲率樣條插值以獲取研究區(qū)域 1° × 1°的熱流值,在沒有SHF的地區(qū)(如松潘甘孜塊體、江南造山帶和右江盆地)插值得到的SHF可信度較低,因此本文通過降低SHF在聯(lián)合反演中所占的權(quán)重(設(shè)置觀測誤差為觀測值的20%)以減小該數(shù)據(jù)對溫度結(jié)構(gòu)的影響,誤差最小值設(shè)置為SHF典型的不確定度10 mW·m-2(Jaupart and Mareschal, 2011).大地水準(zhǔn)面數(shù)據(jù)取自EGM2008全球模型(Pavlis et al., 2012),并濾除了球諧函數(shù)低階部分(2～9階)信號(hào)以消除深部地幔密度異常的影響(Afonso et al., 2019).地形數(shù)據(jù)取自GTOPO30 的地形數(shù)據(jù)(http:/lta.cr.usgs.gov/GTOPO30),并使用一個(gè)低通濾波去除波長小于50 km的高頻成分.每個(gè) 1°×1°網(wǎng)格內(nèi)大地水準(zhǔn)面數(shù)據(jù)和地形數(shù)據(jù)的平均值作為該網(wǎng)格的反演輸入數(shù)據(jù),標(biāo)準(zhǔn)差為觀測數(shù)據(jù)的誤差,為了兼顧實(shí)際觀測數(shù)據(jù)與理論建模數(shù)據(jù)在1°×1°網(wǎng)格內(nèi)的誤差,我們將大地水準(zhǔn)面和高程的最小不確定度分別設(shè)置為1.2 m和150 m(Afonso et al., 2013a,b).

1.2 巖石圈溫度、化學(xué)組分結(jié)構(gòu)的熱動(dòng)力學(xué)正演模擬

為了方便開展研究,我們首先將川滇地區(qū)劃分為多個(gè)1°×1°的單元網(wǎng)格,針對每一個(gè)1°×1°小柱體進(jìn)行一維巖石圈溫度、化學(xué)組分結(jié)構(gòu)的概率密度反演,然后將反演得到的各單元網(wǎng)格模型參數(shù)的后驗(yàn)概率密度分布作為先驗(yàn)信息代入三維巖石圈溫度、巖石圈地幔化學(xué)組分結(jié)構(gòu)反演,進(jìn)一步修正模型參數(shù),最終得到可靠的、高分辨率的川滇地區(qū)三維巖石圈溫度、巖石圈地?；瘜W(xué)組分結(jié)構(gòu).對于一維柱體巖石圈溫度結(jié)構(gòu)反演,我們針對選取的每一個(gè)1°×1°單元的巖石圈、上地幔進(jìn)行模型參數(shù)化設(shè)置.地殼設(shè)置為一層沉積層,兩層結(jié)晶層,每一層由體密度(ρ)、波速比(VP/VS)、生熱率(RHP)、厚度變化范圍(Δh)、熱導(dǎo)率(κ)、熱膨脹系數(shù)(α)、壓縮系數(shù)(β)等物理參數(shù)進(jìn)行定義,前四個(gè)參數(shù)通過反演確定(見表1),后三個(gè)參數(shù)作為反演的輸入數(shù)據(jù),Moho深度數(shù)據(jù)來源于Xu等(2020)基于接收函數(shù)約束的1°×1°反演結(jié)果.地幔設(shè)置為巖石圈地幔和軟流圈地幔(底界面為410 km地幔轉(zhuǎn)換帶的頂部)兩部分,每一層由CFMAS組分系統(tǒng)(CaO-FeO-MgO-Al2O3-SiO2)定義,巖石圈厚度(LAB)作為單獨(dú)變量定義一維巖石圈溫度結(jié)構(gòu)(見表1的先驗(yàn)參數(shù)范圍),lith和sublith上標(biāo)分別表示巖石圈地幔和軟流圈地幔主要氧化物,T-buffer、T-inter、T-Bottom作為隨機(jī)變量定義一維軟流圈溫度結(jié)構(gòu).

表1 模型參數(shù)和先驗(yàn)信息

聯(lián)合巖石學(xué)、地球物理、地球化學(xué)等最新的研究成果(Khan et al., 2007, 2013; Afonso et al., 2008, 2013a,b; Cammarano et al., 2011),Afonso等(2008, 2013a,b)發(fā)展了基于熱動(dòng)力學(xué)模擬研究巖石圈溫度、化學(xué)組分結(jié)構(gòu)的新方法,本方法涉及的主要正演問題如下:(1)由巖石圈溫度、化學(xué)組分模型各參數(shù)的先驗(yàn)信息(見表1)隨機(jī)生成一個(gè)理論模型,該理論模型主要參數(shù)包括LAB深度、不同深度各熱動(dòng)力學(xué)節(jié)點(diǎn)上CaO-FeO-MgO-Al2O3-SiO2等主要化學(xué)氧化物的百分比含量;(2)以模型給出的巖石圈厚度為底部邊界(本文中巖石圈底部為1250 ℃等溫面)求解一維穩(wěn)態(tài)熱傳導(dǎo)方程,得到巖石圈一維溫度分布,計(jì)算理論地表熱流觀測;(3)根據(jù)得到的溫度分布、隨機(jī)生成的模型組分結(jié)構(gòu)以及模型的初始?jí)毫Ψ植记蠼飧鳠釀?dòng)力學(xué)節(jié)點(diǎn)上特定溫度、壓力、化學(xué)組分條件下的吉布斯自由能最小化問題,得到熱平衡狀態(tài)下橄欖巖中各種礦物(石榴石、尖晶石、斜方輝石、單斜輝石、橄欖石等)的含量,并根據(jù)巖石物理性質(zhì)實(shí)驗(yàn)結(jié)果給出各熱動(dòng)力學(xué)節(jié)點(diǎn)上巖石的密度、波速等物理學(xué)特性;(4)由巖石圈一維密度分布求解壓力分布,代入步驟3修正模型的初始?jí)毫Ψ植?重新計(jì)算熱動(dòng)力學(xué)平衡方程并得到新的礦物組成和物理學(xué)特性,通過4～5次迭代計(jì)算得到穩(wěn)定的壓力和密度分布,在此基礎(chǔ)上計(jì)算理論一維海拔和大地水準(zhǔn)面高;(5)考慮溫度-壓力條件下地震波波速的衰減以及可能存在的部分熔融作用后,根據(jù)巖石的波速特性計(jì)算巖石圈的一維波速結(jié)構(gòu),并計(jì)算理論地震面波頻散曲線. 熱動(dòng)力學(xué)正演模擬框架如圖3的流程圖所示.

圖3 熱動(dòng)力學(xué)正演模擬流程圖

由于不同的地球物理觀測數(shù)據(jù)對巖石圈淺部/深部、溫度/化學(xué)組分結(jié)構(gòu)存在不同的靈敏性和解析度,因此當(dāng)隨機(jī)生成的巖石圈溫度、化學(xué)組分模型能夠在一定程度上擬合所有的地球物理觀測數(shù)據(jù)時(shí),該模型可以被認(rèn)為在當(dāng)前地球物理數(shù)據(jù)約束下最接近真實(shí)巖石圈結(jié)構(gòu)的模型.目前,該方法廣泛應(yīng)用于華南地區(qū)(Shan et al., 2014, 2021)、華北地區(qū)(Guo et al., 2016)、蒙古—貝加爾地區(qū)(Fullea et al., 2012)、愛爾蘭島(Jones et al., 2014; Fullea et al., 2014)等復(fù)雜構(gòu)造單元巖石圈結(jié)構(gòu)的研究,并獲得了諸多重要的科學(xué)認(rèn)識(shí).

1.3 概率密度反演

聯(lián)合多種地球物理觀測資料反演巖石圈結(jié)構(gòu)是高度非線性的問題、地球物理觀測資料不可避免地存在觀測誤差、所采用的物理理論(如均衡等)較真實(shí)情況進(jìn)行了簡化、且反演的解可能存在非唯一性(Afonso et al., 2013a; Shan et al., 2014),因此研究采用概率密度反演的方法,基于熱動(dòng)力學(xué)模擬將不同的地球物理觀測融合在一個(gè)體系中.我們采用馬爾可夫鏈蒙特卡洛采樣方法(MCMC)對巖石圈溫度、化學(xué)組分結(jié)構(gòu)的參數(shù)空間進(jìn)行采樣,該方法在每一個(gè)1°×1°的單元內(nèi)共使用250000次抽樣模擬,在非自適應(yīng)階段抽取一半樣本,隨后每25000次抽樣模擬更新一次模型參數(shù)空間的分布,結(jié)果展示反演收斂后的最后3000個(gè)模型.

由于反演中用到的多種地球物理觀測資料對巖石圈不同深度的溫度、化學(xué)組分結(jié)構(gòu)存在不同的靈敏度,不同觀測數(shù)據(jù)的誤差和數(shù)據(jù)點(diǎn)(例如:面波頻散曲線的周期采樣點(diǎn))也不盡相同,因此為了合理地平衡不同數(shù)據(jù)集的權(quán)重,需要進(jìn)行測試以選擇適當(dāng)?shù)哪繕?biāo)函數(shù).本文定義的目標(biāo)函數(shù)如下所示:

2 計(jì)算結(jié)果

2.1 反演結(jié)果可靠性測試

為了測試反演采用的地球物理觀測(地形、大地水準(zhǔn)面高、地表大地?zé)崃鳌?～150 s瑞雷面波相速度頻散曲線)是否能夠有效地約束巖石圈溫度、化學(xué)組分結(jié)構(gòu),本文構(gòu)建四種具有不同溫度和化學(xué)組分結(jié)構(gòu)的巖石圈模型,分別對應(yīng)著薄/厚和飽滿/難熔的巖石圈(如表2所示).其中,巖石圈地?；瘜W(xué)組分結(jié)構(gòu)用橄欖巖的鎂值(Mg#)定義.鎂值(Mg#=100×n(MgO)/(n(MgO)+n(FeO))通常被認(rèn)為是地幔飽滿或難熔的重要指標(biāo),難熔型地幔通過熔體抽取或變質(zhì)作用導(dǎo)致原始橄欖巖化學(xué)組分發(fā)生較大改變,熔融程度較高;飽滿型地?；瘜W(xué)組分接近原始橄欖巖,熔融程度較低.由于MgO是地幔橄欖巖中熔點(diǎn)最高的主要金屬氧化物,所以Mg#越高表明熔體抽取程度越高,通常對應(yīng)更古老的巖石圈地幔.大量研究認(rèn)為Mg#大于92為難熔型地幔,小于90為飽滿型地幔,而90～92之間為過渡性地幔(鄭建平, 1999; Afonso et al., 2013a,b).

表2 可靠性測試中的初始(真實(shí))和反演得到的(平均±標(biāo)準(zhǔn)差)模型參數(shù)

基于表2給出的四種代表性巖石圈結(jié)構(gòu)模型,通過正演獲取相應(yīng)的地球物理觀測值,并賦以真實(shí)的觀測誤差,再以正演得到的地表地形、大地水準(zhǔn)面高、地表大地?zé)崃骱兔娌l散曲線等觀測資料作為約束,反演巖石圈的LAB深度和巖石圈地幔的Mg#.如圖4所示,黑色線段為給定巖石圈結(jié)構(gòu)模型的LAB深度與Mg#值,灰色柱狀圖為地球物理觀測數(shù)據(jù)(模型正演得到)反演得到的LAB深度與巖石圈地幔Mg#的后驗(yàn)概率密度分布.盡管在反演開始時(shí)模型參數(shù)中LAB深度(50～330 km)、巖石圈地幔Mg#(86～94)的初始先驗(yàn)值取值范圍較大,為自然界中該參數(shù)可能取值的最大范圍,但聯(lián)合多種地球物理數(shù)據(jù)反演得到的LAB深度、巖石圈地幔Mg#的后驗(yàn)概率密度分布很好地恢復(fù)了初始模型設(shè)定的LAB深度、巖石圈地幔Mg#.一方面說明反演對模型參數(shù)的先驗(yàn)信息要求不高,同時(shí)表明本文反演所采用的地球物理觀測數(shù)據(jù)及其質(zhì)量能夠有效地約束研究區(qū)域的巖石圈溫度、化學(xué)組分結(jié)構(gòu).對于巖石圈厚度恢復(fù)情況,我們發(fā)現(xiàn)相較于較厚的巖石圈模型(模型3、4),本文所采用的地球物理觀測數(shù)據(jù)對較薄的巖石圈模型(模型1、2)LAB深度約束更好,這是因?yàn)楸疚姆囱莸玫降膸r石圈溫度結(jié)構(gòu)主要由面波數(shù)據(jù)控制(Shan et al., 2014),巖石圈速度結(jié)構(gòu)對溫度結(jié)構(gòu)最為敏感,而面波頻散數(shù)據(jù)(8～150 s)對速度結(jié)構(gòu)的分辨率隨著深度增大逐漸降低,因此對于較深的LAB分辨率較低.先前的研究表明,在反演地?；瘜W(xué)組分結(jié)構(gòu)時(shí),需要聯(lián)合對溫度敏感的地震觀測與對密度敏感的非地震觀測(地形、大地水準(zhǔn)面)(Afonso et al., 2008, 2013a),基于密度或者溫度結(jié)構(gòu)的單方面約束得到的化學(xué)組分結(jié)果誤差較大(Shan et al., 2014, 2021).本文反演得到的巖石圈地幔Mg#較LAB深度的反演結(jié)果存在更大不確定度,這主要是由于反演使用的地球物理觀測可能對化學(xué)組分結(jié)構(gòu)約束較弱,且化學(xué)組分結(jié)構(gòu)反演的多解性較強(qiáng)(Afonso et al., 2013a).但是,可以看到對于給定化學(xué)組分差異較大的巖石圈模型結(jié)構(gòu),反演得到的巖石圈地幔Mg#相對變化可以在一定程度上反映巖石圈地幔的熔融程度.因此,本文所使用的多種地球物理觀測資料能夠有效地獲取川滇地區(qū)巖石圈溫度、化學(xué)組分結(jié)構(gòu).

圖4 觀測數(shù)據(jù)對不同的巖石圈結(jié)構(gòu)模型恢復(fù)程度

川滇地區(qū)地表熱流資料較為稀缺且分布極不均勻,本文收集了川滇地區(qū)及其鄰區(qū)所有地表熱流數(shù)據(jù)并插值到川滇地區(qū),在沒有實(shí)際地表熱流的地區(qū)(如松潘甘孜塊體、江南造山帶和右江盆地)插值得到的地表熱流值不確定度較高,可能會(huì)導(dǎo)致這些地區(qū)反演的巖石圈溫度結(jié)構(gòu)不可靠.本文選取地表熱流分布較好的網(wǎng)格點(diǎn)(位于川滇菱形塊體),對地表熱流數(shù)據(jù)參與反演和地表熱流不參與反演得到的巖石圈厚度結(jié)果進(jìn)行對比,圖5a顯示地表熱流在反演過程中能夠?qū)r石圈厚度(溫度)結(jié)構(gòu)提供一定的約束,但地表熱流數(shù)據(jù)不參與反演約束的巖石圈溫度結(jié)構(gòu)也比較可靠.此外,圖5b顯示地表熱流不參與反演得到的理論地表熱流值與實(shí)際觀測的地表熱流值具有較好的一致性,表明本文反演得到的理論地表熱流數(shù)據(jù)能夠?yàn)闆]有地表熱流觀測的地區(qū)提供一定的參考.

圖5 地表熱流不參與和參與反演對巖石圈溫度結(jié)構(gòu)的影響(a) 地表熱流不參與反演和其參與反演得到的巖石圈厚度模型概率密度分布對比圖; (b) 地表熱流不參與反演和其參與反演得到的理論地表熱流概率密度分布對比圖,紅色線段和黑色虛線分別為觀測值及標(biāo)準(zhǔn)差.

2.2 數(shù)據(jù)擬合程度

如前所述,本文將研究區(qū)域離散為1°×1°的柱體,反演的擬合誤差如圖6b所示.總體趨勢上,松潘甘孜塊體反演誤差大于其他塊體,我們深入分析位于松潘甘孜塊體內(nèi)部誤差相對較大的9號(hào)格點(diǎn)(圖6b中紫色強(qiáng)調(diào)區(qū)域).圖6e為瑞雷面波相速度頻散曲線的擬合程度,灰色曲線代表反演穩(wěn)定后MCMC生成的最后3000個(gè)模型頻散曲線的數(shù)據(jù)集合,黑色誤差棒為不同周期相速度的觀測值及標(biāo)準(zhǔn)差,可以看到模型給出的面波頻散曲線集合與實(shí)際觀測擬合程度較好.圖6f—h中的藍(lán)色柱狀圖分別表示9號(hào)格點(diǎn)反演得到的地表地形、大地水準(zhǔn)面高和地表大地?zé)崃鞯饶Ｐ陀^測數(shù)據(jù)的后驗(yàn)概率密度分布,紅色線段與黑色虛線表示實(shí)際觀測及標(biāo)準(zhǔn)差.如圖6f、6g所示,模型給出的地形、大地水準(zhǔn)面與實(shí)際觀測資料間的擬合程度較高.而圖6h中模型給出的地表大地?zé)崃骱篁?yàn)概率密度分布(68±15 mW·m-2)與插值給出的地表熱流觀測值(78 mW·m-2)存在一定偏差,可能是由于位于松潘甘孜塊體的9號(hào)點(diǎn)缺乏實(shí)際熱流值觀測,僅依靠數(shù)學(xué)方法插值得到的熱流值可靠性較差,而反演輸出的理論地表熱流值能夠?yàn)樵摰貐^(qū)地表熱流值提供一定的參考依據(jù).雖然位于松潘甘孜塊體內(nèi)的9號(hào)網(wǎng)格點(diǎn)在整個(gè)研究區(qū)域中擬合程度相對較低,但反演給出的模型觀測值與實(shí)際觀測資料擬合程度仍然較好.

圖6 模型觀測與實(shí)際觀測的擬合程度(a) 研究區(qū)域網(wǎng)格劃分、地表熱流及地震臺(tái)站分布.灰色網(wǎng)格線為反演劃分的1°×1°單元的邊界,彩色圓點(diǎn)為地表熱流分布,黑色三角形為地震臺(tái)站分布; (b) 正則化后的反演誤差.數(shù)字代表反演過程中1°×1°的網(wǎng)格點(diǎn)編號(hào); (c) 反演得到的LAB深度概率密度分布,9號(hào)網(wǎng)格點(diǎn)LAB深度平均值為230 km; (d) 反演得到的巖石圈地幔鎂值概率密度分布,9號(hào)網(wǎng)格點(diǎn)巖石圈地幔鎂值平均值為90; (e) 面波頻散曲線的擬合程度,紅色線段和黑色誤差棒分別為不同周期的觀測值及誤差,灰色集合為模型給出的頻散曲線集合; (f—h) 絕對高程、大地水準(zhǔn)面高和地表大地?zé)崃鞯臄M合程度,紅色線段和黑色虛線分別為觀測值及標(biāo)準(zhǔn)差,藍(lán)色直方圖為模型觀測值的后驗(yàn)概率密度分布.

2.3 巖石圈厚度分布

巖石圈由地殼和一部分上地幔物質(zhì)組成,是地球表層黏度較大的外殼,在地質(zhì)演化歷史上能夠支撐并傳輸較大的差異應(yīng)力(徐義賢等, 2019).根據(jù)巖石圈的熱定義,巖石圈被認(rèn)為是上地幔熱邊界層以上的部分,巖石圈底界面到地表是一個(gè)相對不對流、溫度梯度較高、內(nèi)部熱傳輸方式主要為熱傳導(dǎo)機(jī)制的區(qū)域,巖石圈在溫度較高、對流作用強(qiáng)烈、黏滯系數(shù)較小的軟流圈上方移動(dòng)(Jordan, 1988),軟流圈內(nèi)部熱傳輸方式以對流為主(Sleep, 2005).小尺度地幔對流的數(shù)值模擬研究結(jié)果(Zlotnik et al., 2008)表明溫度高于1250 ℃時(shí)地幔熱傳遞方式以對流為主,因此本研究中LAB定義為1250 ℃等溫面.位于松潘甘孜塊體9號(hào)網(wǎng)格點(diǎn)的LAB深度如圖6c所示,其中灰色直方圖為反演得到的LAB深度的后驗(yàn)概率密度分布,結(jié)果顯示該區(qū)域LAB深度為230±20 km.相較于松潘甘孜塊體,川滇地區(qū)南部巖石圈厚度較薄(如印支塊體70±10 km).由于反演使用的所有數(shù)據(jù)分辨率都隨深度的增加而降低,因此本文反演對深部結(jié)構(gòu)約束較弱,導(dǎo)致松潘甘孜塊體巖石圈厚度的不確定度高于川滇南部地區(qū).

如圖7、8所示,川滇地區(qū)巖石圈結(jié)構(gòu)復(fù)雜,南北地區(qū)差異顯著.印支塊體巖石圈較為平坦,厚度為70±10 km,而西北部的松潘甘孜塊體巖石圈較厚,厚度為220±20 km.華南塊體西緣巖石圈的起伏較大,其中川滇菱形塊體巖石圈厚度自80 km北向遞增至200 km,標(biāo)準(zhǔn)差自<10 km也呈北向遞增至～20 km.江南造山帶西緣巖石圈厚度自～100 km北向遞增至～220 km,標(biāo)準(zhǔn)差～15 km.比較其他研究者用不同數(shù)據(jù)和方法獲得的結(jié)果,Yang等(2017)利用P波與S波接收函數(shù)進(jìn)行LAB界面成像,得到的印支塊體LAB深度與本文結(jié)果較為一致,為60～80 km;而在松潘甘孜塊體接收函數(shù)得到的LAB深度整體偏薄,～140 km,在華南塊體西緣得到的LAB深度在約27°N時(shí)自80 km北向遞增至140 km,巖石圈厚度的變化趨勢與本文研究的結(jié)果大體一致,但反演得到的巖石圈厚度較本文結(jié)果偏薄.接收函數(shù)反演得到的LAB深度與本文結(jié)果存在差異的原因可能在于川滇地區(qū)存在巖石圈內(nèi)部不連續(xù)界面,接收函數(shù)檢測到的界面可能是該地區(qū)巖石圈內(nèi)部的不連續(xù)界面(Yuan and Romanowicz, 2010; 徐義賢等, 2019).Priestley和McKenzie(2013)根據(jù)S波波速與溫度的經(jīng)驗(yàn)轉(zhuǎn)換關(guān)系給出的川滇地區(qū)LAB深度(如圖7b所示)與本文研究結(jié)果相近,但在華南塊體西緣由南往北巖石圈遞增趨勢不明顯.導(dǎo)致上述差異的原因,一方面本文反演所使用地震面波頻散曲線分辨率較高;另一方面,本文反演聯(lián)合多種地球物理觀測資料,基于熱動(dòng)力學(xué)平衡方程直接反演巖石圈溫度結(jié)構(gòu),降低了波速-溫度經(jīng)驗(yàn)轉(zhuǎn)換關(guān)系造成的誤差.綜上,本文給出的川滇地區(qū)巖石圈厚度模型能夠同時(shí)擬合多種地球物理觀測數(shù)據(jù),并具有較高的分辨率,反演得到的巖石圈厚度分布與地表構(gòu)造、斷層分布等存在較好的相關(guān)性,因此該巖石圈模型更能反映真實(shí)的川滇地區(qū)巖石圈厚度.

圖8 川滇地區(qū)巖石圈溫度結(jié)構(gòu)剖面圖(a) 沿25°N的剖面圖; (b) 沿27°N的剖面圖; (c) 沿102°E的剖面圖; (d) 沿104°E的剖面圖.虛線示意1250 ℃等溫面,即LAB深度分布.

地球化學(xué)研究(Shang et al., 2022)估算的馬關(guān)地區(qū)新生代地幔橄欖巖來源深度為55～70 km,平衡溫度為1016～1148 ℃.而本文利用地球物理資料聯(lián)合反演得到的馬關(guān)地區(qū)巖石圈(底界面溫度為1250 ℃)厚度為60～80 km,導(dǎo)致二者差異的原因可能包括:(1)幔源捕虜體在出露于地表的過程中不可避免地會(huì)遭受不同程度的損壞,這可能會(huì)影響巖石圈溫壓結(jié)果的估計(jì).(2)聯(lián)合反演是高度非線性的問題、所采用的物理理論(如均衡等)較真實(shí)情況進(jìn)行了簡化等,因此與真實(shí)巖石圈存在一定的出入.但在馬關(guān)地區(qū)本文研究與地球化學(xué)研究約束的巖石圈橄欖巖溫壓分布整體趨勢并不矛盾,進(jìn)一步說明本文單方面利用地球物理資料約束的巖石圈溫度結(jié)構(gòu)較為可靠.

雖然本文研究得到的川滇地區(qū)巖石圈溫度結(jié)構(gòu)與其他研究均顯示川滇地區(qū)南部存在巖石圈減薄的現(xiàn)象,但是無法為該地區(qū)巖石圈減薄的動(dòng)力學(xué)機(jī)制提供更多的信息.目前,大多數(shù)研究主張的巖石圈減薄機(jī)制主要集中為交代和拆沉兩種模式,其分別對應(yīng)巖石圈減薄的化學(xué)和物理兩類不同的機(jī)制.交代是指難熔型巖石圈地幔通過巖石圈-軟流圈之間物質(zhì)能量的相互交換作用轉(zhuǎn)變?yōu)轱枬M型巖石圈地幔,這種作用使得巖石圈在減薄的基礎(chǔ)上,化學(xué)組分結(jié)構(gòu)也發(fā)生了改變(吳福元等, 2003; Shan et al., 2021).拆沉是指因板塊擠壓導(dǎo)致巖石圈縮短增厚、等溫面下移,深部巖石圈根由于溫度較低導(dǎo)致其密度高于周邊地幔,深部巖石圈地幔重力失穩(wěn)垮塌到下伏軟流圈地幔中,殘留的巖石圈地幔仍保留了大部分古老難熔的物質(zhì)(Chung et al., 1997; Chen et al., 2017).Feng等(2022)結(jié)合地震接收函數(shù)與地球化學(xué)的研究結(jié)果,認(rèn)為印度板塊與亞歐板塊南北向碰撞導(dǎo)致川滇地區(qū)南部巖石圈增厚,增厚的巖石圈由于重力不穩(wěn)發(fā)生垮塌,造成川滇地區(qū)南部巖石圈發(fā)生減薄.Lei和Zhao(2016)認(rèn)為印度板塊向東的平俯沖作用導(dǎo)致川滇地區(qū)南部巖石圈發(fā)生大規(guī)模減薄,而后印度板塊的轉(zhuǎn)返過程導(dǎo)致軟流圈物質(zhì)上涌,進(jìn)而對殘留的巖石圈進(jìn)行了有效的熱侵蝕,改變了古老巖石圈地幔的化學(xué)組成.這些爭議使得川滇地區(qū)南部巖石圈減薄的動(dòng)力學(xué)機(jī)制更加模糊,要澄清該問題需要結(jié)合化學(xué)組分結(jié)構(gòu)的證據(jù).

2.4 化學(xué)組分結(jié)構(gòu)

如前所述,前人在川滇地區(qū)開展了大量巖石圈波速、溫度/厚度的成像工作,均發(fā)現(xiàn)川滇地區(qū)南部巖石圈發(fā)生了減薄,然而目前關(guān)于該地區(qū)巖石圈減薄的機(jī)制仍存在較大爭議,原因在于巖石圈波速、溫度/厚度成像研究無法提供該地區(qū)巖石圈在減薄的同時(shí)化學(xué)組分是否也發(fā)生改變的證據(jù).地球化學(xué)研究能夠直接提供來自巖石圈深部物質(zhì)的化學(xué)組分信息,然而川滇地區(qū)只有馬關(guān)出露新生代幔源捕虜體(Liu et al., 2013),難以得到對該地區(qū)巖石圈全面、可靠的化學(xué)組分結(jié)構(gòu)的認(rèn)識(shí).本文利用空間覆蓋更好的地球物理觀測能夠反演獲得川滇地區(qū)巖石圈地幔化學(xué)組分結(jié)構(gòu)信息,為研究該地區(qū)的動(dòng)力學(xué)過程提供更多的、獨(dú)立的參考依據(jù).

由于地殼化學(xué)組分結(jié)構(gòu)十分復(fù)雜,本研究只討論組分以橄欖巖為主、相對較為簡單的上地幔組分結(jié)構(gòu).橄欖巖是上地幔最重要的組成部分,橄欖巖的化學(xué)組分影響著人類對于巖石圈起源、演化等過程的認(rèn)識(shí),能夠提供豐富的地球動(dòng)力學(xué)信息.本文采用多觀測概率密度反演得到了川滇地區(qū)巖石圈地幔、軟流圈地幔Mg#的后驗(yàn)概率密度分布,其中反演得到的軟流圈地幔橄欖巖Mg#均低于90,為接近于原始地幔橄欖巖組分的飽滿型橄欖巖,符合之前學(xué)者們對軟流圈地?；瘜W(xué)組分結(jié)構(gòu)的認(rèn)識(shí)(鄭建平, 1999; Afonso et al., 2013a,b),因此本文不討論反演得到的軟流圈地幔化學(xué)組分結(jié)果,主要分析巖石圈地幔的化學(xué)組分結(jié)構(gòu).

本文得到的巖石圈地?；瘜W(xué)組分結(jié)構(gòu)單純來自地球物理觀測數(shù)據(jù)的約束,為了檢測由地球物理觀測約束的巖石圈地?；瘜W(xué)組分結(jié)構(gòu)是否可靠,我們將馬關(guān)地區(qū)出露的幔源捕虜體樣本Mg#分布與本文得到的馬關(guān)地區(qū)巖石圈地幔Mg#分布進(jìn)行對比(如圖9所示).結(jié)果表明,由地球化學(xué)、巖石學(xué)觀測得到的馬關(guān)地區(qū)巖石圈地幔橄欖巖鎂值分布(Liu et al., 2013)與地球物理觀測反演得到的結(jié)果較為吻合,均以飽滿型地幔為主,表明利用地球物理觀測資料能夠在一定程度上約束巖石圈的化學(xué)組分結(jié)構(gòu),為缺乏地球化學(xué)觀測的地區(qū)提供獨(dú)立、重要的化學(xué)組分信息.

圖9 馬關(guān)地區(qū)新生代幔源橄欖巖Mg#分布(紅色柱體)與反演得到的該地區(qū)巖石圈橄欖巖Mg#概率密度分布(灰色柱體)對比圖

圖10為反演得到的川滇地區(qū)巖石圈地幔Mg#及其標(biāo)準(zhǔn)差分布,結(jié)果表明川滇地區(qū)不同塊體之間巖石圈地幔Mg#存在顯著差異.其中,川滇地區(qū)南部反演得到的巖石圈地幔Mg#較其他區(qū)域存在較大的標(biāo)準(zhǔn)差.為了進(jìn)一步顯示不同塊體間巖石圈地幔的組分差異,我們在位于川滇地區(qū)南部的印支塊體、川滇菱形塊體和江南造山帶西緣各選取兩處反演網(wǎng)格點(diǎn),三個(gè)區(qū)域共六處網(wǎng)格點(diǎn)巖石圈地幔Mg#的后驗(yàn)概率密度分布如圖11所示.結(jié)果表明,印支塊體巖石圈地幔以過渡型與難熔型為主,同時(shí)該地區(qū)LAB深度較淺并且水平向起伏不明顯.Feng等(2022)開展的接收函數(shù)界面研究與Chen等(2017)開展的巖石學(xué)、地球化學(xué)研究認(rèn)為:印支塊體約30 Ma發(fā)生過巖石圈增厚過程,之后由于巖石圈重力不穩(wěn)發(fā)生拆沉作用.結(jié)合本文反演的巖石圈地?；瘜W(xué)組分結(jié)構(gòu)和前人的研究,我們認(rèn)為拆沉作用在印支塊體巖石圈減薄過程中占主導(dǎo)位置.反演得到的川滇菱形塊體巖石圈地幔以飽滿型與過渡型為主,此外,川滇菱形塊體南部巖石圈地幔的Mg#小于塊體北部地區(qū),同時(shí)由于川滇菱形塊體巖石圈自南向北逐漸增厚,我們推測川滇菱形塊體巖石圈破壞過程始于南部地區(qū),并逐漸向北推進(jìn),由南往北不均勻的熱侵蝕和改造作用是導(dǎo)致川滇菱形塊體巖石圈減薄的主要?jiǎng)恿W(xué)機(jī)制.反演得到的華南塊體西緣(如江南造山帶西緣)巖石圈物質(zhì)以過渡型與難熔型為主,而Shan等(2021)反演得到的華南塊體東緣巖石圈地幔以飽滿型為主.我們認(rèn)為華南塊體東西兩側(cè)組分變化的原因在于巖石圈受到的力學(xué)機(jī)制不同.華南塊體東側(cè)受到太平洋板塊的平俯沖作用,巖石圈地幔表現(xiàn)為飽滿型地幔;華南塊體西側(cè)受到印支塊體拆沉作用導(dǎo)致的地幔對流的影響,巖石圈地幔組分也發(fā)生了變化,以飽滿型橄欖巖為主,如馬關(guān)地區(qū)等.但馬關(guān)地區(qū)巖石圈地幔的組分結(jié)構(gòu)并不能完全代表整個(gè)華南塊體西緣,我們反演結(jié)果顯示,江南造山帶西緣其主體的巖石圈地?；瘜W(xué)組分仍然以過渡和難熔型地幔為主.巖石圈地?；瘜W(xué)組分結(jié)構(gòu)結(jié)果表明華南塊體東西兩側(cè)巖石圈減薄過程所受到的動(dòng)力機(jī)制不一樣,不同地區(qū)巖石圈地幔化學(xué)組分之間存在較大差異.

圖10 川滇地區(qū)巖石圈地幔Mg#分布圖(a) 反演得到的巖石圈地幔鎂值平均值分布.紅色正方形為出露幔源捕虜體樣本的馬關(guān)地區(qū),紫色正方形為選取的巖石圈地幔鎂值分布的點(diǎn); (b) 反演得到的巖石圈地幔鎂值分布的標(biāo)準(zhǔn)差.

圖11 川滇地區(qū)選取的六個(gè)點(diǎn)的巖石圈地幔Mg#后驗(yàn)概率密度分布(a) 印支塊體選取的兩個(gè)點(diǎn); (b) 川滇菱形塊體選取的兩個(gè)點(diǎn); (c) 江南造山帶西緣選取的兩個(gè)點(diǎn).

本文反演得到的四川盆地巖石圈熱厚度(約220 km)與近期接收函數(shù)(張耀陽等, 2018)得到的巖石圈底部速度不連續(xù)面(約190 km深度)存在一定的差異,導(dǎo)致差異的原因一方面可能是二者巖石圈定義方式不同,另一方面也可能是各自研究中巖石圈厚度不確定度存在差異所導(dǎo)致.此外,本文反演得到的四川盆地巖石圈地幔為飽滿型橄欖巖,這與傳統(tǒng)認(rèn)知的難熔型克拉通巖石圈地幔并不一致.需要指出的是,本文無法提供巖石圈地幔垂向組分分層的可靠信息,而是用整個(gè)巖石圈地幔平均深度的Mg#結(jié)果來代表巖石圈地幔組分結(jié)構(gòu),這樣給出的巖石圈地幔Mg#結(jié)果能夠指示不同塊體化學(xué)組分結(jié)構(gòu)的態(tài)勢,但無法提供具體深度處的組分結(jié)構(gòu).本文反演得到的四川盆地巖石圈地幔為飽滿型橄欖巖,我們推測四川盆地巖石圈地?？赡芙?jīng)歷了重富集作用,由峨眉山大火成巖省引起的鎂鐵質(zhì)巖漿作用導(dǎo)致四川盆地巖石圈下方熔體滲流進(jìn)入原有的難熔型橄欖巖地幔,古老難熔的巖石圈地幔物質(zhì)被新生飽滿的物質(zhì)所替換(Yang et al., 2021),這與最近有關(guān)四川盆地內(nèi)部廣泛的巖漿活動(dòng)報(bào)道一致(Lu et al., 2019; Liu and Qiu, 2020).

2.5 不同觀測資料對反演結(jié)果的敏感性測試

本文利用地表熱流、地形、大地水準(zhǔn)面、瑞雷面波相速度頻散曲線等多種觀測數(shù)據(jù)聯(lián)合反演巖石圈溫度、化學(xué)組分結(jié)構(gòu),由于各種觀測數(shù)據(jù)對巖石圈不同深度的溫度、化學(xué)組分結(jié)構(gòu)存在敏感性差異,聯(lián)合多種觀測資料最終可能導(dǎo)致反演結(jié)果存在爭議.為了認(rèn)識(shí)爭議產(chǎn)生的來源,我們以松潘甘孜塊體為例(圖6b所示的網(wǎng)格點(diǎn)),開展了不同觀測資料對巖石圈溫度、化學(xué)組分結(jié)構(gòu)的敏感性測試.圖12顯示單獨(dú)以地表熱流、地形、大地水準(zhǔn)面作為觀測約束的巖石圈厚度存在較大的標(biāo)準(zhǔn)差,這主要是由于模型中的熱參數(shù)(如地殼各層的厚度、熱產(chǎn)率、熱傳導(dǎo)系數(shù)等)不確定度較大所導(dǎo)致,單獨(dú)以面波頻散曲線(8～150 s周期)作為觀測能夠較好地約束巖石圈溫度結(jié)構(gòu).然而,單一的地球物理觀測資料難以有效約束巖石圈化學(xué)組分結(jié)構(gòu),因?yàn)閱畏矫娴拿芏?控制地形、大地水準(zhǔn)面)或波速(控制面波頻散曲線)結(jié)構(gòu)作為約束得到的化學(xué)組分結(jié)果誤差較大,而聯(lián)合多種觀測資料是獲取巖石圈化學(xué)組分結(jié)構(gòu)的有效途徑.綜上,我們測試發(fā)現(xiàn)不同的觀測數(shù)據(jù)對巖石圈溫度、化學(xué)組分結(jié)構(gòu)存在不同的敏感性,但不同觀測資料反演得到的巖石圈溫度、化學(xué)組分結(jié)構(gòu)整體趨勢并不矛盾,因此聯(lián)合多種觀測資料能夠有效約束巖石圈溫度、化學(xué)組分結(jié)構(gòu).

圖12 不同觀測數(shù)據(jù)對巖石圈溫度和化學(xué)組分結(jié)構(gòu)的敏感性分析(a—c) 反演得到的LAB深度; (d—f) 反演得到的巖石圈地幔Mg#.圖中直方圖分別代表以地表熱流、高程、大地水準(zhǔn)面、面波頻散曲線、聯(lián)合多種觀測數(shù)據(jù)作為約束得到的反演結(jié)果.

3 討論

先前大量研究均認(rèn)為川滇地區(qū)南部存在巖石圈減薄的現(xiàn)象,與本文得到的巖石圈溫度結(jié)構(gòu)較為一致.但對巖石圈減薄過程的動(dòng)力學(xué)認(rèn)識(shí)仍存在兩種爭議:印度板塊與亞歐板塊的南北向碰撞(Chung et al., 1997; Chen et al., 2017; Zhang et al., 2017; Feng et al., 2022)及印度板塊的東向俯沖(Huang et al., 2015; Lei and Zhao, 2016).其中,第一種觀點(diǎn)認(rèn)為印度板塊與歐亞板塊的南北向碰撞導(dǎo)致了川滇地區(qū)巖石圈增厚,而增厚的巖石圈由于重力失穩(wěn)垮塌,使得川滇地區(qū)南部巖石圈減薄,垮塌的巖石圈形成地幔上升流,進(jìn)而對殘留的巖石圈地幔物質(zhì)進(jìn)行改造.而另一種觀點(diǎn)則認(rèn)為印度板塊的東向平俯沖作用導(dǎo)致川滇地區(qū)南部巖石圈發(fā)生大規(guī)模減薄,而后印度板塊的折返過程引起地幔上升流,進(jìn)而對殘留的巖石圈進(jìn)行了有效的熱侵蝕,完全改變了古老巖石圈地幔的化學(xué)組成.可以看到,要澄清川滇地區(qū)南部巖石圈減薄機(jī)制的關(guān)鍵在于認(rèn)識(shí)川滇地區(qū)南部各構(gòu)造單元巖石圈地幔的化學(xué)組分是否發(fā)生改變.

印支塊體地處紅河斷裂以南,受印度板塊東向俯沖和青藏高原物質(zhì)東向擠壓作用的影響,其構(gòu)造活動(dòng)較為強(qiáng)烈,巖石圈較薄.由于缺少地球化學(xué)證據(jù),該地區(qū)巖石圈是否發(fā)生過增厚過程很難確定.但是,本文反演得到的印支塊體巖石圈溫度、化學(xué)組分結(jié)構(gòu)表明:印支塊體LAB深度為70±10 km,遠(yuǎn)小于穩(wěn)定大陸的LAB深度值(～200 km),其巖石圈地幔以過渡型和難熔型橄欖巖為主,表明印支塊體在經(jīng)歷較長的演化歷史后,巖石圈在減薄的同時(shí)仍保留部分古老難熔的巖石圈地幔物質(zhì).如果印度板塊的東向俯沖是導(dǎo)致印支塊體巖石圈減薄的主要因素,那么該區(qū)域巖石圈地幔應(yīng)該以飽滿型橄欖巖為主,然而反演得到的印支塊體巖石圈化學(xué)組分不符合上述特征,說明熱化學(xué)侵蝕作用在該區(qū)域巖石圈減薄過程中不占據(jù)主導(dǎo)地位,而巖石圈增厚拆沉作用可能是導(dǎo)致印支塊體巖石圈減薄的主要?jiǎng)恿W(xué)機(jī)制.我們認(rèn)為印支塊體巖石圈減薄過程受印度板塊與歐亞板塊南北向碰撞影響,導(dǎo)致青藏高原物質(zhì)東向流出,高原物質(zhì)東向擠壓作用于印支塊體導(dǎo)致該地區(qū)巖石圈發(fā)生增厚拆沉,拆沉作用使其巖石圈發(fā)生減薄并形成地幔上升流,導(dǎo)致印支塊體巖石圈地幔整體呈現(xiàn)高溫狀態(tài),但該上升流并未對印支塊體殘留的巖石圈地幔物質(zhì)形成有效的侵蝕和改造.地幔界面成像(Feng et al., 2022)和地震速度成像(Huang et al., 2015)發(fā)現(xiàn)印支塊體深部地幔存在速度不連續(xù)面,高速異常區(qū)域可能對應(yīng)下沉的巖石圈塊體,在深部地幔結(jié)構(gòu)層面支持印支塊體巖石圈增厚拆沉的動(dòng)力學(xué)機(jī)制.

川滇菱形塊體構(gòu)造活動(dòng)強(qiáng)烈,與其他塊體相比,該地區(qū)具有較高的地表熱流值與溫度分布,而且該地區(qū)的地表熱流、溫度結(jié)構(gòu)、LAB深度分布均與斷層的空間分布有較好的吻合.本文反演得到的川滇菱形塊體LAB深度自南(～80 km)往北(～200 km)逐漸增大,巖石圈地幔以飽滿型橄欖巖為主,而且川滇菱形塊體南部巖石圈地幔的飽滿程度高于北部,說明由南往北的熱化學(xué)侵蝕作用在該區(qū)域巖石圈減薄過程中占據(jù)主導(dǎo)地位.結(jié)合印支塊體巖石圈減薄機(jī)制及川滇菱形塊體深部上地幔的溫度結(jié)構(gòu),我們認(rèn)為由于印度板塊與亞歐板塊的碰撞作用,導(dǎo)致印支塊體巖石圈發(fā)生增厚拆沉,拆沉作用形成的地幔上升流由南往北作用于川滇菱形塊體巖石圈,導(dǎo)致川滇菱形塊體巖石圈由南往北遞增,而地幔上升流對巖石圈地幔物質(zhì)形成有效的侵蝕和改造,導(dǎo)致川滇菱形塊體巖石圈地幔以飽滿型橄欖巖為主,并在地表出露大量富鉀鎂鐵質(zhì)巖石(Chen et al., 2017).

華南塊體西緣(如江南造山帶西緣)LAB深度自南(～80 km)往北(～200 km)逐漸增大,巖石圈地幔以過渡型和難熔型橄欖巖為主,而Shan等(2021)反演得到的華南塊體東緣巖石圈地幔以飽滿型為主.華南塊體東西兩側(cè)組分結(jié)構(gòu)存在較大差異,這可能與華南塊體東西兩側(cè)巖石圈受到的力學(xué)機(jī)制不同有關(guān),華南塊體東側(cè)受到太平洋板塊的平俯沖作用,巖石圈地幔表現(xiàn)為飽滿型地幔;華南塊體西側(cè)受到印支塊體拆沉作用導(dǎo)致的地幔對流的影響,部分地區(qū)巖石圈地幔組分也發(fā)生了變化,如馬關(guān)地區(qū),但馬關(guān)地區(qū)巖石圈地幔的組分結(jié)構(gòu)并不能完全代表整個(gè)華南塊體西緣,華南塊體西緣巖石圈地?；瘜W(xué)組分結(jié)構(gòu)結(jié)果顯示,該地區(qū)主體的巖石圈地?；瘜W(xué)組分仍然以過渡型和難熔型橄欖巖為主,加之該地區(qū)巖石圈厚度由南往北逐漸增加,可能是由于印支塊體拆沉引起的地幔對流作用導(dǎo)致華南塊體西南緣(如江南造山帶西緣)巖石圈發(fā)生減薄,但地幔對流作用并未對華南塊體西側(cè)的巖石圈地幔物質(zhì)形成有效的侵蝕和改造.

綜上,我們認(rèn)為約60 Ma印度板塊對歐亞板塊北向俯沖導(dǎo)致青藏高原地殼、巖石圈開始縮短增厚(Royden et al., 1997).約45 Ma印度大陸板塊與亞歐板塊發(fā)生硬碰撞(Chung et al., 2005),加厚了印支塊體的巖石圈,并在30～40 Ma發(fā)生大規(guī)模的巖石圈拆沉作用(Feng et al., 2022),而印支塊體仍保留有古老難熔的巖石圈地幔物質(zhì).印支塊體巖石圈大規(guī)模的拆沉作用導(dǎo)致川滇菱形塊體、江南造山帶西緣巖石圈底部自南向北被加熱,深部上地幔熱結(jié)構(gòu)的突變導(dǎo)致巖石圈發(fā)生部分熔融和富鉀巖漿作用,川滇菱形塊體巖石圈地幔物質(zhì)被新生物質(zhì)所取代,而江南造山帶西緣仍保留有古老難熔的巖石圈地幔物質(zhì).

4 結(jié)論與展望

本文利用地表地形、大地水準(zhǔn)面高、地表大地?zé)崃骱偷卣鹈娌l散曲線等多種地球物理觀測資料,基于熱動(dòng)力學(xué)模擬和概率密度反演,研究了川滇地區(qū)巖石圈溫度、化學(xué)組分結(jié)構(gòu),得到以下主要結(jié)論:

(1)川滇地區(qū)巖石圈結(jié)構(gòu)復(fù)雜,不同構(gòu)造單元之間差異顯著.位于川滇地區(qū)南部的印支塊體巖石圈厚度較薄(<80 km),相對起伏較小,溫度結(jié)構(gòu)的橫向差異不明顯.川滇菱形塊體、華南塊體西緣巖石圈由南向北逐漸增加,導(dǎo)致這一現(xiàn)象的原因可能是深部地幔熱物質(zhì)自南向北侵蝕造成的.

(2)印支塊體巖石圈地幔以難熔型和過渡型的橄欖巖為主,表明仍保存有古老難熔的巖石圈地幔物質(zhì),結(jié)合該地區(qū)巖石圈較薄的特征,我們認(rèn)為新生代巖石圈增厚拆沉是印支塊體巖石圈減薄的主要?jiǎng)恿W(xué)機(jī)制.川滇菱形塊體巖石圈地幔以飽滿型為主,表明古老難熔的巖石圈地幔物質(zhì)已經(jīng)被新的飽滿型橄欖巖物質(zhì)所取代,結(jié)合該地區(qū)自南向北逐漸增加的巖石圈結(jié)構(gòu),我們認(rèn)為印支塊體拆沉作用所導(dǎo)致的地幔對流作用是川滇菱形塊體巖石圈減薄的主要?jiǎng)恿W(xué)機(jī)制,熱侵蝕作用在川滇菱形塊體巖石圈減薄過程中扮演了重要角色.華南塊體西側(cè)巖石圈地幔以過渡型與難熔型為主,結(jié)合該地區(qū)自南向北逐漸遞增的巖石圈結(jié)構(gòu)以及華南塊體東西兩側(cè)的化學(xué)組分結(jié)構(gòu),我們認(rèn)為華南塊體東西兩側(cè)經(jīng)歷的動(dòng)力作用不同,華南塊體西側(cè)巖石圈發(fā)生減薄的原因在于印支塊體拆沉導(dǎo)致的地幔對流作用,但上升的地幔物質(zhì)并未對華南塊體西側(cè)的巖石圈地幔物質(zhì)形成有效的侵蝕和改造.

本文雖然研究的是川滇地區(qū)三維巖石圈溫度、化學(xué)組分結(jié)構(gòu),但在反演時(shí)我們將川滇地區(qū)劃分為多個(gè)1°×1°的區(qū)域,實(shí)際得到的巖石圈、化學(xué)組分結(jié)構(gòu)為偽三維結(jié)構(gòu),在以后的工作中我們會(huì)嘗試開發(fā)真三維反演程序,并引入遠(yuǎn)震體波資料來更好地約束深部結(jié)構(gòu).

致謝本文采用了美國紐約州立大學(xué)石溪分校沈偉森博士以及中國科學(xué)技術(shù)大學(xué)張智奇博士發(fā)布的瑞雷面波相速度頻散曲線,審稿人提出建設(shè)性的修改意見,特此感謝.本文圖件均采用GMT軟件繪制(Wessel and Smith, 1998).