闽西南基性岩脉中捕获锆石SIMS U-Pb年龄及Hf、O同位素特征
张文慧, 王翠芝, 李晓敏, 刘文元
福州大学紫金矿业学院,福建 福州 350108

作者简介:张文慧(1970-),女,山东日照人,讲师,主要从事岩石地球化学方面的研究.E-mail:zhangw5@163.com

摘要

测定了闽西南地区5件基性岩脉的全岩主微量元素,锆石U-Pb年龄和Hf,O同位素。5件岩脉主要为辉绿岩,其全岩SiO2含量为45%~53%,稀土元素显示轻稀土富集的右倾配分模式。基性岩脉中大部分锆石具明显振荡环带和扇状环带,为典型岩浆结晶锆石特征。锆石U-Pb年龄(96~2 400 Ma)分布较为分散。除少量锆石年龄(96~142 Ma)可能指示岩脉形成年龄外,其余均为捕获锆石。捕获锆石年龄主要分布在4个范围:早元古代(2 467~1 796 Ma),中晚元古代—震旦纪(1 343~647 Ma),志留纪—晚三叠世(427~225 Ma),晚侏罗世(159~140 Ma)。Hf-O同位素显示早元古代锆石来源于接近球粒陨石均一储库的地幔。中晚元古代以后年龄的锆石其Hf-O同位素均具有亏损地幔岩浆与地壳组分混合的特征;志留纪—晚侏罗世锆石主要来源于S型壳源花岗岩的重熔;早白垩世晚期的锆石εHf( t) 值与δ18O值清晰地显示出亏损地幔与地壳岩浆混合的趋势。闽西南基性岩脉中锆石的二阶段Hf亏损地幔模式年龄(TDM2)峰值主要分布在1.6~1.9 Ga,说明早元古代晚期幔源岩浆作用形成的基性岩地壳可能是形成后期花岗岩的主要源区。

关键词: SIMS锆石U-Pb测年; Hf-O同位素; 基性岩脉; 闽西南
中图分类号:P588.12+4 文献标志码:A 文章编号:1001-8166(2016)03-0320-15
Zircon SIMS U-Pb Age, Hf and O Isotopes of Mafic Dikes, Southwest Fujian Province
Zhang Wenhui, Wang Cuizhi, Li Xiaomin, Liu Wenyuan
College of Zijin Mining, Fuzhou University, Fuzhou 350108,China

First author:Zhang Wenhui (1970-), female, Rizhao City, Shandong Province, Lecturer. Research area include petro geochemistry.E-mail:zhangw5@163.com

Abstract

The study in this paper determined whole rock major and trace elements, zircon U-Pb age and Hf, O isotopes of 5 mafic dikes in the southwestern Fujian province. The 5 dikes are mainly diabase and the whole rock SiO2 content are between 45%~53%. Most zircons of the mafic dikes display obvious oscillatory zoning and fan-shaped zoning, and have the typical magmatic zircon crystallization characteristics. Zircon U-Pb age is dispersed with 96~2 400 Ma range. In addition to the minimum age (96~142 Ma) which might be the age of the formation of dikes, the remaining are captured zircon. The captured zircon age was mainly distributed in 4 groups: Early Proterozoic (2 467~1 796 Ma); Middle and late Proterozoic (1 343~647 Ma); Silurian to late Triassic Epoch (427~225 Ma); and Late Jurassic (159~140 Ma). Hf-O isotope shows that the early Proterozoic zircon was derived from the mantle of the homogeneous chondrite reservoir, and the others show magmatic mixing characteristics between depleted mantle and crust. Zircon’s εHf( t) and δ18O of the early Late Cretaceous clearly show the mixing trend of depleted mantle and crustal magma. The peak of zircon Hf two-stage depleted mantle model age TDM2 was mainly distributed in the 1.6~1.9 Ga. The Early Proterozoic mafic crust might be the main source for latter granite.

Keyword: SIMS zircon U-Pb age; Hf-O isotopes; Mafic dikes; Southwestern Fujian.

