閆紀(jì)元，李旭平，鄢全樹

1) 山東科技大學(xué)，山東省沉積成礦作用實(shí)驗(yàn)室，山東青島，266590； 2)國家海洋局第一海洋研究所，海洋沉積與環(huán)境地質(zhì)國家海洋局重點(diǎn)實(shí)驗(yàn)室，山東青島，266061

內(nèi)容提要: 南海新生代玄武巖中發(fā)育有大量的單斜輝石斑晶，與橄欖石、角閃石等斑晶共存，多數(shù)被熔蝕，呈渾圓狀，部分發(fā)育有很窄的反應(yīng)邊。研究表明單斜輝石具有巨晶的特征，是在高壓下在與寄主巖同源的巖漿中形成，沒有經(jīng)歷長時間的沉淀生長而直接向上運(yùn)移被帶到地表，因而具有巨晶的主量及微量元素的特征卻不具有巨晶的形態(tài)。單斜輝石稀土元素含量不高，HREE及LREE虧損，而MREE富集；大離子親石元素Rb、Ba、Sr均出現(xiàn)明顯的虧損，高場強(qiáng)元素Nb、Zr出現(xiàn)虧損而Hf略富集。Nb、Ta與Zr、Hf分餾明顯。Th較為富集，而Pb強(qiáng)烈虧損，U的變化范圍較大。研究發(fā)現(xiàn)南海新生代玄武巖中的單斜輝石的來源較為簡單，為地幔柱的直接產(chǎn)物，并沒有受到洋中脊—地幔柱相互作用的影響，由于巖漿上升速度較快壓力迅速下降，橄欖石大量晶出，引起了巖漿成分的變化，致使單斜輝石與寄主巖成分和結(jié)構(gòu)上沒有達(dá)到平衡。

南海位于歐亞板塊、印—澳板塊以及菲律賓板塊三大板塊交匯處，是西太平洋大型的邊緣海之一 (Karig, 1971; Briais et al., 1993; 鄢全樹等，2007；王葉劍等，2009)。南海擴(kuò)張期后(15.5Ma，Briais et al., 1993)，南海海盆及周邊地區(qū)發(fā)生了范圍較廣的板內(nèi)火山作用，影響了中南半島、雷瓊半島、珠江口盆地以及南海洋盆廣大地區(qū)。近年來，國內(nèi)外學(xué)者對南海新生代玄武巖作了大量的研究(王賢覺，1984；Tu et al., 1992; 鄢全樹等，2008；王葉劍等，2009)。但是對南海新生代玄武巖中單礦物的研究工作仍然有限，鄢全樹等(2007，2008)利用電子探針對單斜輝石和橄欖石的主量元素成分進(jìn)行了測定，并分析了巖石形成過程的溫壓條件，氧逸度條件以及潛在的地幔溫度等。而單礦物微量元素及稀土元素的精確測量結(jié)果，以及對微量元素及稀土元素的富集規(guī)律、演化趨勢和特征的研究鮮有報(bào)道。

單斜輝石作為玄武巖類巖石中的主要礦物相，同時也是微量及稀土元素的重要賦存礦物相，其化學(xué)成分記錄了巖漿成因、演化及巖漿作用、物理化學(xué)條件等重要信息 (白志民， 2000；鄢全樹，2007)。產(chǎn)于地幔深部的巨晶礦物更是可以用來研究玄武質(zhì)巖漿在地幔深部的地球化學(xué)特征和演變規(guī)律(于津海等，2001)。本文利用激光剝蝕技術(shù)(LA- ICP- MS)對南海新生代玄武巖中的單斜輝石進(jìn)行了微量、稀土元素的分析測定，重點(diǎn)討論了單斜輝石的微量和稀土元素特征，并探討了單斜輝石及其母巖漿的成因，對深入了解南海底新生代玄武巖的特征及巖漿的形成、演化研究具有重要的啟示意義。

1 地質(zhì)背景及樣品分析方法

圖1 南海及周邊地區(qū)地質(zhì)簡圖及取樣位置圖(據(jù)鄢全樹，2007，有修改)Fig. 1 Sketch geological map of the South China Sea and its adjacent areas and simpling locaties of this study(modified after Yan Quanshu et al., 2007)

南海是西太平洋邊緣海中面積最大的海盆(圖1)，自北向南可分為北部大陸邊緣、洋盆和南部大陸邊緣三個部分，其北部是與大陸的平滑過渡區(qū)，南部和東部是由于印度—澳大利亞板塊與巽他大陸碰撞導(dǎo)致南海擴(kuò)張停止時形成的擠壓俯沖區(qū)，西部是以紅河斷裂帶以及越東斷裂為標(biāo)志的走滑引張區(qū)(Briais et al., 1993; 萬玲等，2005；張漢泉和吳廬山，2005；姚伯初，1996；鄢全樹等，2007)。根據(jù)楊文采和宋海濱(2014)的研究，南海的擴(kuò)張屬于過成熟期的洋陸轉(zhuǎn)換作用。這一時期洋脊已經(jīng)完成俯沖，海盆逐漸縮小而俯沖帶匯聚，構(gòu)造十分復(fù)雜。南海擴(kuò)張開始于印度板塊與歐亞板塊碰撞之后，而同時期太平洋板塊向歐亞板塊之下發(fā)生俯沖(姚伯初，1996)。這可能為南海的擴(kuò)張?zhí)峁┝藬D壓的地應(yīng)力(楊文采和宋海濱，2014)，導(dǎo)致南海近南北方向發(fā)生擴(kuò)張，形成了具有洋殼屬性的中央海盆(趙海玲等，1999；閻全人等，2000)以及中央海盆中早期近東西方向的磁異常帶(Brais et al., 1993;趙海玲等，1999)。而中新世后，由于印度—澳大利亞板塊與巽他大陸碰撞，擴(kuò)張停止(Lee et al., 1994)，開始向菲律賓及巽他大陸下形成俯沖(李昌年等，2005；鄢全樹和石學(xué)法，2007)。擴(kuò)張期后新近系—第四系的板內(nèi)火山活動影響了南海洋盆及周邊的廣大地區(qū)。

