磷灰石单颗粒(U-Th)/He测年实验流程的建立及验证
王英, 郑德文, 武颖, 李又娟, 王一舟
中国地震局地质研究所, 地震动力学国家重点实验室, 北京 100029

〔作者简介〕 王英, 女, 1984年生, 2014年于中国科学院地质与地球物理研究所获地球化学博士学位, 工程师, 主要从事(U-Th)/He、 40Ar/39Ar 年代学和稀有气体同位素研究, 电话: 010-62009091, E-mail: wangying@ies.ac.cn

摘要

(U-Th)/He同位素定年是近几十年快速发展起来的1种能够对岩石低温冷却历史进行有效研究的新型放射性同位素定年方法。由于磷灰石的 He 封闭温度较低, 所以它能反映低温阶段(40~80℃)的热历史信息, 在低温热年代学领域具有很好的应用前景。但是, 由于磷灰石(U-Th)/He 测年的影响因素多且测试流程复杂, 目前国内对于(U-Th)/He 同位素定年方法和流程的研究还处于起步阶段。依托中国地震局地质研究所新建立的(U-Th)/He 同位素年代学实验室, 利用配备的Alphachron He同位素质谱仪对4批共75个国际标样Durango磷灰石单颗粒进行激光熔融和4He的含量测量, 随后应用自动进样安捷伦7900 ICP-MS和同位素稀释剂法测定母体同位素 U、 Th 的含量, 最后计算得到Durango磷灰石的Th/U同位素比值在17.23~23.60之间, 全部年龄分布在28.61~34.51Ma, 平均年龄为(31.71±1.55)Ma(1σ), 与国际标定年龄在误差范围内一致。

关键词: (U-Th)/He定年; 激光熔融法; 同位素稀释剂法; Durango磷灰石
中图分类号:P597 文献标志码:A 文章编号:0253-4967(2017)06-1143-15
MEASUREMENT PROCEDURE OF SINGLE-GRAIN APATITE(U-Th)/He DATING AND ITS VALIDATION BY DURANGO APATITE STANDARD
WANG Ying, ZHENG De-wen, WU Ying, LI You-juan, WANG Yi-zhou
State Key Laboratory of Earthquake Dynamics, Institute of Geology,China Earthquake Administration, Beijing 100029, China
Abstract

(U-Th)/He isotopic dating has been developed very quickly in recent years, due to the recognition that the thermal history of rock at low temperature can be effective revealed by such dating method. In particular, He closure temperature in apatite(40~80℃) is very low, so apatite(U-Th)/He ages can reflect the thermal history information of the low-temperature stage, and have a good application prospect in the field of low-temperature thermal chronology. However, because of many influence factors and complicated measurement procedures, the development of apatite He dating in China remains in its early stage. In this study, a measurement procedure was established at the(U-Th)/He dating laboratory of Institute of Geology, China Earthquake Administration. We measured the daughter isotopic helium by diode laser heating four batches of a total seventy-five grains of Durango apatite in an Alphachron helium mass spectrometry system. Then the apatite grains were dissolved to precisely measure the concentration of parent nuclides(U, Th)using the solution isotope dilution method through an automatic sampling ICP-MS(Agilent 7900). Results show that the Th/U values of Durango apatite grains were in the range of 17.23 to 23.60, while all the 75ages were in the range of 28.61 to 34.51Ma with an average of(31.71±1.55)Ma(1σ), which are consistent with the international calibrated ages.

Keyword: (U-Th)/He dating; laser heating; isotope dilution method; Durango apatite
0 引言

磷灰石广泛存在于不同类型的岩石矿物中。由于矿物富含U和Th元素, 因此被广泛应用于各种地质同位素定年, 如U-Pb定年、 裂变径迹和(U-Th)/He定年。磷灰石的(U-Th)/He定年是近几十年快速发展起来的1种同位素定年方法。Zeitler等详细研究了磷灰石中的He扩散, 提出磷灰石(U-Th)/He是潜在的低温热年代计(Zeitler et al., 1987)。随后(U-Th)/He定年技术得到快速发展, 尤其是美国加利福尼亚州理工学院的Farley 教授自20世纪90年代起一直致力于磷灰石(U-Th)/He定年技术的研究, 做了大量的基础工作(Farley et al., 1996, 2002; Farley, 2000)。由于磷灰石的He 封闭温度较低, 能反映近地表更低温阶段(40~80℃)的热历史信息, 所以磷灰石的(U-Th)/He定年被广泛应用于地质体定年、 造山带演化、 盆地热历史恢复、 构造隆升历史、 古地形、 金属和非金属矿床剥露时间等研究(郑德文等, 2000; Farley et al., 2002; 保增宽等, 2005; 常健等, 2012; 蒋毅等, 2012; 王杰等, 2012; Dai et al., 2013; Tripathy-lang et al., 2013; 郑德文等, 2016)。

(U-Th)/He定年的基本原理与其他放射性同位素定年一样, 利用母体同位素238U、235U、232Th 和 147Sm 随时间的衰变和子体同位素4He随时间的积累来测定矿物的形成或冷却时间。子体同位素4He的累积方程为

4He=8× 238U(eλ 238t-1)+7× 235U( eλ235t-1)+6× 232Th( eλ232t-1)+ 147Sm(eλ 147t-1) (1)