闽西南地区位于南岭成矿带东段与东南沿海岩浆岩带的衔接处, 是中生代太平洋板块与欧亚板块相互作用的重要地区。区内中生代岩浆活动强烈, 但以花岗岩类为主。花岗质岩石的地球化学特征通常受到源岩性质、岩浆的形成和演化等多种因素的影响, 没有明确的构造环境指示意义, 而幔源基性岩浆的地球化学组成则能够较好地反映地幔组成、热结构及其构造背景[1]。基性岩脉是起源于地幔的基性岩浆充填于先存破裂或区域性节理形成的, 是深部热动力作用在地表的重要表现形式[2]。闽西南地区基性岩脉分布较少, 且基性岩脉在成岩后普遍受到了后期蚀变作用的影响。对基性岩传统的全岩主微量及Sr, Nd同位素的研究, 得到的是现今岩浆成分的一个平均值。如果一个岩石由多种组份构成, 则全岩的同位素分析结果是不能对岩浆多源区性质及混合等特征进行分辨的, 而且岩浆演化及后期蚀变作用过程会改变其原始岩浆的地球化学组成。近年来研究显示, 通过识别不同成因锆石的同位素特征可以认识它的演化[3]。近年来迅速发展起来的锆石原位Hf-O同位素分析可以避开岩浆后期演化及蚀变的影响, 能够深入探寻岩浆源区特征及深部演化过程。基于此, 我们采集了闽西南地区几个代表性区域的基性岩脉, 对其中所含的锆石进行了详细的 U-Pb 定年和Hf-O 同位素测试。结果显示, 几件基性岩脉样品中均含有数量较多且不同时期的捕获锆石, 为我们研究闽西南地区的岩浆期次及岩浆源区演化特征提供了线索。同时, 也显示出该地区基性岩脉成因的复杂性, 岩脉全岩的地球化学性质可能受到岩浆多源区混合特征的影响, 已经不能准确反映岩脉形成时的源区特征和演化过程。

1 地质背景

福建省岩浆岩广泛发育。自元古代开始 (可能有太古代), 地壳演化各主要阶段均有岩浆活动。岩石类型齐全, 以酸性和中酸性岩类占绝对优势。闽西南坳陷带位于福建西南部(图1), 北以宁化— 南平一线与闽西北隆起带相接, 东以政和— 大埔断裂带与闽东火山断陷带相邻, 往西南延入广东省。闽西南坳陷带是叠加在加里东褶皱基础上发展起来的, 主要由晚古生代— 三叠纪地层组成, 褶皱、断裂均很发育, 并有华力西— 印支期花岗岩侵入。基底出露震旦系和下古生界浅变质岩, 属复理石沉积。

闽西南地区是华南金、银、铜、铁、钼金属成矿带的重要组成部分, 区内矿产资源丰富。著名的紫金山铜金矿、马坑铁矿等分布其中, 前者是中国东部环太平洋斑岩铜矿带的4个主要斑岩铜(金)矿之一[5]。这些矿床的成矿物质来源与幔源物质存在着一定联系。尽管这一地区中生代花岗岩广泛分布, 但仅通过花岗岩的野外地质、岩相学、元素和Sr-Nd 同位素地球化学证据, 在许多情况下还难以对幔源岩浆是否参与花岗岩的形成提供明确的约束[3], 幔源岩浆的源区性质以及对成矿的影响还存在许多未知的问题。本文采集了紫金山铜金矿附近悦洋 (YY-1)、四坊 (SF-5) 2件样品, 马坑铁矿 (MK-1) 以及位于两矿区之间(SH-1和SH-2)的2件样品, 共计5件基性岩脉作为研究对象。这些基性脉体主要为辉绿岩, 产状较陡, 侵入到晚侏罗世或早白垩世花岗岩中(图1)。脉体未发生变质或变形, 保持了形成时的应力状态。5件基性岩脉样品详细的岩性描述见表1, 部分样品的镜下照片如图2所示。

图1 闽西南地区地质简图及基性岩脉采样位置[4]Fig.1 Geological sketch map of southwestern Fujian and mafic dikes’ sampling positions[4]

图2 闽西基性岩脉岩相学特征
Act.阳起石; Hy.紫苏辉石; Pl.斜长石
Fig.2 The petrographic characteristics of mafic dikes
Act.Actinolite; Hy.Hypersthene; Pl. Plagioclase

2 分析方法

主微量元素测试在河北省区域地质矿产调查研究所实验室完成。主量元素测试采用XRF荧光测试法, 测试误差小于2%; 稀土和微量元素测试采用等离子体质谱法(ICP-MS), 测试精度优于5%。

锆石从约20 kg的样品中选出, 在无污染条件下经过传统的重液和磁法分选, 并在双目镜下手工挑选出来。本次样品挑选出的锆石数量见表1。将锆石颗粒和锆石标样(Qinghu、Plé sovice和91500)用环氧树脂粘在靶上, 抛光至锆石的内部暴露在外面以后, 在表面喷金。对锆石颗粒拍摄反射光、透射光和阴极发光(CL)图像以选择合适点位。锆石CL图像在中国地质科学院大陆构造与动力学国家重点实验室扫描电子显微镜上拍摄。锆石颗粒先进行原位O同位素分析, 然后在相同点位进行U-Pb年龄和Hf同位素测定。

表1 闽西南基性岩脉岩石学特征 Table 1 Petrological characteristics of mafic dikes
2.1 SIMS锆石O同位素分析