本文研究的單斜輝石寄存于南海海底的新生代玄武巖。南海新生代玄武巖在TAS圖解中投點(diǎn)于玄武巖、堿玄巖—碧玄巖、粗面玄武巖等區(qū)域，全巖化學(xué)成分(表1)以及單斜輝石主量元素的研究均表明其屬于堿性玄武巖系列，年齡在7.9～3.8Ma (鄢全樹，2008)，為擴(kuò)張期后板內(nèi)火山作用的產(chǎn)物，地球化學(xué)特征表明其具有OIB型的稀土及微量元素的特點(diǎn)，同位素的研究則表明其源區(qū)具有不均一性，可能與地幔柱相關(guān)(鄢全樹和石學(xué)法，2007；鄢全樹等，2008)。南海新生代玄武巖中同時有斑晶和微晶產(chǎn)出，本次研究的單斜輝石為玄武巖中較大的斑晶，代表了早期結(jié)晶的產(chǎn)物，只有部分尚與寄主巖保持平衡關(guān)系；而微晶單斜輝石則是巖漿正常結(jié)晶的產(chǎn)物，除了淬火導(dǎo)致的假高壓效應(yīng)造成的不平衡以外，其他與寄主巖達(dá)到了平衡。為了研究和敘述的方便，在下文中除了涉及到巖相學(xué)的描述之外，均直接稱“單斜輝石”。

表1 單斜輝石寄主巖的主量元素(%)分析結(jié)果Table 1 The analyses of major elements (%) of the host rocks for the clinopyroxenes

樣品S04- 11、S04- 12、S04- 12- 2、S04- 12- 5、S04- 12- 11、S04- 12- 12采自南海中央海盆的海山中，S04- 16則采自于北大陸坡底。樣品S04- 11、S04- 12- 11、S04- 12- 12中，斑晶含量10%～20%，主要由輝石、角閃石、橄欖石組成；基質(zhì)約占40%，以斜長石微晶、輝石以及火山玻璃為主；氣孔大量發(fā)育(圖2a，2c)。其中S04- 11中斑晶以輝石和角閃石為主，橄欖石次之；粒徑總體在0.2～0.5mm；輝石晶形較好，部分發(fā)育良好的解理，均呈現(xiàn)出不同程度的熔蝕圓化(圖2a)，部分甚至被熔蝕呈港灣狀。S04- 12- 11中斑晶含量不高，輝石斑晶較小，0.2mm左右，多呈聚晶產(chǎn)出；基質(zhì)中針狀的斜長石發(fā)育，隨巖漿流動方向定向排列(圖2c)。S04- 12- 12中發(fā)育有大量的斜長石，0.2～1mm；輝石斑晶不多，且明顯熔蝕。來自北大陸坡底部的S04- 16氣孔發(fā)育，斑晶以輝石和角閃石為主(圖2d)，約占20%，較大的輝石斑晶可大于1mm，且發(fā)育有由較小的輝石斑晶形成的聚晶；基質(zhì)中發(fā)育有斜長石微晶、輝石以及火山玻璃。S04- 12、S04- 12- 2、S04- 12- 5為火山碎屑巖，其巖屑以堿性玄武巖為主，膠結(jié)物為綠泥石、蒙脫石等粘土礦物；受后期地質(zhì)作用的改造，顯微鏡下顯示出明顯的蝕變特征，大部分的斑晶礦物都被蝕變僅殘留假像，只有部分較大的晶體殘留，多呈渾圓狀(圖2b)。對火山碎屑巖中主要礦物的研究表明它們與堿性玄武巖具有同源演化的特征(鄢全樹和石學(xué)法，2009)。

圖2 南海新生代玄武巖中單斜輝石斑晶的光學(xué)特征Fig. 2 Optical images for phenocryst clinopyroxenes (a) S04- 11樣品(單偏光鏡下)，堿玄巖，可見單斜輝石斑晶，具有微弱的反應(yīng)邊；(b) S04- 12- 2樣品(單偏光鏡下)，火山碎屑巖，堿性玄武巖巖屑發(fā)生蝕變，其中的橄欖石僅殘留假像；(c) S04- 12- 11樣品(單偏光鏡下)，堿玄巖，單斜輝石斑晶呈渾圓狀，基質(zhì)中發(fā)育斜長石微晶，半定向排列；(d) S04- 16樣品(單偏光鏡下)，堿玄巖，可見單斜輝石受熔蝕圓化，局部呈港灣狀。礦物代號：Cpx—單斜輝石；Pl—斜長石；Ol—橄欖；石Vs—?dú)饪?(a) Sample S04- 11. Tephrite. Clinopyroxene phenocrysts have indistinct reaction rims. (b) Sample S04- 12- 2. Volcanoclastic rocks. Fragments of alkali basalts are altered and pseudomorph of olivine can be observed. (c) Sample S04- 12- 1. Tephrite. Clinopyroxene phenocrysty have a rounded shape and feldspar microlites in the groundmass orient towards the direction of magma flow. (c)Sample S04- 16. Tephrite. Clinopyroxene phenocryst has been resorbed into rounded shape and has reaction rim

LA- ICP- MS工作在中國地質(zhì)科學(xué)院國家地質(zhì)實(shí)驗(yàn)測試中心完成。儀器采用德國Finnigan公司的Element2等離子體質(zhì)譜儀與美國New Wave公司的UP213，213nm Nd:YAG固體激光器聯(lián)用。等離子質(zhì)譜測試的參數(shù)：冷卻氣流為16 L/min，樣品氣流為0.85 L/min，輔助氣流為0.73L/min，射頻功率為1200W；激光系統(tǒng)工作參數(shù)：采樣錐孔、截取錐孔分別為1.0mm和0.7mm。測試采用氦氣作為載體進(jìn)行激光剝蝕，載氣流為0.61L/min，激光波長為213nm，能量密度達(dá)41.2J/cm2，激光束斑直徑40μm，頻率10Hz，輸出能量在2mJ。測試過程中首先遮擋激光束空白計(jì)數(shù)12s，然后對樣品中的某一個點(diǎn)位連續(xù)剝蝕43s，停止剝蝕，再利用氦氣清洗進(jìn)樣系統(tǒng)15s，單點(diǎn)總分析時間約為70s。測試采用NIST 612和KL2- G兩種作為外標(biāo)，采用Si作為內(nèi)標(biāo)。

2 主量元素特征

