張?zhí)煜? 頡 煒*

IODP385-U1550單斜輝石指示加利福尼亞灣瓜伊馬斯盆地海底擴(kuò)張?jiān)缙贛ORB深部結(jié)晶作用

張?zhí)煜?, 2, 頡 煒1, 2*

(1. 河海大學(xué) 自然資源部海洋災(zāi)害預(yù)報(bào)技術(shù)重點(diǎn)實(shí)驗(yàn)室, 江蘇 南京 210024; 2. 河海大學(xué) 海洋學(xué)院, 江蘇 南京 210024)

洋中脊玄武巖(MORB)分離結(jié)晶深度是控制其成分變化以及洋殼增生的關(guān)鍵因素之一。為了探究海底擴(kuò)張?jiān)缙谘笾屑箮r漿分離結(jié)晶的最深深度, 本文對(duì)國(guó)際大洋發(fā)現(xiàn)計(jì)劃(IODP)385航次在加利福尼亞灣瓜伊馬斯盆地U1550站位獲得的玄武質(zhì)巖漿巖樣品進(jìn)行了詳細(xì)的巖石學(xué)和礦物學(xué)研究。U1550站位巖心主要由拉斑玄武巖、輝綠巖以及輝綠巖捕虜體組成, 其單斜輝石主要為普通輝石, 透輝石次之。單斜輝石?熔體溫壓計(jì)計(jì)算得出U1550玄武巖單斜輝石結(jié)晶溫度(1126～1170 ℃, 誤差±45 ℃)和壓力(–0.5～5.4 kbar, 誤差±1.4 kbar)相對(duì)輝綠巖及輝綠巖捕虜體單斜輝石(1158～1212 ℃, 誤差±45 ℃; 1.6～5.9 kbar, 誤差±1.4 kbar)均較低, 結(jié)晶深度(0～17.8 km)相對(duì)后者(5.3～19.5 km)較淺。與全球其他洋脊段MORB單斜輝石相比, U1550單斜輝石具有更陡峭的冷卻軌跡。綜合研究表明, 在海底擴(kuò)張?jiān)缙? 擴(kuò)張速率很可能均為慢速擴(kuò)張, 其巖漿結(jié)晶深度可超過10 km。這暗示著在海底擴(kuò)張?jiān)缙趲r石圈地幔厚度很厚, 地幔部分熔融產(chǎn)生巖漿的終止深度也相對(duì)更深, 導(dǎo)致巖漿結(jié)晶深度深, 部分巖漿難以抵達(dá)淺部洋殼。但不同于慢速?超慢速擴(kuò)張成熟洋脊段, 海底擴(kuò)張?jiān)缙谏畈繋r漿房巖漿供給可能并不匱乏, 而是這些巖漿房很可能位于更深的部位。

IODP 385; 加利福尼亞灣; 海底擴(kuò)張?jiān)缙? 單斜輝石; 巖漿深部結(jié)晶

0 引 言

洋中脊玄武巖(mid-oceanic ridge basalts, MORB)是地球內(nèi)部向地表輸送的最主要物質(zhì), 是洋殼和大洋巖石圈的重要組成部分, 在全球板塊構(gòu)造體系中扮演著十分重要的角色。大洋地幔在洋中脊下方發(fā)生減壓熔融形成原始巖漿上涌, 并在向上運(yùn)移過程中發(fā)生分離結(jié)晶、巖漿混合等演化作用最終形成MORB, 因此MORB是了解大洋地幔源區(qū)性質(zhì)和大洋巖石圈物質(zhì)組成及演化過程的絕佳研究對(duì)象(Gale et al., 2013; 張國(guó)良等, 2017; 牛耀齡, 2022)。其中, MORB分離結(jié)晶的最深深度是控制其成分變化以及洋殼增生(oceanic crust accretion)的關(guān)鍵因素之一, 代表洋中脊巖漿開始結(jié)晶的深度和地幔巖發(fā)生部分熔融生成洋中脊巖漿的終止深度(兩者往往一樣或相近), 可以限定洋殼增生的最大深度, 因?yàn)橄卵髿さ男纬膳cMORB分離結(jié)晶出的礦物堆晶作用密切相關(guān)(Coogan, 2014; White and Klein, 2014)。同時(shí), 洋中脊巖漿的分離結(jié)晶作用主要由熱傳導(dǎo)冷卻(conductive cooling)和熱損失(heat loss)驅(qū)動(dòng), 而驅(qū)動(dòng)下洋殼結(jié)晶的大部分熱損失又與熱液循環(huán)(hydrothermal circulation)密切相關(guān)。

前人根據(jù)對(duì)巖漿結(jié)晶壓力、洋脊擴(kuò)張速率和巖漿供應(yīng)量三者關(guān)系的研究指出, 控制洋中脊巖漿結(jié)晶深度的一級(jí)因素主要為洋脊擴(kuò)張速率和巖漿供應(yīng)量: 擴(kuò)張速率的快慢控制著深部熱傳導(dǎo)冷卻和熱損失, 以及提供巖漿上升通道的斷裂帶(fracture zone)的發(fā)育; 巖漿供給量控制洋殼內(nèi)部淺部巖漿房(magma chamber)和熔體透鏡體(melt lens)的發(fā)育(Michael and Cornell, 1998; Herzberg, 2004; Wanless and Behn, 2017)。在快速擴(kuò)張洋脊和巖漿供應(yīng)充足的慢速擴(kuò)張洋脊(如大西洋靠近冰島地幔柱的Reykjanes洋脊段), 淺部巖漿房和熔體透鏡體非常發(fā)育, 巖漿分離結(jié)晶作用主要發(fā)生在洋殼內(nèi)部, 巖漿記錄的結(jié)晶深度相對(duì)較淺(往往不超過深約6～7 km的莫霍面); 而在典型的慢速?超慢速擴(kuò)張洋脊及一些洋中脊斷裂帶(fracture zones)附近, 深部巖石圈具有更強(qiáng)的熱傳導(dǎo)冷卻和熱損失, 巖漿供應(yīng)量不足, 淺部巖漿房和熔體透鏡體不發(fā)育或短暫存在, 巖漿分離結(jié)晶作用不局限于下洋殼深度, 可以超過10 km, 甚至可以深達(dá)30 km(Sleep and Barth, 1997; Michael and Cornell, 1998; Herzberg, 2004; Eason and Sinton, 2006)。Bennett et al. (2019)通過斜長(zhǎng)石熔體包裹體研究也指出超慢速Gakkel洋中脊最深的巖漿結(jié)晶深度可高達(dá)16.4 km。但也有一些學(xué)者通過橄欖石熔體包裹體的研究認(rèn)為在慢速?超慢速擴(kuò)張洋脊巖漿結(jié)晶深度應(yīng)該也不超過10 km(Shaw et al., 2010; Colman et al., 2015; Wanless et al., 2015)。這些關(guān)于洋中脊巖漿結(jié)晶深度的認(rèn)識(shí)均聚焦在海底擴(kuò)張成熟期的洋脊段(比如東太平洋洋隆EPR、大西洋洋中脊MAR、西南印度洋洋中脊SWIR和北冰洋Gakkel洋脊等), 但對(duì)于尚處在海底擴(kuò)張?jiān)缙诘难蠹苟窝芯可跎?。海底擴(kuò)張?jiān)缙谘笾屑箮r漿分離結(jié)晶的最深深度是多少？是否可達(dá)巖石圈地幔深度至今尚不清楚。目前, 可以借助地球物理和地球化學(xué)方法研究洋中脊巖漿結(jié)晶深度問題。其中, 地球物理方法主要利用地震波速度變化反映巖漿產(chǎn)生與結(jié)晶的大致位置(Detrick et al., 1990; Singh et al., 2006; Combier et al., 2015), 但無法對(duì)巖漿演化過程進(jìn)行精細(xì)描述。地球化學(xué)方法包括全巖主量元素、礦物熔融包裹體和單礦物分析。MORB中的礦物簡(jiǎn)單直接, 記錄著巖漿分離結(jié)晶和巖漿混合等作用的豐富信息(Streck, 2008; 葛振敏等, 2020; 楊錦等, 2022)。其中, 單斜輝石在MORB中廣泛發(fā)育, 并且可以在地幔到地表的不同深度范圍內(nèi)結(jié)晶, 對(duì)溫度和壓力有很強(qiáng)的敏感性, 因此常常被用來計(jì)算巖漿的結(jié)晶壓力與溫度, 其攜帶的原位信息可以反演巖漿結(jié)晶過程的物理化學(xué)狀態(tài)(Villiger et al., 2007; Putirka, 2008; 羅雕等, 2022)。