锆石微区原位O同位素分析, 在中国科学院地质与地球物理研究所离子探针实验室, Cameca IMS-1280型双离子源多接收器二次离子质谱仪(SIMS)上进行, 分析方法详见参考文献[6]。强度约为2 nA 的一次133Cs+离子束通过10 kV加速电压轰击样品表面。样品表面信号采集大小约为20 μ m。垂直入射的电子枪均匀覆盖于100 μ m范围内, 以中和样品的表面电荷效应。以2个法拉第杯同时接收16O 和18O。单点测量时间约为5 min, 单组18O/16O数据内精度一般优于0.2‰ ~0.3‰ (1σ ), 标样的外部精度为0.5‰ (2SD)。分析过程中的仪器质量分馏用Penglai锆石标准[7]进行校正, 测量的18O/16O比值通过 VSMOW值(18O/16O=0.0020052)校正后, 加上仪器质量分馏校正因子IMF即为该点的 δ 18O值:(δ 18O)M=((18O/16O)M/0.0020052-1)× 1000(‰ ), IMF=(δ 18O)M(standard)-(δ 18O)VSMOW, δ 18O样品=(δ 18O)M+IMF。

2.2 SIMS锆石 U-Pb定年

锆石U/Pb 定年在中国科学院地质与地球物理研究所离子探针实验室, Cameca IMS-1280型二次离子质谱仪(SIMS)上进行, 分析方法详见参考文献[8]。用强度为10 nA的一次 O2-离子束通过-13 kV加速电压轰击样品表面, 束斑大小约为20 μ m× 30 μ m。二次离子经60 eV的能量窗过滤, 其质量分辨率约为5 400, 用以从同量异位素干扰中分离出Pb+的峰值。用氧气流来增加样品室氧气的压力达到约5× 10-6T以提高锆石中Pb+离子的敏感度。每个样品点分析7组数据, 测量时间约为12 min。锆石标样Plé sovice和Qinghu锆石[9]交替测定。以长期监测标准样品获得的标准偏差1.5%(1σ )[10]和单点测试内部精度共同传递得到样品单点误差。锆石样品的Pb/U比值利用标准锆石91500(1 065 Ma)进行校正; Th和U含量用标准锆石91500 ( Th=29 μ g/g, U=81 μ g/g)[11]计算。普通Pb用测量的204Pb进行校正。单点分析的同位素比值及年龄误差均为1σ 。数据结果处理采用ISOPLOT软件[12]

2.3 LA-MC-ICPMS锆石Hf同位素分析

锆石Hf 同位素分析在中国科学院地质与地球物理研究所多接收等离子质谱仪实验室, 配备有Geolas-193型紫外激光剥蚀系统(LA)的Neptune 型多接收电感耦合等离子体质谱仪 (MC-ICP-MS)上进行, 测试方法详见参考文献[13]。样品分析时激光剥蚀的时间为26 s, 激光剥蚀脉冲功率为100 mJ时的脉冲频率是10 Hz, 束斑大小约为60 μ m。176Lu和176Yb对176Hf的同质异位素干扰是通过监测175Lu 和172Yb 信号强度, 采用175Lu/176Lu 比值(0.02655)和176Yb/172Yb比值(0.5886)[14]进行校正。利用标准锆石91500与锆石样品交叉分析对仪器漂移进行外部监控。在测试过程中, 锆石标样Mud Tank的176Hf/177Hf比值测定的平均值为0.282491± 27(2σ ), 锆石标样GJ的176Hf/177Hf比值测定的平均值为0.281971± 26(2σ )。在误差范围内与文献综合的结果一致[9, 15]。ε Hf值与两阶段模式年龄据参考文献[16]中方法计算, 其中(176Lu/177Hf)CHUR=0.0332, (176Hf/177Hf)CHUR, 0=0.282772, (176Lu/177Hf)DM=0.0384, (176Hf/177Hf)DM=0.28325。fcc, fs, fDM 分别为大陆地壳、样品和亏损地幔的fLu/Hft为样品形成时间, λ =1.867× 10-11a-1

3 分析结果
3.1 全岩主微量元素特征

5件基性岩脉的全岩SiO2含量为45%~53% (表2, 图3)。MgO (4.31%~7.49%)含量明显低于原始玄武岩浆成分[17], 表明岩浆经历了不同程度的演化作用。在稀土元素蛛网图中 (图4a), 基性岩脉表现为轻稀土富集的右倾模式。样品SF-5显示轻微的Eu负异常, 与邻近紫金山铜金矿交代蚀变岩的稀土元素Eu负异常[18]相似, 而且样品SF-5在薄片观察中蚀变作用明显, 因此其Eu异常可能是受到蚀变作用的影响; 而同样经历蚀变作用的样品SH-2远离成矿区, 其正Eu异常可能是不同成分蚀变热液反向作用的结果。样品SH-1与MK-1的Gd负异常与闽北赤门角闪辉长岩相似[19]。在微量元素蛛网图中 (图4b), Th, U含量从富集到亏损变化幅度较大, 出现Nb, Sr, P, Ti的亏损。样品SF-5与SH-2在标准化图中表现异常 (Eu, Th, U, P的反向变化等)。