單斜輝石中的Al具有特殊意義，其在6次配位體和4次配位體中的分布與溫壓條件緊密相關(guān)(劉英俊等，1984), 根據(jù)n(AlⅥ)/n(AlⅣ)可以定性衡量單斜輝石的結(jié)晶壓力(Aoki and Kushiro, 1968; Thompson, 1974; Wass, 1979)，劃分不同壓力下形成的單斜輝石 (Aoki and Kushiro, 1973)。南海新生代玄武巖中單斜輝石的n(AlⅥ)/n(AlⅣ)變化范圍為0.60～1.83，平均為1.28，近于匈牙利Nograd火成巖省中單斜輝石巨晶 (Dobosi and Jenner, 1999)。在單斜輝石的n(AlⅣ) —n(AlⅥ)判別圖(圖3a)上，本次研究的單斜輝石大部分落入了“麻粒巖中的單斜輝石”以及“玄武巖包體中的單斜輝石”范圍內(nèi)，且絕大部分落入中國東部新生代玄武巖中的巨晶單斜輝石圈定的范圍，高于太平洋Marshall海山鏈 (Dieu, 1995) 和印度洋Kerguelen熱點(diǎn)地區(qū) (Damasceno et al., 2000) 堿性玄武巖中單斜輝石，進(jìn)一步說明了南海新生代玄武巖中的單斜輝石在成分上大部分都是巨晶，只不過是沒有足夠的時間形成巨晶或者是在形成巨晶之后在上升過程中由于溫壓條件的變化導(dǎo)致碎裂而使得粒度變小(邱家驤，廖群安，1996)。

表2 南海新生代玄武巖中單斜輝石的電子探針分析結(jié)果(%)

注：Si、AlⅥ、AlⅣ均以6個氧原子為基礎(chǔ)計(jì)算，Mg#=n(Mg)/ [n(Mg)+n(Fe2+)]。

在單斜輝石的n(Si)—n(Al)關(guān)系圖上(圖3b)(Kushiro, 1960；鄢全樹等，2007)，南海新生代玄武巖中的單斜輝石落在巨晶單斜輝石的范圍內(nèi)，而太平洋Marshall海山鏈堿性玄武巖中的單斜輝石則大部分落在了堿性玄武巖斑晶的范圍內(nèi) (Dieu, 1995) ，僅有Wodejebato海山的部分樣品落在了巨晶單斜輝石范圍附近。Wodejebato海山中堿性玄武巖中的單斜輝石礦物化學(xué)組成變化范圍較大，可以用分離結(jié)晶程度、源區(qū)富集程度或者結(jié)晶壓力的不同來解釋 (Dieu, 1995)。分析認(rèn)為，落在巨晶單斜輝石范圍這些樣品可能較其他的樣品具有更高的結(jié)晶壓力。

3 微量元素及稀土元素特征

3.1稀土元素特征

分析結(jié)果(表3)表明，南海新生代玄武巖中單斜輝石的稀土元素含量不高，且變化范圍較大。ΣREE變化范圍為47.53×10-6～96.74×10-6，平均為62.73×10-6；LREE/HREE變化范圍為2.31～6.00，平均為3.05，相應(yīng)的LaN/YbN為1.47～7.12，平均為2.30；CeN/YbN變化范圍為1.47～6.70，平均為2.96。