式(1)中, 4He238U235U232Th147Sm 为矿物中对应的同位素原子数, λ 238λ 235λ 232λ 147 分别为238U、235U、232Th和 147Sm 的衰变常数, t为矿物中的He同位素累积时间。由于矿物中的Sm含量通常较低, 半衰期较长, 因此 147Sm 对4He的贡献率通常< 1% (Farley et al., 2002); 只有当矿物中的U含量< 5μ g/g时, 147Sm 对4He的贡献率才会> 5% (Reiners et al., 2006)。所以, 大多数实验室在测量母体同位素浓度时不测量 147Sm 的含量。式(1)中的子体同位素4He可以通过He提取和测量系统进行熔融和测量, 母体同位素238U、235U和232Th可以通过ICP-MS测试得到, 而各母体同位素衰变常数是已知的, 通过计算可以得到子体同位素He的累积时间t

目前磷灰石(U-Th)/He定年已经作为1种严格、 定量的技术广泛应用于低温热年代学和低温热演化的研究。过去这一技术在国外实验室的研究相对较多, 许多相关研究单位都建立了(U-Th)/He年代学实验室。但是由于磷灰石单颗粒(U-Th)/He测年实验流程相对复杂, 需要激光加热、4He气体提取纯化、 颗粒酸溶等一系列过程, 同时还需要配备高灵敏度的质谱仪准确测定子体同位素 4He 和母体同位素 238U、 235U 和 232Th 的含量, 以前国内样品的(U-Th)/He年龄大多在国外测试完成。然而近年来, 随着(U-Th)/He定年技术的发展和广泛应用, 国内多个单位相继建立了(U-Th)/He同位素年代学实验室, 并致力于开展有关的(U-Th)/He 测试技术研究(陈文等, 2010; 张彦等, 2011a, b; 孙敬博等, 2015; 吴林等, 2016)。

产自墨西哥的Durango磷灰石是国际上(U-Th)/He 年代学实验室广泛使用的标准样品。这是由于Durango磷灰石基本没有包裹体, 容易获得粒径较大的颗粒。Zeitler等(1987)、 Wolf等(1996)和Farley(2000)利用Durango磷灰石详细研究了He在磷灰石中的扩散, 发现不同颗粒的He的扩散行为相似, 可作为不同实验室之间的标样。目前不同实验室对于Durango磷灰石测得的(U-Th)/He的平均年龄为31~33Ma(House et al., 2000; Evans et al., 2005; McDowell et al., 2005; Foeken et al., 2006; Reiner et al., 2006; 吴林等, 2016)。

本文详细报道了中国地震局地质研究所2015年新建立的(U-Th)/He 同位素年代学实验室的仪器本底水平、 测试精度和稳定性等指标, 已经达到国际同类实验室的标准, 在此基础上, 首先使用Durango磷灰石单颗粒作为标样, 系统地建立单颗粒样品激光熔融, 子体同位素 4He 的提取, 纯化和测量, 母体同位素 U、 Th 含量的稀释剂法质谱测定, 年龄的计算等实验流程。计算得到75颗Durango磷灰石的平均年龄为(31.71± 1.55)Ma(1σ ), 与国际标定年龄在误差范围内一致, 表明实验流程的准确可行。

1 仪器介绍
1.1 Alphachron He同位素质谱仪

Alphachron He同位素质谱仪是澳大利亚科学仪器公司(Australian Scientific Instruments Pty Limited, 简称ASI)的产品, 专门用于提取、 纯化和测量矿物单颗粒或多颗粒4He同位素的含量。该系统由Labview软件控制, 全自动运行, 主要由气体提取系统、 纯化系统、4He气体瓶、3He稀释剂瓶以及四级杆质量分析器等5部分组成:

(1)气体提取系统。 该系统包括装载在同一X-Y可调轨道底座上的970nm二极管激光器和相应的透镜聚焦系统、 电源、 以及能容纳25个样品颗粒的固定激光样品室和蓝宝石窗口。

(2)高真空气体纯化系统。 最高的真空水平< 2.67× 10-8Pa。由机械泵、 涡轮分子泵、 离子泵3级真空系统组成。2组SAES AP10N Mk3型Zr-Al泵分别在高温(350℃)和室温下工作, 用于去除样品中释放的化学活性气体, 如H2O、 OH、 N2、 H2、 H、 HD、 SO2、 O2、 CO2、 CO和12C、 CH、 CH2、 CH3、 CH4、 C2H6等碳氢化合物。该系统具有很低的本底水平, 4He本底< 6× 10-12ml。

(3)4He气体瓶和3He稀释剂瓶。 该仪器共有2瓶高纯4He标准气体(分别简称Q-tank, D-tank)和1瓶高纯3He稀释剂气体。每个气瓶与真空管线通过1段约0.3ml体积的pipette管线和2个气动阀门进行连接和取样。分别用来量取标准4He气体和3He稀释剂气体。

(4)四级杆质量分析器。分析的质量范围为2~260amu, 分辨率可以调节并且任意质量数10%峰高处的峰宽应为0.65~0.80amu。经2组Zr-Al泵纯化后的气体引入四级杆质量分析器进行同位素分析。

1.2 Agilent 7900 ICP-MS

Agilent 7900系列ICP-MS是安捷伦公司2014年推出的1款质谱仪。它降低了所有原来ICP-MS所面临的分子离子干扰, 主要由以下几个部分组成:

(1)样品引入装置。 样品引入包括样品提升和样品雾化。本台ICP-MS样品提升使用ASX-520自动进样器和蠕动泵相结合的方式完成。使用标准的1.02mm内径的样品管时, 在0.1r/s转速下, 蠕动泵提升样品的能力大约为0.4ml/min。