图3 主量元素硅碱图(据参考文献[20]修改)
Pc.苦橄玄武岩; B.玄武岩; O1.玄武安山岩; O2.安山岩; O3.英安岩; R.流纹岩; S1.粗面玄武岩; S2.玄武质粗面安山岩; S3.粗面安山岩; T.粗面岩、粗面英安岩; F.副长石岩; U1.碱玄岩、碧玄岩; U2.响岩质碱玄岩; U3.碱玄质响岩; Ph.响岩; Ir.Irvine分界线, 上方为碱性, 下方为亚碱性
Fig.3 Silicon vs. alkali diagram of major elements (modified after reference [20])
Pc.Picrite basalt; B.Basalt; O1.Basaltic andesite; O2.Andesite; O3.Dacite; R.Rhyolite; S1.Trachybasalt; S2.Basaltic trachyandesite; S3.Trachyandesite; T.Trachyte and trachydacite; F.Feldspathoidite; U1.Tephrite and basanite; U2.Phonolitic tephrite; U3.Tephriticphonolite; Ph.Phonolite; Ir.Irvine boundary, upper is alkali series, lower is subalkali series

表2 闽西南基性岩脉全岩主量(%)及微量(μ g/g)元素组成 Table 2 Major (%) and trace (μ g/g) elements of mafic dikes

图4 稀土(a)微量元素(b)标准化图解
闽北赤门晚白垩世角闪辉长岩, 强不相容元素及Sr, Nd同位素显示富集岩石圈地幔源区特征[19]; 球粒陨石及原始地幔值据参考文献[17]
Fig.4 Rare earth (a) trace element (b) standardized diagram
Late Cretaceous Chimen bojite in northern Fujian, strong incompatible element and Sr, Nd isotopes show characteristics of enriched lithospheric mantle source[19]; The chondrite meteorites and the primitive mantle after reference[17]

3.2 锆石形态

基性岩脉由于结晶程度较低, 原生岩浆锆石一般数量较少, 但部分基性岩浆在源区混合或上升过程中会捕获其他成因的锆石。这些捕获锆石不仅是我们探索基性岩浆成因的重要线索, 也为研究区域的岩浆活动及岩浆源区演化特征提供依据。本次挑选出的锆石多为无色、透明的晶体。大部分锆石呈半自形到自形晶体, 长100~300 μ m, 长宽比大多为1∶ 1~3∶ 1。CL图像显示锆石内部结构主要为振荡环带和扇状环带 (图5), 为典型的岩浆结晶锆石的特征。样品MK-1中的锆石颗粒MK-1@8及MK-1@11呈浑圆状, 无环带结构, 且具有较低的Th, U含量(表3), 锆石可能受到后期的热动力作用影响。但这些锆石的Th/U比值(0.51及0.11)仍> 0.1; 且在206Pb/238U-207Pb/235U图中表现出很好的谐和性(图5), 显示早元古代的年龄特征; 其Hf, O值也与其他早元古代年龄锆石相近(表4表5); 说明这些锆石尽管受到后期各种地质作用影响, 其Pb/U及Hf, O同位素值可能仍保留了其形成时的特征。

3.3 锆石U/Pb年龄

样品YY-1测定的14颗锆石的U含量为112~963 μ g/g, Th含量为22~527 μ g/g, Th/U比值为0.05~0.73 (表3)。在206Pb/238U-207Pb/235U图解上, 数据显示2组年龄分布(图5, YY-1):206Pb/238U年龄加权平均值分别为 (107.0± 2.0) Ma和 (156.4± 2.3) Ma, 前一组年龄可能代表岩体的结晶年龄, 也可能为捕获锆石; 后一组年龄则应该是岩浆后期捕获的围岩锆石, 与这一地区广泛分布的晚侏罗世早期花岗岩相符合。

样品SF-5测定的14颗锆石的U含量为278~1 649 μ g/g, Th含量为64~933 μ g/g, Th/U比值为0.112~1.473 (表3)。在206Pb/238U-207Pb/235U 图解上, 数据显示的年龄较为分散 (图5, SF-5 ), 除样品SF-5@2显示97.5 Ma的206Pb/238U年龄, 其余锆石均显示震旦纪和中元古代年龄特征。

表3 闽西南基性岩脉锆石 U/Pb 年龄组成 Table 3 Zircon U/Pb age of mafic dikes

样品SH-1测定了13颗锆石的U/Pb年龄。锆石的U含量为159~1 552 μ g/g, Th含量为53~1 842 μ g/g, Th/U比值为0.16~1.34 (表3)。 206Pb/238U-207Pb/235U 图解上, 数据显示中生代— 中元古代的年龄特征, 年龄分布较为分散 (图5, SH-1)。

样品SH-2完成8颗锆石的U/Pb年龄测定。锆石的U含量为281~1 034 μ g/g, Th含量为179~2 087 μ g/g, Th/U比值为0.45~2.02(表3)。除1件锆石显示1 311 Ma年龄外, 在206Pb/238U-207Pb/235U 图解上, 数据主要显示2组年龄分布 (图5, SH-2):一组206Pb/238U年龄加权平均值为 (100.3± 5.7) Ma (n=4) 和一组由2个158 Ma组成的年龄值。其年龄分布与样品YY-1相似。