圖3 南海新生代玄武巖中單斜輝石的n(AlⅣ) —n (AlⅥ)圖解(a) 和n(Si) —n(Al)圖解(b)(據(jù)Aoki and Shiba, 1973)Fig. 3 The n(AlⅣ)—n(AlⅥ)diagram and n(Si)—n(Al) diagram of the clinopyroxenes in Cenozoic basalts from the South China Sea (after Aoki and Shiba, 1973) (a) A以下，火成巖中單斜輝石；A—B，玄武巖包體中單斜輝石；B—C之間，麻粒巖中單斜輝石；C以上，榴輝巖中單斜輝石. (b)A—巨晶單斜輝石；B—堆積巖中單斜輝石；C—堿性玄武巖中斑晶；D—拉斑玄武巖中斑晶(數(shù)據(jù)區(qū)域根據(jù)邱家驤和廖群安，1996) (a) Below A , the clinopyroxene in igneous rocks; A—B, theclinopyroxene in xenoliths of basalts; B—C, the clinopyroxene in granulites; Above C, the clinopyroxene in eclogites. (b) A—the clinopyroxene megacrysts; B—the clinopyroxene in cumulates; C—the phenocryst in alkali basalts; D—the phenocryst in tholeiitic basalts (the areas are all from Qiu Jiaxiang and Liao Qun'an, 1996)

圖4 南海新生代玄武巖中單斜輝石的球粒隕石標(biāo)準(zhǔn)化蛛網(wǎng)圖(a)和原始地幔標(biāo)準(zhǔn)化蛛網(wǎng)圖(b); 球粒隕石數(shù)據(jù)來自于 McDonough and Sun (1995), 原始地幔數(shù)據(jù)來自于Sun and McDonough (1989)Fig.4 REE patterns normalized to CI chondrite (McDonough and Sun, 1995) and multi- trace element patterns normalized to Primitive mantle (Sun and McDonough ,1989) for the clinopyroxenes in Cenozoic basalts from the South China Sea

單斜輝石中HREE與LREE均有不同程度的虧損，HREE的虧損程度高于LREE，而MREE富集。LaN/SmN在0.43～1.68，平均為0.68，而GdN/LuN在2.28～4.00，平均為2.88，高于LaN/SmN。表現(xiàn)在球粒隕石標(biāo)準(zhǔn)化曲線(圖4a)上為兩端下凹，中間凸起的配分形式，最高點(diǎn)為Sm。δEu平均為1.00，無Eu異常。δPr平均為1.02，不存在陳道公等(1997)中推測的巨晶中Pr的正異常。

樣品S04- 11中，04- 11- 11_1的LREE和LILE總體較高，LaN/YbN值為7.12，高于總體的平均值，而不同于其他的樣品，推測是在上升過程中受到巖漿的同化。對寄主玄武巖中斜長石的研究(鄢全樹等，2008)表明，S04- 11樣品噴出的位置距離南海擴(kuò)張后的洋中脊斷裂區(qū)較遠(yuǎn)，巖漿上升速度慢，這使得巖漿對輝石的同化有了更多的反應(yīng)時間，同化作用較強(qiáng)而引起了微量元素含量的變化。樣品S04- 12- 11中的單斜輝石出現(xiàn)了標(biāo)準(zhǔn)化曲線相互交叉的現(xiàn)象，這與陳道公等(1997)，黃婉康等(1991)觀察到的現(xiàn)象相似，但該樣品中的配分曲線呈現(xiàn)出多變的趨勢，考慮到鏡下觀察該樣品中有斜長石的晶出，很有可能是斜長石的同時晶出導(dǎo)致了單斜輝石稀土分配規(guī)律的變化(夏群科等，1998)，而非巖漿演化的標(biāo)志。

表3 南海新生代玄武巖中單斜輝石的微量元素離子探針分析結(jié)果 (×10-6) Table 3 The LA- ICP- MS analyses of the clinopyroxenes in Cenozoic basalts from the South China Sea (×10-6)

續(xù)表3

3.2 微量元素特征

除一個特殊的樣品之外，本次研究中大離子親石元素(LILE)較低，Rb變化范圍在0.00×10-6～0.33×10-6，平均為0.09×10-6；Ba變化范圍在0.02×10-6～4.79×10-6，平均為0.62×10-6；Sr變化范圍在48.61×10-6～111.60×10-6，平均為74.23×10-6。而樣品04- 11- 11_1受到巖漿同化作用，Rb、Ba、Sr的含量分別為10.80×10-6,165.90×10-6和104.40×10-6，高于原始的單斜輝石。與寄主玄武巖中Rb、Ba、Sr的含量相比，單斜輝石中LILE的含量較低。

高場強(qiáng)元素(HFSE)中Zr的變化范圍為50.53×10-6～163.40×10-6，平均值為82.88×10-6；Hf的變化范圍為1.98×10-6～6.42×10-6，平均為3.36×10-6；Nb的變化范圍為0.08×10-6～1.35×10-6，平均為14.11×10-6；Ta的變化范圍為0.08×10-6～1.35×10-6,平均為0.58×10-6。均低于寄主玄武巖中HFSE的含量。

微量元素原始地幔標(biāo)準(zhǔn)化蛛網(wǎng)圖上(圖4b)，LILE變化范圍較大，出現(xiàn)虧損，HFSE呈現(xiàn)明顯的分餾。LILE在蛛網(wǎng)圖上出現(xiàn)Rb、Ba、Sr明顯的負(fù)異常；HFSE中Nb、Ta出現(xiàn)明顯的負(fù)異常，Hf具有弱的正異常。Nb/Ta與Zr/Hf的值變化較大，Nb/Ta為2.02～11.08，Zr/Hf為20.22～33.03，平均值分別為4.75和24.44，這與單斜輝石在玄武巖中的分配系數(shù)相符，是單斜輝石結(jié)晶分異的表現(xiàn)(劉艷榮等，2012)。Th、U、Pb的平均值分別為0.19×10-6,0.04×10-6以及0.25×10-6，但U的含量相對變化較大。在原始地幔標(biāo)準(zhǔn)化圖解上，Th較為富集，而Pb強(qiáng)烈虧損。單斜輝石根據(jù)Nb/U的比值可分為兩類，一類為堿性玄武巖中的單斜輝石，Nb/U在26.43～123.67；另一類為火山碎屑巖中的單斜輝石，Nb/U在3.39～23.11，它們很有可能是在形成后固結(jié)成巖過程中受到了地殼物質(zhì)的影響。

4 討論

4.1 南海新生代玄武巖單斜輝石地球化學(xué)特征

主量元素的研究已表明南海新生代玄武巖中的單斜輝石具有巨晶的特征，且大部分落在巨晶單斜輝石 (Irving and Frey，1984) 的稀土元素的范圍之內(nèi)，并與之具有相似的球粒隕石標(biāo)準(zhǔn)化曲線(圖5a)。稀土元素在輝石中的分布主要受其晶體化學(xué)性質(zhì)的控制(肖森宏，1990)，稀土元素標(biāo)準(zhǔn)化曲線的相似性，表明了其在輝石中具有相似的賦存規(guī)律，但是本次研究的單斜輝石稀土元素總量總體較高，則可能暗示了南海新生代玄武巖源區(qū)的特殊性。

與中國東部新生代堿性玄武巖中的巨晶單斜輝石(肖森宏，1990；陳道公等，1997；肖燕等，2008)相比(圖5a)，南海新生代玄武巖中的單斜輝石具有稀土元素總量高，重稀土分異弱的特點(diǎn)。中國東部新生代堿性玄武巖中的巨晶單斜輝石ΣREE在19.42×10-6～41.39×10-6，LREE/HREE平均為3.14，略高于本次研究，而其GdN/LuN遠(yuǎn)高于本次研究，且Mg#平均為0.74，低于本次研究，說明中國東部新生代玄武巖中的巨晶單斜輝石演化程度總體上要高于本次研究。中國東部新生代玄武巖的單斜輝石隨巖漿演化，形成了巨晶單斜輝石曲線相互交叉的現(xiàn)象(陳道公等，1997；黃婉康等，1991)，這是本次研究中的單斜輝石所不具有的。且本次研究中稀土配分曲線同中國東部新生代玄武巖巖漿演化早期形成的單斜輝石配分曲線極為相似，暗示了南海新生代玄武巖中的單斜輝石形成于巖漿演化較為早期的階段，單斜輝石的Mg#較高也說明了這一點(diǎn)。但南海新生代玄武巖中單斜輝石稀土元素總體仍偏高，這表明它們源區(qū)性質(zhì)與中國東部(大陸板內(nèi)堿性玄武巖)的不同，源區(qū)應(yīng)該更加富集稀土元素。

圖5 南海新生代玄武巖中單斜輝石的微量元素與巨晶單斜輝石(a)、太平洋Marshall海山鏈玄武巖中的單斜輝石(b) 以及印度洋Kerguelen Archipelago和Kerguelen Plateau玄武巖中單斜輝石(b、c)微量元素含量的比較Fig. 5 Comparison of trace element patterns between clinopyroxenes in Cenozoic basalts from the South China Sea and megacryst clinopyroxenes in basalts throughout the world , clinopyroxenes in basats from the Marshall Seamount Chain and clinopyroxenes in Kerguelen Archipelago and Kerguelen Plateau 圖5 (a)球粒隕石數(shù)據(jù)來自于 McDonough and Sun (1995)；圖5(b)、(c)球粒隕石數(shù)據(jù)來自Anders and Grevesse (1989) CI chondrite values inFig. 5 (a) from McDonough and Sun (1995) and CI chondrite values inFig. 5 (b), (c) from McDonough and Sun (1995)

與產(chǎn)自相似構(gòu)造環(huán)境的太平洋Marshall海山鏈(Dieu, 1995)和印度洋Kergulen地區(qū)(Damasceno et al., 2000)相比，具有相似的REE擴(kuò)展標(biāo)準(zhǔn)化曲線，三者均具有Sr、Zr的負(fù)異常，且南海新生代玄武巖中單斜輝石的REE擴(kuò)展曲線絕大部分落在Mashall海山鏈、Kergulen群島玄武巖中的單斜輝石圈定的范圍內(nèi)(圖5b)。

Sr—Nd—Pb同位素以及Hf同位素的研究表明，Kerguelen Plateau 北部邊緣1140站位的拉斑玄武巖受到了地幔柱—洋中脊相互作用的影響，顯示出了Kerguelen 地幔柱與SEIR組分的混合(Yang et al., 1998；Weis and Frey, 2000； Damasceno et al., 2000)。玄武巖中的單斜輝石化學(xué)成分也表現(xiàn)出了OIB與MORB過渡的特征，隨著受洋中脊影響程度的增強(qiáng)，LREE含量升高， REE擴(kuò)展配分曲線上LREE虧損程度減小(圖5c)(Damasceno et al., 2000)。而Kerguelen 群島Mont Crozier處的弱堿性玄武巖則代表了Kerguelen 地幔柱的特征，沒有地殼的混染，也沒有受到虧損源的影響，單斜輝石的成分相對一致(Damasceno et al., 2000)。

鄢全樹和石學(xué)法(2006，2007)認(rèn)為南海新生代玄武巖的形成受到了洋中脊—地幔柱相互作用的影響。其同位素組成類似于Kerguelen Plateau(Frey et al., 2000；鄢全樹等，2008)，其源區(qū)由兩個端元混合而成，一個是似印度洋MORB型的軟流圈地幔，另一個則是EM2端元，可能是來自于海南地幔柱(Zou Haibo et al., 2000；鄢全樹等，2008)。但南海新生代玄武巖中的單斜輝石中卻沒有出現(xiàn)類似于1140站位玄武巖單斜輝石中REE規(guī)律變化的現(xiàn)象(圖5c)，反而與Mont Crozier弱堿性玄武巖中的單斜輝石擴(kuò)展配分曲線頗為一致。且大部分落在Mont Crozier弱堿性玄武巖范圍之內(nèi)(圖5b、c)。這表明本文研究的這些單斜輝石可能并未受到洋中脊—地幔柱相互作用的影響，而是地幔柱的直接產(chǎn)物，是在地幔柱上升過程中分離結(jié)晶形成的(鄢全樹等，2007)。但更精確的結(jié)論需要必要的同位素分析以及更大范圍的采樣。

4.2南海新生代玄武巖中單斜輝石的成因

多數(shù)研究認(rèn)為巨晶單斜輝石是堿性玄武巖漿在地幔高壓環(huán)境下形成的(Irving, 1984; 陳道公等，1997; Woodland and Jugo，2007)，是巖漿在特殊的位置適當(dāng)?shù)臏貕簵l件下長時間結(jié)晶的產(chǎn)物，生長速度很慢(10-10～10-11cm/s)(Cashman and Marsh, 1988; Marsh et al., 1991)，生長1cm的晶體需要上千年的時間(Richter and Carmichael，1993),后來被巖漿攜帶經(jīng)過噴發(fā)至地表環(huán)境。