大陸裂解形成的年輕擴(kuò)張盆地(如加利福尼亞灣、紅海、亞丁灣等)很可能記錄著海底擴(kuò)張?jiān)缙诘难髿そM分和洋中脊巖漿作用的關(guān)鍵信息。紅海(Red Sea)和亞丁灣(Gulf of Aden)是典型的海底擴(kuò)張?jiān)缙诨蜓髿づ咛テ诘漠a(chǎn)物, 分別于約5 Ma和17 Ma開始裂解并形成新的洋殼(Garfunkel and Beyth, 2006; Bonatti et al., 2015), 而加利福尼亞灣在10 Ma左右開始裂解(Sumy et al., 2013; 李三忠等, 2019), 灣內(nèi)的瓜伊馬斯盆地則在約7 Ma開始裂解并形成新的洋殼(Lizarralde et al., 2007, 2011; Miller and Lizarralde, 2013; Teske et al., 2021)。瓜伊馬斯盆地位于加利福尼亞灣南北部的交界處, 是海灣南部最北端的擴(kuò)張海盆。盆地中發(fā)育廣泛且又劇烈的洋中脊巖漿作用以及熱液活動(dòng)(Curray and Moore, 1982; Lizarralde et al., 2007, 2011; Teske et al., 2021)。國(guó)際大洋發(fā)現(xiàn)計(jì)劃于2019年在瓜伊馬斯盆地進(jìn)行科考, 即IODP385航次?？瓶即诖似陂g共進(jìn)行了8個(gè)站位的鉆探工作, 其中的U1550屬于在軸洋中脊站位, 取得了海底擴(kuò)張?jiān)缙贛ORB巖心樣品(Teske et al., 2021)。本文通過對(duì)U1550巖心樣品中單斜輝石進(jìn)行詳細(xì)的微區(qū)原位主量元素分析, 計(jì)算單斜輝石結(jié)晶的溫度與壓力條件, 為探討海底擴(kuò)張?jiān)缙谘笾屑箮r漿分離結(jié)晶的物理化學(xué)特征及最深結(jié)晶深度提供重要線索, 對(duì)認(rèn)識(shí)并理解尚處于擴(kuò)張?jiān)缙诘难笾屑箮r漿作用具有重要意義。

1 地質(zhì)背景與樣品描述

加利福尼亞灣是東太平洋北部的一個(gè)狹長(zhǎng)的邊緣海, 位于北美大陸和下加利福尼亞半島之間, 屬于非常年輕的(轉(zhuǎn)換型)被動(dòng)大陸邊緣裂谷(Curray and Moore, 1982; Mercier de Lépinay et al., 2016; 汪品先, 2018; 孫珍等, 2021)。海灣呈NW-SE走向, 長(zhǎng)約1300 km, 寬100～150 km, 呈北窄南寬喇叭狀。灣內(nèi)由一系列狹長(zhǎng)的右旋轉(zhuǎn)換斷層和短小的新生洋脊組成, 隸屬于東太平洋海隆(EPR)最北端(圖1a)?；贕ripp and Gordon (2002)開發(fā)的洋中脊板塊邊界模型(HS3-NUVEL1A)計(jì)算, 灣內(nèi)洋脊的擴(kuò)張速率明顯低于EPR主體洋中脊的擴(kuò)張速率, 按Dick et al. (2003)的分類應(yīng)屬于慢速擴(kuò)張洋中脊(全擴(kuò)張速率約 48 mm/a)。加利福尼亞灣北部仍以陸殼為主, 主要構(gòu)造作用是大陸裂解, 南部以洋殼為主, 主要構(gòu)造作用是洋脊?轉(zhuǎn)換斷層共同主導(dǎo)的海盆擴(kuò)張(Atwater and Stock, 1998; Dixon et al., 2000; Plattner et al., 2007; Zhang et al., 2009; Sumy et al., 2013)。

圖中各洋脊段的全擴(kuò)張速率基于 Gripp and Gordon (2002)開發(fā)的洋中脊板塊邊界模型(HS3-NUVEL1A)計(jì)算所得。圖(a)和(b)的海洋測(cè)深底圖引自GeoMapAPP(Ryan et al., 2009); 圖(c)據(jù)IODP 385航次報(bào)告(Teske et al., 2021)修改, 圖1c中 U1545-1552 是 IODP385 的航次站位, 477～481是DSDP64的航次站位(DSDP64 航次報(bào)告 Curray and Moore, 1982)。

30 Ma之前整個(gè)加利福尼亞海岸是一條統(tǒng)一的俯沖帶。這條俯沖帶將北美板塊與一個(gè)現(xiàn)今已經(jīng)破碎的法拉隆(Farallon)板塊分隔開。法拉隆板塊在30 Ma左右沿著北美海岸向北美板塊之下發(fā)生向東的斜向俯沖, 分隔太平洋板塊和法拉隆板塊的洋脊也發(fā)生俯沖, 導(dǎo)致太平洋板塊開始直接與北美板塊接觸(圖 2a)。20 Ma前, 由于太平洋板塊以較高的速率向NNE向運(yùn)動(dòng), 從而在與北美板塊相接觸的區(qū)域演化成圣安德烈斯右旋走滑斷層。而原先的法拉隆板塊則分裂成兩個(gè)次級(jí)微板塊, 即胡安德富卡板塊和科科斯板塊。伴隨著法拉隆板塊的俯沖, 門多西諾(Mendocino)和里維拉(Riviera)三節(jié)點(diǎn)分別向北和向南躍遷, 圣安德烈斯走滑斷層也多次向東躍遷, 其長(zhǎng)度向北和向南發(fā)生了拓展(圖 2b)。20～10 Ma期間, 加利福尼亞一帶的俯沖作用完全停止, 被圣安德烈斯走滑斷層所替代(圖 2c)。在 10 Ma左右, 圣安德烈斯走滑斷層向東躍遷進(jìn)入墨西哥西側(cè)內(nèi)陸, 導(dǎo)致斷層縮短, 形成圣安德烈斯走滑斷層大拐彎(big bend)。屬于島弧型陸殼的下加利福尼亞(Baja California)開始以約 43～48 mm/a 的速度從北美板塊分離, 向北漂移并入太平洋板塊體系內(nèi), 從而開啟了現(xiàn)今的加利福尼亞灣裂谷和轉(zhuǎn)換斷層系統(tǒng)(圖 2d; Atwater and Stock, 1998; Dixon et al., 2000; Plattner et al., 2007; 李三忠等, 2019)。

MTJ. 門多西諾三節(jié)點(diǎn); RTJ. 里維拉三節(jié)點(diǎn); JF. 胡安德富卡板塊; LA. 洛杉磯。

瓜伊馬斯盆地位于加利福尼亞灣南部和北部的交界處, 是海灣南部最北端的擴(kuò)張海盆, 記錄著從大陸裂解伸展過渡到海底擴(kuò)張的地質(zhì)演化信息(圖1b)。盆地發(fā)育厚約 1～2 km富含有機(jī)質(zhì)的沉積物, 且發(fā)育大規(guī)模的巖漿作用和熱液活動(dòng)。相對(duì)于典型的洋中脊多發(fā)育在軸(on-axis; 距脊軸＜5 km)的噴出玄武巖, 瓜伊馬斯盆地同時(shí)發(fā)育豐富的離軸(off-axis; 0至＞50 km)玄武巖和淺層侵入巖——玄武質(zhì)巖床和巖墻(Curray and Moore, 1982; Lizarralde et al., 2007, 2011; Miller and Lizarralde, 2013; Teske et al., 2021)。

國(guó)際大洋鉆探發(fā)現(xiàn)計(jì)劃IODP385航次U1550站位位于瓜伊馬斯盆地內(nèi)最長(zhǎng)洋中脊的軸部(圖1b)。本次研究的4個(gè)巖心樣品取自該站位的A孔和B孔(圖3a), 兩孔所在位置的海水深度約2000 m, 與DSDP64航次的481站位鉆孔相鄰。A孔和B孔均打到了玄武巖層, 并打到玄武巖層約3 m長(zhǎng)就終孔。A孔玄武巖層位于海底以下204 m(～204 mbsf)處, 回收到約1.2 m玄武質(zhì)巖心; B孔玄武巖層位于～170 mbsf處, 獲得約1.6 m玄武質(zhì)巖心。所獲巖心主要呈細(xì)粒顯晶質(zhì)?隱晶質(zhì)結(jié)構(gòu)(圖3b), 礦物粒度大多部分＜2 mm,少數(shù)粒度可高達(dá)5 mm。在U1550A鉆孔中, 玄武質(zhì)巖從上到下從隱晶質(zhì)玄武巖過渡到粒度較粗的輝綠巖(或稱粒玄巖, 圖3f), 礦物粒度逐漸變大, 在幾厘米手標(biāo)本尺度中呈短程漸變關(guān)系, 指示較緩慢的巖漿結(jié)晶過程。在U1550B鉆孔中, 輝綠巖呈捕虜體的形式被隱晶質(zhì)玄武巖捕獲, 暗示至少有兩期巖漿脈期次作用(magma pulses; 圖3b、c)。玄武質(zhì)巖心相對(duì)新鮮, 海水蝕變作用較弱, 不發(fā)育或發(fā)育很少的氣孔或杏仁, 但發(fā)育很多細(xì)小的方解石脈(圖3b)。

為了便于后續(xù)討論, 本文將所獲巖心分為玄武巖、U1550A輝綠巖以及輝綠巖捕虜體三類。玄武巖主要呈灰黑色, 塊狀構(gòu)造, 發(fā)育斑狀和填間結(jié)構(gòu)等, 斑晶以斜長(zhǎng)石和單斜輝石為主(圖3c、d、e)。斜長(zhǎng)石斑晶約占2%～20%, 粒徑為0.5～3 mm, 多為針狀, 呈中空骸晶結(jié)構(gòu), 少量呈柱狀或者板狀, 可見聚晶結(jié)構(gòu)。單斜輝石斑晶約占1%～8%, 粒徑為0.2～1 mm, 多呈不規(guī)則它形?半自形粒狀, 也發(fā)育骸晶結(jié)構(gòu)?；|(zhì)為隱晶質(zhì)?玻璃質(zhì)結(jié)構(gòu), 主要由針狀斜長(zhǎng)石(含量5%～50%)、粒狀單斜輝石(含量2%～30%)微晶(＜0.2 mm)和玻璃(含量5%～50%)組成, 含少量的磁鐵礦、黃鐵礦和碳酸鹽礦物等。U1550A輝綠巖及捕虜體呈灰色?灰黑色, 塊狀構(gòu)造, 發(fā)育輝綠結(jié)構(gòu)或粒玄結(jié)構(gòu)、含長(zhǎng)嵌晶結(jié)構(gòu)(圖3b、c、f)。主要礦物為斜長(zhǎng)石(含量50%～70%)和單斜輝石(含量25%～45%), 其中斜長(zhǎng)石呈自形?半自形長(zhǎng)條狀, 粒徑為1～5 mm, 單斜輝石多為它形?半自形板狀, 粒度為0.5～3 mm。副礦物主要為磁鐵礦和少量黃鐵礦等。