MK-1测定了15颗锆石的U/Pb年龄。锆石的U含量为14~794 μ g/g, Th含量为2~916 μ g/g, Th/U比值为0.11~1.42 (表3)。在206Pb/238U-207Pb/235U图解上, 大部分数据显示为早元古代的年龄特征 (图5, MK-1), 还包括一个140.4 Ma(MK-1@2)的最小年龄值和2个晚二叠世年龄值 (257.5 Ma, 267.6 Ma)。

图5 锆石的U-Pb年龄谱图和阴极发光(CL)图像
CL图像中白色圆代表U/Pb定年和Hf-O同位素分析的位置, 数字表示U/Pb年龄/ε Hf(t)值/δ 18O值
Fig.5 Zircon U-Pb age spectra and cathode luminescence (CL) images
White circle represents the analysis location of Hf-O and U/Pb isotopes in the CL image, the digital means U/Pb age/ε Hf(t)/δ 18O values

图6 基性岩脉锆石中δ 18O值正态分布图Fig.6 Normal distribution of δ 18O values in zircon

3.4 锆石O同位素组成

锆石δ 18O值 (表4) 及分布频数见图6。样品YY-1测定了35颗锆石的δ 18O值, 3颗锆石的δ 18O值为11.2~11.8, 其余为5.80‰ ~9.95‰ , 8‰ ~9‰ 有一个明显的峰值。SF-5测定的34颗锆石, 除一颗锆石的δ 18O值为11.1, 其余为5.65‰ ~9.03‰ 。SH-1测定的32颗锆石中, 除去2颗高δ 18O值锆石(10.7和12.2)和2颗低δ 18O值锆石(4.11和4.81)外, 其余为5.61‰ ~9.95‰ 。SH-2测定了22颗锆石的δ 18O值, 其值为5.79‰ ~9.37‰ 。MK-1测定了23颗锆石的δ 18O值, 除1颗锆石的δ 18O值为11.0, 其余为4.03‰ ~8.90‰ 。

表4 闽西南基性岩脉锆石O同位素组成* Table 4 Zircon O isotopes of mafic dikes

在年龄vs. δ 18O图中 (图7), 锆石的δ 18O值显示为3组数据:早元古代锆石δ 18O值与地幔值(5.3± 0.3)‰ [21, 22]相近, 显示幔源特征; 中元古代— 震旦纪的锆石δ 18O值主要为6‰ ~9‰ , 与古生代以后的锆石δ 18O值均位于地幔值之上, 只是后者有少数锆石显示出更高的锆石δ 18O值, 显示更多壳源的特征。

图7 闽西南基性岩脉锆石δ 18O值与年龄相关图
锆石δ 18O地幔值范围据参考文献[21, 22]
Fig.7 Zircon δ 18O vs. age diagram
Mantle δ 18O value after the references [21, 22]

3.5 锆石Lu-Hf 同位素组成

锆石176Lu/177Hf比值为0.0002~0.0029, 176Hf/177Hf比值为0.280798~0.282662, 以测得的锆石年龄进行计算, 初始176Hf/177Hfi比值为0.280745~0.282660, ε Hf(t)值为-16.5~+4.5 (表5)。锆石的ε Hf(t)变化随年龄趋势线不明显(图8), 在古生代之前分布于球粒陨石均一储库(CHUR) 左右及以下(偏负值)。古生代之后锆石的ε Hf(t)值都在CHUR以下, 且变化范围相对较小, 主要位于基性地壳演化线与亏损地幔之间; 古生代以前的捕获锆石位于亏损地幔与平均的中酸性地壳之间, 说明这一地区的下地壳形成基本都受到亏损幔源岩浆的影响。

表5 闽西南基性岩脉锆石Hf同位素组成 Table 5 Zircon Hf isotopes of mafic dikes

图8 闽西南基性岩脉锆石年龄与ε Hf(t)值相关图
底图据参考文献[18]
Fig.8 Zircon age vs. ε Hf(t) diagram
Base map after the reference [18]

4 讨 论

来自幔源的基性岩浆不仅携带了岩浆源区岩石的地球化学特征, 在岩浆上升过程中还会发生岩浆混合、围岩混染以及后期的蚀变作用, 使原生幔源岩浆成分发生改变。传统的全岩同位素组分研究只能对过程的结果给出评价, 而对演化过程本身的性质较难判断。在微量元素蜘蛛图中, 全岩标准化值虽然保留了基性岩浆轻稀土富集, 重稀土亏损的基本特征, 但样品之间的元素相对变化却两两不同(图4)。由于本次研究的基性岩脉均形成于晚中生代或以后, 推测不完全是地幔源区成分差异的结果。通常, 幔源岩浆在上升过程中与壳源岩浆混合或混染地壳围岩也被用来解释微量元素之间的差异。例如, MK-1捕获了最多数量的早元古代锆石, 显示最低的δ 18O值(图6), 其不相容元素标准化值与其他均有不同(图4)。但是根据以下的观察:样品YY-1与SH-2有相似的2组锆石年龄组成, 只是2件样品的后期蚀变程度不同; SH-1与SH-2岩性较为相似(表1), 但SH-2的蚀变更为强烈; SF-5显示了与SH-1, SH-2相似的锆石ε Hf(t)值与δ 18O值变化范围, 但SF-5为强蚀变辉绿岩(表1); 因此, 导致闽西南基性岩脉微量元素标准化值两两差异显著的原因, 可能主要是基性岩脉成岩后的蚀变作用。