(2)离子源。 本仪器中使用的离子源是感应耦合等离子体(Inductively Coupled Plasma, 简称ICP)。高纯Ar气由多个单独管线提供, 其中连接到炬管的等离子体气和辅助气流量一般为15ml/min、 1ml/min, 用来维持等离子体的形成和形状。

(3)采样锥和截取锥。 本台仪器的采样锥底座由黄铜制成, 锥尖为镍合金, 孔径为1mm; 截取锥材质为镍合金, 孔径为0.45mm。

(4)离子传输透镜。 包括提取透镜, 轴向偏转透镜和八极杆反应系统。离子束在该部分被提取, 并轴向偏转与中性粒子分开。八极杆碰撞反应池系统(ORS)充入高纯He气可极大地消除分子、 离子的干扰, 而加在八极杆上的电场使得离子束沿中心轴飞行时收缩变窄, 极大地提高了离子传输效率。

(5)双曲面四级杆质量分析器。 采集过程中四级杆连续扫描质量范围, 大约0.1s内, 可采集0~260amu之间全部的质量数。

(6)电子倍增器(EM)。 主要用来放大(增益)到达检测器的离子数量, 增益倍数最高可达108

(7)真空系统。 系统配备了1个机械泵和1个分子涡轮泵, 机械泵除了作为分子涡轮泵的前级真空泵, 还维持接口区的真空; 分子涡轮泵用于中间区和分析器区, 即从离子传输透镜到四级杆质量分析器之间的区域。

1.3 仪器的检测限、 本底水平和稳定性分析

ICP-MS的仪器检测限DL=3× sd/a, 其中sd代表连续11次测量MilliQ超纯水(电阻率> 18.2MΩ )所得到的238U或者232Th计数值的标准偏差; a代表相应同位素的灵敏度。a值通过以下方法获得: 配置7个不同浓度的标准溶液, 浓度分别是0μ g/L、 0.5μ g/L、 1μ g/L、 5μ g/L、 10μ g/L、 20μ g/L、 50μ g/L, 利用ICP-MS获得计数值cps和浓度μ g/L 的关系式, 即标准曲线的关系式y=a× x+b, 其中a值即为灵敏度值, 不同元素的a值不同, 所以不同元素的检测限DL亦不同。

为了更好地了解本实验室Alphachron He同位素质谱仪的本底水平和稳定性, 我们对其本底进行了1.5a以上的长期监测。在此期间共测试27批样品(1批装载25个样品, 一般情况下含有1个空Pt囊, 3个Durango磷灰石和21个地质样品), 包括地质样品和本文中介绍的Durango磷灰石。这里所涉及的本底主要有2种, 分别是冷本底(Cold Blank简称CB)和热本底(Pt Hot Blank, 简称Pt HB 或 HB)。冷、 热本底与磷灰石样品测试时的流程完全一致, 冷本底是激光输出能量为0时, 样品盘和管线的4He含量; 热本底是采用和样品测试一样的激光能量即8A对空Pt囊进行5min持续加热所获得的4He含量。1批样品一般可获得至少11个冷本底值和1个热本底值。

2 实验流程
2.1 样品的准备

实验室测试用的标准样品Durango磷灰石是较大的晶体去除了20μ m的外表面破碎得到的碎片。在双目正交偏光显微镜下挑选形状较好、 不含包裹体和裂隙以及晶体宽度> 80μ m的磷灰石单颗粒。挑选出的磷灰石颗粒用纯度99.96%的Pt囊包裹(长度:(0.99± 0.099)mm; 内径:(0.64± 0.01)mm; 外径:(0.74± 0.01)mm)。Pt囊包裹的优势在于: 1)熔点高, 在激光加热磷灰石提取4He的过程中不会变形; 2)良好的热传导性, Pt囊能快速传递激光能量有助于磷灰石中4He的萃取; 3)便于样品转移, 被提取4He的磷灰石颗粒需要从激光样品盘转移至聚四氟乙烯(PFA)材质的溶样小瓶中, Pt包裹可有效避免颗粒因激光加热碎裂而导致的不完整转移。为了保证激光束能作用在Pt囊平面上提高加热效率, 需要将装入样品的Pt囊两侧开口平行捏紧且留少许缝隙(便于4He气体扩散以及后续酸溶提取U、 Th过程中酸的渗入), 整体呈稍扁平的方块型。

2.2 4He的提取、 纯化与分析

将Pt囊包裹的磷灰石单颗粒放入Alphachron He同位素质谱仪的激光样品室, 设置激光的工作电流为8A(约1, 000℃), 加热时间为5min。每个样品都进行2次4He提取实验, Durango磷灰石的多次测量结果显示, 第1次加热5min即可释放 > 99% 的4He气体, 第2次加热5min释放出的4He气体含量几乎等同于本底水平。提取的4He与稀释剂3He混合后, 经过SAES AP10N 锆铝泵纯化2min, 以除去H2、 CO、 CO2、 SO2和H2O等活性气体。纯化后的气体引入四级杆质量分析器进行同位素分析。 4He含量的测定使用同位素稀释法。稀释剂(Spike)为纯净的3He, Q-tank中的标准气体(Qstandard)为纯净的4He, 2种气体均取1 pipette(约0.3cc)体积混合后, 测试得到(4He/3He)Spiked Qstandard 比值, 而1 pipette中4HeQstandard 绝对含量的计算公式为

4HeQstandard=3HeSpike × (4He/3He)Spiked Qstandard (2)

当样品中的4He与1 pipette3He混合后, 样品中4Hesample 绝对含量的计算公式为