2 分析方法與結(jié)果

2.1 分析方法

選取最新鮮的樣品開展全巖地球化學(xué)分析。先用砂紙對(duì)巖心樣品表面進(jìn)行打磨去皮, 再用去離子水洗凈, 晾干后使用剛玉碎樣機(jī)破碎成粒徑為0.5 cm左右的小碎粒, 再使用1%稀鹽酸浸泡以去除樣品中的碳酸鹽雜質(zhì), 去離子水洗凈后再次晾干, 最后用瑪瑙缽研磨至 200目以下粉末以備全巖化學(xué)分析。全巖主量元素分析在武漢上譜分析科技有限責(zé)任公司利用日本理學(xué) PrimusⅡ X射線熒光光譜儀(XRF, ZSXPrimusⅡ)分析完成。實(shí)驗(yàn)流程如下: ①將200目樣品置于105 ℃烘箱中烘干12小時(shí); ②稱取～1.0 g烘干樣品置于恒重陶瓷坩堝中, 在1000 ℃馬弗爐中灼燒2小時(shí), 取出待冷卻至室溫再稱量, 計(jì)算燒失量(LOI); ③分別稱取6.0 g助熔劑(Li2B4O7︰LiBO2︰LiF=9︰2︰1)、0.6 g樣品、0.3 g氧化劑(NH4NO3)置于鉑金坩堝中, 在1150 ℃熔樣爐中熔融14分鐘, 取出坩堝轉(zhuǎn)移到耐火磚上冷卻, 然后將玻璃片取出以備XRF測(cè)試, 分析精度優(yōu)于5%。

稀土元素含量在武漢上譜分析科技有限責(zé)任公司利用Agilent 7700e ICP-MS分析完成。實(shí)驗(yàn)流程如下: ①將200目樣品置于105 ℃烘箱中烘干12小時(shí); ②準(zhǔn)確稱取粉末樣品50 mg置于Teflon溶樣彈中; ③先后依次緩慢加入1 mL高純HNO3和1 mL高純HF; ④將Teflon溶樣彈放入鋼套, 擰緊后置于190 ℃烘箱中加熱24小時(shí)以上; ⑤待溶樣彈冷卻, 開蓋后置于140 ℃電熱板上蒸干, 然后加入1 mL HNO3并再次蒸干; ⑥加入1 mL高純HNO3、1 mL MQ水和1 mL內(nèi)標(biāo)In(濃度為1 μg/g), 再次將Teflon溶樣彈放入鋼套, 擰緊后置于190 ℃烘箱中加熱12小時(shí)以上; ⑦將溶液轉(zhuǎn)入聚乙烯料瓶中, 并用2% HNO3稀釋至100 g以備ICP-MS測(cè)試。使用GB/T14506.30-2010《硅酸鹽巖石化學(xué)分析方法》來測(cè)試樣品的元素含量。對(duì)于＞10 μg/g的稀土元素, 分析精度一般優(yōu)于5%, 對(duì)于＜10 μg/g的稀土元素, 分析精度一般優(yōu)于8%。

單斜輝石主量元素分析測(cè)試在兩個(gè)實(shí)驗(yàn)室完成。其中一部分?jǐn)?shù)據(jù)在中國(guó)科學(xué)院地球化學(xué)研究所礦床地球化學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室JEOL JXA-8530F電子探針(EPMA)分析儀上完成。儀器分析條件為: 加速電壓25 kV, 測(cè)試電流為10 nA。根據(jù)點(diǎn)位的情況, 分析選用“spot”模式(束斑直徑1 μm)或“circle”模式(束斑直徑10 μm)。選取標(biāo)樣為: 正長(zhǎng)石(K)、鉻透輝石(Ca、Cr、Mg、Si)、斜長(zhǎng)石(Na)、鎂鋁榴石(Ti、Al、Fe、Mn)、橄欖石(Ni)。另一部分?jǐn)?shù)據(jù)在中國(guó)科學(xué)院廣州地球化學(xué)研究所礦物學(xué)與成礦學(xué)重點(diǎn)實(shí)驗(yàn)室JEOL JXA-8230電子探針分析儀上完成。儀器分析條件為: 加速電壓15 kV, 測(cè)試電流為20 nA, 分析選用“spot”模式(束斑直徑1 μm)。選取標(biāo)樣為: 正長(zhǎng)石(K、Na)、透輝石(Ca、Si)、鎂鋁榴石(Fe)、方鎂石(Mg)、鐵鋁榴石(Al)、金紅石(Ti)、薔薇輝石(Mn)。

2.2 全巖主量元素特征

