LA-MC-ICP-MS方解石U-Pb定年技术 &nbsh1; LA-MC-ICP-MS方解石U-Pb定年技术* 谢博航 吴石头** 杨岳衡 王浩 赵子福 黄超 谢烈文 1. 中国科学院地质与地球物理研究所，岩石圈演化国家重点实验室，北京 100029 2. 中国科学院地球科学创新研究院，北京 100029 3. 中国科学院大学地球与行星科学学院，北京 100049 4. 中国科学院壳幔物质与环境重点实验室，中国科学技术...
谢博航 吴石头** 杨岳衡 王浩 赵子福 黄超 谢烈文
1. 中国科学院地质与地球物理研究所，岩石圈演化国家重点实验室，北京 100029 2. 中国科学院地球科学创新研究院，北京 100029 3. 中国科学院大学地球与行星科学学院，北京 100049 4. 中国科学院壳幔物质与环境重点实验室，中国科学技术大学地球和空间科学学院，合肥 230026 5. 中国科学院比较行星学卓越创新中心，合肥 230026
同位素地质年代学是地球科学研究中最为基础的方向之一，为厘定深时地质过程发生和持续时间提供了定量制约，主要包括U-Pb、Lu-Hf、Sm-Nd、Rb-Sr、Re-Os和Ar-Ar等体系(李献华等, 2022)。近年来激光剥蚀电感耦合等离子体质谱技术(LA-ICP-MS)快速发展，其具有空间分辨率高(微米级)、分析速度快(～2分钟/剥蚀点)和普及率高等优点，已广泛应用于地球科学领域。目前基于LA-ICP-MS分析技术已开发了多种副矿物U-Pb定年方法，包括U含量较高(>50×10-6)的锆石(Yuanetal., 2004; 谢烈文等, 2008)、独居石(Liuetal., 2012)、磷钇矿(刘志超等, 2011)、钙钛矿(Yangetal., 2009)、榍石(孙金凤等, 2012; Maetal., 2019)、磷灰石(Chewetal., 2011)和U含量较低(
方解石U含量低(一般低于5×10-6)，并且普通铅占比高，因此采用LA-ICP-MS进行方解石U-Pb定年更具有挑战性，通常需要高灵敏度的仪器，如扇形磁场(SF)-ICP-MS(Cooganetal., 2016; Robertsetal., 2017; Wuetal., 2022)和多接收(MC)-ICP-MS(Yokoyamaetal., 2018; 程婷等, 2020; Zhangetal., 2021)。目前也有报道采用四级杆(Q)-ICP-MS进行方解石U-Pb定年，但往往需要非常大的激光束斑(>200μm)(Yangetal., 2021)，这使得空间分辨率变差，在进行期次多变的样品分析时受限较大。相比于SF-ICP-MS来说，MC-ICP-MS有一定的优势，如其灵敏度更高，并且同位素同时分析，使得分析精度更优，这使得其在做低U含量样品或高空间分辨率下(～85μm)有一定的优势。程婷等(2020)采用Nu Plasma II MC-ICP-MS建立了方解石U-Pb定年方法，并报道了AHX-1a标样的多次分析数据。Yokoyamaetal. (2018)采用Neptune MC-ICP-MS建立了方解石U-Pb定年方法，作者采用12次方高阻接收206Pb和207Pb信号，11次方正常电阻接收其它信号，显著提高了灵敏度，但法拉第杯具有较高的检测背景(～8000cps或～0.1mv)，这使得不能准确测定低U含量的样品，并且12次方高阻响应时间长，需要对响应时间进行校正。Zhangetal. (2021)进一步表明Neptune Plus MC-ICP-MS可采用混合接收器模式(即法拉第杯接收U和Th，离子计数器接收Pb)进行方解石U-Pb定年，拓宽了可检测的范围，特别是低U含量的样品。
虽然当前LA-MC-ICP-MS方解石U-Pb定年技术已相对成熟，但MC-ICP-MS仪器，如Neptune Plus型号的MC-ICP-MS，可配置不同型号的锥组合(Standard采样锥、Jet采样锥、H截取锥和X截取锥)。不同的锥组合具有特定的优势，Yokoyamaetal. (2018)和Zhangetal. (2021) 证明 住所证明下载场所使用证明下载诊断证明下载住所证明下载爱问住所证明下载爱问 了对于U-Pb定年来说，Jet+X锥要比S+H锥灵敏度更高，但不同锥组合的优势未进行详细的研究。此外，目前LA-MC-ICP-MS方解石U-Pb定年都是在100μm左右的束斑下进行的，LA-MC-ICP-MS在高空间分辨率(～85μm)模式下尚未得到验证。本文系统探究了四种锥组合(S+H、Jet+H、S+X和Jet+X)的灵敏度情况及其应用于方解石U-Pb定年的优势。基于此开发了高空间分辨率下(～85μm)LA-MC-ICP-MS方解石U-Pb定年技术，为期次多变的脉体方解石应用提供有力技术保证。我们采用LA-MC-ICP-MS初步测定了两个实验室内部方解石U-Pb年龄标样(TLM和LSJ07)。方解石U-Pb体系易开放，为了获得可靠的年龄结果，我们讨论了方解石U-Pb体系的封闭性问题，基于理论计算初步探讨了方解石U-Pb定年封闭温度情况。
本研究共用到8个标准物质，分别为2个硅酸盐玻璃标准物质和6个方解石U-Pb年龄标准物质，先前采用多种分析技术系统的表征了NIST SRM 614和ARM-3的元素含量和同位素比值(Jochumetal., 2011; Wuetal., 2019, 2021)。WC-1、Duff Brown Tank、ASH-15D、JT四个方解石年龄标准物质已经进行了ID-TIMS分析；TLM和LSJ07是实验室内部质量监控标准物质。表1汇总了这8个标准物质信息，以下对这8个标准物质进行简单介绍。
表1 本文使用的硅酸盐玻璃Pb-Pb标准物质和方解石U-Pb定年标准物质
NIST SRM 614：该样品为钠钙质人工合成的玻璃标准物质(Jochumetal., 2011)。虽然NIST SRM 614玻璃的基体不同于方解石，但在测试过程中207Pb/206Pb比值未见明显分馏(Guillongetal., 2020)。目前在进行方解石207Pb/206Pb仪器分馏校正时，绝大多实验室采用的是NIST SRM 614。该样品的微量元素含量范围为0.8×10-6～3.0×10-6。
ARM-3：该样品为安山质玻璃标准物质。Wuetal. (2019, 2021)采用多种分析技术对ARM-3玻璃的主量元素、微量元素和多个同位素比值(Li、B、Si、O、Mg、Sr、Nd、Hf和Pb同位素)进行了表征。该样品的微量元素含量为～5×10-6。前期研究表明ARM-3是良好的207Pb/206Pb校正标准物质(Wuetal., 2022)。
WC-1：该样品为1kg左右海相方解石胶结物。WC-1是目前唯一一个已系统表征，且广泛分发于LA-ICP-MS实验室的方解石U-Pb定年标准物质。该方解石胶结物充填于二叠纪(Capitanian阶)礁杂岩内与断层相关的不整合海相岩脉，出露在位于西德克萨斯州特拉华盆地西侧的瓜达卢佩山脉上。Robertsetal. (2017)采用ID-TIMS方法得到WC-1的U-Pb年龄为254.4±6.4Ma(2s)。该样品的U含量为～5×10-6，放射性成因Pb占比为85%～98%。
Duff Brow Tank：该样品为石灰岩手标本样品，其采自于Duff Brown Tank地层。该地层是富黏土粉砂岩中粗粒长石砂岩透镜体与Long Point石灰岩的互层。Hilletal. (2016)采用ID-MC-ICP-MS方法得到Duff Brown Tank的U-Pb年龄为64.04±0.67Ma(2s)，该样品的U含量为5×10-6～20×10-6。
ASH-15D：该样品采自于以色列南部内盖夫沙漠中部的喀斯特洞穴中，Nurieletal. (2021)采用ID-TIMS方法得到ASH-15D的U-Pb年龄为2.965±0.011Ma(2s)(未进行U不平衡校正)。其放射性成因Pb占比低，需采用SF-ICP-MS或MC-ICP-MS测试，灵敏度才能达到 要求 对教师党员的评价套管和固井爆破片与爆破装置仓库管理基本要求三甲医院都需要复审吗 。该样品的U含量为2×10-6～5×10-6。
JT：该样品为方解石脉，其采自于瑞士北部磨拉石盆地中三叠统Mushelkalk群的泥晶灰岩中。JT是方解石脉网的一部分，Guillongetal. (2020)采用ID-TIMS方法得到JT的U-Pb年龄为13.797±0.031Ma(2s)，该样品的U含量为～0.5×10-6。
TLM：该样品为白色方解石晶体，其采自于塔里木盆地TS2井奥陶系鹰山组二段，采样深度为6106.4m，该地层岩性整体为黄灰色泥晶灰岩、含沥青质泥晶灰岩及沥青质泥晶灰岩。鹰山组二段中上部晶洞充填了白色巨晶方解石。TLM的U含量是在0.1×10-6～0.5×10-6之间。
LSJ07：该样品从网络经销商购买得到。这是个与辉锑矿伴生方解石晶体，质量约15g。购买信息显示其来自于湖南冷水江辉锑矿山，但并没有详细的采集信息。LSJ07的U含量是在1×10-6～16×10-6之间。
LA-MC-ICP-MS原位方解石U-Pb定年分析在中国科学院地质与地球物理研究所进行，采用的仪器设备是Neptune Plus MC-ICP-MS和PhotoMachine Analyte G2型号193nm准分子激光剥蚀系统。仪器的基本情况参考Xieetal. (2017)。
将方解石样品颗粒整齐排列在双面胶带上，放上PVC环，将环氧树脂和固化剂混合均匀后注入PVC环，待树脂固化后对其进行抛光至样品表面露出。制备好的环氧树脂靶放置于HelEx-Ⅱ激光剥蚀池进行测试，所采用的仪器参数和杯结构见表2。对于实际样品，可直接采用光学薄片进行分析。
表2 LA-MC-ICP-MS仪器参数和杯结构信息
Neptune Plus在低分辨模式下同时接收202Hg、204Pb、206Pb、207Pb、208Pb、232Th和238U离子信号强度，中心杯虚拟质量数设置为223.95，H3和H4分别接收232Th和238U，IC4、IC5、IC3A、IC2、IC1B1分别接收202Hg、204Pb、206Pb、207Pb、208Pb。IC3、IC2、IC1B1是传统的二次电子倍增器(SEM)，IC4、IC5是小型电子倍增器(CDD)，离子流通过导向装置被导入SEM或CDD。图1为溶液模式和激光模式下的套峰图。
图1 溶液模式和激光模式下的U-Pb定年的套峰图Fig.1 Peak plots of U-Pb dating in solution mode and laser mode
IC1B1、IC2、IC3、IC4、IC5的暗噪声和工作电压是采用仪器操作软件中的PCL程序确定的，单个IC的暗噪声均不超过0.013cps，这5个IC的工作电压均在2500V到2800V之间。各IC产率是在溶液进样模式下通过动态跳峰方法测定得出。使这5个IC和中心法拉第杯依次接收约12万cps的信号，积分时间为33.554s，采集5个cycle。最后根据其信号强度的平均值与中心法拉第杯信号强度的平均值标准化得到(1mv=62500cps)每个IC的产率。5个IC的产率在80%到90%之间。
激光剥蚀时采用的束斑直径为20～85μm(其中ARM-3采用35μm)，剥蚀频率为10Hz，单点分析包括5s的背景信号采集、23s的剥蚀信号采集和5s的清洗。Duff Brown Tank、ASH-15D、JT、TLM、LSJ07作为未知样品。对于未知样品，每间隔15个分析点插入两个NIST SRM 614、两个ARM-3和3个方解石标准物质WC-1，每个未知样品进行30～50次点分析。
数据处理采取了Robertsetal. (2017)和Wuetal. (2022)描述的两步校正法。首先将Neptune Plus的原始数据导入到Iolite 3.7数据处理软件中，NIST SRM 614作为标准物质，采用X_U-Pb_Geochron4结合VisualAge数据处理策略(Patonetal., 2010; Petrus and Kamber, 2012)计算得到207Pb/206Pb和206Pb/238U的比值及误差。基于NIST SRM 614分馏曲线对WC-1和其他样品的206Pb/238U分馏进行校正(Patonetal., 2010)。然后在Microsoft Excel中，基于WC-1标准物质校正计算出238U/206Pb的偏差系数，并用该系数校正NIST SRM 614与方解石之间的基体效应(由元素分馏程度不同造成的偏差)，校正后的数据在Isoplot软件中作图(Ludwig, 2003)。误差传递参照于Horstwoodetal. (2016)所描述的模型。
图2 四种锥组合(S+H、Jet+H、S+X和Jet+X)在引入不同量氮气情况下的238U和206Pb信号强度和Th+/U+比值Fig.2 238U and 206Pb signal intensities and Th+/U+ ratios of four cone combinations (S+H, Jet+H, S+X and Jet+X) with different amounts of nitrogen
方解石U含量低，通常需要采用高灵敏度仪器进行测试分析。热电公司Neptune Plus型号的MC-ICP-MS可配置不同的锥组(如Jet+X锥)。前人研究表明对于U-Pb定年来说，Jet+X锥要比S+H锥灵敏度更高(Yokoyamaetal., 2018; Zhangetal., 2021)，但并未系统探究其它锥组的灵敏度情况(如S+X和Jet+H)。本文系统探究了四种锥组合(S+H、Jet+H、S+X和Jet+X)在引入不同量氮气情况下的U和Pb灵敏度，结果如图2所示。Jet采样锥的孔径要大于S采样锥；X截取锥和H截取锥的孔径虽然都是0.8mm，但是H锥入口为圆柱形，出口为喇叭形，而X锥的出入口均为喇叭形。因此与其他锥组合相比，Jet+X锥能为离子束传输提供更大空间，提高了灵敏度。
从图2中可以看出，不同锥组的U和Pb的灵敏度具有显著差别。在不引入氮气的情况下，对于Pb元素来说，不同锥组的灵敏度为Jet+X>S+X≈Jet+H>S+H，对于U元素来说，Jet+X≈S+X≈Jet+H>S+H。四种锥组下，氮气的引入均使得最高灵敏度对应的载气流速降低，并且使得峰的平顶区明显变宽。对于S+H锥和Jet+H锥来说，氮气的引入并没有使得绝对灵敏度提高，而对于S+X锥和Jet+X锥来说，氮气的加入对U和Pb的最高灵敏度有影响。氮气的引入使得U的最高信号明显增强，但一定程度的抑制了Pb的最高灵敏度，这其中的原因尚不清楚，可能与元素离子化效率有关。Heetal. (2018)和Wuetal. (2020, 2022)基于Element XR SF-ICP-MS仪器表明氮气对信号的增敏可能与X锥有关，本研究结果进一步表明这种现象同样适用于Neptune Plus型号的MC-ICP-MS仪器。
本研究结果表明，不同锥组结合氮气增敏效果对不同的元素具有不同的效果，在采用这些锥组合时，应该先进行条件实验。对于U-Pb定年来说，综合各项参数，Jet+X锥结合4mL·min-1氮气时条件最优。图3绘制了不同锥组的U和Pb的信号强度。Element XR SF ICP-MS的灵敏度也绘制于图3以做对比(数据引自Wuetal., 2022)。从图3中可以看出在仪器最佳灵敏度模式下(Jet+X锥结合氮气增敏，85μm束斑)，信号提高了3～5倍，灵敏度达到：238U>55万cps/10-6，206Pb>18万cps/10-6。这与前人采用Nu Plasma Ⅱ MC-ICP-MS在100μm束斑下的灵敏度相近(U>50万cps/10-6，Pb>40万cps/10-6，程婷等, 2020)。虽然这个灵敏度下Th+/U+比值仅有0.4左右(可能与氧化物产率升高有关)，但我们发现这种条件并不影响方解石U-Pb定年精确度和准确度(见讨论2.2部分)。
Duff Brown Tank：该样品进行了两次分析，一次采用S+H锥结合4mL·min-1氮气，一次为Jet+X锥结合4mL·min-1氮气。图4绘制了Tera-Wasserburg图，上交点固定于0.738(Hilletal., 2016)。该样品的普通铅占比为18%～40%，238U/206Pb单点分析的精度优于3%(2SE)，207Pb/206Pb单点分析的精度优于2%(2SE)(见电子版附表1)。固定上交点得到的年龄结果分别为66.4±0.8Ma(2s，n=24/24)和66.0±0.4Ma(2s，n=50/50)，略高于之前报道的ID-TIMS结果(64.04±0.67Ma，2s)。这个结果与我们之前采用Element XR测出的数据(66.2±0.5Ma，2s，Wuetal., 2022)一致。考虑到WC-1的ID-TIMS定值误差(～2.5%)，目前LA-ICP-MS方解石U-Pb定年不确定度为3.0%～3.5%，因此，我们所得的Duff Brown Tank结果与ID-TIMS结果一致。
图4 Duff Brown Tank的LA-MC-ICP-MS数据Tera-Wasserburg图上交点固定为0.738. 共进行了两次实验，分别采用S+H锥和引入少量N2以及Jet+X锥和引入少量N2的仪器条件. 图中的n=a/b，其中“a”表示剔除异常值的数据总数，“b”表示分析总次数. 由Isoplot软件进行画图(Ludwig, 2003)Fig.4 Tera-Wasserburg diagrams of LA-MC-ICP-MS data of Duff Brown TankThe upper intersection is anchored to 0.738. Two experiments were carried out, using the instrument conditions of S+H cone with a small amount of N2, Jet+X cone with a small amount of N2. n=a/b in the figure, where “a” represents the total number of data excluding outliers, and “b” represents the total number of analyses. These data were plotted by isoplot software (Ludwig, 2003)
图5 ASH-15D的LA-MC-ICP-MS数据Tera-Wasserburg图上交点固定为0.832. 共进行了两次实验，均采用Jet+X锥和引入少量N2的仪器条件. 图中的n=a/b，其中“a”表示剔除异常值的数据总数，“b”表示分析总次数. 由Isoplot软件进行画图(Ludwig, 2003)Fig.5 Tera-Wasserburg diagram of LA-MC-ICP-MS data of ASH-15DThe upper intersection is anchored to 0.832. Two experiments were carried out, both using Jet+X cone with a small amount of N2. n=a/b in the figure, where “a” represents the total number of data excluding outliers, and “b” represents the total number of analyses. These data were plotted by isoplot software (Ludwig, 2003)
图6 JT的LA-MC-ICP-MS数据Tera-Wasserburg图上交点固定为0.839. 共进行了三次实验，第一次采用S+H锥和引入少量N2仪器条件，第二次和第三次采用Jet+X锥和引入少量N2的仪器条件. 图中的n=a/b，其中“a”表示剔除异常值的数据总数，“b”表示分析总次数. 由Isoplot软件进行画图(Ludwig, 2003)Fig.6 Tera-Wasserburg diagrams of LA-MC-ICP-MS data of JTThe upper intersection is anchored to 0.839. A total of three experiments were conducted. The first experiment, the S+H cone with a small amount of N2 were used, and the last two experiments, the Jet+X cone with a small amount of N2 were used. n=a/b in the figure, where “a” represents the total number of data excluding outliers, and “b” represents the total number of analyses. These data were plotted by isoplot software (Ludwig, 2003)
ASH-15D：该样品进行了两次分析，均采用Jet+X锥结合4mL·min-1氮气仪器参数。图5绘制了Tera-Wasserburg图，上交点固定于0.832(Nurieletal., 2021)。结果显示该样品的普通铅占比变化非常大(5%～90%)，238U/206Pb的单点分析的精度为4%～10%(2SE)，207Pb/206Pb单点分析的精度为2%～20%(2SE)(见附表2)。两次分析固定上交点(207Pb/206Pb0=0.832，Nurieletal., 2021)后的下交点年龄分别为：2.94±0.07Ma(2s，n=30/30)和3.03±0.08Ma(2s，n=30/30)。考虑到LA-ICP-MS方解石U-Pb定年不确定度为3.0%～3.5
浙公网安备 33021202002483