4Hesample=3HeSpike × (4He/3He)Spiked Sample (3)

联合方程(2)和(3)可以得到样品中的4Hesample 绝对含量:

4Hesample=4HeQstandard× [(4He/3He)Spiked Sample/(4He/3He)Spiked Qstandard] (4)

式(2)— (4)中, 4He/3He比值可通过仪器测量获得, 4HeQstandard 为第N次充入pipette中4He气体的绝对含量, 简称QN。从标准4He气体瓶(VQtank≈ 3, 300ml)每取出1 pipette体积(Vpipette≈ 0.3ml)的气体, 瓶中剩余的4He绝对含量相应减少, 减少程度与气体瓶和pipette体积有关, 定义减少系数为DF(Depletion Factor), 则DF=VQtank/(VQtank+Vpipette), 因此QN=QI× DFN-I, 其中NI均是Q的测量次数, DF的准确数值通过2个标准4He气体瓶即Q-tank和T-tank相互校正得来。最终用理想气体状态方程(PV=nRT)将获得的气体绝对体积转化成物质的量, 其中T取室温20℃(293K), 压强P用1个标准大气压, R为理想气体常数(8.314J· mol-1· K-1)。

2.3 U、 Th分析

磷灰石颗粒的母体同位素U、 Th的分析和计算亦采用稀释剂法(Evans et al., 2005)。我们所使用的稀释剂U、 Th浓度分别约为15ng/mL和5ng/mL, 溶解于7mol/L的硝酸溶液中, 稀释剂中的同位素比值 235U/238U=(838± 7), 230Th/232Th=(10.45± 0.5)。为了避免PFA溶样瓶引入杂质, 需对其进行彻底清洗。在化学超净实验室的通风橱中先后使用50%硝酸、 MilliQ超纯水、 50%盐酸、 MilliQ超纯水于100℃加热板上蒸煮6h, 最后进行烘干。每次酸蒸煮之后均先用超纯水超声清洗15min, 接着用流动的超纯水对每个小瓶进行冲洗, 再换上MilliQ超纯水蒸煮。将提取完4He的磷灰石颗粒连同Pt囊一起放入溶样瓶中, 加入25μ L稀释剂, 在超声波中超声震荡15min, 促进稀释剂中的硝酸溶液通过Pt囊缝隙溶解磷灰石颗粒中的U、 Th。超声后的溶液随后放置在恒温100℃的加热板上1h, 再在室温下静置4h, 保证磷灰石颗粒中的U、 Th完全溶解。最后加入MilliQ超纯水稀释到2.5mL。通过自动进样器和蠕动泵使待测溶液进入雾化器, 雾化后形成的气溶胶通过炬管进入等离子体蒸发、 解离和原子化, 最终电离成待测离子, 电子倍增器测定230Th、232Th、234U、235U 和238U等同位素离子的计数, 从而获得 235U/238U和 230Th/232Th 的比值。本台ICP-MS在测试时的参数设置见表1。在测量时, 一般会选取He测量模式即在八极杆碰撞反应池中充入高纯He气, 尽可能减少产生的 195Pt40Ar 和 198Pt40Ar 等分子离子对同质量数235U和238U测试的干扰。每2个样品之间均用5%体积的HNO3和超纯水先后清洗进样管90s。

表1 ICP-MS7900测试参数 Table1 Measurement parameters of ICP-MS7900

每批进行ICP-MS测定的样品需要加入2个试剂空白(Blank, 7M HNO3), 2个MilliQ超纯水, 1、 2个Pt囊空白(7M HNO3+空Pt囊+水), 3个加入稀释剂的标准溶液(标准溶液U、 Th浓度均为24.08 ng/mL)。通过测试标准溶液计算稀释剂中的238U和232Th的原子数来计算未知样品中的U和Th的含量。

测试样品中的238U和232Th的原子数能够根据式(5)、(6)计算:

238Uatmos=238Usp× (RUsa-RUsp)/(RUmin-RUsa) (5)

232Thatmos=232Thsp× (RThsa-RUsp)/RThsa (6)

式(5)、(6)中, RUsa=(235U/238U)sa, RThsa=(230Th/232Th)sa, 是加入稀释剂的样品溶液中的同位素的原子数之比, 此值可以通过ICP-MS测量得到; RUsp 是稀释剂中的同位素原子数之比, 235U/238U=(838± 7), 230Th/232Th=(10.45± 0.5); RUmin=(235U/238U)min=1/137.88, 是自然界中矿物的同位素比值。

稀释剂中的238U和232Th原子数依据式(7)、(8)计算:

238Usp=238Ustd× (RUstd-RUsst)/(RUsst-RUsp) (7)

232Thsp=232Thstd× RThsst/(RThsp-RUsst) (8)

式(7)、(8)中, RUstd=(235U/238U)std, 标准溶液中的同位素比值, 由于Th是1种单同位素元素, 自然成因的体系中Th的同位素比值为零, 因此假定标准溶液中的(230Th/232Th)std 为0; RUsst=(235U/238U)sstRThsst=(230Th/232Th)sst, 是加入稀释剂的标准溶液中的同位素比值。

标准溶液中的238U和232Th的原子数可以根据238U和232Th在其中的物质的量乘以阿伏伽德罗常数求得。那么标准溶液中的238U摩尔数和232Th可以依据式(9)、(10)计算:

238Ustd=(UCstd/238)× UVstd× 103× (1/RUstd)/((1/RUstd)+1) (9)