U1550鉆孔樣品全巖主量和微量元素組成見表1。U1550鉆孔樣品的燒失量(LOI)為0.37%～0.98%, 暗示其或多或少受到了海水或熱液蝕變的影響, 與巖相學(xué)觀察一致。在下文的圖解和討論中, 所有主要氧化物均已扣除燒失量, 換算成100%硅酸鹽成分, 并且全鐵均換算為FeOT。玄武巖主量元素成分變化較小(SiO2=49.36%～50.22%; MgO=6.88%～7.29%; FeOT=9.16%～9.93%; Mg#=56～59), 具有相對(duì)較低的CaO(10.66%～11.75%)和較高的Al2O3(15.94%～16.83%)含量。與玄武巖相比, U1550A輝綠巖含有相等的MgO(6.92%)和FeOT(9.93%), 具有相似的Mg#(56), 但具有較低的SiO2(49.36%), 稍高的CaO(11.75%)。輝綠巖捕虜體暫無全巖主量元素?cái)?shù)據(jù)?？傮w來說, 本次研究的樣品落在前人已發(fā)表的瓜伊馬斯盆地N-MORB數(shù)據(jù)范圍內(nèi)。在TAS分類圖解中, 玄武巖和U1550A輝綠巖均落在亞堿性巖漿范圍內(nèi)(圖4a), 在AFM巖漿分類圖解中, 玄武巖和輝綠巖均為拉斑玄武質(zhì)巖漿(圖4b)。結(jié)合樣品的稀土元素組成(全巖(La/Sm)N＜1), U1550鉆孔玄武巖為典型的在軸N-MORB。

表1 瓜伊馬斯盆地IODP385航次U1550鉆孔玄武巖和輝綠巖的全巖主量(%)和稀土元素(μg/g)組成

瓜伊馬斯N-MORB數(shù)據(jù)引自PetDB數(shù)據(jù)庫(kù)。

2.3 單斜輝石的成分特征

玄武巖中細(xì)粒單斜輝石斑晶SiO2=45.40%～50.66%, FeOT=9.16%～13.36%, MgO=10.56%～14.70%, CaO= 18.96%～21.96%, Mg#范圍為61～73(表2)。在輝石分類圖解中, 樣品測(cè)點(diǎn)主要落在透輝石和普通輝石區(qū)域(圖5a), 端元為Wo40.9～47.3En31.5～41.6Fs15.4～22.5。U1550A輝綠巖中的單斜輝石SiO2=49.31%～50.80%, FeOT=7.66%～9.74%, MgO=14.17%～15.12%, CaO=20.44%～21.81%, Mg#范圍為72～78。在輝石分類圖解中, 樣品的測(cè)點(diǎn)主要落在普通輝石區(qū)域(圖5a), 端元為Wo41.8～44.2En40.0～43.2Fs12.5～15.3。輝綠巖捕虜體中單斜輝石SiO2=48.24%～50.14%, FeOT=6.82%～7.36%, MgO=14.42%～15.13%, CaO=20.03%～22.14%, Mg#為78～79, 相對(duì)富Mg, 主要為普通輝石(圖5a), 端元為Wo42.0～45.3En40.9～43.9Fs11.1～15.7。與玄武巖中的單斜輝石斑晶相比, 輝綠巖中的單斜輝石含有較低的TiO2(1.20%～1.97%)和Na2O(0.36%～0.48%); 與U1550A輝綠巖中單斜輝石相比, 捕虜體中的單斜輝石則有較高的TiO2和相似的Na2O (圖4b、c)。

表2 瓜伊馬斯盆地IODP385航次U1550鉆孔玄武巖和輝綠巖中單斜輝石主量元素組成(%)

圖5 瓜伊馬斯盆地U1550鉆孔單斜輝石Wo-En-Fs分類圖解(a;底圖引自Morimoto et al., 1988)、單斜輝石Mg#與TiO2(b)及Na2O(c)的相關(guān)圖

3 討 論

3.1 U1550鉆孔玄武巖和輝綠巖的差異

U1550鉆孔玄武巖和輝綠巖均呈現(xiàn)N-MORB (全巖(La/Sm)N＜1)的特征, 指示其應(yīng)該來自于同一巖漿源區(qū)。U1550A輝綠巖和玄武巖的全巖成分相似, 但單斜輝石則呈現(xiàn)出一定的成分差異。U1550A輝綠巖的單斜輝石Mg#較高(72～78), 相對(duì)富Mg, 而玄武巖的單斜輝石Mg#較低(61～73), 相對(duì)富Fe, 暗示與U1550A輝綠巖單斜輝石平衡的母巖漿更接近原始巖漿成分。這與輝綠巖具有較低的全巖SiO2含量相一致。同時(shí), 從U1550A輝綠巖到玄武巖, 隨著Mg#降低, 單斜輝石的TiO2和Na2O含量大致呈線性升高趨勢(shì)(圖5b、c), 也表明U1550玄武巖相對(duì)于輝綠巖更演化, 經(jīng)歷了更高程度的結(jié)晶分異作用, 而輝綠巖則應(yīng)形成于Mg#較高的近原始巖漿。在U1550B鉆孔中, 輝綠巖呈捕虜體的形式被隱晶質(zhì)玄武巖捕獲(圖3b、c)。玄武巖發(fā)育斑狀和填間結(jié)構(gòu), 斜長(zhǎng)石斑晶多呈骸晶結(jié)構(gòu), 單斜輝石斑晶多為細(xì)粒不規(guī)則狀, 且多發(fā)育玻璃基質(zhì), 而輝綠巖捕虜體中礦物粒度較粗, 且多發(fā)育輝綠結(jié)構(gòu)或粒玄結(jié)構(gòu)、含長(zhǎng)嵌晶結(jié)構(gòu)(圖3c), 不發(fā)育玻璃基質(zhì)。綜上所述, 我們認(rèn)為U1550玄武巖與U1550A輝綠巖應(yīng)代表上洋殼的玄武巖層(Layer 2A), 而輝綠巖捕虜體則來源于上洋殼輝綠巖巖墻(Layer 2B), 這些樣品均適合單斜輝石溫壓計(jì)(Putirka et al., 1996; Neave and Putirka, 2017)的計(jì)算。但由于U1550鉆孔不發(fā)育玻璃冷凝邊, 且斑晶多為粒度長(zhǎng)但非常窄的中空骸晶結(jié)構(gòu)的針狀斜長(zhǎng)石斑晶, 以及粒度小的中空骸晶結(jié)構(gòu)的它形?半自形輝石斑晶, 暗示其結(jié)晶快速且與熔體未完全發(fā)生分異, 甚至可能是原位結(jié)晶的產(chǎn)物。故在后續(xù)的討論中, 我們用全巖成分(未扣除斑晶成分)來代替與單斜輝石平衡的熔體成分來計(jì)算。由于缺少輝綠巖捕虜體的全巖成分, 但U1550A輝綠巖的單斜輝石與捕虜體單斜輝石具有相似的主量成分, 故本文將U1550A輝綠巖的全巖成分作為捕虜體的平衡熔體進(jìn)行溫壓計(jì)算。

3.2 單斜輝石?熔體平衡關(guān)系判別

單斜輝石是火成巖中廣泛發(fā)育的造巖礦物, 對(duì)溫度和壓力有很強(qiáng)的敏感性, 因此常常被用來計(jì)算巖漿的結(jié)晶壓力與溫度, 進(jìn)而推算巖漿在地殼中的演化過程(Putirka, 2008)。Putirka et al. (1996)基于在玄武質(zhì)巖漿中單斜輝石?熔體間硬玉(NaAlSi2O6; Jd)組分的交換平衡、透輝石(CaMgSi2O6; Di)和鈣鐵輝石(CaFeSi2O6; Hd)組分與硬玉Jd組分的交換平衡、透輝石和鈣鐵輝石DiHd與鈣契爾馬克組分(CaAl2SiO6; CaTs)的交換平衡, 通過實(shí)驗(yàn)巖石學(xué)方法和地球化學(xué)熱力學(xué)理論中的吉布斯自由能的變化?r0與平衡常數(shù)eq的關(guān)系, 即–lneq=?r0, 將溫度和壓力與熱力學(xué)相關(guān)的平衡常數(shù)聯(lián)系起來, 從而建立起適用于玄武質(zhì)巖漿中單斜輝石?熔體平衡的溫壓計(jì)。后續(xù)的研究在此基礎(chǔ)上不斷利用新的相平衡數(shù)據(jù)或提高計(jì)算精度, 或提出特定巖漿體系下的計(jì)算方法(Putirka and Condit, 2003; Putirka, 2008; Neave and Putirka, 2017; Mollo et al., 2018)。

使用單斜輝石?熔體溫壓計(jì)之前往往要判斷單斜輝石與熔體是否平衡。根據(jù)溫壓計(jì)原理, 通過單斜輝石原位主量元素, 可以獲得單斜輝石中“透輝石+鈣鐵輝石(DiHd)”、“頑輝石+斜鐵輝石(EnFs)”和鈣契爾馬克組分(CaTs)含量; 而通過熔體成分回歸分析可以獲得體系內(nèi)預(yù)測(cè)的單斜輝石的DiHd+EnFs+ CaTs組分。對(duì)比預(yù)測(cè)和實(shí)測(cè)組分的切合度(1∶1線; 圖6a), 含量相差在20%誤差范圍內(nèi)可視為平衡(Putirka et al., 1996; Neave and Putirka, 2017)。此外, 還可以通過測(cè)算體系內(nèi)Fe-Mg分配系數(shù)來判斷平衡。Roeder and Emslie (1970)發(fā)現(xiàn)橄欖石與熔體之間的Fe-Mg交換系數(shù)D與溫度、氧逸度和熔體成分無關(guān), 且數(shù)值區(qū)間變化不大(0.3±0.03), 可以用作判斷橄欖石?熔體平衡的標(biāo)準(zhǔn)。后來, 巖石學(xué)家將這個(gè)方法應(yīng)用到輝石領(lǐng)域(Topliss and Carroll, 1995; Putirka, 2008), 如單斜輝石D(Fe-Mg)Cpx-liq(0.28±0.08),斜方輝石K(Fe-Mg)Opx-liq(0.29±0.06)。通過測(cè)定單斜輝石原位主量元素和玄武巖全巖主量元素, 可以計(jì)算出對(duì)應(yīng)的D值(圖6b)。但是, 使用D值來判斷單斜輝石?熔體是否平衡也存在一些局限性。例如,K受到系統(tǒng)冷卻速率的影響(Hammer, 2008), 當(dāng)冷卻速率較大時(shí), 實(shí)驗(yàn)觀測(cè)到的K值往往偏大。相關(guān)實(shí)驗(yàn)發(fā)現(xiàn)當(dāng)單斜輝石的生長(zhǎng)速率遠(yuǎn)大于熔體中Fe和Mg的擴(kuò)散速率時(shí),Fe和Mg都趨向于幾乎恒定的Fe/Mg值(Mollo and Hammer, 2017; 羅雕等, 2022)。同時(shí), 冷卻速率的變化也會(huì)影響單斜輝石及其共生成分的變化(Mollo et al., 2010)?；谶@些原因, Mollo et al. (2013)根據(jù)平衡實(shí)驗(yàn)提出, 將回歸分析獲得的預(yù)測(cè)單斜輝石端元與實(shí)際測(cè)試出的單斜輝石端元做差, 例如ΔDiHd=DiHd預(yù)測(cè)–DiHd實(shí)驗(yàn), 當(dāng)差值的絕對(duì)值≤0.1時(shí)可被認(rèn)為處于平衡狀態(tài)。通過上述三個(gè)平衡判定標(biāo)準(zhǔn), 本次研究共獲得15組平衡數(shù)據(jù)對(duì)用于以下計(jì)算。由于U1550樣品的含水量未知, 基于全球MORB含水量的平均值(0.1%～0.3%; Xia et al., 2013),本次的溫壓計(jì)算將樣品含水量設(shè)定為0.1%。

圖b陰影區(qū)為平衡區(qū)間(0.28±0.08; Putirka, 2008)。

3.3 單斜輝石結(jié)晶溫壓計(jì)算結(jié)果

將篩選后的15組有效數(shù)據(jù)對(duì)(其中玄武巖5組, 輝綠巖10組)進(jìn)行溫壓計(jì)算。本次研究使用多種溫壓計(jì)的計(jì)算結(jié)果進(jìn)行橫向?qū)Ρ? 擬選取最佳的溫度和壓力計(jì)。其中, 壓力計(jì)有以下幾種: ①Putirka et al. (1996)的公式P1, 是Putirka等人提出的第一個(gè)適用于玄武質(zhì)巖漿基于硬玉交換體系的單斜輝石?熔體壓力計(jì)(為了后續(xù)討論, 將此方法簡(jiǎn)寫為Putirka96P1); ②Putirka (2008)的公式30, 是Putirka等人改良后基于硬玉交換體系的單斜輝石?熔體壓力計(jì)(Putirka08Eqn30);③Putirka (2008)的公式32a, 是基于Nimis (1995)模型改進(jìn)的基于單斜輝石成分的壓力計(jì)(Putirka08Eqn32a);④Neave and Putirka (2017)中提出的基于硬玉交換體系的壓力計(jì)(N&P2017)。Mollo et al. (2018)也提出了一個(gè)新的高精度單斜輝石?熔體壓力計(jì), 但其方法主要運(yùn)用于堿性玄武巖的計(jì)算, 所以本次研究不予考慮。計(jì)算結(jié)果顯示, Putirka96P1在幾乎所有的平衡數(shù)據(jù)對(duì)中都得出了最大的壓力值, Putirka08Eqn30的計(jì)算結(jié)果約比Putirka96P1小1.0～1.6 kbar(圖7a)。Eqn32a返回的計(jì)算結(jié)果在輝綠巖中通常最低, 在玄武巖中規(guī)律性相對(duì)較差(圖7a)。前人通過對(duì)比實(shí)驗(yàn)數(shù)據(jù)與壓力計(jì)計(jì)算的結(jié)果, 發(fā)現(xiàn)Putirka96P1和Putirka08Eqn30擬合的回歸線不會(huì)經(jīng)過原點(diǎn), 意味著這兩種方法在高壓范圍(≥8 kbar)獲得的計(jì)算結(jié)果更可靠, 低壓范圍內(nèi)的計(jì)算結(jié)果往往偏大。Neave and Putirka (2017)指出Putirka08Eqn32a擬合的回歸線相對(duì)前兩個(gè)方法較好, 但實(shí)驗(yàn)結(jié)果顯示這種方法得出的壓力往往會(huì)比實(shí)際值高1～2 kbar, 對(duì)于＜10 kbar的貧水樣品尤其如此。因此, 他們提出了新的單斜輝石?熔體壓力計(jì)(N&P2017)。他們搜集了1 atm～20 kbar的數(shù)據(jù)集進(jìn)行測(cè)試, 發(fā)現(xiàn)此方法對(duì)拉斑玄武巖的計(jì)算精度相對(duì)較高, 擬合的回歸線接近零點(diǎn), 適用于更加寬泛的壓力范圍。雖然本次研究的輝綠巖樣品通過N&P2017計(jì)算得出的壓力相對(duì)Putirka08Eqn32a的計(jì)算結(jié)果偏大或相似, 但玄武巖的計(jì)算結(jié)果與前人的結(jié)果(Neave and Putirka, 2017)一樣(圖7a)。鑒于本次研究的樣品均為拉斑質(zhì)的N-MORB, 含水量一般都很低, 所以最終選取N&P2017作為本文計(jì)算U1550玄武巖與輝綠巖中的單斜輝石結(jié)晶壓力的最佳壓力計(jì)(誤差為±1.4 kbar)。