本次研究的單斜輝石在玄武巖中以斑晶形式產(chǎn)出，從粒度大小來判斷并非巨晶，但是其化學(xué)成分上顯示出了巨晶單斜輝石的特點(diǎn)，在n(AlⅣ) —n(AlⅥ)及n(Si)—n(Al)判別圖上也落在了巨晶單斜輝石的范圍內(nèi)(圖3a、b)。這表明了本次研究的單斜輝石是在與巨晶單斜輝石相近的形成條件下形成的，只是由于沒有與巨晶相似的生長條件(不斷的巖漿補(bǔ)充，適當(dāng)位置長時間結(jié)晶)，而沒有形成巨晶，或者如鄢全樹等(2007)解釋的，由于溫壓條件的改變，形成的巨晶發(fā)生碎裂而使得粒度變小。

Wass (1979) 指出單斜輝石的晶出將導(dǎo)致巖漿成分的改變，會使得堿質(zhì)的增高以及巖漿Mg#的降低。因而如果堿性單斜輝石的母巖漿Mg#低于現(xiàn)在的寄主巖，就說明了單斜輝石與現(xiàn)有的寄主玄武巖不存在成分上的平衡性，單斜輝石是其他巖漿形成的，只是由現(xiàn)有的寄主玄武巖帶到地表。反之，如果母巖漿中Mg#高于寄主玄武巖或與寄主玄武巖近似，則表明單斜輝石與寄主玄武巖是平衡的或者同源的。

表4 單斜輝石與熔體之間的分配系數(shù)Table 4 Partitioning coefficients for cpx—melt

注：分配系數(shù)來自Hart and Dunn(1993)和Zack等(1997)。

Kinzler (1997) 的實(shí)驗(yàn)獲得了單斜輝石與熔體的FeO/MgO比值平均分配系數(shù)(0.33)，本文據(jù)此計(jì)算出了與南海新生代玄武巖中的單斜輝石相平衡的母巖漿Mg#的范圍，為0.65～0.79，高于寄主玄武巖的Mg#(0.34～0.59)。 Green(1976)認(rèn)為如果火山巖的Mg#值在0.68～0.72范圍內(nèi)，則代表了該火山巖代表了幔源原始巖漿成分。說明本次研究單斜輝石的母巖漿為相當(dāng)原始的巖漿。鄢全樹(2007)的研究表明只有部分單斜輝石與寄主玄武巖達(dá)到了平衡。鏡下的觀察中，單斜輝石大部分具有很窄的反應(yīng)邊或是被巖漿熔蝕形成渾圓狀的晶形(圖2)，也說明了單斜輝石并未與寄主巖達(dá)到平衡。單斜輝石從母巖漿中的晶出以及之后橄欖石的分離結(jié)晶作用導(dǎo)致了巖漿成分的改變，寄主玄武巖中Mg#急劇下降的同時，單斜輝石與寄主玄武巖的平衡性也被破壞，出現(xiàn)了反應(yīng)邊等現(xiàn)象(圖2a)。

