矿物岩石地球化学通报  2017, Vol. 36 Issue (6): 905-911   PDF    
全球大陆板内玄武岩数据挖掘
杜雪亮 , 李玉琼 , 杜君 , 马蓁 , 王金荣     
兰州大学 地质科学与矿产资源学院 甘肃省西部矿产资源重点实验室, 兰州 730000
摘要: 大陆板内玄武岩(ICB)是一个广义的术语,包括大陆溢流玄武岩(CFB)、大陆裂谷玄武岩(CRB)以及所有出露在大陆板内构造环境的玄武岩。本文研究的大陆板内玄武岩不包括大陆溢流玄武岩和大陆裂谷玄武岩,主要是一些规模比较小的玄武岩,是狭义的ICB。本文详细介绍了研究的方法、数据清洗过程以及数据处理过程。在上述研究的基础上探讨了新生代ICB的时空分布,着重对不同时代的ICB进行了对比,发现中新生代、古生代、元古代的ICB的地球化学特征总体上一致,其中,古生代和元古代之间的相似性更加明显,而中新生代则显示V、Sc、Yb亏损的特征,具体如何解释还需要今后进一步的研究。本文的研究表明,全球不同时代ICB的地球化学特征大体相似,地质理论中的"将今论古"原则适用于ICB,只是在应用时尽量回避Al2O3、Zn、V、Sc、Yb等元素就更加万无一失了。
关键词: 大陆板内玄武岩      大数据      置信椭圆      时代      对比     
A Mining of Data of the Global Intra-continental Basalt
DU Xue-liang, LI Yu-qiong, DU Jun, MA Zhen, WANG Jin-rong     
Key Laboratory of Mineral Resources in Western China(Gansu Province), School of Earth Sciences, Lanzhou University, Lanzhou 730000, China
Abstract: Intra-continental basalt (ICB) is a broad term which is used to describe various basalts, including continental flood basalt (CFB) and continental rift basalt (CRB), and all other basalts exposed in the intra-continental tectonic environment. The intra-continental basalt studied in this paper is only the basalt limitedly distributed within the continent. It does not include continental flood basalt and continental rift basalt. In this paper, the research methods, processes of data cleaning and data processing have been described in details. Based on the above work, we have discussed the temporal and spatial distribution of the Cenozoic ICB. Especially, we have compared the ICB in various ages and found that the Mesozoic-Cenozoic, Paleozoic and Proterozoic ICB have generally similar geochemical characteristics, with more obvious similarities between the Paleozoic and the Proterozoic ICB. The Mesozoic-Cenozoic ICB has characteristics of relative depletion of V, Sc, Yb. How to explain this phenomenon needs further research in the future. The results of this paper show that the global ICB in various ages have generally similar geochemical characteristics. The principle of "The present is the key to the past" could be applied perfectly to comparing geochemical characteristics of the ICB in various ages if some elements, such as Al, Zn, V, Sc, and Yb, can be not considered.
Key words: intra-continental basalt     big data     confidence ellipsoid     time     comparison    

大陆板内玄武岩是一个含糊的说法,有广义和狭义之分。大陆溢流玄武岩和大陆裂谷玄武岩也是大陆板内玄武岩。学术界通常将上述二者分别予以研究和讨论。实际上大陆不止上述两类玄武岩,除了上述两类而外,还有许许多多规模比较小的玄武岩(以及基性侵入岩),它们既不是大陆溢流玄武岩(因为规模不大),也不是大陆裂谷玄武岩(因为没有证据是裂谷),因此,就笼统地归为大陆板内玄武岩了。其实,学术界所谓的大陆板内玄武岩(ICB)是包括大陆溢流玄武岩(CFB)和大陆裂谷玄武岩(CRB)的。而数据库所收集的大体是除了大陆溢流玄武岩和大陆裂谷玄武岩之外的所有的大陆板内玄武岩,因此,本文研究的是狭义的大陆板内玄武岩。