圖7 不同單斜輝石?熔體壓力計(jì)(a)和溫度計(jì)(b)計(jì)算結(jié)果對(duì)比柱狀圖

溫度計(jì)有兩種: ①Putirka et al. (1996)的公式T1, 是Putirka等人提出的第一個(gè)適用于玄武質(zhì)巖漿基于硬玉交換體系的單斜輝石?熔體溫度計(jì)(Putirka96T1); ②Putirka (2008)的公式33, 是Putirka等人改進(jìn)后的單斜輝石?熔體溫度計(jì)(Putirka08Eqn33), 有效降低了約10～20 ℃的計(jì)算誤差, 且在后續(xù)研究中得到廣泛運(yùn)用。計(jì)算結(jié)果顯示, 玄武巖、單斜輝石和輝綠巖捕虜體三個(gè)單斜輝石數(shù)據(jù)對(duì)用Putirka96T1得出的溫度普遍偏大, 比Putirka08Eqn33計(jì)算出的值高約20～40 ℃(圖7b)。因此, 本次研究選取Putirka等人改良后的Putirka08Eqn33作為計(jì)算U1550玄武巖與輝綠巖捕虜體中的單斜輝石結(jié)晶溫度的最佳溫度計(jì)(誤差為±45 ℃)。

溫壓計(jì)算結(jié)果顯示: 玄武巖中單斜輝石結(jié)晶的壓力范圍在–0.5～5.4 kbar(誤差±1.4 kbar), 相當(dāng)于0～17.8 km(誤差±4.6 km)的深度(假定1 kbar≈3.3 km; Putirka and Condit, 2003); 結(jié)晶溫度主要集中于1126～1170 ℃, 誤差±45 ℃(圖7c、d)。U1550A輝綠巖以及輝綠巖捕虜體單斜輝石的結(jié)晶壓力范圍在1.6～5.9 kbar(誤差±1.4 kbar), 相當(dāng)于5.3～19.5 km(誤差±4.6 km)的深度; 結(jié)晶溫度主要集中在1158～1212 ℃(誤差±45 ℃)。與U1550玄武巖中單斜輝石相比, 輝綠巖中單斜輝石結(jié)晶溫度與壓力更大, 尤其是U1550B孔的輝綠巖捕虜體中的單斜輝石顯示最大的結(jié)晶溫度和壓力(圖8), 并且玄武巖的結(jié)晶壓力范圍寬, 而輝綠巖相對(duì)集中于3～6 kbar。Whattam et al. (2020)同樣運(yùn)用N&P2017壓力計(jì)和Putirka08Eqn33溫度計(jì)計(jì)算了全球其他洋脊段MORB中單斜輝石(數(shù)據(jù)引自PetDB數(shù)據(jù)庫(kù), http: //www.earthchem.org/petdb)的結(jié)晶溫度和壓力, 結(jié)果顯示全球MORB中單斜輝石的結(jié)晶壓力范圍在–1.8～5.4 kbar之間, 約0～17.8 km深, 結(jié)晶溫度范圍在1133～1240 ℃之間。Lissenberg and Dick (2008)研究表明一些富鎂單斜輝石也可以通過洋殼內(nèi)的熔體?巖石反應(yīng)(melt-rock reaction)形成“高壓”特征, 并且認(rèn)為單斜輝石Mg#越大, 越“高壓”。相對(duì)于全球其他洋脊段MORB, U1550玄武巖和輝綠巖(捕虜體)單斜輝石具有較低的Mg#(均＜80), 但反而計(jì)算得出較高的結(jié)晶壓力。另外, 在-Mg#圖(圖8b)中, 本次研究的樣品并沒有呈現(xiàn)單一的線性關(guān)系, 因此我們認(rèn)為U1550玄武巖和輝綠巖單斜輝石沒有發(fā)生熔體?巖石反應(yīng), 其計(jì)算所得的壓力值可以代表其真實(shí)的結(jié)晶壓力。相對(duì)于全球慢速?超慢速洋脊段MORB(如大西洋洋中脊MAR、西南印度洋洋中脊SWIR和北冰洋Gakkel洋脊等), U1550樣品中的單斜輝石具有相似的溫度和壓力的變化范圍, 但在給定的溫度下, U1550樣品中的單斜輝石具有更高的結(jié)晶壓力或深度, 或者在給定的壓力或深度下, U1550樣品中的單斜輝石具有更低的結(jié)晶溫度(圖8c), 暗示瓜伊馬斯盆地在軸洋脊段MORB具有更陡峭的冷卻軌跡。

3.4 單斜輝石結(jié)晶溫壓的指示意義

早期研究認(rèn)為, 洋中脊下方會(huì)發(fā)育一個(gè)巨大的在軸巖漿房(AMC; axial magma chamber), 巖漿房壁發(fā)育均勻的輝長(zhǎng)巖, 底部則由一系列層狀鎂鐵?超鎂鐵質(zhì)巖組成。它為上洋殼的巖墻(Layer 2B)和噴出巖(Layer 2A)提供了物質(zhì)來源, 下洋殼(Layer 3)則從此處結(jié)晶增生, 這是經(jīng)典的無限洋蔥洋殼增生模型(infinite onion; Cann, 1974)。一些學(xué)者認(rèn)為這種模型存在一些缺陷, 因?yàn)樵诼?超慢速擴(kuò)張洋脊, 較低的巖漿供應(yīng)率并不能支撐大型巖漿房的存在(Sleep, 1975; Kuznir and Bott, 1976)。因此提出無限韭菜洋殼增生模型(infinite leek), 即小批量的巖漿以侵入體的形式供應(yīng)洋殼增生(Nisbet and Fowler, 1978)。這些模型顯示出一個(gè)共同的特點(diǎn): 即巖漿結(jié)晶的深度仍位于平均6～7 km厚的洋殼內(nèi)。地球物理觀測(cè)指示未蝕變的地幔橄欖巖大規(guī)模就位于海底以下10 km左右, 故他們認(rèn)為即使是在慢速擴(kuò)張洋中脊, 巖漿的結(jié)晶深度和海底熱液循環(huán)的最深溫度應(yīng)該不超過10 km (Sleep and Barth, 1997)。前人也通過橄欖石熔體包裹體的揮發(fā)分(如CO2、H2O等)研究認(rèn)為洋中脊巖漿結(jié)晶作用只局限在下洋殼, 深度一般不超過殼幔邊界, 最深的巖漿結(jié)晶深度不超過10 km (Colman et al., 2015; Wanless et al., 2015), 即使在超慢速擴(kuò)張的北冰洋Gakkel洋中脊, 橄欖石記錄的最深的巖漿結(jié)晶深度不超過9 km(Shaw et al., 2010)。這與地震反射成像探測(cè)出的熔體透鏡體(melt lens)一般發(fā)育在0.76～4.5 km深(Sinton and Detrick, 1992; Van Ark et al., 2007), 最深達(dá)到6 km的觀測(cè)相一致(Chen and Lin, 2004; Wanless and Shaw, 2012)。在快速擴(kuò)張洋脊, 因巖漿房可以長(zhǎng)期穩(wěn)定存在, 巖漿供給量充足, 下洋殼增生后產(chǎn)生的熱量需要很長(zhǎng)的時(shí)間才可散去, 故快速擴(kuò)張洋脊背景下的巖漿結(jié)晶模型常常帶有巖漿房這一特征要素。例如, 前人通過研究阿曼蛇綠巖提出“大巖漿房(large magma chamber)”洋殼增生模型, 源區(qū)洋脊下方可發(fā)育一個(gè)寬約20 km和高4 km的以熔體為主的巖漿房(Pallister and Hopson, 1981; Smewing, 1981)。然而, 相關(guān)地球物理資料表明這個(gè)模型過于理想化, 實(shí)際的以熔體為主的巖漿房(相當(dāng)于熔體透鏡體)可能很小(大致寬約1 km, 高50 m), 但其下可發(fā)育規(guī)模很大的以晶粥為主的巖漿房。這類模型與早期的“無限洋蔥”洋殼增生模型不同之處在于提出大部分巖漿結(jié)晶于以熔體為主的淺層巖漿房?jī)?nèi), 晶粥向下沉降或向外擴(kuò)散形成下洋殼, 被命名為“輝長(zhǎng)巖冰川(gabbro glacier)”模型(Henstock et al., 1993; Henstock, 2002)。后續(xù)一些學(xué)者意識(shí)到下洋殼的輝長(zhǎng)巖(Layer 3)很可能是原位結(jié)晶的產(chǎn)物, 故“席狀巖床”(‘sheeted sill’ model)模型應(yīng)運(yùn)而生(Boudier et al., 1996; Kelemen et al., 1997)。