利用Zack 等(1997)以及Hart and Dunn (1993) 獲得的單斜輝石與熔體間的分配系數(shù)(表4)，計(jì)算出了單斜輝石母巖漿的稀土元素組成。結(jié)果顯示母巖漿變化范圍大部分落在單斜輝石的寄主玄武巖范圍內(nèi)，且具有同寄主玄武巖相似的REE標(biāo)準(zhǔn)化曲線，輕重稀土分異明顯，呈輕稀土富集的型式(圖6a)。這進(jìn)一步說明了單斜輝石與寄主玄武巖的同源性 (Woodland and Jugo, 2007; Liotard et al., 1988)，或者是其母巖漿與寄主玄武巖具有演化關(guān)系。

Liotard等(1988)在研究法國中央地體的巨晶單斜輝石時，提出兩種成因，一是巨晶在高壓環(huán)境下從母巖漿中晶出，后被更原始的巖漿帶出地表，故而與現(xiàn)有的寄主巖漿并沒有達(dá)到平衡，二是在巨晶在不同壓力下持續(xù)結(jié)晶，形成了環(huán)帶結(jié)構(gòu)，輝石晶體邊部與與寄主巖達(dá)到平衡。而南海新生代玄武巖中的單斜輝石的寄主巖石Mg#較低，顯然不是原始巖漿，單斜輝石也不存在環(huán)帶現(xiàn)象，與寄主巖石沒有達(dá)到平衡。這說明了本次研究的單斜輝石成因與兩者均不相同，它是在高壓下在與寄主巖同源的巖漿中形成，沒有經(jīng)歷長時間的沉淀生長而直接向上運(yùn)移被帶到地表。因而具有巨晶的主量及微量元素的特征卻不具有巨晶的形態(tài)。上升的過程中由于上升速度較快壓力迅速下降，橄欖石大量晶出 (Green, 1968)，引起了巖漿成分的變化，致使單斜輝石與寄主巖成分和結(jié)構(gòu)上沒有達(dá)到平衡。

圖6 單斜輝石母巖漿的球粒隕石標(biāo)準(zhǔn)化蛛網(wǎng)圖(a)和單斜輝石/寄主巖REE豐度比(b) 球粒隕石數(shù)據(jù)來自 McDonough and Sun(1995)Fig. 6 Calculated CI chondrite- normalized (McDonough and Sun, 1995) REE patterns of parent magma of the clinopyroxenes in the basalts from the South China Sea and REE abundance ratios for clinopyroxene/host pairs

4.3 南海新生代玄武巖中角閃石與單斜輝石的相對分配系數(shù)

Irving and Frey (1984) 統(tǒng)計(jì)了世界范圍內(nèi)的多個巨晶，并計(jì)算了單斜輝石/寄主巖的REE豐度比，可以近似為單斜輝石的分配系數(shù)。但是這項(xiàng)實(shí)驗(yàn)的前提是單斜輝石同寄主巖達(dá)到平衡。鄢全樹等(2007)的研究表明本次研究中的單斜輝石僅有部分達(dá)到了平衡，而單斜輝石/寄主巖的REE豐度比僅有部分落于Irving and Frey (1984) 劃定的范圍內(nèi)(圖6b)，部分比值甚至超過了1，這在玄武巖漿的分離結(jié)晶過程中是不可能出現(xiàn)的，進(jìn)一步說明了單斜輝石同寄主巖之間存在著成分上的不平衡。巖相學(xué)中的觀察則表明著寄主巖同單斜輝石之間存在著結(jié)構(gòu)上的不平衡(圖2a)。故而不能簡單地將單斜輝石/寄主巖石的REE豐度比作為單斜輝石在堿性玄武巖漿中的分配系數(shù)。

S04- 16的單斜輝石/寄主巖REE豐度比明顯不同于其他的幾組樣品，這是由于S04- 16中廣泛發(fā)育有角閃石，角閃石與輝石的共同結(jié)晶，打亂了稀土元素在輝石中正常的賦存規(guī)律(夏群科等，1998)，使得單斜輝石/寄主巖的稀土比值出現(xiàn)了異常。但角閃石與輝石的共同結(jié)晶為我們提供了計(jì)算稀土元素在角閃石和輝石之間的相對分配系數(shù)(表5)的機(jī)會。

樣品S04- 16中角閃石的微量元素組成較為均勻，沒有明顯的環(huán)帶，且大多數(shù)呈渾圓狀，說明在上升過程中受到了熔蝕，而不與寄主玄武巖處于平衡，這也說明了角閃石和輝石形成于巖漿演化的早期。部分單斜輝石中甚至包裹有角閃石，表明單斜輝石與角閃石在同一時期結(jié)晶，可用來計(jì)算微量元素在兩者之間的相對分配系數(shù)。

計(jì)算出的DHb/Cpx中，DLREE較高，DLa可達(dá)2.14，而DHREE較低，DLu僅為0.87。與法國中央地體(Woodland and Jugo，2007)以及West Eifel 地區(qū)(Witt- Eickschen and Harte, 1994)較為一致(圖7)，但是其LREE的分配系數(shù)高而HREE分配系數(shù)低，這可能是源區(qū)成分差異性所致。

表 5 角閃石與單斜輝石的相對分配系數(shù)Table 5 REE partitioning coefficients for amphibole—clinopyroxene

注：*角閃石與單斜輝石均取S04- 16所測得的樣品的平均值。

圖7 S04- 16樣品中角閃石與單斜輝石的相對分配系數(shù)Fig. 7 REE partitioning coefficients for amphibole—clinopyeoxene in sample S04- 16

5 結(jié)論

(1) 南海新生代玄武巖中的單斜輝石具有相似的微量元素組成。稀土元素含量不高，HREE及LREE虧損，而MREE富集；大離子親石元素Rb、Ba、Sr均出現(xiàn)明顯的虧損，高場強(qiáng)元素中Nb、Zr出現(xiàn)虧損而Hf略富集，Nb、Ta與Zr、Hf分餾明顯；Th較為富集，而Pb強(qiáng)烈虧損，U的變化范圍較大。

(2) 南海新生代玄武巖中的單斜輝石較為原始，源區(qū)物質(zhì)成分也具有特殊性。其來源較簡單，為地幔柱的直接產(chǎn)物，沒有受到洋中脊—地幔柱相互作用的影響。

(3) 南海新生代玄武巖中的單斜輝石具有巨晶的特征，與寄主玄武巖同源，沒有經(jīng)歷長時間的沉淀生長而直接向上運(yùn)移被帶到地表，因而具有巨晶的主量及微量元素的特征，卻不具有巨晶的形態(tài)，上升的過程中巖漿成分發(fā)生變化，致使單斜輝石與寄主玄武巖成分和結(jié)構(gòu)上沒有達(dá)到平衡。