得益于锆石原位同位素测试技术的发展, 我们可以利用岩浆中不同成因的锆石特征, 反演其演化过程。更为重要的是, 基性岩脉中大量捕获锆石的存在, 为研究该区域岩浆活动及岩浆源区演化提供了线索。

锆石对于U-Pb和Hf有较高的封闭温度, 且锆石中Pb/U及Lu/Hf比值极低, 是年龄测定及成因研究的理想矿物。但由于不同地幔储库的Hf同位素组成变化范围很大, 仅用Hf同位素组成有时不能有效限定地幔岩浆是否参与了花岗岩的形成。而锆石对O有非常好的封闭性, O同位素的研究在鉴别幔源岩浆在花岗岩形成过程中的作用提供了更多有效的制约[3]。本次研究的闽西南基性岩脉在野外观察中, 均侵入到晚侏罗世或早白垩世地层。因此, 除每件样品中测得的最小年龄(96~142 Ma)锆石可能为原生锆石外, 其余均为捕获锆石。

本次测得锆石年龄可以分为5组, 各组表现出不同的Hf-O同位素特征(图9):①2 467~1 796 Ma, 为早元古代时期形成锆石, δ 18O值(4.42~8.09, 平均为5.81)与幔源岩浆结晶锆石的δ 18O 值(5.3± 0.3)‰ [21, 22]较为接近, 显示主要为地幔来源, 而锆石的ε Hf(t)值主要分布在-7.9~+2.2, 接近于球粒陨石均一储库。说明闽西南早元古代时期地幔可能并不具亏损地幔特征。②1 343~647 Ma, 为中晚元古代— 震旦纪时期。这一时期锆石的ε Hf(t)值(-15.5~+4.5)变化较大, δ 18O值(6.2~8.8)的变化范围则相对较小, 显示亏损地幔岩浆与地壳组分混合的特征。③427~225 Ma, 志留纪— 晚三叠世时期。较高的δ 18O值, 显示典型S型花岗岩的特征。④159~140 Ma, 晚侏罗世花岗岩, ε Hf(t)值与δ 18O值变化范围较小, 显示源区或围岩类型相对简单。155~160 Ma也是与紫金山铜金矿成矿有关的两期中生代花岗岩中较早的一期[23]。⑤96~112 Ma, 早白垩世晚期的锆石ε Hf(t)值与δ 18O值清晰地显示出亏损地幔与地壳岩浆混合的趋势, 而可能不是岩石圈地幔的直接熔融[24]

图9 ε Hf(t) vs. δ 18O图解
S型花岗岩锆石的Hf-O同位素范围根据广西印支期大容山堇青石花岗岩的数据计算获得[3]
Fig.9 ε Hf(t) vs. δ 18O diagram
Zircon Hf-O isotope value of S type granite according to data calculation of indosinian cordierite granite from Darongshan, Guangxi[3]

1.8~1.9 Ga和0.96~1.1 Ga是华夏块体2个明显的新生地壳增长时期[1]。闽西南基性岩脉中, 锆石Hf的二阶段亏损地幔模式年龄TDM2的峰值主要分布在1.6~1.9 Ga(图10), 说明早元古代晚期幔源岩浆作用形成的基性岩地壳可能是后期形成花岗岩的主要源区。基性岩脉中早元古代锆石年龄的出现, 也说明早元古代晚期基性地壳对后期基性岩浆具有影响。Zhang等[25, 26]通过综合地球物理资料解译华夏块体的地壳结构, 发现华夏中地壳下部(约20 km)及下地壳底部均有基性岩层(2~5 km)存在, 说明大规模幔源岩浆底侵和侵入的存在, 并对该区花岗岩的形成一定具有热源和物质参与的影响。

图10 Hf二阶段模式年龄分布直方图Fig.10 Hf two stage model age distribution histogram

5 结论

闽西南中生代基性岩浆在岩浆上升过程中与壳源岩浆发生混合或围岩混染作用, 捕获丰富的地壳来源锆石。不同程度的蚀变作用可能对基性岩全岩成分造成明显的影响。

闽西南基性岩脉中锆石可分为5组:① 2 467~1 796 Ma, 早元古代时期; ② 1 343~647 Ma, 中晚元古代— 震旦纪时期; ③ 427~225 Ma, 志留纪— 晚三叠世时期; ④ 159~140 Ma, 晚侏罗世时期; ⑤ 112~96 Ma, 早白垩世晚期, 可能为基性岩脉的形成时期。