然而, 一些學(xué)者針對(duì)慢速?超慢速擴(kuò)張洋中脊存在高壓結(jié)晶作用(Bender et al, 1978; Green et al, 1979), 提出了巖漿作用超過10 km的深部結(jié)晶的洋殼增生模型(Elthon et al., 1982; Cannat, 1993, 1996)。近20年, 許多學(xué)者通過地球物理觀測(cè)和玄武巖主量元素壓力計(jì)研究認(rèn)為洋中脊巖漿分離結(jié)晶作用可以超過10 km, 甚至可以深達(dá)30 km(Herzberg, 2004; Villiger et al., 2007; Schlindwein and Schmid, 2016; Wanless and Behn, 2017)。Bennett et al. (2019)通過斜長(zhǎng)石熔體包裹體研究也指出Gakkel洋中脊最深的巖漿結(jié)晶深度可高達(dá)16.4 km, 指示洋中脊巖漿分離結(jié)晶作用可以深達(dá)巖石圈地幔, 不止局限于下洋殼深度。Michael and Cornell (1998)研究表明洋中脊巖漿結(jié)晶作用的壓力與擴(kuò)散速率、巖漿供應(yīng)成反比。在快速擴(kuò)張洋中脊和一些慢速擴(kuò)張洋中脊(如大西洋靠近冰島地幔柱的Reykjanes洋中脊), 巖漿供應(yīng)充足, 巖漿結(jié)晶作用一般只局限在洋殼內(nèi), 結(jié)晶深度淺; 然而在典型的慢速?超慢速擴(kuò)張洋中脊以及一些洋中脊斷裂帶(fracture zones)附近, 巖石圈厚且冷, 地幔部分熔融程度小, 產(chǎn)生的熔體相對(duì)少且熔體相對(duì)富Na, 巖漿供應(yīng)缺乏, 同時(shí), 地幔部分熔融產(chǎn)生巖漿的終止深度也相對(duì)更深, 往往在深部地幔中就開始結(jié)晶, 由于結(jié)晶深度深, 巖漿難以抵達(dá)淺部洋殼。后來, 基于對(duì)慢速擴(kuò)張洋脊產(chǎn)出的橄欖巖體中的輝長(zhǎng)巖以及洋脊自身熱條件的研究, 提出一種介于前兩者之間的模型, 稱為“李子布丁(plum pudding)”模型(Cannat, 1993, 1996)。在這個(gè)模型中, 慢速擴(kuò)張洋脊的洋殼由許多小型的巖漿侵入呈成布丁狀嵌入在不同深度的地?；蜓髿?nèi)。

加利福尼亞灣是東太平洋北部一個(gè)狹長(zhǎng)的邊緣海, 屬于非常年輕的(轉(zhuǎn)換型)被動(dòng)大陸邊緣裂谷, 灣內(nèi)由一系列狹長(zhǎng)的右旋轉(zhuǎn)換斷層和短小的新生洋脊組成, 而瓜伊馬斯盆地位于加利福尼亞灣南北部的交界處, 是海灣南部最北端的擴(kuò)張海盆, 記錄著洋脊擴(kuò)張?jiān)缙诘牡刭|(zhì)演化信息(圖1; Curray and Moore, 1982; Mercier de Lépinay et al., 2016;汪品先, 2018; 孫珍等, 2021)。IODP385航次U1550站位位于瓜伊馬斯盆地內(nèi)最長(zhǎng)洋中脊的軸部, 其擴(kuò)張速率約為 48 mm/a, 屬于慢速擴(kuò)張洋中脊。U1550玄武巖單斜輝石結(jié)晶溫度和壓力的計(jì)算結(jié)果表明, 玄武巖結(jié)晶深度相對(duì)分散, 既存在接近洋中脊巖漿噴口的淺層上洋殼區(qū)域, 也存在更深的下洋殼和巖石圈地幔源區(qū)(深達(dá)17.8±4.6 km)。輝綠巖及輝綠巖捕虜體的結(jié)晶壓力較集中且很大, 最大結(jié)晶深度達(dá)到約19.5±4.6 km。同時(shí), 相對(duì)于全球其他成熟洋脊段MORB, 瓜伊馬斯盆地在軸洋脊段MORB具有更陡峭的冷卻軌跡(圖8c), 暗示其巖漿結(jié)晶深度可能超過了20 km。地球物理學(xué)家通過解譯加利福尼亞灣周邊布置的 NARS-Baja、RESBAN等陸地地震臺(tái)網(wǎng)數(shù)據(jù), 顯示沿加利福尼亞灣走向的二維S波速度異常剖面在 60～70 km深處分布有一系列代表熔體發(fā)育的巖漿房的低速異常體, 瓜伊馬斯盆地下方的低速體從75 km一直延伸到20～25 km, 佐證了以上推測(cè)(Wang et al., 2009; Zhang et al., 2009; Zhang and Paulssen, 2012)。Bonatti and Seyler (1987)研究認(rèn)為紅海北部(全擴(kuò)張速率＜15 mm/a)Zabargad輝長(zhǎng)巖也形成于海底擴(kuò)張?jiān)缙? 其結(jié)晶壓力可高達(dá)8.7 kbar, 即約29 km深。以上討論表明, 在海底擴(kuò)張?jiān)缙? 擴(kuò)張速率很可能均為慢速擴(kuò)張, 其巖漿結(jié)晶深度與前人研究慢速?超慢速擴(kuò)張成熟洋脊段的巖漿結(jié)晶深度可超過10 km的觀點(diǎn)不謀而合。但不同于慢速?超慢速擴(kuò)張成熟洋脊段, 海底擴(kuò)張?jiān)缙谘笾屑箮r漿分離結(jié)晶作用具有更陡峭的冷卻軌跡。這暗示著在海底擴(kuò)張?jiān)缙趲r石圈地幔也很厚, 地幔部分熔融產(chǎn)生巖漿的終止深度也相對(duì)更深, 導(dǎo)致巖漿結(jié)晶深度深, 部分巖漿難以抵達(dá)淺部洋殼。但同時(shí)瓜伊馬斯盆地廣泛發(fā)育在軸和離軸巖漿作用, 結(jié)合地球物理觀測(cè)資料, 指示海底擴(kuò)張?jiān)缙? 深部巖漿房巖漿供給并不匱乏, 這些巖漿房很可能位于更深的部位。

4 結(jié) 論

(1) U1550巖心樣品玄武巖和輝綠巖及捕虜體均屬于拉斑系列。玄武巖和U1550A輝綠巖代表上洋殼的玄武巖層(Layer 2A), 而U1550B輝綠巖捕虜體很可能來自上洋殼輝綠巖巖墻(Layer 2B)。玄武巖中細(xì)粒單斜輝石斑晶主要為透輝石和普通輝石, 而U1550A輝綠巖和U1550B輝綠巖捕虜體中的單斜輝石主要為普通輝石, 相對(duì)富Mg, 低Ti和Na。

(2) 單斜輝石?熔體溫壓計(jì)計(jì)算結(jié)果顯示: 瓜伊馬斯盆地U1550玄武巖的單斜輝石的結(jié)晶溫度(1126～1170 ℃, 誤差±45 ℃)和壓力(–0.5～5.4 kbar, 誤差±1.4 kbar)相對(duì)輝綠巖及輝綠巖捕虜體中單斜輝石(=1158～1212 ℃, 誤差±45 ℃;=1.6～5.9 kbar, 誤差±1.4 kbar)均較低, 結(jié)晶深度(0～17.8 km, 誤差±4.6 km)相對(duì)后者(5.3～19.5 km, 誤差±4.6 km)較淺。與全球其他洋脊段MORB單斜輝石相比, U1550單斜輝石具有更陡峭的冷卻軌跡。

(3) 在海底擴(kuò)張?jiān)缙? 擴(kuò)張速率很可能均為慢速擴(kuò)張, 其巖漿結(jié)晶深度可超過10 km, 深達(dá)巖石圈地幔深處。這暗示著在海底擴(kuò)張?jiān)缙趲r石圈地幔很厚, 地幔部分熔融產(chǎn)生巖漿的終止深度也相對(duì)更深, 導(dǎo)致巖漿結(jié)晶深度深, 部分巖漿難以抵達(dá)淺部洋殼。但不同于慢速?超慢速擴(kuò)張成熟洋脊段, 海底擴(kuò)張?jiān)缙谏畈繋r漿房巖漿供給可能并不匱乏, 這些巖漿房很可能位于更深的部位。