參 考 文 獻(xiàn) /References

陳道公, 夏群科, 支霞臣. 1997. 中國東部新生代玄武巖中巨晶礦物的地球化學(xué). 地球?qū)W報(bào), 18(3): 76～82.

白志民. 2000. 北京西山中生代火山巖中單斜輝石礦物化學(xué)及成因意義. 巖石礦物學(xué)雜志, 19(2): 174～187.

黃婉康, 王俊文, 解廣轟, Basu A R. 1991. 玄武巖中普通輝石、石榴石巨晶的研究(Ⅰ)——REE演化. 地球化學(xué), 2: 153～163.

李昌年, 王方正, 鐘稱生. 2005. 廣西北海潿洲島 (含斜陽島) 第四紀(jì)玄武質(zhì)火山巖的地球化學(xué)性質(zhì)及其源區(qū)特征. 巖石礦物學(xué)雜志, 24(1): 1～11.

劉艷榮, 呂新彪, 梅微, 代玉才. 2012. 新疆北山地區(qū)坡北鎂鐵—超鎂鐵巖體單斜輝石的礦物學(xué)特征及其地質(zhì)意義. 巖石礦物學(xué)雜志, 31(2): 212～224.

劉英俊，曹勵明，李兆麟，王鶴年，儲同慶，張景榮. 1984. 元素地球化學(xué). 北京：科學(xué)出版社，50～372.

邱家驤, 廖群安. 1987. 中國東部新生代玄武巖中單斜輝石巨晶的主要特征及成因信息. 巖石礦物學(xué)雜志, 6(1): 56～64.

邱家驤, 廖群安. 1996. 浙閩新生代玄武巖的巖石成因?qū)W與Cpx礦物化學(xué). 火山地質(zhì)與礦產(chǎn), 17(1): 16～25.

邱家驤, 曾廣策. 1987. 中國東部新生代玄武巖中低壓單斜輝石的礦物化學(xué)及巖石學(xué)意義. 巖石學(xué)報(bào), 4: 1～9.

萬玲, 姚伯初, 吳能友, 馬開義, 董英君. 2005. 南海西部海域新生代地質(zhì)構(gòu)造. 海洋地質(zhì)與第四紀(jì)地質(zhì), 25(2): 45～52.

王賢覺,吳明清,梁德華,印愛武. 1984. 南海玄武巖的某些地球化學(xué)特征. 地球化學(xué), 4: 332～340.

王葉劍, 韓喜球, 羅照華, 邱中炎, 丁巍偉, 李家彪, 高水土, 陳榮華. 2009. 晚中新世南海珍貝—黃巖海山巖漿活動及其演化:巖石地球化學(xué)和年代學(xué)證據(jù). 海洋學(xué)報(bào)(中文版), 31(4): 93～102.

夏群科, 陳道公, 王五一. 1998. 巖漿過程中影響礦物稀土元素質(zhì)量分?jǐn)?shù)的因素──來自女山幔源單斜輝石巨晶的啟示. 礦物巖石, 18(2):9～11.

肖森宏. 1990. 巨晶及其寄主玄武巖中稀土元素的豐度——探討巖漿在高壓下稀土元素的晶液分配及演化. 地球化學(xué), 2: 97～107.

肖燕, 張宏福, 范蔚茗. 2008. 華北晚中生代和新生代玄武巖中單斜輝石巨晶來源及其對寄主巖巖漿過程的制約：以莒南和鶴壁為例. 巖石學(xué)報(bào), 24(1): 65～76.

鄢全樹, 石學(xué)法. 2007. 海南地幔柱與南海形成演化. 高校地質(zhì)學(xué)報(bào), 13(2): 311～322.

鄢全樹, 石學(xué)法, 王昆山, 卜文瑞. 2007. 南海新生代玄武巖中單斜輝石礦物化學(xué)及成因意義.巖石學(xué)報(bào), 23(11): 2981～2989.

鄢全樹, 石學(xué)法. 2008. 南海新生代堿性玄武巖中橄欖石的礦物化學(xué)及南海的地幔潛在溫度. 巖石學(xué)報(bào), 24(1): 176～184.

鄢全樹, 石學(xué)法, 王昆山, 卜文瑞, 肖龍. 2008. 南海新生代堿性玄武巖主量、微量元素及Sr—Nd—Pb同位素研究. 中國科學(xué)(D輯), 38(1): 56～71.

鄢全樹, 石學(xué)法. 2009. 南海盆海山火山碎屑巖的發(fā)現(xiàn)及其地質(zhì)意義. 巖石學(xué)報(bào), 25(12):3327～334.

閻全人, 王宗起, 李贈悅, 李繼亮. 2000. 南海北部陸緣地殼結(jié)構(gòu)特征及其構(gòu)造過程. 地質(zhì)論評, 46(4): 417～423.

楊文采, 宋海濱. 2014. 從地殼上地幔構(gòu)造看洋陸轉(zhuǎn)換作用. 地質(zhì)論評, 60(1): 1～21.

姚伯初. 1996. 南海海盆新生代的構(gòu)造演化史. 海洋地質(zhì)與第四紀(jì)地質(zhì), 16(2): 1～13.

于津海, O'Reilly S Y. 2001. 雷州半島英峰嶺玄武巖中的鐵鋁榴石巨晶及母巖漿成因. 科學(xué)通報(bào), 46(1): 6492～497

張漢泉, 吳廬山. 2005. 南海南部海域構(gòu)造地貌. 海洋地質(zhì)與第四紀(jì)地質(zhì), 25(1): 63～70.

趙海玲, 鄧晉福, 李凱明, 李永華. 1999. 東南沿海及南海新生代火山作用與南海的形成演化.地質(zhì)論評, 45: 809～816.

Anders E, Grevesse N. 1989. Abundances of the elements: Meteoritic and solar. Geochimica et Cosmochimica Acta, 53(1): 197～214.

Aoki K I, Kushiro I. 1968. Some clinopyroxenes from ultramafic inclusions in Dreiser Weiher, Eifel. Contributions to Mineralogy and Petrology, 18(4): 326～337.

Aoki K I, Shiba I. 1973. Pyroxenes from lherzolite inclusions of Itinome- gata, Japan. Lithos, 6(1): 41～51.

Briais A, Patriat P, Tapponnier P. 1993. Updated interpretation of magnetic anomalies and seafloor spreading stages in the South China Sea: Implications for the Tertiary tectonics of Southeast Asi. Journal of Geophysical Research: Solid Earth (1978 2012), 98(B4) : 6299～6328.