232Thstd=(ThCstd/232)× ThVstd× 103 (10)

UCstThCstd 是经过校正后的标准溶液中238U和232Th的浓度, 单位ng/mL; UVstThVstd 是标准溶液的测试体积, 单位μ L。

2.4 年龄的计算

最后将计算得到的4He、238U和232Th的物质的量代入式1就可以计算出测定年龄。计算软件为HelioPlot version 2.0。Durango磷灰石因为不需要α 粒子射出效应的校正, 因此测定年龄为其实际年龄。

3 结果与讨论
3.1 ICP-MS的检测限和灵敏度

在碰撞反应池充入高纯He气时, 测量模式为He气模式。 He模式测量条件下, 本实验室用10mg· L-1 的混标溶液(内含238U和232Th)和超纯水分别配置了上述7种浓度的标准溶液, 获得的标准曲线如图1所示。通过标准曲线可以得到238U的灵敏度> 150, 000cps/μ g· L-1, 232Th的灵敏度> 140, 000cps/μ g· L-1。11个超纯水样品的238U和232Th的测量值的平均标准偏差为8.37cps和9.61cps。经过计算得到本台ICP-MS238U和232Th的仪器检测限分别为0.163ng· L-1 和0.203ng· L-1

图 1 He气模式下ICP-MS的灵敏度测试Fig. 1 Sensitivity test of ICP-MS with He mode.

3.2 4He的本底分析

图2是本实验室自建立以来记录的近1.5a时间内仪器真空管线中4He冷本底和热本底值。4He冷本底为2.6~5.9× 10-12ml。第1— 13批样品测试期间, 冷热本底值的波动主要由真空管线内壁不同程度的释气作用导致; 随着仪器运行时间的延长, 第14批样品测试之后, 释气作用不再明显, 真空度逐渐下降, 最后趋于稳定。相应的冷、 热本底值也趋于稳定, 其中4He的静态冷本底值基本稳定在3× 10-12ml(16批次以后, 仪器运行时间约1a), 4He的热本底基本稳定在2.8× 10-12ml。本次实验采用的Durango磷灰石单颗粒都是在实验仪器的冷本底和热本底趋向于稳定时进行测量的, 测量得到的4He含量都> 0.56× 10-9ml, 约是热本底的200~2, 878倍, 因此本底校正对于测试结果的影响可以忽略不计。

图 2 Alphachron He同位素质谱仪的27批次样品冷热本底测试Fig. 2 Measurements of cold and hot backgrounds for twenty-seven batches of samples in Alphachron helium mass spectrometry system.

3.3 U和Th的本底分析

U、 Th的本底主要由存放稀释剂的硝酸本底和Pt囊的本底2部分组成。本次实验测得存放稀释剂的硝酸中238U平均计数值约32cps, 232Th平均计数值约11cps(表2)。而Durango磷灰石和稀释剂混合后的238U平均计数值> 13, 000cps, 232Th平均计数值> 260, 000cps, 因此试剂本底和稀释样品用的MilliQ超纯水对最终的U、 Th含量测试几乎没有影响。

表2 水和HNO3试剂本底的平均计数值 Table2 Average counts per second(cps) for water and HNO3 reagent blanks

表3是4个批次的Pt囊的238U和232Th的本底分析结果。Pt囊的238U值大部分超出了仪器检出的最低限值, 可以忽略不计。Pt囊的232Th的本底值在10-14 级别, 约占样品的0.1%~1%, 保证本底的校正对于测试精度的影响非常小。

表3 不同批次Pt囊的U、 Th本底值测试结果 Table3 Measurements of backgrounds for 238U, 232Th and 4He in empty Pt foil
3.4 Durango磷灰石单颗粒的年龄

表4是75个Durango磷灰石颗粒的238U、232Th、 232Th/238U和4He同位素以及年龄的测试结果。由于不同颗粒的U、 Th含量存在差异, 而且质量太小(几μ g), 难以计算出U、 Th的绝对含量, 因此常以Th/U同位素比值来判断母体同位素在加热过程中的丢失程度。Th/U的概率分布曲线(图3)显示Th/U同位素比值分布在17.23~23.60之间, 符合正态分布, 平均值为 20.85± 1.27(1σ )。

图 3 75个Durango磷灰石颗粒的 Th/U同位素比值概率分布图Fig. 3 Probability distribution for Th/U value of 75 Durango apatite grains.

表4 中国地震局地质研究所(U-Th)/He实验室Durango磷灰石定年结果 Table4 Ages of Durango apatite measured at(U-Th)/He laboratory of Institute of Geology, China Earthquake Administration

图4是Durango磷灰石颗粒的年龄与Th/U的关系图。4个批次的Th/U比值虽然存在一定的差异, 但是所对应的年龄分布相差不大, 表明Durango磷灰石颗粒Th/U比值差异是由于样品本身的差异造成的。4个批次Durango磷灰石颗粒的平均年龄分别为(31.74± 1.58)Ma(1σ , n=20)、 (31.17± 1.32)Ma(n=18)、 (31.84± 1.69)Ma(n=23)和(32.16± 1.52)Ma(n=14)。每个批次的年龄分布都很接近, 在误差范围内基本一致, 代表批次的可重复性。最后将75个Durango磷灰石颗粒进行统一计算, 得到平均年龄为(31.71± 1.55)Ma(1σ )。图5是75个Durango磷灰石颗粒的概率分布图, 年龄结果分布在28~35之间, 其中60个颗粒分布在30~34之间。正态分布拟合曲线显示基本服从正态分布, 稍微向更老的年龄偏离。