致謝:感謝中國(guó)IODP辦公室以及IODP385航次全體船員及全體科考人員。感謝中國(guó)科學(xué)院地球化學(xué)研究所李響老師, 中國(guó)科學(xué)院廣州地球化學(xué)研究所曹永華老師在實(shí)驗(yàn)中給予的指導(dǎo)。評(píng)審專家中國(guó)科學(xué)院海洋研究所張國(guó)良研究員和中國(guó)科學(xué)院廣州地球化學(xué)研究所楊陽(yáng)副研究員提出了非常寶貴的修改意見, 在此一并致以衷心的感謝！

葛振敏, 鄢全樹, 趙仁杰, 施美娟. 2020. 科科斯脊玄武巖斜長(zhǎng)石礦物化學(xué)及地質(zhì)意義. 海洋學(xué)報(bào), 42(7): 93–107.

李三忠, 曹現(xiàn)志, 王光增, 劉博, 李璽瑤, 索艷慧, 姜兆霞, 郭玲莉, 周潔, 王鵬程, 朱俊江, 汪剛, 趙淑娟, 劉永江, 張國(guó)偉. 2019. 太平洋板塊中?新生代構(gòu)造演化及板塊重建. 地質(zhì)力學(xué)學(xué)報(bào), 25(5): 642–677.

羅雕, 侯通, 王旭東, 田野, 張招崇. 2022. 冷卻過程對(duì)玄武巖體系結(jié)構(gòu)、礦物成分以及熱力學(xué)參數(shù)估算的影響. 礦物巖石地球化學(xué)通報(bào), 41(4): 759–775, 693.

牛耀齡. 2022. 范式革命: 玄武巖記錄有噴發(fā)時(shí)巖石圈厚度的信息, 沒有地幔潛在溫度的記憶. 科學(xué)通報(bào), 67(3): 301–306.

孫珍, 李付成, 林間, 孫龍濤, 龐雄, 鄭金云. 2021. 被動(dòng)大陸邊緣張?破裂過程與巖漿活動(dòng): 南海的歸屬. 地球科學(xué), 46(3): 770–789.

汪品先. 2018. 被動(dòng)大陸邊緣//中國(guó)大洋發(fā)現(xiàn)計(jì)劃辦公室,海洋地質(zhì)重點(diǎn)實(shí)驗(yàn)室(同濟(jì)大學(xué)). 大洋鉆探五十年. 上海: 同濟(jì)大學(xué)出版社: 114–128.

楊錦, 楊帆, 黃小龍, 朱圣柱, 苗秀全, 賀鵬麗. 2022. 南海擴(kuò)張前序巖漿活動(dòng): 解譯華南三水盆地古近紀(jì)玄武質(zhì)巖漿作用過程. 大地構(gòu)造與成礦學(xué), 46(3): 530–551.

張國(guó)良, 羅青, 陳立輝. 2017. 大洋地?；瘜W(xué)組成不均一性成因研究回顧及展望. 海洋地質(zhì)與第四紀(jì)地質(zhì), 37(1): 1–13.

Atwater T, Stock J. 1998. Pacific-North America plate tectonics of the Neogene southwestern United States: An update., 40: 375–402.

Bender J F, Hodges F N, Bence A E. 1978. Petrogenesis of basalts from the project FAMOUS area: Experimental study from 0 to 15 kbar., 41: 277–302.

Bennett E N, Jenner F E, Millet M A. Cashman K V, Lissenberg C J. 2019. Deep roots for mid-ocean-ridge volcanoes revealed by plagioclase-hosted melt inclusions., 572: 235–239.

Bonatti E, Cipriani A, Lupi L. 2015. The Red Sea: Birth of an ocean // Rasul N M A, Stewart I C F. The Red Sea: The Formation, Morphology, Oceanography and Environmentof A Young Ocean Basin. Springer Earth System Sciences, Berlin Heidelberg: 29–44.

Bonatti E, Seyler M. 1987. Crustal underplating and evolutionin the Red Sea Rift: Uplifted gabbro/gneiss crustal complexeson Zabargad and Brothers Islands.:, 92(B12): 12803–12821.

Boudier F, Nicolas A, Ildefonse B. 1996. Magma chambers in the Oman ophiolite: Fed from the top and the bottom., 144(1–2): 239–250.

Cann J R. 1974. A model for oceanic crustal structure developed., 39: 169–187.

Cannat M. 1993. Emplacement of mantle rocks in the seafloor at mid-ocean ridges., 98(B3): 4163–4172.

Cannat M. 1996. How thick is the magmatic crust at slow spreading ridges?, 101(B2): 2847–2857.

Chen Y J, Lin J. 2004. High sensitivity of ocean ridge thermal structure to changes in magma supply: The Galápagos spreading center., 221: 263–273.

Colman A, Sinton J M, Wanless V D. 2015. Constraints from meltinclusions on depths of magma residence at intermediate magma supply along the Galápagos spreading center., 412: 122–131.

Combier V, Seher T, Singh S C, Crawford W C, Cannat M, Escartín J, Dusunur D. 2015. Three-dimensional geometry of axial magma chamber roof and faults at Lucky Strike volcano on the Mid-Atlantic Ridge.:, 120: 5379–5400.

Coogan L. 2014. The lower oceanic crust // Holland H, Turekian K K. Treatise on Geochemistry. Elsevier: 497–541.

Curray J R, Moore D G. 1982. Initial Reports of the Deep Sea Drilling Project, 64: Washington, D C (US Government Printing Office).

Detrick R S, Mutter J C, Buhl P, Kim I I. 1990. No evidence from multichannel reflection data for a crustal magma chamber in the MARK area on the Mid-Atlantic Ridge., 347(6288): 61–64.

Dick H J B, Lin J, Schouten H. 2003. An ultraslow-spreading class of ocean ridge., 426(6955): 405–412.

Dixon T, Farina F, Demets C, Suarez-Vidal F, Fletcher J, Marquez-Azua B, Miller M, Sanchez O, Umhoefer P. 2000. New kinematic models for Pacific-North America Motion from 3 Ma to Present,Ⅱ: Evidence for a “Baja California Shear Zone”., 27: 3961–2964.

Eason D, Sinton J. 2006. Origin of high-Al N-MORB by fractional crystallization in the upper mantle beneath the Galápagos Spreading Center., 252: 423–436.

Elthon D, Casey J F, Komor S. 1982. Mineral chemistry of ultramafic cumulates from the North Arm Massif of the Bay of Islands ophiolite: Evidence for high-pressure crystal fractionation of oceanic basalts., 87(B10): 8717–8734.

Gale A, Dalton C A, Langmuir C H, Su Y, Schilling J G. 2013. The mean composition of ocean ridge basalts.,,, 14: 489–518.

Garfunkel Z, Beyth M. 2006. Constraints on the structural development of Afar imposed by the kinematics of the major surrounding plates.,,, 259(1): 23–42.

Green D H, Hibberson W O, Jaques A L. 1979. Petrogenesis of mid-ocean ridge basalts // McElhinney M W. The Earth: Its Origin, Structure, and Evolution. Elsevier: 265–299.

Gripp A E, Gordon R G. 2002. Young tracks of hotspots and current plate velocities., 150: 321–364.

Hammer J E. 2008. Experimental studies of the kinetics and energetic of magma crystallization., 69: 9–59.

Henstock T. 2002.Compaction control of melt distribution at fast-spreading mid-ocean ridges., 29(7): 411–414.

Henstock T, Woods A, White R. 1993. The accretion of oceanic crystal by episodic sill intrusion., 98(B3): 4143–4161.

Herzberg C. 2004. Partial crystallization of mid-ocean ridge basalts in the crust and mantle., 45(12): 2389–2405.

Kelemen P, Koga K, Shimizu N. 1997. Geochemistry of gabbro sills in the crust-mantle transition zone of the Oman ophiolite: Implications for the origin of the oceanic lower crust., 146: 475–488.

Kuznir N, Bott M H P. 1976. A thermal study of oceanic crust., 47: 83–95.

Lissenberg C J, Dick H J B. 2008. Melt-rock reaction in the lower oceanic crust and its implications for the genesis of mid-ocean ridge basalt., 271: 311–325.

Lizarralde D, Axen G J, Brown H E, Fletcher J M, González-Fernández A, Harding A J, Holbrook W S, Kent G M, Paramo P, Sutherland F, Umhoefer P J. 2007. Variation in styles of rifting in the Gulf of California., 448(7152): 466–469.

Lizarralde D, Soule S A, Seewald J S, Proskurowski G. 2011. Carbon release by off-axis magmatism in a young, sedimented spreading centre.: 50–54.

Mercier de Lépinay M, Loncke L, Basile C, Roest W R, Patriat M, Maillard A, De Clarens P. 2016. Transform continental margins — Part 2: A worldwide review.s, 693: 96–115.

Michael P J, Cornell W C. 1998. Influence of spreading rate and magma supply on crystallization and assimilation beneath mid-ocean ridges: Evidence from chlorine and major element chemistry of mid-ocean ridge basalts., 103: 18325–18356.

Miller N C, Lizarralde D. 2013. Thick evaporites and early rifting in the Guaymas Basin, Gulf of California., 41(2): 283–286.

Mollo S, Blundy J, Scarlato P, De Cristofaro S P, TecchiatoV, Di Stefano F, Vetere F, Holtz F, Bachmann O. 2018. An integrated--H2O-lattice strain model to quantify the role of clinopyroxene fractionation on REE+Y and HFSE patterns of mafic alkaline magmas: Application to eruptions at Mt. Etna., 185: 32–56.

Mollo S, Del Gaudio P, Ventura G, Iezzi G, Scarlato P. 2010. Dependence of clinopyroxene composition on cooling rate in basaltic magmas: Implications for thermobarometry., 118(3–4): 302–312.