狭义的大陆板内玄武岩是最常见的、规模较小的玄武岩。大陆溢流玄武岩非常受学术界关注,中国充其量只有峨眉山玄武岩和塔里木大火成岩省。大陆裂谷玄武岩非常有名,如东非裂谷玄武岩、美国西部盆岭省玄武岩等,而中国并没有这个类型的玄武岩出露。早先常说的攀西裂谷至今也未见有裂谷玄武岩的报道(沈发奎,1989何斌等,2003钟宏等,2009),与其有关的恰恰是峨眉山玄武岩,可见攀西裂谷并非裂谷。目前所出现的汉诺坝玄武岩、长白山玄武岩、五大连池玄武岩、海南玄武岩、女山玄武岩、方山玄武岩、河北义县组玄武岩、山西繁峙玄武岩、碧口玄武岩、塔里木二叠纪玄武岩等等,都可以将其归于狭义的大陆板内玄武岩。

按照早先的认识,狭义的大陆板内玄武岩的地球化学特征通常与OIB类似,由于ICB的规模小,岩浆上升穿过陆壳时不可避免地会与陆壳发生物质交换,因此,是否有混染作用的痕迹是判断是否ICB的重要指标(Arndt and Christensen, 1992; Gibson et al., 2006; Heinonen et al., 2010)。至于ICB的成因,由于CFB通常与地幔柱有关(Wilson,1997),有人认为小规模的ICB也可能是这种成因,而多数学者并不认同这种见解,因为,地球内部不可能有这么多的地幔柱或地幔羽(plume)。因此,ICB很可能是软流圈不均一性引起的,是软流圈局部热散失的一种方式(Courtillot et al., 2003; Conrad et al., 2010Conrad et al., 2011)。至于为什么会出现这种情况,仍然是需要认真研究的。很少有人从全球角度对狭义的ICB进行研究,王金荣等(2016)对此从判别标志出发进行了初步的探讨,得出了很有意义的成果。本文重新收集了数据库的全球ICB数据,对新生代ICB的时空分布以及不同时代的ICB进行了对比研究,取得的认识汇报如下。

1 研究方法

本文主要利用GEOROC数据库资料,对全球大陆板内玄武岩的构造环境、经纬度及地球化学元素信息进行数据挖掘,旨在从数据分析中探究大陆板内玄武岩的源区性质以及不同地区、不同时代之间的关系,试图获得一些有意义的信息和相应的结论。研究总共分为4个步骤进行:数据清洗→提取目标数据→数据处理及数据成图→结果分析。

数据清洗主要参考王金荣等(2016)的方法。下载关于大陆板内玄武岩的全部数据,共计86 271个;清洗后最终获得全球大陆板内玄武岩数据23 622个。

为了寻求地球化学元素与构造环境之间的关系,本次计算依据元素包含的样品数量大于总样品数量10%的原则,选定了31个元素,其中包括10个主量元素、21个微量元素。对样品地球化学数据的处理过程主要包括数据二次清洗、比值计算、置信椭圆方程以及交叠率计算3个部分。置信椭圆方程以及交叠率计算方法和注意事项刘欣雨等(2017)已作了详细的介绍,本文不再赘述。

2 新生代大陆板内玄武岩时空分布

时间上,全球新生代大陆板内玄武岩主要集中分布于中新世,其次分布于更新世、上新世和始新世(图 1);而在地域空间上(图 2),整个新生代的大陆板内玄武岩主要集中分布于亚洲板块。

数据来自PetDB和GEOROC数据库 图 1 新生代ICB时代分布直方图 Figure 1 The histogram of ages of the Cenozoic ICB

数据来自PetDB和GEOROC数据库 图 2 新生代不同时代ICB全球分布图 Figure 2 Map of global distribution of the Cenozoic ICB

在GOEROC数据库中,新生代大陆板内玄武岩分布于不同的地区。本文利用样品所分布的区域,对数据进行了初步的分类,并统计了不同地区大陆板内玄武岩的频率分布情况。其中(图 23),全球新生代ICB主要分布于中亚造山带(Central Asian Orogenic Belt,CAB)、安纳托利亚-朗造山带(Anatolia-Iran belt,AIB)、欧洲造山带(European orogenic belt,EOB),次要分布于阿拉伯努比亚地盾(Arabian-Nubian Shield)、喀尔巴阡山和潘诺尼亚盆地(Carpathian belt and Pannonian Basin)、澳大利亚(Australia)、新西兰岛(New Zealand)等地。