Hf-O同位素显示早元古代锆石来源于接近球粒陨石均一储库的地幔。中晚元古代以后年龄的锆石其Hf-O同位素均显示亏损地幔岩浆与地壳组分混合的特征; 志留纪— 晚侏罗世锆石主要来源于S型壳源花岗岩的重熔; 早白垩世晚期的锆石ε Hf(t)值与δ 18O值清晰地显示出亏损地幔与地壳岩浆混合的趋势。

闽西南基性岩脉中锆石Hf的二阶段亏损地幔模式年龄TDM2的峰值主要分布在1.6~1.9 Ga, 说明早元古代晚期幔源岩浆作用形成的基性岩地壳可能是形成后期花岗岩的主要源区。

致 谢:评审专家对本文提出了宝贵建议。样品在野外采集过程中得到福建紫金矿业集团股份有限公司祁进平工程师的帮助。Hf, O同位素测试与数据处理得到中国科学院地质与地球物理研究所杨岳衡工程师、凌潇潇工程师的帮助, 在此一并感谢!

The authors have declared that no competing interests exist.

参考文献
[1] Li Xianhua, Li Wuxian, He Bin. Buiding of the South China Block and its relevance to assembly and breakup of Rodinia supercontinent: Observations, interpretations and tests[J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2012, 31(6): 543-559.
[李献华, 李武显, 何斌. 华南陆块的形成与Rodinia超大陆聚合—裂解—观察、解释与检验[J]. 矿物岩石地球化学通报, 2012, 31(6): 543-559. ] [本文引用:2]
[2] Shao Ji’an, Zhang Yongbei, Zhang Lüqiao, et al. Early Mesozoic dike swarms of carbonatites and lamprophyres in Datong area[J]. Acta Petrologica Sinica, 2003, 19(1): 93-104.
[邵济安, 张永北, 张履桥, . 大同地区早中生代煌斑岩—碳酸岩岩墙群[J]. 岩石学报, 2003, 19(1): 93-104. ] [本文引用:1]
[3] Li Xianhua, Li Wuxian, Wang Xuance, et al. Role of mantle-derived magma in genesis of early Yanshanian granites in the Nanling Range, South China: In situ zircon Hf-O isotopic constraints[J]. Science in China(Series D), 2009, 39(7): 872-887.
[李献华, 李武显, 王选策, . 幔源岩浆在南岭燕山早期花岗岩形成中的作用: 锆石原位Hf-O同位素制约[J]. 中国科学: D辑, 2009, 39(7): 872-887. ] [本文引用:3]
[4] Fujian Province Geology and Mineral Exploration and Development Bureau. 1∶ 500000 Geological Map of Fujian Province[Z]. Fuzhou: Fujian Province Map Publishing House, 1998.
[福建地质矿产勘查开发局. 福建1∶ 50万地质图(1998年版)[Z]. 福州: 福建省地图出版社, 1998. ] [本文引用:1]
[5] Xiong Xin, Xu Wenyi, Jia Liqiong, et al. Reviews of structural settings of porphyry copper deposits in China[J]. Advances in Earth Science, 2014, 29(2): 250-264.
[熊欣, 徐文艺, 贾丽琼, . 斑岩铜矿成矿构造背景研究进展[J]. 地球科学进展, 2014, 29(2): 250-264. ] [本文引用:1]
[6] Li X, Li Z, He B, et al. The Early Permian active continental margin and crustal growth of the Cathaysia Block: In situ U-Pb, Lu-Hf and O isotope analyses of detrital zircons[J]. Chemical Geology, 2012, 328: 195-207. [本文引用:1]
[7] Li X H, Long W G, Li Q L, et al. Penglai zircon megacrysts: A potential new working reference material for microbeam determination of Hf-O Isotopes and U-Pb age[J]. Geostand ards and Geoanalytical Research, 2010, 34: 117-134. [本文引用:1]
[8] Li X H, Liu Y, Li Q, et al. Precise determination of Phanerozoic zircon Pb/Pb age by multicollector SIMS without external stand ardization[J]. Geochemistry, Geophysics, Geosystems, 2009, 10(4): 1-21. [本文引用:1]
[9] Li Xianhua, Tang Guoqiang, Gong Bing, et al. Qinghu zircon: A working reference for microbeam analysis of U-Pb age and Hf and O isotopes[J]. Chinese Science Bulletin, 2013, 58(20): 1 954-1 961.
[李献华, 唐国强, 龚冰, . Qinghu(清湖)锆石: 一个新的U-Pb年龄和O-Hf同位素微区分析工作标样[J]. 科学通报, 2013, 58(20): 1 954-1 961. ] [本文引用:2]
[10] Li Q, Li X, Liu Y, et al. Precise U-Pb and Pb-Pb dating of Phanerozoic baddeleyite by SIMS with oxygen flooding technique[J]. Journal of Analytical Atomic Spectrometry, 2010, 25(7): 1 107-1 113. [本文引用:1]
[11] Wiedenbeck M, Allé P, Corfu F, et al. Three natural zircon stand ards for U-Th-Pb, Lu-Hf, trace element and REE analyses[J]. Geostand ards and Geoanalytical Research, 1995, 19(1): 1-23. [本文引用:1]
[12] Ludwig K R. User’s Manual for Isoplot/Ex Rev. 3. 75. A Geochronological Toolkit for Microsoft Excel[Z]. Berkeley: Berkeley Geochronological Center, 2012. [本文引用:1]
[13] Wu F, Yang Y, Xie L, et al. Hf isotopic compositions of the stand ard zircons and baddeleyites used in U-Pb geochronology[J]. Chemical Geology, 2006, 234(1/2): 105-126. [本文引用:1]
[14] Chu N, Taylor R N, Chavagnac V R, et al. Hf isotope ratio analysis using multi-collector inductively coupled plasma mass spectrometry: An evaluation of isobaric interference corrections[J]. Journal of Analytical Atomic Spectrometry, 2002, 17(12): 1 567-1 574. [本文引用:1]
[15] Geng Jianzhen, Li Huaikun, Zhang Jian, et al. Zircon Hf isotope analysis by means of LA-MC-ICP-MS[J]. Geological Buletin of China, 2011, 30(10): 1 508-1 513.
[耿建珍, 李怀坤, 张健, . 锆石Hf同位素组成的LA-MC-ICP-MS测定[J]. 地质通报, 2011, 30(10): 1 508-1 513. ] [本文引用:1]
[16] Wu Fuyuan, Li Xianhua, Zheng Yongfei, et al. Lu-Hf isotopic systematics and their applications in prtrology[J]. Acta Petrologica Sinica, 2007, 2(23): 185-220.
[吴福元, 李献华, 郑永飞, . Lu-Hf同位素体系及其岩石学应用[J]. 岩石学报, 2007, 2(23): 185-220. ] [本文引用:1]
[17] McDonough W F, Sun S S. The composition of the Earth[J]. Chemical Geology, 1995, 120: 223-253. [本文引用:1]
[18] Zhang Feifei, Wang Yuejun, Fan Weiming, et al. Zircon U-Pb geochronology and Hf isotopes of the neoproterozoic granites in the central of Jiangnan Uplift[J]. Geotectonica et Metallogenia, 2011, 35(1): 73-84.
[张菲菲, 王岳军, 范蔚茗, . 江南隆起带中段新元古代花岗岩锆石U-Pb年代学和Hf同位素组成研究[J]. 大地构造与成矿学, 2011, 35(1): 73-84. ] [本文引用:1]
[19] Zhang Guishan. Chronology, Geochemistry and Geodynamic Significance of the Mafic-ultramafic Rocks in Fujian Province Since Late Mesozoic[D]. Beijing: Institute of Geochemistry, Chinese Academy of Sciences, 2006.
[张贵山. 福建晚中生代以来基性—超基性岩的年代学、地球化学及其地球动力学意义[D]. 北京: 中国科学院地球化学研究所, 2006. ] [本文引用:1]
[20] Le Maitre R W, Bateman P, Dudek A, et al. Classification of Igneous Rocks and Glossary of Terms[M]. Oxford: Blackwell, 1989. [本文引用:1]
[21] Valley J W, Kinny P D, Schulze D J, et al. Zircon megacrysts from kimberlite: Oxygen isotope variability among mantle melts[J]. Contributions to Mineralogy and Petrology, 1998, 133(1/2): 1-11. [本文引用:2]
[22] Valley J W, Lackey J S, Cavosie A J, et al. 4. 4 billion years of crustal maturation: Oxygen isotope ratios of magmatic zircon[J]. Contributions to Mineralogy and Petrology, 2005, 150(6): 561-580. [本文引用:2]
[23] Yu Bo, Pei Rongfu, Qiu Xiaoping, et al. The evolution series of mesozoic magmatic rocks in the Zijinshan orefield, Fujian Province[J]. Acta Geoscientica Sinica, 2013, 34(4): 437-446.
[于波, 裴荣富, 邱小平, . 福建紫金山矿田中生代岩浆岩演化序列研究[J]. 地球学报, 2013, 34(4): 437-446. ] [本文引用:1]
[24] Li B, Jiang S, Zhang Q, et al. Cretaceous crust-mantle interaction and tectonic evolution of Cathaysia Block in South China: Evidence from pulsed mafic rocks and related magmatism[J]. Tectonophysics, 2015, 661: 136-155. [本文引用:1]
[25] Zhang Z, Badal J, Li Y, et al. Crust-upper mantle seismic velocity structure across Southeastern China[J]. Tectonophysics, 2005, 395(1/2): 137-157. [本文引用:1]
[26] Zhang Z, Zhang X, Badal J. Composition of the crust beneath southeastern China derived from an integrated geophysical data set[J]. Journal of Geophysical Research, 2008, 113(B4): 1-25. [本文引用:1]