广西涠洲岛晚新生代玄武岩地幔源区及岩浆成因

【摘要】 涠洲岛作为我国最年轻的第四纪火山岩岛,其火山活动表现出多期、多旋回和多喷发中心的特征,但其地幔源区特征和岩浆成因依然存在争议。本文对涠洲岛玄武岩开展了详细的矿物学和全岩主、微量元素及SrNdPb同位素研究,以揭示其地幔源区特征和岩浆成因。涠洲岛玄武岩主要为碱性玄武岩,在岩浆上升过程,几乎未受到地壳物质的混染,经历了橄榄石和单斜辉石的分离结晶作用。轻稀土(LREE)富集、重稀土(HREE)亏损,轻、重稀土强烈分馏((La/Yb)N=1442~2864),Nb、Ta明显正异常,显示出与洋岛玄武岩(OIB)相似的微量元素和SrNdPb同位素特征。SrNdPb同位素比值变化较均一,且呈现出亏损地幔端元(DM)与富集地幔端元(EM2)的二元混合趋势。其中,EM2端元可能源于海南地幔柱。Sr/Sr(101~111)正异常,指示源区存在再循环辉长岩洋壳组分。结合已有的地震层析成像结果和岩石地球化学数据,得出南海及周缘地区的晚新生代玄武岩的形成受控于海南地幔柱。伴随着海南地幔柱的上升,再循环的(121~236)和Eu/Eu本文受中国地震局地质研究所基本科研业务专项(IGCEA1712)资助.第一作者简介:杨文健,男,1994年生,硕士生,矿物学、岩石学、矿床学专业,Email:yangwenjian@ies.ac.cn通讯作者:于红梅,女,1981年生,博士,副研究员,从事火山岩显微构造研究,Email:yuhongmei@ies.ac.cn杨文健等:广西涠洲岛晚新生代玄武岩地幔源区及岩浆成因3902辉长岩洋壳经部分熔融与地幔橄榄岩反应生成石榴石辉石岩(贫硅辉石岩),石榴石辉石岩和未反应的地幔橄榄岩混合部分熔融形成涠洲岛玄武岩。关键词晚新生代玄武岩;地幔源区;岩浆成因;海南地幔柱;涠洲岛中图法分类号P588145玄武岩被称为探测地球深部的“岩石探针”,能够有效的揭示深部地幔信息。南海及周缘地区(包括雷琼、北部湾以及中南半岛等地区)广泛出露的晚新生代玄武岩对揭示整个区域的深部地幔物质组分、地幔热状态以及壳幔相互作用本文结合钻孔(ZK、S2和S3)岩芯样和露头样对涠洲岛玄武岩做了详细的全岩主量元素、微量元素、SrNdPb同位素以及单矿物成分测试,以揭示其地幔源区特征和岩浆成因机制,试图对南海以及周缘地区的晚新生代玄武岩成因提供等具有重要意义。重要信息和约束条件。雷琼及北部湾周边地区为我国南方最大的第四纪火山岩分布区(樊祺诚等,2006),处于南海西北缘(图1a),产生于南海扩张后期。涠洲岛位于北部湾海域内,为我国最年轻的第四纪火山岩岛(图1b)。从1963年广东水文工程地质队首次对全岛的火山岩研究至今,有关该岛火山地质、年代学和岩石地球化学等方面的成果丰硕(刘传章,1986;黎希明和刘传章,1991;卢进林,1993;贾大成等,2003;李昌年和王方正,2004;李昌年等,2005;樊祺诚等,2006,2008a,b;黄林培和李昌年,2007;Lietal,2013)。樊祺诚等(2006)结合火山地质特征和年代学结果将全岛火山活动划分为早期(142~049Ma)溢流式喷发和晚期(36~33Ka)射汽岩浆喷发2个阶段。黄林培和李昌年(2007)对火山碎屑岩中的橄榄岩包体进行温度和压力计算指出寄主碧玄岩的起源深度远远大于包体的稳定深度(~40km)。但是,对于涠洲岛玄武岩的地幔源区特征以及岩浆成因机制仍然存在争议(贾大成等,2003;李昌年等,2005;鄢全树和石学法,2007;樊祺诚等,2008a,b),其核心问题在于对其动力学机制认识的不同。SrNdPb同位素特征指示南海以及周缘地区的晚新生代玄武岩地幔端元组成普遍呈现亏损地幔(DM)与富集地幔(EM2)的二元混合趋势(Yanetal,2014)。但是,对于EM2的来源仍然有下部大陆岩石圈地幔(Tuetal,1991;ZhouandMukasa,1997;樊祺诚等,2008a;韩江伟等,2009)和地幔柱(鄢全树等,2008;ZouandFan,2010;Yanetal,2018,2019)之争。此外,贾大成等(2003)认为北部湾及邻区岩浆活动可能源于在5Ma左右红河断裂由左行走滑转变为右行走滑的构造性质转变。樊祺诚等(2008a,b)通过橄榄岩包体的ReOs同位素和熔体包裹体研究指出雷琼及北部湾周边地区的晚新生代玄武岩起源于岩石圈软流圈地幔相互作用的结果。鄢全树和石学法(2007)认为海南地幔柱主导了南海扩张停止以后南海海盆以及周缘地区的岩浆活动,且整个区域玄武岩的微量元素和同位素呈现与洋岛玄武岩(OIB)相似1地质背景南海处于印度澳大利亚板块、欧亚板块、太平洋板块相互作用叠加区,构造岩浆作用广泛发育,被誉为天然的岩浆构造作用试验场(TaylorandHayes,1983;Briaisetal,1993)。海底磁异常条带揭示南海扩张始于32Ma,止于16Ma(Briaisetal,1993)。由于洋脊抽汲作用在南海海盆周缘地区少见扩张期(32~16Ma)玄武岩(Huangetal,2013),而在南海北缘(海南岛、涠洲岛及雷州半岛等)以及中南半岛地区出现大面积的扩张期后(<16Ma)玄武岩,并显示与OIB相似的地球化学特征(Yanetal,2014)。涠洲岛屹立于南海西北缘的北部湾海域内(图1b),全岛面积约25km22km2(樊祺诚等,2006)。第四纪火山活动贯穿整个造岛过程,形的斜阳岛,共同构成北部湾海域内的一对姊妹火山岛,地貌上北低南高。其与东南角面积不足成一套连续的火山岩地层。下更新统火山岩层由玄武岩、火山集块岩和凝灰岩构成,呈夹层产于湛江组地层中,最大厚度达40m;中更新统火山岩层以玄武岩为主、底部见凝灰角砾岩,厚度介于21~127m之间,与湛江组地层呈不整合接触;上更新统火山岩层顶部为沉凝灰岩、底部为玄武岩,最大厚度达70m,与石峁岭组地层呈平行不整合接触(黎希明和刘传章,1991;卢进林,1993)。李昌年和王方正(2004)根据地层切割关系识别出一套由沉凝灰岩、火山角砾岩和集块岩组成的全新统火山岩地层,与湖光岩组地层呈角度不整合接触。樊祺诚等(2006)综合年代学与火山地质特征,将全岛火山活动划分为早中更新世(142~049Ma)溢流玄武岩喷发和晚更新世末期(36~33ka)射汽岩浆喷发。早中更新世溢流玄武岩构成本岛的底座,主体淹没在海平面以下,在潮间带可见呈球状风化玄武岩出露;晚更新世末期射汽岩浆喷发产物遍及全岛,以火山基浪堆积物为主夹杂不同粒度火山碎的地球化学特征。同时,越来越多的地球物理学资料显示在屑岩,巨厚的基浪堆积层发育交错层理、陷落构造、增生火山海南岛下方存在低速体,并认为海南地幔柱制约着海南岛及邻区的火山活动(LebedevandNolet,2003;Leietal,2009;Weietal,2012;Xiaetal,2016;WeiandZhao,2020),这些认识为重新厘定涠洲岛玄武岩的地幔源区特征以及岩浆成因机制提供了重要科学依据。砾以及爬升层理。已发现有南湾、横路山、大岭、横岭山等多个火口(刘敬合等,1991;卢进林,1993;亓发庆等,2003)。本次取样位置见图1b。样品包括:直接出露的球状风化玄武岩(样品号:NW2、NW5、HLS2、HLS3、18WZ1、18WZ2、18WZ3、18WZ4、18WZ5、NW7、NW15)和火山碎屑4902ActaPetrologicaSinica岩石学报2020,36(7)图1南海及周缘地区晚新生代玄武岩分布图(a,据Yanetal,2018)和涠洲岛地质简图(b,据樊祺诚等,2006修改)ZK、S2和S3分别为钻孔取样位置,其它均为露头取样位置Fig.1ThedistributionoftheLateCenozoicbasaltsintheSouthChinaSeaandsurroundingareas(a,afterYanetal,2018),andthebriefgeologicalmapofWeizhouIsland(b,modifiedafterFanetal,2006)岩中的同源角砾(样品号:WZ1、WZ10、XYD3、XYD4、ZZL3、ZZL4、);钻孔岩芯中的玄武岩(样品号:S33、S37、S39、S311、S24、S214、S220、S230)和火山碎屑岩中的同源角砾(样品号:ZK14、ZK23)。2分析方法MS(NeptunePlus)完成测试分析,具体实验操作流程见文献(Liuetal,2008)。SrNd同位素比值分别采用88Sr/86Sr=146Nd/144Nd=07129进行校正(Linetal,2016)。8375209、其中,Pb同位素比值采用205Tl/203Tl=238714进行校正,同时由于Tl和Pb分馏行为存在一定差异,选择标样NBSSRM981作为外标进行数据二次校正。矿物电子探针分析在中国地质大学(北京)电子探针实验室采用EPMA1600型电子探针完成测试。其中,测试电压15kV,电流1×10-8A,束斑直径1μm,分析精度优于5%。3结果31岩相学和矿物学特征挑选新鲜样品,洗净并在烘箱中烘干,再用玛瑙研磨粉碎至200目粉末。测试全岩主量元素26件、微量元素16件、SrNdPb同位素6件。样品NW7和NW15的主量元素分析在武汉上谱分析科技有限责任公司采用日本理学PrimusⅡX射线荧光光谱仪(XRF)完成测试,其余样品在河北省区域地质矿产调查研究所实验室采用AxiosmaxX射线荧光光谱仪(XRF)完成,分析精度优于5%。全岩微量元素在武汉上谱分析科技有限责任公司采用Agilent7700eICPMS分析完成,分析精度优于5%,具体测试过程见文献(Liuetal,2008)。全岩SrNdPb同位素在武汉上谱分析科技有限责任公司采用美国ThermoFisherScientific公司的MCICP玄武岩整体呈灰黑色,隐晶质结构,以块状构造为主、少数见气孔构造。显微镜下玄武岩呈斑状结构,斑晶主要为橄榄石和少量单斜辉石(图2a,b)。橄榄石斑晶含量约5%~10%,呈粒状,半自形,大小主要为50~250μm;单斜辉石斑晶含量少于5%,大小为50~200μm;基质为间粒结构,不规则排列的长条状针状斜长石中间充填有橄榄石、辉石微晶和钛铁氧化物,斜长石呈长条状、针状,自形半自形,大小介于30~300μm。火山碎屑岩中同源角砾的显微结构与前者存在差异,显微镜下呈玻基斑状结构,斑晶主要为橄榄石(图2c,d),含量约5%~10%,呈粒状,半自形自形,大小主要为50~300μm;少量单斜辉石斑晶,含量1%~3%,呈粒状,杨文健等:广西涠洲岛晚新生代玄武岩地幔源区及岩浆成因5902表1涠洲岛玄武岩橄榄石电子探针分析结果(wt%)Table1ElectronicalmicroprobeanalysesofolivinesfrombasaltsinWeizhouIsland(wt%)NiObdl44820260071315023TiO2Al2O3Cr2O3MgOSiO2MnOFeO测点号CaOTotalbdl00345660270221332024001Na2OK2OFoS230114072bdl006bdl00299338586S2301239280030010064607027027132802600200299588608S230213977bdl0010104500027024139802600200199678515S230413970014bdl99588594S230424070bdl0090014484029020134502700300699938559S23043392900300600144060290161463022002bdl98778429S241139360010060074474023020144702000500199408464bdl44420240191445027003002100808457S24124117001002S24213984bdl002bdl4526023018137302400200199538546bdl100648459S24224046bdl00200244850280221456022001bdl99227949bdlS2423401100700200739840280271832026bdl100238271S24313984001004bdl43470230201619024bdlS2432403201200701043160280231564028001bdl100208310S24333929bdl00400239820290312048012003001100417761S2461410400600600744030250181429032002001100348459bdl10bdl8506S2471409400400200344460220161392022bdlS24724102bdl00300744400240191440024001bdl100598461bdl00299747808S24733997bdl00601239190300271960021S24813842005005bdl37990260272159014bdl0019877758218WZ2114066003003006445302601813490250050039955854718WZ2123962bdl005007411202203217600180060029925806318WZ2213986bdl0030024005022039195502300200110037785018WZ2514045bdl006009448202301513340260030039946856918WZ252413400100200445080210141316032001bdl10033859218WZ261407800400500844050230211392036002bdl99748494bdl4056024023186801800900210031794618WZ262401800201218WZ2714113bdl00700443260210291531025bdl00110056834318WZ27241250030080014112023026169203000100210023812418WZ2733974bdl00600744200210281464031bdl0019952843218WZ1214043005010001440702202313650270020029906852018WZ1224028bdl006001436402502414330230010029906844518WZ131404000100200445080230191407033003bdl10041851018WZ13241200050080054001030020173701800100299478041bdl00110026855918WZ1814101bdl0020044502022019135002518WZ1824113005bdl100258536bdl432202202015370280010019955833718WZ1834019003002bdl00110002853118WZ3113964bdl0050104559021019139902518WZ312407000500300644640230171413030bdl00110032849218WZ3214122bdl003009432102101214520270020029972841318WZ3224132bdl01500443450260131458022007001100248415bdl18WZ3313897005005bdl44740230151429030bdl9880848018WZ3414119bdl00300545440240151373033bdl00110114855118WZ3514137bdl00200345270240191355035001bdl10103856218WZ3524010009004007411802201916780240040039896813918WZ3614006bdl00300542590220241665029bdl001100158201bdl002447102001513670310026902ActaPetrologicaSinica岩石学报2020,36(7)续表1ContinuedTable1SiO2测点号CaOMnOFeOTotalTiO2Al2O3Cr2O3MgObdl44720260181385031Na2OK2ONiOFo18WZ381399000500300939660250301876012bdl00199177903XYD4113999bdl00400144750250201367025003bdl99188537XYD4124051bdl004bdl00199868520XYD41340570030080044362027018141402100200199168461bdl99468780XYD4214079bdlbdl00546500190091152032001XYD4224151bdl00400245050250181321025005bdl100578587XYD423407500300300345180260231347030002002100318567XYD461405600200400944960220151374022001bdl100028536S391140570060110014331028025146501800200299468404S39124042bdl00300144700290111396020bdl99728509S39314132bdl00600644830260241355034001001100678549bdl001100797881S39334083bdl005bdl99608559S39514060bdl00400244760270171342031001S3952402000400100644010260171431023bdlbdl99298457S396140350090100054471027015138002800100299828523bdl4472027014137502600100199668529S39624044002005bdl39930280321913024bdl注:bdl(belowthedetectionlimits),表示低于检测限;表2、表3同图2涠洲岛玄武岩显微镜照片(a、b)分别为钻孔玄武岩(S214)正交偏光和单偏光下的显微照片,基质为间粒结构;(c、d)分别为钻孔火山碎屑岩中同源角砾(ZK14)正交偏光和单偏光下的显微照片,呈玻基斑状结构Ol橄榄石;Cpx单斜辉石;Pl斜长石;G玻璃Fig.2MicrophotographsofbasaltsinWeizhouIsland杨文健等:广西涠洲岛晚新生代玄武岩地幔源区及岩浆成因7902表2涠洲岛玄武岩单斜辉石电子探针分析结果(wt%)Table2ElectronicalmicroprobeanalysesofclinopyroxenesofthebasaltsinWeizhouIsland(wt%)测点号SiO2TiO2Al2O3Cr2O3MgOCaOMnOFeONiONa2OK2OTotalWoEnFsS23024514818438901414682188010641007038bdl10087458342761141S2302552931373140291485219901359500103300310102453042561214S23026517219639501514362213012627005037bdl10106462041711209S2451505820449806814462089010603007038bdl10022447043041226S24525306140268024155420780126390020320011005542654438129718WZ283513118042603014611938013658004041bdl988041454346150918WZ284525713928301115291952015720bdl035001994240664431150318WZ151532213431304215221991021603bdl033001998141274385148818WZ411538511624302415692086011611008031bdl1008442484443130918WZ412537211125802415472095010626bdl033bdl1007642794395132618WZ413536513922500815731989016696013028bdl1005340534458148918WZ421533211125305116302099012573004030bdl1009442814622109718WZ4224987274449bdl14012040023775bdl045bdl999443774180144318WZ423539111922701915902099010656007033bdl1015242524480126718WZ431533811526102215682068016613007033bdl1004242454476127818WZ432526315223600815422021018703012034bdl998941824436138218WZ433533713828402415442082014620004034bdl10080427044041327表3涠洲岛玄武岩斜长石电子探针分析结果(wt%)Table3ElectronicalmicroprobeanalysesofplagioclasesofthebasaltsinWeizhouIsland(wt%)测点号AbSiO2TiO2Al2O3Cr2O3MgOCaOMnOFeONiONa2OK2OTotalAnOr18WZ23157850182625bdl018882002049bdl53607599904542499746218WZ24154170152899006012110800403900744003399795702409420418WZ242544100729040010121105005049bdl43003999935728403423818WZ243558800128190020131038005037bdl468046100165353436428318WZ28156090132793bdl0121016005046bdl484046100235223449628118WZ28254210102907bdl0091102bdl04600742503699635758401722618WZ11155440112893bdl0151047bdl056001452050100695439425131118WZ11256580132764003012898bdl05500254404499914641509126818WZ11355730162814007010982bdl04800348104599795155456628018WZ14154660152878bdl0181101003053bdl432040100055705405024618WZ16154370152906bdl0111118bdl03100542403799825798397522718WZ16253020132880bdl0201137001038bdl42803598555819396521618WZ163549101228760010111098bdl046003450040100275600415824218WZ171545700328800050121088004038bdl44504099725608414724518WZ44153880092932bdl0101136bdl054bdl41303099715922389518418WZ442541101729170030141158003049006405034100165999379420718WZ443544301428890030121090bdl0490094490379995560041752248902ActaPetrologicaSinica岩石学报2020,36(7)表4涠洲岛玄武岩主量元素分析结果(wt%)Table4Majorelement(wt%)analysesofbasaltsinWeizhouIsland样品号SiO2TiO2Al2O3Fe2O3FeOMnOMgOCaONa2OK2OP2O5LOIH2O+H2OTotalMg#NW24749240120553458401810209432851830781340920849973688NW54652223116553352802011108632921960783112071459973720HLS24479259127055060601696810322051800843052491029955669HLS34585260126750563501795610163660870871901450639970664ZZL34581256119344966801810219973172050841801301329970677ZZL444532661247371785017105010153822040870920750409970663XYD343172271135557578018114210203711481033512882979968711XYD4467220511204085700169841017353173088367128080997270018WZ149862551386374678015738885354191062050018021997560818WZ249322211373405635015843942316159052082051036997564818WZ349682151381417610015807955326156051073027038997464318WZ449862241455534433012633955305177054209141156997762618WZ54551247122644068901810169803541900811781361179971671WZ1462924712434137160189979713582000810990810239971664WZ104867223133870433202372410372591710542441661769977665ZK144461238119261248501797611231591660824553452469965696S334628257132783232201490410292061700722111371809972698S37481222413173776840169839252491450571891401319976672S394641249130572240901497110322081730701781481269972702S31146202411268340738016102610052511830672201471189972672S24482822513154466190189019562681520581891461109975662S2144678249129036872501610139673641660720650110249972670S2204561246123865045501610809992371420722762371289973720S2304731241129547262001699097931317606507404702299736799995689NW74673233119411895900191031974316158074134NW1546532421240123038201897595925016608021110023700注:样品NW7和NW15中的Fe2O3表示Fe2O3T,FeO含量通过滴定法测得半自形自形,大小为50~200μm;发育气孔;基质主要为黑品(18WZ2、18WZ3、S24、S230、XYD4)中的橄榄石均落在色玻璃质。橄榄石成分见表1,其Fo[Fo=100Mg/(Mg+Fe)]值为758~878,NiO含量为012%~036%,且Fo值与NiO呈正相关趋势,明显区别于地幔橄榄石演化趋势(Sato,1977)(图3a)。所有样品的CaO含量为019%~030%,都高于地幔捕虏晶(CaO<01%)范围,应属岩浆成因(ThompsonandGibson,2000)(图3b)。MnO含量为007%~039%,随着橄榄石Fo值降低,MnO含量升高,呈负相关趋势(图3c)。与橄榄石Fo显示(图3d),当橄榄石与熔体达到平全岩Mg#衡时,数据点将集中落在平衡曲线上;若数据点落在平衡曲线的下方或上方,则显示橄榄石来源于富铁的基质或岩浆演化过程中早期结晶的富镁斑晶。在图3d中,样品S39中的橄榄石均落在平衡曲线的下方,指示富铁的基质来源;而样品18WZ1中的橄榄石同时落在平衡曲线上方、中间和下方,分别显示富镁斑晶、平衡结晶和富铁基质来源特征;其他样平衡曲线中间和下方,显示平衡结晶和富铁基质来源特征。单斜辉石成分见表2,由次透辉石和普通辉石构成(Wo405462En4174622Fs109151),其中TiO2含量为111%~274%,Al2O3含量为225%~498%。斜长石成分见表3,主要由中拉长石(An455600)组成。32全岩主、微量元素特征全岩主、微量元素成分见表4和表5。其中,SiO2含量为4317%~4986%,MgO含量为633%~1142%,Mg#606~720。在TAS图中(图4),几乎所有样品都集中落于碱性玄武岩范围内。在哈克图解中(图5),MgO与SiO2、Al2O3呈良好负相关,与Fe2O3、MnO、CaO、CaO/Al2O3具正T相关关系,与TiO2、Na2O、K2O相关性不明显。在球粒陨石标准化稀土元素配分曲线上(图6b),呈现值为杨文健等:广西涠洲岛晚新生代玄武岩地幔源区及岩浆成因9902表5涠洲岛玄武岩微元素分析结果(×10-6)Table5Traceelement(×10-6)analysesofbasaltsinWeizhouIsland18WZXYD3518WZ3XYD4样品号S39S24WZ1WZ1ZK14S37NW2S214ZZL4HLS318WZ218WZ1S230Sc180184186160152168189191183186181166189207188190196177192194186153178174168183184182168169186172182183VCr264287259312243206297291269271270218278329281306315Co524541556532473451517486532533535507506522503528532Ni229234222283213160210194225222224207211246210235239Rb38434639233235937132429532842341942727432128633735191210581021114010317967457509909939921843729930840871912SrY221222225228209324235255223229224211196215200213204Zr215228231215190190160155219218216206163194171200187Nb684781792813636527438428721728727703471605502636574Ba493606548516511501440580532534528509413486423495474Hf495538514472442435381378488505506480396466405482461Ta386439451429345301252245415405404394263342283355320Pb476372340502501298243235348465440263232284245284263Th640679663747669453366370641653653619401498427531480U152157170162160078091060159166157132103081103124124La430477474547450349279284443447442447293372312387352Ce816901889103854632531522843858851853572718600744676Pr957105104117980784659638987101983100683857733895807Nd379409412460391325272263390392389387276343286353320Sm755821864919796714599564793796748771597700624730673Eu245265265294246234204194253263259247198250210239229Gd675713738787694701576552732707745700567679589675642Tb095102102111095100084079099104100099083092088093092Dy490500486505460519414390485483467449407443409456414Ho081081081080073093074070078083082076068075070076072Er194195196207192239181186198201199188178189182188185Tm024024024026022031024023027024025023021022023023023Yb148131127130119164123114142130130123125127131131130Lu018019021017017021018017019021019018018017017017018∑REE19932377216924622064166613781352205720792058205614361778150618361677∑LREE18212201199222751897147912281209187919041881188912891614135516701519∑HREE17251765177518631672186814941431178017531767167614671644150916591576∑LREE/∑HREE1055113311221221113579282284510561086106511278799828981007964(La/Yb)N19722465253528642570144215421697211423422308246215851995161920061839(Tb/Yb)N294354367390363276310316319367352366301333307323322Eu/Eu105106102106102101106106102107106103104111106104107Sr/Sr121128124122131129144149127126127236135138149125140注:WZ1为WZ1样品的重复测量结果;N表示原始地幔标准化,标准化数值据SunandMcDonough(1989)0012ActaPetrologicaSinica岩石学报2020,36(7)图3橄榄石斑晶的成分变化特征(a)橄榄石Fo与NiO含量关系图解,地幔橄榄石趋势线据Sato(1977),夏威夷橄榄石和普通橄榄石组分变化区域(据Sobolevetal,2005),海南岛玄武岩中的橄榄石成分变化区域(据Wangetal,2012);(b)橄榄石Fo与CaO含量关系图解(底图据ThompsonandGibson,2000);(c)橄榄石Fo与MnO含量关系图解(底图据Wangetal,2012);(d)橄榄石Fo与全岩Mg#关系图解,FeMg分配系数KD据RoederandEmsile(1970)Fig.3Variationsinthecompositionsofolivinephenocrysts轻稀土(LREE)富集、重稀土(HREE)亏损的右倾模式,(La/Yb)N=1442~2864,不存在Eu明显负异常。在原始地幔标准化微量元素蛛网图上(图6a)具有明显的NbTa正异常。同时,与南海海山玄武岩、越南玄武岩以及海南岛玄武岩相比,均显示与SunandMcDonough(1989)提出的洋岛玄武岩(OIB)相似的微量元素特征。33全岩同位素特征全岩SrNdPb同位素组成见表6。其中,位素组成变化特征。在0703492~0703741,208Pb/204Pb=385574~389179,156525,87Sr/86Sr=143Nd/144Nd=0512914~0512977,207Pb/204Pb=1555871~206Pb/204Pb=184361~187327,显示相对均一的同87Sr/86Sr143Nd/144Nd图中(图7a),所有的样品均落在Staudigeletal(1984)定义的洋岛玄武岩(OIB)范围内,并与中南半岛、雷琼半岛以及南海海山玄武岩呈现相似的同位素组成,呈现出DM与EM2的二元混合趋势(鄢全树等,2008;石学法和鄢全树,2011;徐义刚等,2012;Yanetal,2014)。在143Nd/144Nd206Pb/204Pb图中(图7b),图4涠洲岛玄武岩TAS分类图解(底图据LeBasetal,1986)文献数据:李昌年等,2005;樊祺诚等,2008b;Lietal,2013后图同Fig.4TASdiagramofbasaltsinWeizhouIsland(basemapafterLeBasetal,1986)杨文健等:广西涠洲岛晚新生代玄武岩地幔源区及岩浆成因1012图5涠洲岛玄武岩哈克图解Fig.5HarkerdiagramsforbasaltsinWeizhouIsland表6涠洲岛玄武岩SrNdPb素分析结果Table6SrNdPbisotopicelementcompositionsofbasaltsinWeizhouIsland2σ208Pb/204Pb87Sr/86Sr2σ样品号143Nd/144Nd2σ2σNW20703667000001105129140000011388614000181563240000718686600007ZZL40703534000000605129540000012387123000151558710000618609700007XYD30703492000000905129590000008385574000151559360000618436100006ZK230703637000001205129380000012389179000181565250000718732700008S370703545000000705129770000014387696000221560060000818654700010S2300703741000000605129440000004387338000221559120000918621900009207Pb/204Pb2σ206Pb/204Pb也显现出DM与EM2的混合趋势。在206Pb/204Pb207Pb/204Pb206Pb/204Pb208Pb/204Pb图中(图7c,d),所有的样品均位于北半球参考线(NHRL)上方,显示Hart(1984)提出的Dupal和同位素异常特征。4讨论41岩浆演化过程411地壳混染情况产于大陆板内的幔源岩浆在上升到地表的过程中,可能会受到地壳物质不同程度的混染。涠洲岛玄武岩中地幔橄榄岩包体的存在指示岩浆快速上升且几乎不与地壳发生同化混染作用(黄林培和李昌年,2007;樊祺诚等,2008a)。在微量元素蛛网图上(图6a),涠洲岛玄武岩显示NbTa正异常,暗示不存在地壳物质的混染。同时,在NbNb/U图上(图8a),大多数样品均落于OIB(47±10;Hofmannetal,1986)范围,远高于平均大陆地壳(Nb/U=615;RudnickandGao,2003);在(La/Nb)N(Th/Nb)N图上(图8b),所有样品(La/Nb)N和(Th/Nb)N比值都明显低于下地壳组成((La/Nb)N=17,(Th/Nb)N=20,RudnickandGao,2003),但极2012ActaPetrologicaSinica岩石学报2020,36(7)图6涠洲岛玄武岩原始地幔标准化微量元素蛛网图(a,标准化值据SunandMcDonough,1989)和球粒陨石标准化稀土元素配分图(b,标准化值据McDonoughandSun,1995)洋岛玄武岩(OIB)和洋中脊玄武岩(EMORB和NMORB)组分据SunandMcDonough(1989);越南玄武岩据Anetal(2017);海南岛玄武岩据Wangetal(2012);南海海山玄武岩据鄢全树等(2008)Fig.6Primitivemantlenormalizedtraceelementdiagrams(a,normalizedvaluesafterSunandMcDonough,1989)andchondritenormalizedrareearthelementdiagrams(b,normalizedvaluesafterMcDonoughandSun,1995)ofbasaltsinWeizhouIsland其接近OIB范围,进一步排除了涠洲岛玄武岩受地壳混染的影响(Daietal,2011)。此外,与整个南海及周缘地区的晚新生代玄武岩几乎不受地壳物质的混染具有一致性(Yanetal,2014;Anetal,2017;Hoangetal,2018)。综上所述,涠洲玄武岩几乎不受地壳物质的混染影响,其微量元素、同位素特征代表了其地幔源区特征。412岩浆结晶分异涠洲岛玄武岩Mg#=606~720,Ni=160×10-6~283×10-6(表4、表5),与原始岩浆,Cr=206×10-6~329×10-6(Mg#>70,Ni>400×10-6~500×10-6,Cr>1000×10-6(Freyetal,1978;WilkinsonandLeMaitre,1987)相比,仅部分样品的Mg#(>70)接近原始岩浆,而大部分玄武岩Mg#低于原始岩浆成分,说明部分岩浆经历了结晶分异过程。在哈克图解中(图5),MgO与SiO2、Al2O3呈良好负相关,与Fe2O3、MnO、CaO、CaO/Al2O3具正相关关系,与TiO2、Na2O、TK2O相关性不明显,说明可能发生了橄榄石和单斜辉石结晶分异。在稀土元素配分曲线上(图6b)不存在明显Eu负异常,暗示不存在斜长石的分离结晶。同时,在相容元素(Cr、Co、Ni、Sc)与MgO关系图中(图9),MgO与Ni、Co呈良好正相关,与Cr、Sc不显示相关性,也说明岩浆经历了橄榄石的结晶分异。岩石地球化学特征和岩相学观察到橄榄石和单斜辉石斑晶(图2)具有一致性,综合表明涠洲岛玄武岩经历了橄榄石和少量单斜辉石的结晶分异作用。42地幔源区特征421源区岩石学特征涠洲岛玄武岩呈现强烈的轻、重稀土分馏,ΣLREE/ΣHREE=792~1221,(La/Yb)N=1442~2864,暗示其源区可能存在石榴石。Wangetal(2002)指出(Tb/Yb)N比值可以约束地幔源区特征,涠洲岛玄武岩均显示高(Tb/Yb)N比值(276~390),暗示其源区存在石榴石((Tb/Yb)N>18;Wangetal,2002)。此外,在熔体/石榴石中DDy/Yb=026,而在熔体/尖晶石中DDy/Yb=1(MckenzieandONions,1991),因此Dy/Yb比值可以用来约束地幔源区矿物。微量元素模拟结果(图10)显示尖晶石二辉橄榄岩和尖晶石石榴石二辉橄榄岩具有低Dy/Yb比值(前者小于15,后者小于27),而涠洲岛玄武岩具有高Dy/Yb(312~390)比值,指示源区存在石榴石。同时,本区玄武岩并不落在单一的模拟的部分熔融曲线上,而是介于模拟的石榴石辉石岩和石榴石二辉橄榄岩熔融曲线之间(图10),暗示其源区岩石可能由石榴石橄榄岩和石榴石辉石岩的混合组成(Anetal,2017;Hoangetal,2018;Kimetal,2019)。实验岩石学表明橄榄岩部分熔融产生的熔体与很多具有OIB特征的玄武岩出现解耦(HiroseandKushiro,1993;HiroseandKawamoto,1995;Dasguptaetal,2007;Davisetal,2011),而辉石岩、角闪石岩、碳酸盐化橄榄岩、榴辉岩、碳酸盐化榴辉岩等均可以作为碱性玄武岩的源区母岩(Hirschmannetal,2003;Sobolevetal,2005,2007;Dasguptaetal,2007;Piletetal,2008;Chenetal,2009;Zengetal,2010)。全岩和橄榄石的主、微量元素以及非传统稳定同位素(如Mg同位素等)研究暗示,南海以及周缘地区的晚新生代玄武岩的源区存在辉石岩和碳酸盐化橄榄岩(Wangetal,2012;Liuetal,2015;Yanetal,2015,2018;Anetal,2017;Lietal,2017;Zhangetal,2017;Hoangetal,2018)。杨文健等:广西涠洲岛晚新生代玄武岩地幔源区及岩浆成因3012图7涠洲岛玄武岩SrNdPb同位素比值图解不同地幔端元(DM、EM1、EM2、HIMU)据ZindlerandHart(1986);洋岛玄武岩(OIB)范围据Staudigeletal(1984);北半球参考线(NHRL)据Hart(1984);中南半岛玄武岩据Hoangetal(1996),ZhouandMukasa(1997),Anetal(2017),Hoangetal(2018),Yanetal(2018);南海玄武岩据Tuetal(1992),鄢全树等(2008);雷琼半岛玄武岩据Tuetal(1991),ZouandFan(2010),Wangetal(2013),韩江伟等(2009),朱炳泉和王慧芬(1989)Fig.7SrNdPbisotoperatiodiagramsofbasaltinWeizhouIsland通常碳酸盐化地幔橄榄岩部分熔融形成的岩浆具有明显Ti、Zr、Hf负异常(Zengetal,2010),但在微量元素蛛网图(图7a)上涠洲岛玄武岩呈现Zr、Hf正异常及Ti弱负常,故排除源区为碳酸盐化地幔橄榄岩。Herzberg(2011)指出全岩CaO含量能够有效的识别源区岩石。为了消除单斜辉石分离结晶的影响,剔除MgO<9%的样品。在CaOMgO图(图11a)中,所有的样品均落在辉石岩熔体区域内。Herzberg(2011)认为辉石岩熔融产生的熔体具有高FeOT/MnO(>60)比值,而橄榄岩熔融产生的熔体具有低FeOT/MnO(<60)比值。在图FeOT/MnOMgO图(图11b)中绝大多数样品集中落在辉石岩熔体区域内,暗示其源区存在辉石岩。此外,与Gaffneyetal(2005)提出混合的橄榄岩辉石岩源区形成的夏威夷玄武岩相比显示相似的CaO含量和FeOT/MnO比值(图11),推测涠洲岛玄武岩也可能具有混合的源区特征。Sobolevetal(2005)认为辉石岩和地幔橄榄岩的混合部分熔融可以形成夏威夷玄武岩中高Ni橄榄石斑晶,并强调再循环的榴辉岩洋壳经部分熔融后交代周围的橄榄岩生成二阶段辉石岩,所形成的的二阶段辉石岩与残存的地幔橄榄岩混合发生部分熔融形成我们所观察到的夏威夷玄武岩(Sobolevetal,2007)。在图3a中,绝大多数橄榄石也都落在夏威夷玄武岩中的橄榄石成分区域内,并与海南玄武岩中橄榄石显现很高的相似性(Wangetal,2012),进一步说明涠洲岛玄武岩的源区母岩具有石榴石辉石岩和地幔橄榄岩混合特征。此外,Strackeetal(2003)指出深海辉长岩作为再循环洋壳的独特组分具有典型的Sr正异常特征,故Sr正异常可以识别再循环(辉长岩)洋壳。Yuetal(2010)对蛟河辉石岩包体研究指出Sr、Eu正异常反映了源区存在富斜长石堆积的辉长岩洋壳。在微量元素蛛网图(图6a)中,涠洲岛玄武岩也显现出Sr(Sr/Sr=121~236)和Eu(Eu/Eu=101~111)正异常,说明了其源区存在再循环洋壳(辉长岩)组分。4012ActaPetrologicaSinica岩石学报2020,36(7)图8涠洲岛玄武岩NbNb/U(a,底图据Hofmannetal,1986)和(La/Nb)N(Th/Nb)N(b,底图据Daietal,2011)关系图解大陆地壳Nb/U比值、LC(下地壳)、MC(中地壳)和UC(上地壳)的(Th/Nb)N及(La/Nb)N比值据RudnickandGao(2003);N表示原始地幔标准化,标准化值据SunandMcDonough(1989)Fig.8DiagramsofNbvsNb/U(a,afterHofmannetal,1986)and(La/Nb)Nvs(Th/Nb)NforWeizhouIslandbasalts(b,afterDaietal,2011)HofmannandWhite(1982)指出在地幔柱中普遍存在再循环洋壳,伴随着地幔柱上升携带的再循环洋壳经部分熔融产生的熔体与橄榄岩反应生成辉石岩(Sobolevetal,2005)。综上所述,我们认为涠洲岛玄武岩源自石榴石辉石岩和地幔橄榄岩混合的部分熔融。422地幔端元组成南海以及周缘地区广泛分布的晚新生代玄武岩具有与OIB相似的微量元素和同位素地球化学特征(鄢全树等,2008;石学法和鄢全树,2011;徐义刚等,2012;Yanetal,2014),暗示着它们可能源于同一个源区(鄢全树和石学法,2007)。涠洲岛玄武的SrNdPb同位素变化范围较窄(图7),表明其具有较为均一的地幔源区特征。此外,前人研究指出中国东部新生代玄武岩SrNdPb同位素组成具有区域性变化,在东南沿海地区普遍显示DM与EM2混合的特征(Liuetal,1994;Zouetal,2000)。在SrNdPb同位素图中(图7),涠洲岛玄武岩也呈现出DM与EM2的二元混合趋势,与南海及周缘地区的晚新生代玄武岩具有一致性(鄢全树等,2008;石学法和鄢全树,2011;徐义刚等,2012;Yanetal,2014)。在微量元素蛛网图和稀土元素球粒陨石标准化图上(图6),呈现出明显NbTa正异常和LREE富集、HREE亏损的OIB地球化学特征,而明显区别来源于软流圈地幔的洋中脊玄武岩(EMORB和NMORB),推测涠洲岛玄武岩中的DM组分代表了海南地幔柱自身的特征。然而,对于EM2al(1992)认为在缺乏地幔柱的作用下南海及其周缘玄武岩中EM2组分来源于下部大陆岩石圈地幔(SCLM)。然而,层析成像技术揭示了海南地幔柱的存在(LebedevandNolet,2003;Montellietal,2006;Zhao,2007;Leietal,2009;Weietal,2012;Xiaetal,2016),海南地幔柱不仅对南海及其周缘的岩浆活动提供了热量,还提供了物质组分。Anetal(2017)基于全球代表SCLM的地幔包体统计发现其NdHf同位素具有明显的解耦现象,不呈现线性正相关趋势。而南海及周缘的晚新生代玄武岩的NdHf同位素呈正相关趋势(Yanetal,2015),暗示EM2可能不是来源于SCLM。在微量元素蛛网图上(图6a),涠洲岛玄武岩显示NbTa正异常,也排除了EM2源于SCLM的可能性,说明了EM2最可能源于海南地幔柱。43岩浆成因模型从化学成分来讲,地幔柱成因的火山岩具有OIB的特点,但显示OIB特征的岩浆未必源自地幔柱活动的产物(徐义刚等,2012)。中国东部新生代玄武岩与南海及其周缘地区的晚新生代玄武岩均显示OIB特征,但二者的成因机制存在显著差异。徐义刚等(2018)指出在地幔过渡带中滞留的太平洋板块脱碳、脱水以及熔融和交代作用导致了中国东部新生代玄武岩显示OIB特征,而南海以及周缘地区的晚新生的来源认识存在争议,一部分学者认为源于下部大陆岩石圈地幔或岩石圈地幔(Tuetal,1991;韩江伟等,2009;Anetal,2017),另一部分学者认为源于下地幔或地幔柱(Yanetal,2008,2015;ZouandFan,2010;Wangetal,2013)。Tuet代玄武岩受控于海南地幔柱的活动(鄢全树和石学法,2007)。从大地构造背景来看,涠洲岛玄武岩与南海及周缘扩张期后(<16Ma)玄武岩具有一致的构造背景特征,均集中落于地幔柱玄武岩区域内(徐义刚等,2012)。大量的地球杨文健等:广西涠洲岛晚新生代玄武岩地幔源区及岩浆成因5012图9涠洲岛玄武岩相容元素(Ni、Cr、Co、Sc)与氧化物(MgO)关系图解Fig.9Diagramsoftherelationshipbetweencompatibleelements(Ni,Cr,Co,andSc)andoxide(MgO)forWeizhouIslandbasalts物理学资料表明在海南岛的下方存在低速异常,指示海南地幔柱的存在(LebedevandNolet,2003;Montellietal,2006;Zhao,2007;Leietal,2009;Xiaetal,2016)。海南地幔柱活动不仅导致了南海及周缘地区的地幔显示高温热异常(Wangetal,2012;Anetal,2017;Yangetal,2019),而且还形成了整个区域分布面积超过4×106km2具有OIB特征的晚新生代玄武岩(鄢全树和石学法,2007)。区域和全球层析成像显示海南地幔柱起源于下地幔(LebedevandNolet,2003;Montellietal,2006;Zhao,2007;Weietal,2012;Xiaetal,2016),与夏威夷地幔柱成因玄武岩相比,其具有异常高230Th含量,表明它呈慢速上升(ZouandFan,2010)。Xiaetal(2016)根据层析成像结果提出海南地幔柱呈双层上升模式(Weietal,2012;WeiandZhao,2020),即在地幔过渡带地的2009)。Zhangetal(2006)指出位于峨眉山地幔柱轴部的丽江地区显示最高的热异常(T=1630~1680℃)。进一步研究认为,随着距离轴部的距离增加,呈现温度递减的趋势(李永生,2012;汪云峰等,2013)。原始岩浆组分被广泛用来估算地幔潜在温度(Tp)及熔融条件(McKenzieandBickle,1988;Putirka,2005;Zhangetal,2006;Putirkaetal,2007;HerzbergandAsimow,2008,2015)。根据估算的原始岩浆计算获得处于海南地幔柱中心的海南岛显示高地幔潜在温度(Tp=1500~1580℃,Wangetal,2012),而处于地幔柱边部的越南地区显示较低地幔潜在温度(Tp=1468~1490℃,Anetal,2017),但是由于它们的岩浆源区存在辉石岩(Wangetal,2012;Anetal,2017),故计算值将比实际地幔潜在温度偏高(HerzbergandAsimow,2008)。但以上这些计算结果在幔柱头发生横向扩展呈蘑菇状,随着柱头的持续上升到达岩一定程度依然可以反映海南地幔柱热结构的差异,且与峨眉石圈底部时分散呈布丁状,形成独立的岩浆房或地幔柱分支(Xiaetal,2016;Fanetal,2017)。山地幔柱热结构相似,均显示轴部温度高、边部温度低的特征(Zhangetal,2006;李永生,2012;汪云峰等,2013),由此由于地幔柱热结构存在差异,核部温度远高于边缘位置(Renetal,2005;Zhangetal,2006;HerzbergandGazel,可以推测出涠洲岛的地幔潜在温度介于海南岛和越南地区之间。此外,涠洲岛与南海、雷琼半岛以及中南半岛等地区6012图10涠洲岛玄武岩YbDy/Yb图解(底图据Anetal,2017;Yanetal,2018)石榴石二辉橄榄岩、尖晶石石榴石二辉橄榄岩和尖晶石二辉橄榄岩非实比熔融曲线,以及榴辉岩(Cpx82Grt18和Cpx75Grt25)实比熔融曲线据Anetal(2017);石榴石辉石岩(Opx5Cpx45Grt50)实比熔融曲线据Yanetal(2018).Grt石榴石;Ol橄榄石;Opx斜方辉石;Cpx单斜辉石Fig.10DiagramofYbvsDy/YbforWeizhouIslandbasalts(afterAnetal,2017;Yanetal,2018)的晚新生代玄武岩均显示与OIB相似的地球化学特征(Yanetal,2014),但涠洲岛玄武岩具有明显窄的SrNdPb同位素变化范围(图7),综合说明了海南地幔柱的物质组成和温ActaPetrologicaSinica岩石学报2020,36(7)升过程中携带的再循环物质在地幔过渡带中的分层(Weietal,2012;Ballmeretal,2013;Xiaetal,2016;WeiandZhao,2020)。因此,基于海南地幔柱的双层上升模型(Weietal,2012;Xiaetal,2016;WeiandZhao,2020),我们认为起源于下地幔的海南地幔柱上升驱使再循环辉长岩洋壳上浮,当柱头达到地幔过渡带时几乎不受浮力影响而发生横向扩展形成一个化学成分和温度结构极其不均一的热化学堆积层(Ballmeretal,2013),随着源自柱尾持续的热和物质供给使得这个高度不均一的热化学堆积层重新上升(Xiaetal,2016),当进一步上升抵达地幔浅部时,由于过高热异常而使得再循环辉长岩洋壳优先达到其固相线温度经部分熔融交代周围的地幔橄榄岩生成石榴石辉石岩(贫硅辉石岩),混合的石榴石辉石岩和未反应的地幔橄榄岩再次发生部分熔融形成我们所观察到的涠洲岛碱性玄武岩。5结论(1)涠洲岛玄武岩显示与OIB相似的微量元素和同位素特征,在岩浆上升过程,几乎未受到地壳物质的混染,经历了橄榄石和单斜辉石的分离结晶作用。(2)SrNdPb同位素特征显示地幔端元组成呈DM与EM2二元混合趋势,其中DM和EM2组分均来源于海南地幔柱。(3)主、微量元素地球化学特征显示源区母岩由地幔橄榄岩和石榴石辉石岩(贫硅辉石岩)组成。其中,石榴石辉石岩与再循环辉长岩洋壳熔体交代地幔橄榄岩有关。(4)南海及周缘地区的晚新生代岩浆活动受控于海南地度结构具有不均一性的特点。海南地幔柱的不均一性可能幔柱。伴随着海南地幔柱的上升,再循环的辉长岩洋壳经部诱发了海南地区的火山活动,且这种不均一性源自地幔柱上分熔融交代地幔橄榄岩生成石榴石辉石岩,石榴石辉石岩和图11涠洲岛玄武岩MgOCaO(a,底图据HerzbergandAsimow,2008)和MgOFeOT/MnO(b,底图据Herzberg,2011)图解夏威夷玄武岩据Gaffneyetal(2005);FeOT=FeO+08998×Fe2O3Fig.11DiagramsofMgOvsCaO(a,afterHerzbergandAsimow,2008)andMgOvsFeOT/MnO(b,afterHerzberg,2008)forWeizhouIslandbasalts杨文健等:广西涠洲岛晚新生代玄武岩地幔源区及岩浆成因未反应的地幔橄榄岩混合部分熔融形成涠洲岛玄武岩。致谢感谢中国地质大学(北京)郝金华老师在电子探针分析工作中的热情支持;感谢中国地震局地质研究活动火山研究室魏费翔、韦伟和赵勇伟等老师在文章撰写和修改过程中的热情指导;感谢广西壮族自治区地震局和广西北海水文工程矿产地质勘查研究院提供钻孔岩芯样品;衷心感谢审稿人对本文提出的宝贵意见和建议。ReferencesAnAR,ChoiSH,YuYandLeeDC.2017.PetrogenesisofLateCenozoicbasalticrocksfromsouthernVietnam.Lithos,272-273:192-204BallmerMD,ItoG,WolfeCJandSolomonSC.2013.DoublelayeringofathermochemicalplumeintheuppermantlebeneathHawaii.EarthandPlanetaryScienceLetters,376:155-164BriaisA,PatriatPandTapponnierP.1993.UpdatedinterpretationofmagneticanomaliesandseafloorspreadingstagesintheSouthChinaSea:ImplicationsfortheTertiarytectonicsofSoutheastAsia.JournalofGeophysicalResearch:SolidEarth,98(B4):6299-6328ChenLH,ZengG,JiangSY,HofmannAW,XuXSandPanMB.2009.SourcesofAnfengshanbasalts:SubductedlowercrustintheSuluUHPbelt,China.EarthandPlanetaryScienceLetters,286(3-4):426-435DaiJG,WangCS,HébertR,LiYL,ZhongHT,GuillaumeR,BezardRandWeiYS.2011.LateDevonianOIBalkalinegabbrointheYarlungZangboSutureZone:RemnantsofthePaleoTethys?GondwanaResearch,19(1):232-243DasguptaR,HirschmannMMandSmithND.2007.Partialmeltingexperimentsofperidotite+CO2at3GPaandgenesisofalkalicoceanislandbasalts.JournalofPetrology,48(11):2093-2124DavisFA,HirschmannMMandHumayunM.2011.Thecompositionoftheincipientpartialmeltofgarnetperidotiteat3GPaandtheoriginofOIB.EarthandPlanetaryScienceLetters,308(3-4):380-390FanCY,XiaSH,ZhaoF,SunJL,CaoJH,XuHLandWanKY.2017.NewinsightsintothemagmatisminthenorthernmarginoftheSouthChinaSea:Spatialfeaturesandvolumeofintraplateseamounts.Geochemistry,Geophysics,Geosystems,18(6):2216-2239FanQC,SunQ,LongAM,YinKJ,SuiJL,LiNandWangTH.2006.GeologyanderuptionhistoryofvolcanoesinWeizhouIslandandXieyangIsland,NorthernBay.ActaPetrologicaSinica,22(6):1529-1537(inChinesewithEnglishabstract)FanQC,SuiJL,SunQ,LiN,ZhaoYWandDuXX.2008a.MeltinclusionsofperidotitesinWeizhouIsland,NorthernBay:Newevidenceofmantlemetasomatisminsubcontinentlithosphericmantle.ActaPetrologicaSinica,24(11):2495-2500(inChinesewithEnglishabstract)FanQC,SunQ,SuiJLandLiN.2008b.TraceelementandisotopicgeochemistryofvolcanicrocksanditstectonicimplicationsinWeizhouIslandandXieyangIsland,NorthernBay.ActaPetrologicaSinica,24(6):1323-1332(inChinesewithEnglishabstract)FreyFA,GreenDHandRoySD.1978.Integratedmodelsofbasaltpetrogenesis:AstudyofquartztholeiitestoolivinemelilititesfromsoutheasternAustraliautilizinggeochemicalandexperimentalpetrologicaldata.JournalofPetrology,19(3):463-513GaffneyAM,NelsonBKandBlichertToftJ.2005.MeltingintheHawaiianplumeat1~2MaasrecordedatMauiNui:Theroleofeclogite,peridotite,andsourcemixing.Geochemistry,Geophysics,Geosystems,6(10):Q10L11HanJW,XiongXLandZhuZY.2009.GeochemistryofLateCenozoic7012basaltsfromLeiqiongarea:TheoriginofEM2andthecontributionfromsubcontinentallithospheremantle.ActaPetrologicaSinica,25(12):3208-3220(inChinesewithEnglishabstract)HartSR.1984.AlargescaleisotopeanomalyintheSouthernHemispheremantle.Nature,309(5971):753-757HerzbergCandAsimowPD.2008.Petrologyofsomeoceanicislandbasalts:PRIMELT2.XLSsoftwareforprimarymagmacalculation.Geochemistry,Geophysics,Geosystems,9(9):Q09001HerzbergCandGazelE.2009.Petrologicalevidenceforsecularcoolinginmantleplumes.Nature,458(7238):619-622HerzbergC.2011.IdentificationofsourcelithologyintheHawaiianandCanaryIslands:Implicationsfororigins.JournalofPetrology,52(1):113-146HerzbergCandAsimowPD.2015.PRIMELT3MEGA.XLSMsoftwareforprimarymagmacalculation:PeridotiteprimarymagmaMgOcontentsfromtheliquidustothesolidus.Geochemistry,Geophysics,Geosystems,16(2):563-578HiroseKandKushiroI.1993.Partialmeltingofdryperidotitesathighpressures:Determinationofcompositionsofmeltssegregatedfromperidotiteusingaggregatesofdiamond.EarthandPlanetaryScienceLetters,114(4):477-489HiroseKandKawamotoT.1995.Hydrouspartialmeltingoflherzoliteat1GPa:TheeffectofH2Oonthegenesisofbasalticmagmas.EarthandPlanetaryScienceLetters,133(3-4):463-473HirschmannMM,KogisoT,BakerMBandStolperEM.2003.Alkalicmagmasgeneratedbypartialmeltingofgarnetpyroxenite.Geology,31(6):481-484HoangN,FlowerMFJandCarlsonRW.1996.Major,traceelement,andisotopiccompositionsofVietnamesebasalts:InteractionofhydrousEM1richasthenospherewiththinnedEurasianlithosphere.GeochimicaetCosmochimicaActa,60(22):4329-4351HoangTHA,ChoiSH,YuY,PhamTH,NguyenKHandRyuJS.2018.GeochemicalconstraintsonthespatialdistributionofrecycledoceaniccrustinthemantlesourceoflateCenozoicbasalts,Vietnam.Lithos,296-299:382-395HofmannAWandWhiteWM.1982.Mantleplumesfromancientoceaniccrust.EarthandPlanetaryScienceLetters,57(2):421-436HofmannAW,JochumKP,SeufertMandWhiteWM.1986.NbandPbinoceanicbasalts:Newconstraintsonmantleevolution.EarthandPlanetaryScienceLetters,79(1-2):33-45HuangLPandLiCN.2007.PeridotitexenolithfragmentsfromLateQuaternarypyroclasticflowinWeizhouIsland,Beihai,Guangxi.ActaPetrologicaetMineralogica,26(4):321-328(inChinesewithEnglishabstract)HuangXL,NiuYL,XuYG,MaJL,QiuHNandZhongJW.2013.GeochronologyandgeochemistryofCenozoicbasaltsfromeasternGuangdong,SEChina:ConstraintsonthelithosphereevolutionbeneaththenorthernmarginoftheSouthChinaSea.ContributionstoMineralogyandPetrology,165(3):437-455JiaDC,QiuXL,HuRZandLuY.2003.GeochemicalnatureofmantlereservoirsandtectonicsettingofbasaltsinBeibuGulfanditsadjacentregion.JournalofTropicalOceanography,22(2):30-39(inChinesewithEnglishabstract)KimJI,ChoiSH,KohGW,ParkJBandRyuJS.2019.PetrogenesisandmantlesourcecharacteristicsofvolcanicrocksonJejuIsland,SouthKorea.Lithos,326-327:476-490LeBasMJ,LeMaitreRWStreckeisenAandZanettinB.1986.Achemicalclassificationofvolcanicrocksbasedonthetotalalkalisilicadiagram.JournalofPetrology,27(3):745-750LebedevSandNoletG.2003.UppermantlebeneathSoutheastAsiafromSvelocitytomography.JournalofGeophysicalResearch,108(B1):2048LeiJS,ZhaoDP,SteinbergerB,WuB,ShenFLandLiZX.2009.NewseismicconstraintsontheuppermantlestructureoftheHainanplume.PhysicsoftheEarthandPlanetaryInteriors,173(1-2):33-50LiCNandWangFZ.2004.HolocenevolcaniceffusioninWeizhouIslandanditsgeologicalsignificance.JournalofMineralogyandPetrology,801224(4):28-34(inChinesewithEnglishabstract)LiCN,WangFZandZhongCS.2005.GeochemistryofQuaternarybasalticvolcanicrocksofWeizhouIslandinBeihaiCityofGuangxiandadiscussiononcharacteristicsoftheirsource.ActaPetrologicaetMineralogica,24(1):1-11(inChinesewithEnglishabstract)LiNS,YanQS,ChenZHandShiXF.2013.GeochemistryandpetrogenesisofQuaternaryvolcanismfromtheisletsintheeasternBeibuGulf:EvidenceforHainanplume.ActaOceanologicaSinica,32(12):40-49LiSG,YangW,KeS,MengXN,TianHC,XuLJ,HeYS,HuangJ,WangXC,XiaQK,SunWD,YangXY,RenZY,WeiHQ,LiuYS,MengFCandYanJ.2017.DeepcarboncyclesconstrainedbyalargescalemantleMgisotopeanomalyineasternChina.NationalScienceReview,4(1):111-120LiXMandLiuCZ.1991.PreliminarystudyonthepetrologicalcharacteristicsofalkalibasaltsinWeizhouIsland,Guangxi.GeologyofGuangxi,(1):43-51(inChinesewithEnglishabstract)LiYS.2012.QuantitativemodelingoftheEmeishanlargeigneousprovincemagmatismprocesses.Ph.D.Dissertation.Beijing:ChinaUniversityofGeosciences,1-163(inChinese)LinJ,LiuYS,YangYHandHuZC.2016.CalibrationandcorrectionofLAICPMSandLAMCICPMSanalysesforelementcontentsandisotopicratios.SolidEarthScience,1(1):5-27LiuCQ,MasudaAandXieGH.1994.MajorandtraceelementcompositionsofCenozoicbasaltsineasternChina:Petrogenesisandmantlesource.ChemicalGeology,114(1-2):19-42LiuCZ.1986.GeologicalandgeomorphologicalcharacteristicsofthecraterinsouthbayofWeizhouIsland.GeologyofGuangxi,(1):57-64(inChinesewithEnglishabstract)LiuJH,LiGZandNongHQ.1991.FeaturesofgeomorphyandQuaternarygeologyoftheWeizhouIsland.JournaloftheGuangxiAcademyofScience,7(1):27-36(inChinesewithEnglishabstract)LiuJQ,RenZY,NicholsARL,SongMS,QianSP,ZhangYandZhaoPP.2015.PetrogenesisofLateCenozoicbasaltsfromNorthHainanIsland:Constraintsfrommeltinclusionsandtheirhostolivines.GeochimicaetCosmochimicaActa,152:89-121LiuYS,HuZC,GaoS,GütherD,XuJ,GaoCGandChenHH.2008.InsituanalysisofmajorandtraceelementsofanhydrousmineralsbyLAICPMSwithoutapplyinganinternalstandard.ChemicalGeology,257(1-2):34-43LuJL.1993.ThecharacteristicsofvolcaniceruptioninWeizhouIsland,BeihaiCity,Guangxi.GuangxiGeology,6(3):53-58(inChinesewithEnglishabstract)McDonoughWFandSunSS.1995.ThecompositionoftheEarth.ChemicalGeology,120(3-4):223-253McKenzieDandBickleMJ.1988.Thevolumeandcompositionofmeltgeneratedbyextensionofthelithosphere.JournalofPetrology,29(3):625-679MckenzieDandONionsRK.1991.Partialmeltdistributionsfrominversionofrareearthelementconcentrations.JournalofPetrology,32(5):1021-1091MontelliR,NoletG,DahlenFAandMastersG.2006.Acatalogueofdeepmantleplumes:Newresultsfromfinitefrequencytomography.Geochemistry,Geophysics,Geosystems,7(11):Q11007PiletS,BakerMBandStolperEM.2008.Metasomatizedlithosphereandtheoriginofalkalinelavas.Science,320(5878):916-919PutirkaKD.2005.MantlepotentialtemperaturesatHawaii,Iceland,andthemidoceanridgesystem,asinferredfromolivinephenocrysts:Evidenceforthermallydrivenmantleplumes.Geochemistry,Geophysics,Geosystems,6(5):Q05L08PutirkaKD,PerfitM,RyersonFJandJacksonMG.2007.Ambientandexcessmantletemperatures,olivinethermometry,andactivevs.passiveupwelling.ChemicalGeology,241(3-4):177-206QiFQ,LiGZ,SunYF,LiangWandDuJ.2003.BasicgeomorphologicfeatureoftheWeizhouIslandoftheBeibuBay.AdvanceinMarineScience,21(1):41-50RenZY,IngleS,TakahashiE,HiraboNandHirataT.2005.TheActaPetrologicaSinica岩石学报2020,36(7)chemicalstructureoftheHawaiianmantleplume.Nature,436(7052):837-840RoederPLandEmslieRF.1970.Olivineliquidequilibrium.ContributionstoMineralogyandPetrology,29(4):275-289RudnickRLandGaoS.2003.Compositionofthecontinentalcrust.In:HollandHDandTurekianKK(eds.).TreatiseonGeochemistry.Amsterdam:Elsevier,1-64SatoH.1977.Nickelcontentofbasalticmagmas:Identificationofprimarymagmasandameasureofthedegreeofolivinefractionation.Lithos,10(2):113-120ShiXFandYanQS.2011.GeochemistryofCenozoicmagmatismintheSouthChinaSeaanditstectonicimplications.MarineGeology&QuaternaryGeology,31(2):59-72(inChinesewithEnglishabstract)SobolevAV,HofmannAW,SobolevSVandNikogosianIK.2005.AnolivinefreemantlesourceofHawaiianshieldbasalts.Nature,434(7033):590-597SobolevAV,HofmannAW,KuzminDV,YaxleyGM,ArndtNT,ChungSL,DanyushevskyLV,ElliottT,FreyFA,GarciaMO,GurenkoAA,KamenetskyVS,KerrAC,KrivolutskayaNA,MatvienkovVV,NikogosianLK,RochollA,SigurdssonLA,SushchevskayaNMandTeklayM.2007.Theamountofrecycledcrustinsourcesofmantlederivedmelts.Science,316(5823):412-417StaudigelH,ZindlerA,HartSR,LeslieT,ChenCYandClagueD.1984.Theisotopesystematicsofajuvenileintraplatevolcano:Pb,Nd,andSrisotoperatiosofbasaltsfromLoihiSeamount,Hawaii.EarthandPlanetaryScienceLetters,69(1):13-29StrackeA,BizimisMandSaltersVJM.2003.Recyclingoceaniccrust:Quantitativeconstraints.Geochemistry,Geophysics,Geosystems,4(3):8003SunSSandMcDonoughWF.1989.Chemicalandisotopicsystematicsofoceanicbasalts:Implicationsformantlecompositionandprocesses.In:SaundersADandNorryMJ(eds.).MagmatismintheOceanBasins.GeologicalSociety,London,SpecialPublications,42(1):313-345TaylorBandHayesDE.1983.OriginandhistoryoftheSouthChinaSeaBasin.In:HayesDE(ed.).TheTectonicandGeologicEvolutionofSoutheastAsianSeasandIslands:Part2.WashingtonDC:AmericanGeophysicalUnion,27:23-56ThompsonRNandGibsonSA.2000.TransienthightemperaturesinmantleplumeheadsinferredfrommagnesianolivinesinPhanerozoicpicrites.Nature,407(6803):502-506TuK,FlowerMFJ,CarlsonRW,ZhangMandXieGH.1991.Sr,Nd,andPbisotopiccompositionsofHainanbasalts(SouthChina):ImplicationsforasubcontinentallithosphereDupalsource.Geology,19(6):567-569TuK,FlowerMFJ,CarlsonRW,XieGH,ChenCYandZhangM.1992.MagmatismintheSouthChinaBasin:1.IsotopicandtraceelementevidenceforanendogenousDupalmantlecomponent.ChemicalGeology,97(1-2):47-63WangK,PlankT,WalkerJDandSmithEI.2002.AmantlemeltingprofileacrosstheBasinandRange,SWUSA.JournalofGeophysicalResearch,107(B1):ECV51ECV521WangXC