数据来自PetDB和GEOROC数据库 图 3 不同时代ICB全球分布图 Figure 3 Map of global distribution for the ICB in various epochs
3 大陆板内玄武岩年代对比

在全球ICB时代分布直方图(图 4)中,全球大陆板内玄武岩主要集中分布于中新生代,元古宙次之,古生代最末。利用时代的分类(图 3),对中新生代、古生代、元古宙的相关数据进行研究,进而从图中得到一些规律和发现。

数据来自PetDB和GEOROC数据库 图 4 全球不同时代ICB时代分布直方图 Figure 4 The histogram of ages of global ICB in various epochs

不同时代大陆板内玄武岩数据成图共计533640个。研究表明,中新生代、古生代、元古宙的大陆板内玄武岩样品数据对应的置信椭圆区域大体是重叠的,意味着不同时代的大陆板内玄武岩在地球化学特征上保持一致性。但是,有一些元素却呈现出一定的规律性差异。在图 5中,古生代和元古宙的大陆板内玄武岩样品的置信椭圆分布区域相似;而中新生代的大陆板内玄武岩样品相对与前两者而言,明显地呈现出贫V、Sc、Yb的特征。

图 5 不同时代ICB置信椭圆分布图(图a: Sc; 图b: Yb; 图c:V) Figure 5 The Confidence profiles of the ICB in various epochs(Fig.a: Sc; Fig.b: Yb; Fig.c: V)
4 讨论

(1) 人们通常认为,产于板内构造环境的板内玄武岩(ICB)与洋岛玄武岩(OIB)具有类似的地球化学特征,与地幔柱的活动有关,是来自下地幔的岩浆活动(Morgan and Morgan, 1999; 王金荣等,2016)。虽然这种说法也存在争议(Anderson,1994; Meibom and Ander-son, 2003)。

认同地幔柱假说的学者在讨论地幔柱假说时,总会要讨论为OIB在地球化学上高度富集的特点。然而其物源至今并不十分清楚,且学术界争议很大。Hofmann和White等(Hofmann and White, 1892; Sobolev et al., 2000; 2005; 2007)认为OIB源区物质是循环的古洋壳(ROC),而许多学者(Woodhead and McCulloch, 1989; Weaver,1991; Chauvel et al., 1992; Woodhead and Devey, 1993; White and Duncan, 1996; Hofmann,1997; Eisele et al., 2002; Murphy et al., 2002)认为OIB中富集的组分来自循环陆壳(RCC,陆源沉积物)。牛耀龄(2010)认为ROC模型的诸多问题之一在于,OIB高度富集(如PM$ \gg $1),而洋壳则相当亏损(PM < 1),故循环洋壳不可能是OIB的源区物质。而RCC模型中,尽管陆壳(CC)和OIB富集很多不相容元素(如PM$ \gg $1),但是CC极度富Pb而贫Nb、Ta、P和Ti。CC的这一独特的地球化学性质强度很大,若的确是OIB源区重要组成的话,那它一定会在OIB有所显示,然而OIB并不具有这种CC特征(图 7),因此CC不会是OIB的富集组分。

CFB数据王金荣等(2016);标准化值据OIB数据安屹等(2017) 图 7 全球CFB、ICB与OIB均值微量元素蛛网图(a)和稀土元素配分图(b) Figure 7 Primitive mantle-normalized trace element spidergrams(a) and chondrite-normalized REE patterns (b)of CFB, ICB and OIB