Cashman K V, Ferry J M. 1988. Crystal size distribution (CSD) in rocks and the kinetics and dynamics of crystallization. Contributions to Mineralogy and Petrology, 99(4) : 401～415.

Damasceno D, Scoates J S, Weis D, Frey F A, Giret A. 2002. Mineral chemistry of mildly alkalic basalts from the 25 Ma Mont Crozier section, Kerguelen Archipelago: constraints on phenocryst crystallization environments. J. Petrol., 43 : 1389～1413.

Damasceno D, Scoates J S, Weis D, Shimizu N, Frey F A, Wallace P. 2000. Implications of trace element abundances in clinopyroxene phenocrysts and basaltic lavas from the Kerguelen Archipelago and Plateau. Spring AGU Meeting, Washington DC. EOS Transactions, American Geophysical Union 81 , S430.

Dieu J J. 1995. 30. Mineral compositions in Leg 144 lavas and ultramafic xenoliths: the roles of cumulates and carbonatite metasomatism in magama petrogenesis[OL]. [2014- 6- 15] http://www-odp.tamu.edu/publications/144_SR/144TOC.HTM

Dobosi G, Jenner G A. 1999. Petrologic implications of trace element variation in clinopyroxene megacrysts from the Nograd volcanic province, north Hungary: a study by laser ablation microprobe- inductively coupled plasma- mass spectrometry. Lithos, 46(4) : 731～749.

Frey F A, Coffin M F, Wallace P J, Weis D, Zhao X, Weis S W, Wahnert V, Teagle D A H, Saccocia P J, Reusch D N, Pringle M S, Nicolaysen K E, Neal C R, Muller R D, Moore C L, Mahoney J J, Keszthelyi L, Inokuchi H, Duncan R A, Delius H, Damuth, J E, Damasceno D, Coxall H K, Borre M K, Boehm F, Barling J, Arndt N T, Antretter M. 2000. Origin and evolution of a submarine large igneous province: the Kerguelen Plateau and Broken Ridge, southern Indian Ocean. Earth and Planetary Science Letters, 176(1): 73～89.

Green D H, Ringwood A E. 1968. The stability fields of aluminous pyroxene peridotite and garnet peridotite and their relevance in upper mantle structure. Earth and Planetary Science Letters, 3(2) : 151～160.

Green D H. 1976. Experimental testing of equilibrium partial melting of peridotite under water- saturated, high pressure conditions . Can. Mineral., 14(3) : 255～268.

Hart S R, Dunn T. 1993. Experimental cpx/melt partitioning of 24 trace elements. Contributions to Mineralogy and Petrology, 113(1) : 1～8.

Irving A J, Frey F A. 1984. Trace element abundances in megacrysts and their host basalts: constraints on partition coefficients and megacryst genesis. Geochimica et Cosmochimica Acta, 48(6): 1201～1221.

Karig D E. 1971. Origin and development of marginal basins in the western Pacific. Journal of geophysical research, 76(11): 2542～2561.

Kinzler R J. 1997. Melting of mantle peridotite at pressures approaching the spinel to garnet transition: Application to mid- ocean ridge basalt petrogenesis. Journal of Geophysical Research, 102(B1): 853～874.

Kushiro I. 1960. Si—Al relation in clinopyroxenes from igneous rocks. American journal of science, 258(8): 548～554.

Lee T Y, Lawver L A. 1994. Cenozoic plate reconstruction of the South China Sea region. Tectonophysics, 235(1～2): 353～388.

Liotard J M, Briot D, Boivin P. 1988. Petrological and geochemical relationships between pyroxene megacrysts and associated alkali- basalts from Massif Central (France). Contributions to Mineralogy and Petrology, 98(1): 81～90.

Marsh B D, Gunnarsson B, Congdon R, Carmody, R. 1991. Hawaiian basalt and Icelandic rhyolite: indicators of differentiation and partial melting. Geologische Rundschau, 80(2): 481～510.

McDonough W F, Sun S S. 1995. The composition of the Earth. Chemical Geology, 120(3): 223～253.

Richter K, Carmichael I S E. 1993. Mega- xenocrysts in alkali olivine basalts; fragments of disrupted mantle assemblages. American Mineralogist, 78(11～12): 1230～1245.

Sun S S, McDonough W F. 1989. Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. Geological Society, London, Special Publications, 42(1): 313～345.

Thompson R N. 1974. Some high- pressure pyroxenes. Mineral. Mag., 39(307): 768～787.

Tu K, Flower M F, Carlson R W, Xie Guanghong, Chen Chuyung, Zhang Ming. 1992. Magmatism in the South China Basin: 1. Isotopic and trace- element evidence for an endogenous Dupal mantle component. Chemical Geology, 97(1): 47～63.

Wass S Y. 1979. Multiple origins of clinopyroxenes in alkali basaltic rocks. Lithos, 12(2): 115～132.

Weis D, Frey F A. 2000. ODP Site 1140 on the Northern Kerguelen Plateau: 34 Ma Pillow Basalts Related to the Kerguelen Plume and the Southeast Indian Ridge, EOS.

Witt- Eickschen G, Harte B. 1994. Distribution of trace elements between amphibole and clinopyroxene from mantle peridotites of the Eifel (western Germany): an ion- microprobe study. Chemical Geology, 117(1): 235～250.

Woodland A B, Jugo P J. 2007. A complex magmatic system beneath the Deves volcanic field, Massif Central, France: evidence from clinopyroxene megacrysts. Contributions to Mineralogy and Petrology, 153(6): 719～731.

Yang Huaijen, Frey F A, Weis D, Giret A, Pyle D, Michon G. 1998. Petrogenesis of the Flood Basalts Forming the Northern Kerguelen Archipelago: Implications for the Kerguelen Plume. Journal of Petrology, 39(4): 711～748.

Zack T, Foley S F, Jenner G A. 1997. A consistent partition coefficient set for clinopyroxene, amphibole and garnet from laser ablation microprobe analysis of garnet pyroxenites from Kakanui, New Zealand. Neues Jahrbuch Fur Mineralogie—Abhandlungen, 172(1): 23～42.

Zou Haibo, Zindler A, Xu Xisheng, Qu Qi. 2000. Major, trace element, and Nd, Sr and Pb isotope studies of Cenozoic basalts in SE China: mantle sources, regional variations, and tectonic significance. Chemical Geology, 171(1): 33～47.