图 5 75个Dugango磷灰石颗粒的年龄概率分布图Fig. 5 Probability distribution for ages of 75 Durango apatite grains.

图 4 4个批次的75个Durango磷灰石颗粒的He年龄与Th/U比值的关系图Fig. 4 He age versus Th/U for 4 batches of 75 Durango apatite grains.

表5是本实验室和其他实验室测得的Durango磷灰石的年龄及Th/U值的对比。House等(2000)首先使用Nd-YAG激光器间接加热11个Pt囊包裹的Durango磷灰石单颗粒, 测得Th/U分布在15~23之间, 平均的(U-Th)/He年龄为(32.1± 1.7)Ma(1σ )。与传统加热炉相比, 激光间接加热法更加迅速, 产生更低的He背景值。McDowell等(2005)测定了24个Durango磷灰石单颗粒, 得到Th/U分布在16~20之间, 平均的(U-Th)/He年龄为(31.13± 1.01)Ma(1σ )。他们同时还测试了Durango磷灰石产出层位相邻地层中透长石的 40Ar/39Ar 年龄, 将Durango磷灰石的年龄限定在(31.44± 0.18)Ma(2σ )。Evans等(2005)测试了70个Durango磷灰石颗粒, 得到平均年龄为(31.5± 1.6)Ma(2σ )。Foeken et al.(2006)利用二极管激光加热来提取He, 测试了8个阿姆斯特丹自由大学(VU)来源的Durango和3个加利福尼亚州理工学院(CIT)来源的 Durango磷灰石颗粒, 得到Th/U分布在17~25之间, 平均年龄为(32.8± 1.8)Ma(1σ )。Reiners 等(2006)测试了169颗Durango磷灰石颗粒, 发现所有颗粒的Th/U分布在15~22之间, 平均的(U-Th)/He年龄为(31.9± 2.2)Ma(2σ )。最近, 中国科学院地质与地球物理研究所的吴林等(2016)测试了40个Durango磷灰石颗粒, 得到Th/U分布在16.43~23.72之间, 平均年龄为(31.61± 2.7)Ma(2σ )。这些数据表明, 本实验室测得的Durango磷灰石不仅Th/U分布范围与其他实验室具有可比性, 而且(U-Th)/He年龄在误差范围内一致, 表明新建立实验流程的可靠性。

表5 Durango磷灰石颗粒年龄与Th/U同位素比值的国内外对比 Table5 Comparison of(U-Th)/He age and Th/U of Durango apatite from different institutions in the world

虽然Durango磷灰石颗粒不含包裹体, 容易获得较大颗粒, 避免α 粒子射出效应校正的误差, 具有可重复的He扩散行为, 是1个公认的国际标样(Zeitler et al., 1987, Wolf et al., 1996; Farley, 2000)。但是我们的定年实验表明, 不同颗粒年龄的2个标准偏差在9%。House等(2000)对比了传统炉子加热和激光加热的(U-Th)/He定年结果, 发现这2个结果具有8%的偏差。Reiners等(2006)的结果显示169个Durango磷灰石单颗粒年龄的2个标准偏差在6%。

造成(U-Th)/He年龄分散的可能原因有以下几种(Reiners et al., 2006): 1)颗粒的母体同位素U、 Th分布不均一。Boyce等(2005)利用激光剥蚀ICP-MS对不同的Durango磷灰石颗粒的切面进行测定, 发现有些颗粒内部存在明显的U、 Th分带。通常常规的地质样品磷灰石颗粒也存在明显的母体同位素U、 Th分带。这些分带不仅影响颗粒的α 粒子射出效应校正, 而且还影响4He的浓度、 扩散和丢失(Farley et al., 2011)。2)不同颗粒的He扩散性质不同。3)普通双目镜无法识别的晶体内部的包裹体阻碍提取4He之后的U、 Th、 Sm的溶解和完整测量。4)周围富含U、 Th的矿物向磷灰石颗粒注入4He。5)不同的He扩散性或周围环境中不同的He分配系数导致He在磷灰石颗粒中的保存存在差异。6)α 粒子射出效应校正的不确定性导致定年重复性差。除了上述原因之外, 颗粒的粒径大小、 冷却速率、 辐射等因素都会影响(U-Th)/He年龄的分散性(Fitzgerald et al., 2006; Shuster et al., 2006; Flowers et al., 2007, 2009; Flowers, 2009; Brown et al., 2013)。因此, 要得到标准样品和地质样品的准确年龄, 首先需要使用激光剥蚀ICP-MS或扫描电镜评估样品的U、 Th分布。另外, 近年来发展的激光原位微区(U-Th)/He定年方法能够识别颗粒间的U、 Th、 Sm分布不均和定量分析母体同位素分布不均对年龄造成的影响, 是未来发展的趋势。中国地震局地质研究所(U-Th)/He年代学实验室新安装的准分子激光剥蚀系统RESOlution M-50不仅可与ICP-MS进行连接, 还可用真空管与He萃取和测量系统(Alphachron)进行连接, 为未来发展激光微区(U-Th)/He定年方法奠定了基础。