ICB同OIB一样,源自于下地幔岩浆活动。从143 Nd/144 Nd -87 Sr/86 Sr图(图 6)看,全球新生代的大陆板内玄武岩样品集中分布于原始地幔(PM)和OIB的近乎全部范围;在207 Pb/204 Pb -206 Pb/204 Pb图(图 6)中,中亚造山带(CAB)的新生代大陆板内玄武岩样品主要分布在富集端元EMI和EMII、原始地幔(PM)以及MORB范围内,其中安纳托利亚伊朗造山带(AIB)的样品集中分布于富集端元EMII,欧洲造山带(EOB)主要分布于MORB范围。中亚造山带(CAB)覆盖范围广,故其样品分布范围具有多样性;而欧洲造山带(EOB)分布于MORB范围,可能是由于大西洋洋脊活动充当了其源区物质的缘故。相对于207 Pb/204 Pb来说,206 Pb/204 Pb的比值略高,也指示了源区地幔富集206 Pb/204 Pb的特征。在图143 Nd/144 Nd-207 Pb/204 Pb和图87 Sr/86 Sr-207 Pb/204 Pb图(图 6)中,富集的程度更加明显。大部分样品投到了原始地幔(PM)和全硅酸盐地球(BSE)的范围内,部分投入EMII的范围内。而在图87 Sr/86 Sr -207 Pb/204 Pb图(图 6)中,欧洲造山带(EOB)依旧主要分布于MORB范围。总体上,新生代的大陆板内玄武岩的同位素地球化学特征总体表现为富集,且集中分布范围具有地域差异性。

Rollinson(2014); Zindler和Hart(1986) 图 6 全球新生代不同地区大陆板内玄武岩Pb-Sr-Nd相关图 Figure 6 Nd-Sr-Pb isotopic diagrams of the global Cenozoic ICB in different regions

虽然ICB与OIB在物源上具有相似性,但是在岩浆的上升过程中,OIB处于洋壳内部,不受陆壳混染的影响。而形成ICB的来自下地幔的岩浆经上、下地幔过渡带、上地幔以及陆壳,将会有一些物质的加入。因此与地幔柱有关的大陆板内玄武岩同位素组成不仅具有似原始地幔的特征,同时表现出岩浆上升过程中的“多源”同位素特征(朱炳泉,2007)。这也可能是上述置信椭圆分布图中,不管按照地区划分和按照年代划分,一些元素出现一定程度的富集或者贫乏的原因。具体如何解释还需要今后进一步的研究。元素Pb就是一个很好的例子,ICB与洋岛玄武岩OIB比较,明显富集Pb,这是由于OIB不受陆壳的影响,因而Pb含量较低;而大ICB很难避免陆壳混染的影响,造成其最明显的特征之一就是富集Pb(图 7)。

表 1 全部的CFB、ICB、OIB样品的主元素、微量元素、稀土元素和同位素含量表 Table 1 A table of major, trace, REE elements contents and isotope composition for the CFB, ICB and OIB

(2) 将今论古法是由赫顿(Hutton,1788)提出的一个原则。其经典的概念是“今天是过去的钥匙”。本文研究表明,在不同时代大陆板内玄武岩的置信椭圆分布图中,中新生代、古生代、元古代的置信椭圆区域大体是重叠的,虽然在V、Sc、Yb元素上,中新生代相对于古生代和元古代呈现亏损的特征,但并不影响整体的分布,这也就意味着不同时代的大陆板内玄武岩在地球化学特征上保持一致性。因此,利用地质理论中的“将今论古”原则来研究大陆板内玄武岩是完全可行,只是在应用时尽量回避V、Sc、Yb等元素就更加万无一失了。

5 结论

(1) 大陆板内玄武岩同洋岛玄武岩一样,源自于下地幔岩浆活动。与其不同的是,形成大陆板内玄武岩的来自下地幔的岩浆经上、下地幔过渡带、上地幔以及陆壳,将会有一些物质的加入。因此与地幔柱有关的大陆板内玄武岩同位素组成不仅具有似原始地幔的特征,同时表现出岩浆上升过程中的“多源”同位素特征。

(2) 对比不同时代的大陆板内玄武岩,发现中新生代、古生代、元古代的ICB的地球化学特征总体上一致。其中,古生代和元古代之间的相似性更加明显,而中新生代则呈现出亏损V、Sc、Yb的特征。暗示地质理论中的“将今论古”原则适用于ICB,只是在应用时尽量回避Al2O3、Zn、V、Sc、Yb等元素就更加万无一失了。