Mollo S, Hammer J E. 2017. Dynamic crystallization in magmas // Heinrich W, Abart R. Mineral Reaction Kinetics: Microstructures, Textures, Chemical and Isotopic Signatures. EMU Notes in Mineralogy, European Mineralogical Union and Mineralogical Society of Great Britain and Ireland, 373–418.

Mollo S, Putirka K, Misiti V, Soligo M, Scarlato P. 2013. A new test for equilibrium based on clinopyroxene-melt pairs: Clues on the solidification temperatures of' Etne analkaline melts at post-eruptive conditions., 352: 92–100.

Morimoto N. 1988. Nomenclature of pyroxenes., 39(1): 55–76.

Neave D A, Putirka K D. 2017. A new clinopyroxene-liquid barometer, and implications for magma storage pressures under Icelandic rift zones., 102: 777–794.

Nimis P. 1995. A clinopyroxene geobarometer for basaltic systems based on crystals-structure modeling., 121: 115–125.

Nisbet E G, Fowler C M R. 1978. The Mid-Atlantic Ridge at 37 and 45°N: Some geophysical and petrological constraints., 54: 631–660.

Pallister J, Hopson C. 1981. Samail Ophiolite plutonic suite: Field relations, phase variation, cryptic variation and layering, and a model of a spreading ridge magma chamber.:, 86(B4): 2593–2644.

Plattner C, Malservisi R, Dixon T H, Lafemina P, Sella G F, Fletcher J, Suarez-Vidal F. 2007. New constraints on relative motion between the Pacific Plate and Baja California microplate (Mexico) from GPS measurements., 170: 1373–1380.

Putirka K, Condit C. 2003. Cross section of a magma conduitsystem at the margins of the Colorado Plateau., 31: 701–704.

Putirka K D. 2008. Thermometers and barometers for volcanic systems., 69: 61–120.

Putirka K D, Johnson M, Kinzler R J, Longhi J, Walker D. 1996. Thermobarometry of mafic igneous rocks based on clinopyroxene-liquid equilibria, 0–30 kbar., 123: 92–108.

Roeder P L, Emslie R F. 1970. Olivine-liquid equilibrium., 29: 275– 289.

Ryan W B F, Carbotte S M, Coplan J O, O’Hara S, Melkonian A, Arko R, Weissel R A, Ferrini V, Goodwillie A, Nitsche F, Bonczkowski J, Zemsky R. 2009. Global multi-resolution topography synthesis.,,,10, Q03014.

Schlindwein V, Schmid F. 2016. Mid-ocean-ridge seismicity reveals extreme types of ocean lithosphere., 535: 276–279.

Shaw A M, Behn M D, Humphris S E, Sohn R A, Gregg P M. 2010. Deep pooling of low degree melts and volatile fluxes at the 85°E segment of the Gakkel Ridge: Evidence from olivine-hosted melt inclusions and glasses., 289: 311–322.

Singh S C, Crawford W C, Carton H, Seher T, Combier V, Cannat M, Pablo C J, Dusunur D, Escartin J, Miranda J M. 2006. Discovery of a magma chamber and faults beneath a Mid-Atlantic Ridge hydrothermal field., 442(7106): 1029–1032.

Sinton J M, Detrick R S. 1992. Mid-ocean ridge magma chambers., 97(B1): 197–216.

Sleep N H. 1975. Formation of oceanic crust: Some thermal constraints., 80: 4037– 4042.

Sleep N H, Barth G A. 1997. The nature of the lower crust and shallow mantle emplaced at low spreading rates., 279: 181–191.

Smewing J. 1981. Mixing characteristics and compositional differences in mantle-derived melts beneath spreading axes: Evidence from cyclically layered rocks in the ophiolite of North Oman.:, 86(B4): 2645–2659.

Streck M J. 2008. Mineral textures and zoning as evidence for open system processes., 69(1): 595–622.

Sumy D F, Gaherty J B, Kim W Y, Diehl T, Collins J A. 2013. The mechanisms of earthquakes and faulting in the southern Gulf of California., 103: 487–506.

Teske A, Lizarralde D, H?fig T W, and the Expedition 385 Scientists. 2021. International Ocean Discovery Program Expedition 385 Scientific Report.

Topliss M J, Carroll M R. 1995. An experimental study of the influence of oxygen fugacity on Fe-Ti oxide stability, phase relations and mineral-melt equilibria in ferro- basaltic systems., 36: 1137–1170.

Van Ark E, Detrick R S, Canales J P, Carbotte S M, Harding A J, Kent G M, Nedimovic M R, Wilcock W S D, Diebold J B, Babcock J M. 2007. Seismic structure of the Endeavour Segment, Juan de Fuca Ridge: Correlations with seismicity and hydrothermal activity., 112, B02401.

Villiger S, Müntener O, Ulmer P. 2007. Crystallization pressures of mid-ocean ridge basalts derived from major element variations of glasses from equilibrium and fractional crystallization experiments., 112, B01202.

Wang Y, Forsyth D W, Savage B. 2009. Convective upwelling in the mantle beneath the Gulf of California., 462(7272): 499–501.

Wanless V D, Behn M D. 2017. Spreading rate-dependent variations in crystallization along the global mid-ocean ridge system.,,, 18: 3016–3033.

Wanless V D, Shaw A M. 2012. Lower crustal crystallization and melt evolution at mid-ocean ridges., 5: 651–655.

Wanless V D, Shaw A M, Behn M D, Soule S A, Escartin J, Hamelin C. 2015. Magmatic plumbing at Lucky Strike volcano based on olivine-hosted melt inclusion compositions.,,, 16: 126–147.

Whattam S A, Shervais J W, Reagan M K, Coulthard D A, Stern R J, Pearce J A, Jones P, Seo J, Putirka K, Chapman T, Heaton D, Li H, Nelson W R, Shimizu K, Stern R J. 2020. Mineral compositions and thermobarometry of basalts and boninites recovered during IODP Expedition352 to the Bonin forearc., 105(10): 1490–1507.

White W M, Klein E M. 2014. Composition of the oceanic crust // Holland H, Turekian K K. Treatise on Geoche-mistry. Elsevier: 497–541.

Xia Q K, Liu J, Liu S C, Kovács I, Feng M, Dang L. 2013. High water content in Mesozoic primitive basalts of the North China Craton and implications on the destruction of cratonic mantle lithosphere., 361: 85–97.

Zhang X, Paulssen H. 2012. Geodynamics of the Gulf of California from surface wave tomography., 192–193: 59–67.

Zhang X, Paulssen H, Lebedev S, Meier T. 2009. 3D shear velocity structure beneath the Gulf of California from Rayleigh wave dispersion., 279: 255–262.

Deep Crystallization Processes of MORB During the Early Stage of Seafloor Spreading: Implication from Clinopyroxene of IODP385-U1550 in Guaymas Basin, Gulf of California

ZHANG Tianxiang1, 2, XIE wei1, 2*

(1.MNR Key Laboratory of Marine Hazards Forecasting, Hohai University, Nanjing 210024, Jiangsu, China; 2. College of Oceanography, Hohai University, Nanjing 210024, Jiangsu, China)

The depth of MORB fractional crystallization is one of the key factors to control the MORB composition and oceanic crustal accretion.In order to determine the deepest crystallization depth associated with mid ocean ridge magma process in the early stage of seafloor spreading, we take detailed petrological and mineralogical study on the basaltic rock samples in the Guaymas Basin which were drilled in the site U1550 during the International Ocean Drilling Program (IODP) Expeditions 385. U1550 samples belong to tholeiitic magma, mainly composed of basalt, dolerite and dolerite xenoliths.The clinopyroxene in our samples is predominantly augite and diopside.The crystallization temperature (1126 – 1170 ℃, ±45 ℃) and pressure (–0.5 – 5.4 kbar, ±1.4 kbar) of clinopyroxene in basalt calculated by the clinopyroxene-liquid thermobarometers are both lower than those of dolerite and dolerite xenoliths (1158 – 1212 ℃, ±45 ℃; 1.6 – 5.9 kbar, ±1.4 kbar), and the crystallization depth (0 – 17.8 km) is shallower than the latter (5.3 – 19.5 km). Compared to clinopyroxenes from the other global spreading centers, U1550 clinopyroxenes have a steeper cooling trajectory. Comprehensive studies suggest that in the early stage of seafloor spreading, the spreading rate is likely to be slow, and the depth of magma crystallization can exceed 10 km. It implies thatthe mantle lithosphere was very thick in the early stage of seafloor spreading, and the terminate depth of magma produced the mantle partial melting was relatively deeper, resulting in deep crystallization of magma and difficulty to upwell shallow oceanic crust. However, unlike mature oceanic ridges spread at low speed, the magma supply of deep magma chambers may not be scarce in the early stage of seafloor spreading, instead the magma chamber could be located in deeper depth than we thought.

IODP 385; Gulf of California; early stage of seafloor spreading; clinopyroxene; deep crystallization of magma

2023-03-12;

2023-05-12;

2023-05-25

國(guó)家自然科學(xué)基金項(xiàng)目(42276067、41873031)和南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室(廣州)人才團(tuán)隊(duì)引進(jìn)重大專項(xiàng)(GML2019ZD0202)聯(lián)合資助。

張?zhí)煜?1998–), 男, 碩士研究生, 海洋科學(xué)專業(yè)。E-mail: qifeix@163.com

頡煒(1985–), 男, 副教授, 從事洋中脊巖漿巖地球化學(xué)研究。E-mail: weixie@hhu.edu.cn

P67; P595

A

1001-1552(2023)06-1345-018

10.16539/j.ddgzyckx.2023.01.302