本研究采用ASX-520型自动进样器完成样品U、 Th同位素比值的自动测试工作。首先编写样品自动测试列表, 然后将待测溶液放入相应的样品架孔位上, 最后自动进样器通过进样针按顺序提取待测溶液。为了避免样品间的污染, 进样针在提取样品之后均先用超纯水冲洗至少10s, 接着先后提取稀HNO3和超纯水对管路进行彻底清洗。先后步骤及各步的时间一旦被设定, 自动进样器便根据设置的参数稳定有序地提取样品溶液并完成自身清洗。因此, 自动进样器不仅增加了ICP-MS的测试效率, 还在很大程度上降低了人为误差。然而, 自动进样器也有自身的局限性, 其配备的进样管路较长, 在一定程度上增加了进样溶液的需求量。文中溶解磷灰石中的U和Th之后均会加超纯水至2.5ml, 这个体积对标样Drango年龄的准确测定无影响, 因为Durango Th/U比值较大, 利用高灵敏的ICP-MS7900仍然可以准确进行测定。但是, 对Th/U比值较小或Th、 U含量很低的样品来说, 稀释倍数越多, 在现有灵敏度水平的条件下, 测量精度受到的影响越大。通过条件实验, 我们可将配置的溶液体积缩减至1ml, 从而大大提高元素分析时的灵敏度。

The authors have declared that no competing interests exist.