参考文献
[] Anderson D L. 1994. The sublithospheric mantle as the source of continental flood basalts; the case against the continental lithosphere and plume head reservoirs. Earth and Planetary Science Letters, 123(1-3): 269–280. DOI:10.1016/0012-821X(94)90273-9
[] Arndt N, Christensen U. 1992. The role of lithospheric mantle in continental flood volcanism; thermal and geochemical constraints. Journal of Geophysical Research, 97(B7): 10967–10981. DOI:10.1029/92JB00564
[] Chauvel C, Hofmann A W, Vidal P. 1992. HIMU-EM:The French Polynesian connection. Earth and Planetary Science Letters, 110(1-4): 99–119. DOI:10.1016/0012-821X(92)90042-T
[] Conrad C P, Bianco T A, Smith E I, Wessel P. 2011. Patterns of intraplate volcanism controlled by asthenospheric shear. Nature Geoscience, 4(5): 317–321. DOI:10.1038/ngeo1111
[] Conrad C P, Wu B J, Smith E I, Bianco T A, Tibbetts A. 2010. Shear-driven upwelling induced by lateral viscosity variations and asthenospheric shear:A mechanism for intraplate volcanism. Physics of the Earth and Planetary Interiors, 178(3-4): 162–175. DOI:10.1016/j.pepi.2009.10.001
[] Courtillot V, Davaille A, Besse J, Stock J. 2003. Three distinct types of hotspots in the Earth's mantle. Earth and Planetary Science Letters, 205(3): 295–308.
[] Eisele J, Sharma M, Galer S J, Blichert-Toft J, Devey C W, Hofmann A. 2002. The role of sediment recycling in EM-1 inferred from Os, Pb, Hf, Nd, Sr isotope and trace element systematics of the Pitcairn hotspot. Earth and Planetary Science Letters, 196(3-4): 197–212. DOI:10.1016/S0012-821X(01)00601-X
[] Gibson S A, Thompson R N, Day J A. 2006. Timescales and mechanisms of plume-lithosphere interactions:40Ar/39Ar geochronology and geochemistry of alkaline igneous rocks from the Paraná-Etendeka large igneous province. Earth and Planetary Science Letters, 251(1-2): 1–17. DOI:10.1016/j.epsl.2006.08.004
[] Heinonen J S, Carlson R W, Luttinen A V. 2010. Isotopic(Sr, Nd, Pb and Os)composition of highly magnesian dykes of Vestfjella, western Dronning Maud Land, Antarctica:A key to the origins of the Jurassic Karoo large igneous province. Chemical Geology, 277(3-4): 227–244. DOI:10.1016/j.chemgeo.2010.08.004
[] Hofmann A W. 1997. Mantle geochemistry:The message from oceanic volcanism. Nature, 385(6613): 219–229. DOI:10.1038/385219a0
[] Hofmann A W, White W M. 1982. Mantle plumes from ancient oceanic crust. Earth and Planetary Science Letters, 57(2): 421–436. DOI:10.1016/0012-821X(82)90161-3
[] Hutton J. 1788. Theory of the Earth; or an Investigation of the Laws observable in the Composition, Dissolution, and Restoration of Land upon the Globe. Earth and Environmental Science Transactions of The Royal Society of Edinburgh, 1(2): 209–304. DOI:10.1017/S0080456800029227
[] Meibom A, Anderson D L. 2004. The statistical upper mantle assemblage. Earth and Planetary Science Letters, 217(1-2): 123–139. DOI:10.1016/S0012-821X(03)00573-9
[] Morgan J P, Morgan W J. 1999. Two-stage melting and the geochemical evolution of the mantle:A recipe for mantle plum-pudding. Earth and Planetary Science Letters, 170(3): 215–239. DOI:10.1016/S0012-821X(99)00114-4
[] Murphy D T, Collerson K D, Kamber B S. 2002. Lamproites from Gaussberg, Antarctica:Possible transition zone melts of Archaean subducted sediments. Journal of Petrology, 43(6): 981–1001. DOI:10.1093/petrology/43.6.981
[] Sobolev A V, Hofmann A W, Kuzmin D V, Yaxley G M, Arndt N T, Chung S L, Danyushevsky L V, Elliott T, Frey F A, Garcia M O, Gurenko A A, Kamenetsky V S, Kerr A C, Krivolutskaya N A, Matvienkov V V, Nikogosian I K, Rocholl A, Sigurdsson I A, Sushchevskaya N M, Teklay M. 2007. The amount of recycled crust in sources of mantle-derived melts. Science, 316(5823): 412–417. DOI:10.1126/science.1138113
[] Sobolev A V, Hofmann A W, Nikogosian I K. 2000. Recycled organic crust observed in 'ghost plagioclase' within the source of Mauna Loa lavas. Nature, 404(6781): 986–990. DOI:10.1038/35010098
[] Sobolev A V, Hofmann A W, Sobolev S V, Nikogosian I K. 2005. An olivine-free mantle source of Hawaiian shield basalts. Nature, 434(7033): 590–597. DOI:10.1038/nature03411
[] Weaver B L. 1991. The origin of ocean island basalt end-member compositions:Trace element and isotopic constraints. Earth and Planetary Science Letters, 104(2-4): 381–397. DOI:10.1016/0012-821X(91)90217-6
[] White W M, Duncan R A. 1996. Geochemistry and geochronology of the Society Islands:New evidence for deep mantle recycling. In:BASU A, HART S, Eds. Earth Processes:Reading the Isotopic Code. Washington, DC:American Geophysical Union: 183–206.
[] Wilson M. 1997. Thermal evolution of the Central Atlantic passive margins:Continental break-up above a Mesozoic super-plume. Journal of Geological Society, 154: 491–495. DOI:10.1144/gsjgs.154.3.0491
[] Woodhead J D, Devey C W. 1993. Geochemistry of the Pitcairn seamounts, Ⅰ:Source character and temporal trends. Earth and Planetary Science Letters, 116(1-4): 81–99. DOI:10.1016/0012-821X(93)90046-C
[] Woodhead J D, Mcculloch M T. 1989. Ancient seafloor signals in Pitcairn Island lavas and evidence for large amplitude, small length-scale mantle heterogeneities. Earth and Planetary Science Letters, 94(3-4): 257–273. DOI:10.1016/0012-821X(89)90145-3
[] Zindler A, Hart S R. 1986. Chemical geodynamics. Annual Review of Earth and Planetary Sciences, 14: 493–571. DOI:10.1146/annurev.ea.14.050186.002425
[] 安屹, 杨婧, 陈万峰, 王金荣, 张旗, 潘振杰, 焦守涛. 2017. N-MORB、E-MORB和OIB的区别及其可能的原因:大数据的启示. 地质科学, 52(3): 727–742. DOI:10.12017/dzkx.2017.046
[] 何斌, 徐义刚, 肖龙, 王雅玫. 2003. 攀西裂谷存在吗. 地质论评, 49(6): 572–582.
[] 刘欣雨, 张旗, 张成立. 2017. 全球新生代安山岩构造环境有关问题探讨. 地质科学, 52(3): 649–667. DOI:10.12017/dzkx.2017.042
[] 牛耀龄. 2010. 板内洋岛玄武岩(OIB)成因的一些基本概念和存在的问题. 科学通报, 55(2): 103–114.
[] 沈发奎. 1989. 攀西裂谷构造-岩浆活动及其成因探讨. 岩石矿物学杂志, 8(4): 311–320.
[] 王金荣, 潘振杰, 张旗, 陈万峰, 杨婧, 焦守涛, 王淑华. 2016. 大陆板内玄武岩数据挖掘:成分多样性及在判别图中的表现. 岩石学报, 32(7): 1919–1933.
[] 钟宏, 徐桂文, 朱维光, 胡瑞忠, 何德锋. 2009. 峨眉山大火成岩省太和花岗岩的成因及构造意义. 矿物岩石地球化学通报, 28(2): 99–110.
[] 朱炳泉. 2007. 全球幔源岩Pb-Sr-Nd同位素体系. 地学前緣, 14(2): 24–36.