参考文献
[1] 保增宽, 袁万明, 王世成, . 2005. 磷灰石(U-Th)/He定年技术及应用简介[J]. 岩石矿物学杂志, 24(2): 126132.
BAO Zeng-kuan, YUAN Wan-ming, WANG Shi-cheng, et al. 2005. Apatite(U-Th)/He dating and its application[J]. Acta Petrologica et Mineralogica, 24(2): 126132(in Chinese). [本文引用:1]
[2] 常健, 邱楠生. 2012. 塔里木盆地自然样品磷灰石(U-Th)/He封闭温度及其地质意义[J]. 科学通报, 57(23): 22112219.
CHANG Jian, QIU Nan-sheng. 2012. Closure temperature of(U-Th)/He system in apatite obtained from natural drillhole samples in the Tarim Basin and its geological significance[J]. Chinese Science Bulletin, 57(26): 34823490. [本文引用:1]
[3] 陈文, 张彦. 2010. ( U-Th)/He同位素定年实验室在中国地质调查局系统内建成[J]. 中国地质, 37(3): 840.
CHEN Wen, ZHANG Yan. 2010. Establishement of(U-Th)/He isotope dating laboratory in the China Geological Survey System[J]. Geology in China, 37(3): 840(in Chinese). [本文引用:1]
[4] 蒋毅, 常宏. 2012. 磷灰石(U-Th)/He定年方法综述[J]. 岩石矿物学杂志, 31(5): 757766.
JIANG Yi, CHANG Hong. 2012. Apatite(U-Th)/He dating: A review[J]. Acta Petrologica et Mineralogica, 31(5): 757766(in Chinese). [本文引用:1]
[5] 孙敬博, 孙腾飞, 陈文, . 2015. 新疆东天山红云滩地区构造-热演化探讨: 来自Ar-Ar和(U-Th)/He热年代学的约束[J]. 岩石学报, 31(12): 37323742.
SUN Jing-bo, SUN Teng-fei, CHEN Wen, et al. 2015. Thermo-tectonic evolution history of Hongyuntan area, eastern Tianshan, Xinjiang: Constrained from Ar-Ar and (U-Th)/He dating[J]. Acta Petrologica Sinica, 31(12): 37323742(in Chinese). [本文引用:1]
[6] 王杰, 秦建中, 刘文汇, . 2012. 川东北元坝地区中生代构造与动态热演化史-磷灰石、 锆石(U-Th)/He定年分析[J]. 石油实验地质, 34(1): 1924.
WANG Jie, QIN Jian-zhong, LIU Wen-hui, et al. 2012. Mesozoic tectonics and dynamic thermal history in Yuanba area of northeastern Sichuan Basin: Application of(U-Th)/He dating of apatite and zircon[J]. Petroleum Geology & Experiment, 34(1): 1924(in Chinese). [本文引用:1]
[7] 吴林, 王非, 单竞男, . 2016. 国际标样Durango磷灰石(U-Th)/He年龄测定[J]. 岩石学报, 32(6): 18911900.
WU Lin, WANG Fei, SHAN Jing-nan, et al. 2016. (U-Th)/He dating of international stand ard Durango apatite[J]. Acta Petrologica Sinica, 32(6): 18911900(in Chinese). [本文引用:3]
[8] 张彦, 陈文. 2011a. 4He同位素含量测试技术研究[J]. 地质论评, 57(2): 300304.
ZHANG Yan, CHEN Wen. 2011a. Study on the4He content measurement[J]. Geological Review, 57(2): 300304(in Chinese). [本文引用:1]
[9] 张彦, 陈文, 胡明月. 2011b. 非稀释剂法电感耦合等离子体质谱测定磷灰石中铀钍钐同位素的含量[J]. 岩矿测试, 30(6): 727731.
ZHANG Yan, CHEN Wen, HU Ming-yue. 2011b. Determination of238U, 232Th and 147Sm isotopes in apatites by inductively coupled plasma-mass spectrometry with non-spike method[J]. Rock and Mineral Analysis, 30(6): 727731(in Chinese). [本文引用:1]
[10] 郑德文, 王非, 张培震, . 2000. 磷灰石U-Th/He法: 一种低温热年代计[J]. 地震地质, 22(4): 427435. doi: DOI: 103969/j. issn. 0253-4967. 2000. 04. 012.
ZHENG De-wen, WANG Fei, ZHANG Pei-zhen, et al. 2000. Apatite(U-Th)/He dating method-a low temperature thermochronometer[J]. Seismology and Geology, 22(4): 427435(in Chinese). [本文引用:1]
[11] 郑德文, 武颖, 庞建章, . 2016. U-Th/He热年代学原理、 测试及应用[J]. 第四纪研究, 36(5): 10271036.
ZHENG De-wen, WU Ying, Pang Jian-zhang, et al. 2016. Fundamentals, dating and application of U-Th/He dating thermochronology[J]. Quaternary Sciences, 36(5): 10271036(in Chinese). [本文引用:1]
[12] Boyce J W, Hodges K V. 2005. U and Th zoning in Cerro de Mercado(Durango, Mexico)fluorapatite: Insights regarding the impact of recoil redistribution of radiogenic4He on(U-Th)/He thermochronology[J]. Chemical Geology, 219(1-4): 261274. [本文引用:1]
[13] Brown R W, Beucher R, Roper S, et al. 2013. Natural age dispersion arising from the analysis of broken crystals. Part I: Theoretical basis and implications for the apatite(U-Th)/He thermochronometer[J]. Geochimica et Cosmochimica Acta, 122: 478497. [本文引用:1]
[14] Dai J G, Wang C S, Hourigan J, et al. 2013. Insights into the early Tibetan plateau from(U-Th)/He thermochronology[J]. Journal of the Geological Society, 170(6): 917927. [本文引用:1]
[15] Evans N J, Byrne J P, Keegan J T, et al. 2005. Determination of uranium and thorium in zircon, apatite, and fluorite: Application to laser(U-Th)/He thermochronology[J]. Journal of Analytical Chemistry, 60(12): 11591165. [本文引用:2]
[16] Farley K A. 2000. Helium diffusion from apatite: General behavior as illustrated by Durango fluorapatite[J]. Journal of Geophysical Research, 105(B2): 29032914. [本文引用:2]
[17] Farley K A, Shuster D L, Ketcham R A. 2011. U and Th zonation in apatite observed by laser ablation ICPMS, and implications for the(U-Th)/He system[J]. Geochimica et Cosmochimica Acta, 75(16): 45154530. [本文引用:1]
[18] Farley K A, Stockli D F. 2002. (U-Th)/He dating of phosphates: Apatite, monazite, and xenotime[J]. Reviews in Mineralogy and Geochemistry, 48(1): 559577. [本文引用:3]
[19] Farley K A, Wolf R A, Silver L T. 1996. The effects of long alpha-stopping distances on(U-Th)/He ages[J]. Geochimica et cosmochimica acta, 60(21): 42234229. [本文引用:1]
[20] Fitzgerald P G, Baldwin S L, Webb L E, et al. 2006. Interpretation of(U-Th)/He single grain ages from slowly cooled crustal terranes: A case study from the Transantarctic Mountains of southern Victoria Land [J]. Chemical Geology, 225(1-2): 91120. [本文引用:1]
[21] Flowers R M. 2009. Exploiting radiation damage control on apatite(U-Th)/He dates in cratonic regions[J]. Earth and Planetary Science Letters, 277(1-2): 148155. [本文引用:2]
[22] Flowers R M, Ketcham R A, Shuster D L, et al. 2009. Apatite(U-Th)/He thermochronometry using a radiation damage accumulation and annealing model[J]. Geochimica et Cosmochimica Acta, 73(8): 23472365. [本文引用:1]
[23] Flowers R M, Shuster D L, Wernicke B P, et al. 2007. Radiation damage control on apatite(U-Th)/He dates from the Grand Canyon region, Colorado Plateau[J]. Geology, 35(5): 447450. [本文引用:1]
[24] Foeken J P T, Stuart F M, Dobson K J, et al. 2006. A diode laser system for heating minerals for(U-Th)/He chronometry[J]. Geochemistry, Geophysics, Geosystems, 7(4). doi: DOI:10.1029/2005GC001190. [本文引用:1]
[25] House M A, Farley K A, Stockli D. 2000. Helium chronometry of apatite and titanite using Nd-YAG laser heating[J]. Earth and Planetary Science Letters, 183(3-4): 365368. [本文引用:1]
[26] McDowell F W, McIntosh W C, Farley K A. 2005. A precise 40Ar- 39Ar reference age for the Durango apatite(U-Th)/He and fission-track dating stand ard[J]. Chemical Geology, 214(3-4): 249263. [本文引用:1]
[27] Reiners P W, Nicolescu S. 2006. Measurement of parent nuclides for(U-Th)/He chronometry by solution sector ICP-MS [R]. ARHDL Report 1. Tucson, AZ: University of Arizona. [本文引用:4]
[28] Shuster D L, Flowers R M, Farley K A. 2006. The influence of natural radiation damage on helium diffusion kinetics in apatite[J]. Earth and Planetary Science Letters, 249(3-4): 148161. [本文引用:1]
[29] Tripathy-Lang A, Hodges K V, Monteleone B D, et al. 2013. Laser(U-Th)/He thermochronology of detrital zircons as a tool for studying surface processes in modern catchments[J]. Journal of Geophysical Research: Earth Surface, 118(3): 13331341. [本文引用:1]
[30] Wolf R A, Farley K A, Silver L T. 1996. Helium diffusion and low-temperature thermochronometry of apatite[J]. Geochimica et Cosmochimica Acta, 60(21): 42314240. [本文引用:1]
[31] Zeitler P K, Herczeg A L, McDougall I, et al. 1987. U-Th-He dating of apatite: A potential thermochronometer[J]. Geochimica et Cosmochimica Acta, 51(10): 28652868. [本文引用:2]