矿物岩石地球化学通报  2016, Vol. 35 Issue (3): 441-447   PDF    
高Cr地质样品的Mg同位素分析方法
李瑞瑛, 柯珊 , 何永胜, 高庭, 孟旭男    
中国地质大学(北京), 地质过程与矿产资源国家重点实验室, 地球科学与资源学院, 北京 100083
摘要: 使用多接收等离子体质谱仪测定Mg同位素比值时,Cr元素会干扰测定结果。高Cr地质样品如铬尖晶石和铬铁矿中的Cr/Mg值可高达5:1。针对这些样品,本研究检测了Cr对Mg同位素测试的干扰程度,评估了实验室流程中Cr和Mg的分离效果,探讨了Cr在淋洗过程中的行为。实验采用1 mol/L HNO3及2.3 mL AG50W-X8型阳离子树脂,仅通过一次分离纯化便能将Mg与Cr等其他基质元素有效分离,且Mg的回收率可达99.8%,测试结果可靠。此外,Cr在分离过程中分为两阶段被洗脱,这与其在样品中不同的赋存价态有关。
关键词: Mg同位素     MC-ICP-MS     高Cr样品     基质元素    
High Precision Magnesium Isotope Measurement for High-Cr Samples
LI Rui-ying, KE Shan, HE Yong-sheng, GAO Ting, MENG Xu-nan    
China University of Geosciences(Beijing), State Key Lab of Geological Processes and Mineral Resources, School of Earth Sciences and Resources, Beijing 100083, China
Abstract: High chromium contents of samples will significantly interfere results ofthe Mg isotopic analysis using a multi-collector inductively coupled plasma mass spectrometry(MC-ICP-MS). There are high-Cr samples, such as chrome spinel and chromite with Cr/Mg ratios up to 5:1. This study investigates the intensity of interference of Cr on Mg isotopic analyses, evaluates the separation efficiency of Cr from Mg in laboratory process, and discusses the Cr behavior in eluting process. Experiment results show that Mg can be purified efficiently from Cr and other matrix elements through cluting process once with AG50W-X8 cation resin in 1 mol/L HNO3 medium. The recovery of yielded Mg is higher than 99.8%, ensuring the accuracy and precision of Mg isotope measurement. In addition, the behavior of Cr during elution process is also discussed. During the separation process, chromium was eluted from the column in early and late two stages respectively, owing to different valence states of Cr in solutions.
Key words: Mg isotopic analyses     MC-ICP-MS     high-Cr samples     matrix elements    

Mg作为主要的造岩元素之一,广泛分布在地球的各个圈层中,参与了大部分地球化学过程。它有3个稳定同位素:24Mg(78.99%),25Mg(10.00%)和26Mg(11.01%),其中24Mg和26Mg相对质量差高达8%,使得Mg同位素在自然样品中易于产生显著的分馏。因此,利用Mg同位素可以示踪物质源区和地质过程。近年来,随着多接收电感耦合等离子体质谱仪(MC-ICP-MS)技术的发展,Mg同位素的测试精度显著提高,δ26Mg的分析精度可达0.05‰(2SD)(An et al.,2014),使Mg同位素在地学领域的广泛应用成为可能。

最近的研究表明,高Cr样品的Mg同位素特征在地质温度计和研究铬铁矿成因方面可能具有良好的应用前景。根据地幔橄榄岩中尖晶石-橄榄石间的Mg同位素分馏系数,Young等(2009)Liu等(2011)提出Mg同位素有望成为一种新的适用于地幔橄榄岩的地质温度计。此外,Xiao等(2016)尝试利用西藏罗布莎铬铁矿的Mg-Fe同位素特征探讨铬铁矿的成因,发现橄榄石和铬铁矿之间发生了Mg-Fe同位素交换,且不同产状的铬铁矿同位素交换程度不同。这些研究中涉及到的尖晶石和铬铁矿等都是极富Cr的矿物,其Cr/Mg可高达5 ︰ 1。理论上,在MC-ICP-MS测试过程中,Cr会产生50 Cr2+52 Cr2+等双电荷电子,影响25Mg和26Mg的测试结果。何学贤等(2008)研究发现,当Cr/Mg浓度比超过0.05时,δxMg值(x=25,26)会受到Cr基体效应的影响而发生不规律的变化,严重偏离其真值。因此,为保证数据的准确性,高Cr样品在测试前必须对Mg和Cr进行有效的化学分离。

为将元素Cr和Mg分离,前人针对高Cr样品主要采用了多次化学纯化的办法。如Young等(2009)先后采用AG50W-X12阳离子树脂和AG1-X8阴离子树脂,以及1 mol/L HNO3、不同浓度的HCl和丙酮,通过三柱法提纯了San Carlos富铬尖晶石(Cr/Mg=2 ︰ 1)中的Mg。Liu等(2011)利用1 mL AG50W-X8阳离子树脂和1 mol/L HNO3淋洗酸经过4次或更多次化学分离,将尖晶石(Cr/Mg=1.5 ︰ 1)中的Mg进行了纯化。但是,这些方法均未评估Cr与Mg的分离效果,实际数据的准确性未知,并且化学流程较为复杂和耗时。

本文通过基质元素添加实验,测试了元素Cr对Mg同位素的干扰效应。参考前人实验方法,评估了Mg与Cr及其他元素的分离效果,并探讨了元素Cr在淋洗过程中的化学行为。

1  实验材料和质谱分析 1.1  实验材料

本次实验在中国地质大学(北京)同位素地球化学千级净化实验室内完成。实验所使用的溶液由GSB Mg、Na、Al、K、Ca、Ti、Cr、Mn和Fe单元素标准溶液(北京钢铁研究总院,1000 μg/g)配制,地质标样为美国地质调查局(USGS)的玄武岩标样BCR-2和BHVO-2。实验使用树脂为AG50W-X8阳离子交换树脂(Bio-Rad,200~400目),离子交换柱和器皿均为聚四氟乙烯(PFA)材质。

实验所用的HNO3、HF、HCl是由BV-Ⅲ级原酸(北京化学试剂研究所)经过1次亚沸蒸馏(Savillex DST-1000)获得的超纯试剂,本底分别为82 pg/g、258 pg/g、164 pg/g。H2O为经Millipore系统纯化的超纯水,电阻率为18.2 MΩ。

1.2  质谱分析

样品的Cr各元素浓度和Mg同位素组成的测试,均在本实验室Neptune Plus型多接收器电感耦合等离子体质谱仪完成,并采用低分辨模式。样品和标样都以3% HNO3为介质,以高纯氩气作为进样载气,通过自动进样器(Cetac ASX-110FR)湿法进样。

淋洗曲线中各元素浓度使用中心杯动态跳扫24Mg、23Na、27Al、39K、44 Ca、48Ti、52 Cr、55Mn、56Fe的信号强度获得。每次测量采集5次跳扫数据点。浓度计算时,仪器灵敏度由一已知各元素浓度的混合标准溶液的信号获得,混合标液的组成为:Mg、Na、Al、K、Cr、Mn各0.4 μg/g,Fe为0.8 μg/g,Ca和Ti为1.6 μg/g。

Mg同位素比值采用样品-标样间插法进行测试。样品和标样都使用同批同瓶3% HNO3进行配制,且Mg元素浓度差严格控制在10%以内,以避免酸度和浓度不匹配带来的影响,保证测试的精准度(An et al.,2014Teng and Yang,2014)。本实验中所用人工模拟溶液由GSB Mg标准溶液配制,测试时将GSB Mg原溶液作为报道和间插标样,则模拟溶液的δxMg(x=25,26)理论值为0。在样品和标样的测试之间,用3% HNO3清洗120 s,样品提取时间为90 s。为获得更好的外部重现性,每个样品-标样序列重复测量4次,每次采集30个数据点(cycle),每点的积分时间为4.194 s。此外,在样品-标样测试序列的首尾都会测试一个纯3% HNO3,作为仪器的在线背景值进行扣除。实验室的长期外精度为0.06‰(2SD),实验的全流程空白低于10 ng,平均在7ng左右。

2  基质元素Cr对Mg同位素测试的干扰效应

在Mg同位素分析测试过程中,基质元素的存在会严重影响Mg同位素的测试(Galy et al.,2001Chang et al.,2003何学贤等,2008)。一方面,因为MC-ICP-MS采用的是质谱原理,当基质离子与Mg离子(24Mg+25Mg+26Mg+)的质荷比相同时,便会影响25Mg/24Mg、26Mg/24Mg值的测试,造成同质异位素干扰。另一方面,基质元素(如Na、Al、Fe等)的存在会降低Mg的离子化效率,导致不可预期的质量分馏,产生基体效应。

何学贤等(2008)研究发现,在Nu Plasma HR型质谱仪和膜去溶进样条件下,当样品中的Cr/Mg(浓度)小于0.05时,δ25Mg、δ26Mg值约为0。随着Cr浓度的增加,δ25Mg、δ26Mg值出现波动,相对于理论值或重或轻,不规律地协同变化,并且δ25Mg和δ26Mg值的线性拟合结果(δ25Mg=0.53δ26Mg)服从质量分馏线的斜率。因此,认为是基体效应的影响。

本次实验用GSB Mg、Cr单元素标准溶液,配置了Cr/Mg浓度比不同的若干样品,直接测试Mg同位素的比值。根据实验结果,当Cr/Mg浓度比不超过0.5时,δ25Mg、δ26Mg值在误差范围内与理论值一致(图 1a),未见显著的同质异位素干扰和基体效应影响。随着Cr/Mg浓度比增加,δ26Mg值明显变大,δ25Mg值却逐渐变小;当Cr/Mg浓度比达到5 ︰ 1时,δ26Mg值可达0.32‰,δ25Mg值低至-0.09‰(图 1b)。理论上,随Cr浓度增大,电离产生的50 Cr2+52 Cr2+同质异位素会越多,δ25Mg和δ26Mg的测量值都越重。但实验观察与理论推测不符,说明Cr元素对Mg同位素测试的影响是同质异位素干扰和基体效应的综合体现,只是这两种效应对25Mg和26Mg的影响程度不同。原因在于50 Cr和52 Cr丰度相差很大,50 Cr相对丰度很低(4.35%),产生的50 Cr2+较少,对25Mg的同质异位素干扰较小,而基体效应更加显著;52 Cr丰度高(83.79%),会产生大量的52 Cr2+,使得26Mg受到的同质异位素干扰更显著,基体效应的表现相对较弱。总之,当Cr/Mg>0.5时,Mg同位素值会同时受到同质异位素和基体效应的影响,同质异位素的干扰使其变重,基体效应的影响使其变轻。

灰色区域代表实验室长期外精度(±0.06‰) 图 1 Cr含量对Mg同位素测定值的影响 Fig. 1 Effects of Cr concentrations on δxMg values
3  化学分离流程 3.1  样品准备

本次实验用GSB Mg和其他单元素标准溶液配制了2种人工混合溶液A和B,具体的元素组成比例参见 表 1。 为了便于同时检测其他主要基质元素与Mg的分离效果,溶液A不仅Cr含量很高,还加入了硅酸岩样品中的常见主量元素。 溶液B为模拟的天然铬铁矿,其元素组成参考Zhou等(2014)罗布莎蛇绿岩套铬铁矿样品L16(TiO2=0 . 09%,Al2O3=4 . 30%,Cr2O3=64 . 2%,FeO=23 . 8%,MgO=6 . 21%,均为质量分数)配制。 由于各单元素标准溶液的保存介质不同,因此进行了转换介质。 将混合溶液于100~130℃蒸干后,加入5~10滴浓HNO3,溶解后再蒸干,此过程重复2遍,以确保完全转换为化学分离所需要的HNO3介质。 然后将样品溶解在1 mol/L HNO3中,待分离提纯。

表 1 混合溶液A、B元素含量 Table 1 Compositions of mixing solutions A and B
3.2  淋洗流程

实验采用的离子交换柱内径为4.0 mm,高度为20.5 cm,树脂量为2.3 mL,淋洗介质为1 mol/L HNO3和6 mol/L HCl。

树脂在使用前分别用6 mol/L HCl、1 mol/L HNO3与Milli-Q H2O交替清洗若干次,并封存在Milli-Q H2O中。洗干净的树脂装柱后,加入6 mL 1 mol/L HNO3平衡树脂,取100 μL混合溶液上柱,溶液A第1次淋洗过程详见 表 2。淋洗结果表明Ca、Cr、Al、Fe等元素在54 mL 1 mol/L HNO3中不能被完全洗脱(图 2a)。为此,基于第一次淋洗方案,第2次淋洗过程延长了HNO3淋洗体积,并在淋洗的后期使用了浓度高的6 mol/L HCl(表 2)。淋洗结果如 图 2b,残余离子在后期HCl中全部被洗脱,为树脂的重复使用奠定了基础。

表 2 混合溶液A淋洗流程 Table 2 Eluting processes of the mixing solution A

图 2 溶液A淋洗曲线 Fig. 2 Eluting profiles of the solution A
4  结果与讨论 4.1  淋洗曲线测定结果及分离流程确立

由淋洗结果(图 2),在1 mol/L HNO3中,少量的Cr(~10%)在0~4 mL被洗脱,Na、K、Ti分别在0~12 mL、12~18 mL、12~20 mL之间被洗脱,Mg的洗脱期间为24~36 mL,Mn在32~54 mL被洗脱,Ca在40 mL后被HNO3洗出;样品中的Al、Fe和大部分Cr在6 mol/L HCl中才能被洗脱。 其中约6%的Mn在32~36 mL被洗脱,与Mg的淋洗区间相交,这对于Mn/Mg浓度比为1 ︰ 1的样品,不会影响Mg同位素的测试结果(见本文5.2)。 且大多数岩石样品Mn的含量比Mg低很多(Mn/Mg<0.1)。 而高Mn的样品,或许可以通过2遍甚至更多遍的分离纯化,除去Mn元素。

为了尽可能减少基质元素的影响,并保证样品中Mg的高回收率,根据溶液A的2次淋洗实验结果,最终确定了如下分离纯化流程(表 3)。

表 3 离子交换淋洗流程 Table 3 Mg purification scheme
4.2  Mg同位素回收率及测试结果

前人研究表明,Mg在经过离子交换色谱柱时会产生明显的分馏(Chang et al.,2003Teng et al.,2007Huang et al.,2011),其分馏值可高达2.5‰(Chang et al.,2003),因此,必须尽可能将Mg完全回收才能保证测试值可靠。本文使用以下方法计算了该分离流程的回收率:

回收率=(Mg区间接样量 / 全流程Mg接样量)×100%

回收率=(Mg区间接样量 / Mg上样量)×100%

经计算,溶液A的回收率大于99.8%。模拟的铬铁矿溶液B按 表 3中淋洗流程进行了化学分离,并计算得Mg的回收率为99.8%,能够确保该分离流程不会导致Mg同位素发生分馏。

另外,前人研究发现当Mg的进样量不同时,会导致淋洗曲线漂移(An et al.,2014)。为进一步验证本流程分离效果是否受Mg上样量的影响,分别将Mg含量为20 μg和10 μg 的GSB Mg单元素标准溶液进行了提纯和测试,并与溶液A和溶液B的测试结果(图 3)对比,结果表明,GSB Mg单元素溶液和多元素混合溶液(溶液A和B)的测试结果均在误差范围内与理论值吻合。如,上样量为20 μg和10 μg的GSB Mg单元素溶液测得的δ26Mg值分别为0.01‰±0.08‰和-0.04‰±0.04‰,与理论值吻合;混合溶液A(Mg含量为20 μg)和B(Mg含量为10 μg)的δ26Mg值分别是0.03‰±0.02‰和-0.04‰±0.07‰,在允许的误差范围内也与理论值一致(图 3)。

灰色区域为实验室长期外精度±0.06‰;红色圆圈代表未经过分离流程的GSB Mg原溶液;方形为经过分离流程的GSB Mg单元素溶液,蓝色和橙色分别代 表 10μg和20 μg的上样量;绿色三角为B溶液, 上样量为10 μg,紫色三角为A溶液,上样量为20 μg 图 3 同进样量,不同溶液的δ26Mg值 Fig. 3 δ26Mg values of different mixing solutions and sampling amounts

为进一步验证流程的可靠性,本次实验选取了地质标样进行基质添加前后的测试对比。 在玄武岩标样BCR-2和BHVO-2中加入GSB Cr,使Cr/Mg浓度比为5 ︰ 1。 经化学分离和测试,BCR-2的δ26Mg值为-0.16‰±0.02‰,BHVO-2的δ26Mg值为-0.26‰±0.04‰,分别与本实验中未添加GSB Cr的结果,以及已经发表的BCR-2、BHVO-2的δ26Mg值相吻合(表 4图 4)。

表 4 玄武岩标样BCR-2和BHVO-2的δ26Mg值 Table 4 δ26Mg values of BCR-2 and BHVO-2

红色实心圆圈代表添加了GSB Cr的地质标样δ26Mg值;黑色实心圆圈代表未添加GSB Cr的地质标样结果;灰色实心圆圈为已发表数据 图 4 玄武岩标样BCR-2、BHVO-2的δ26Mg值 Fig. 4 δ26Mg values of reference materials BCR-2 and BHVO-2
4.3  Cr的淋洗行为

在混合溶液淋洗过程中,一个有趣的现象是Cr在整个淋洗过程中分2个阶段被淋洗出来:约10%的Cr在前4 mL就被HNO3洗脱,而90%在后期用较浓的HCl才可洗脱。 这一行为应与Cr的价态密切相关。 Cr是多价态元素,在水溶液中主要以3价和6价形式存在,其中Cr(Ⅲ)呈绿色,在溶液中多以Cr3+,CrOH2+,Cr(OH)+2,Cr(OH)3等阳离子的形式存在;Cr(Ⅵ)呈黄色,以H2CrO4,HCrO-4,CrO2-4,HCr2O-7,Cr2O2-7等阴离子形式为主。 GSB Cr单元素标准溶液本身呈绿色,说明主要以Cr(Ⅲ)存在。 淋洗时被吸附在树脂中的样品也呈暗绿色,表明样品中Cr的主要赋存形态仍为Cr(Ⅲ)阳离子,并在淋洗初期被树脂吸附,在后期6 mol/L HCl中被脱。 由此,可以推测在前期洗脱的少量Cr为非阳离子形式的Cr(Ⅵ)。 那么,从Cr的原始标准溶液到混合溶液的分离过程中,应有少量的Cr(Ⅲ)被氧化至Cr(Ⅵ)。

在强酸环境下(pH=0),Cr3+是Cr的主要存在形式,Cr6+几乎不存在;但随着温度的升高,Cr6+的赋存几率逐渐增大,并以H2CrO4的形式出现(图 5)(Beverskog and Puigdomenech,1997)。 因此,在介质转换过程中,由于高温和浓硝酸的作用,少量的Cr(Ⅲ)可能会被氧化至Cr(Ⅵ)赋存在样品溶液中。 除了高温和浓硝酸的氧化作用外,还可能还与溶液中的其他变价离子的氧化性有关,如Mn(Ⅳ)可作为催化剂将Cr(Ⅲ)氧化至Cr(Ⅵ)。

图 5 20℃、100℃和200℃时,不同pH-Eh条件下Cr在溶液中的赋存状态(据Beverskog and Puigdomenech,1997) Fig. 5 h-pH diagrams of various valence states of Cr in aqueous solutions at 25℃,100℃ and 200℃(modified from Beverskog and Puigdomenech,1997)
5  结论

(1)Cr元素的存在会对Mg同位素的测试同时造成同质异位素干扰和基体效应影响。当Cr/Mg浓度比小于0.5时,δ25Mg、δ26Mg值在误差范围内与理论值保持一致,Cr的存在不会影响Mg同位素比值的测试;而当Cr/Mg浓度比超过0.5时,δ26Mg值明显变重,δ25Mg值却逐渐变轻。

(2)本实验流程利用2.3 mL AG50W-X8型阳离子树脂和1 mol/L HNO3的淋洗介质,仅通过一次化学分离便能够将高Cr样品中的Mg与Cr等基质元素有效分离,同时保证高的Mg回收率和准确的测试结果。

(3)在样品处理和分离过程中,由于元素Cr赋存价态不同,使得Cr以阳离子和阴离子形式分为2个阶段洗脱。

致谢: 研究过程中得到李曙光老师和朱建明老师的指导,在此致以诚挚的感谢。

参考文献
[1] An Y J, Wu F, Xiang Y X, Nan X Y, Yu X, Yang J H, Yu H M, Xie L W, Huang F. 2014. High-precision Mg isotope analyses of low-Mg rocks by MC-ICP-MS. Chemical Geology, 390: 9-21
[2] Beverskog B, Puigdomenech I. 1997. Revised pourbaix diagrams for chromium at 25-300℃. Corrosion Science, 39(1): 43-57
[3] Bizzarro M, Baker J A, Haack H, Lundgaard K L. 2005. Rapid timescales for accretion and melting of differentiated planetesimals inferred from 26Al-26Mg chronometry. The Astrophysical JournalLetters, 632(1): L41-L44
[4] Bizzarro M, Paton C, Larsen K, Schiller M, Trinquier A, Ulfbeck D. 2011. High-precision Mg-isotope measurements of terrestrial and extraterrestrial material by HR-MC-ICPMS: Implications for the relative and absolute Mg-isotope composition of the bulk silicate Earth. Journal of Analytical Atomic Spectrometry, 26(3): 565-577
[5] Chang V T C, Makishima A, Belshaw N S, O'Nions R K. 2003. Purification of Mg from low-Mg biogenic carbonates for isotope ratio determination using multiple collector ICP-MS. Journal of Analytical Atomic Spectrometry, 18(4): 296-301
[6] Galy A, Belshaw N S, Halicz L, O'Nions R K. 2001. High-precision measurement of magnesium isotopes by multiple-collector inductively coupled plasma mass spectrometry. International Journal of Mass Spectrometry, 208(1): 89-98
[7] Huang F, Zhang Z F, Lundstrom C C, Zhi X C. 2011. Iron and magnesium isotopic compositions of peridotite xenoliths from Eastern China. Geochimica et Cosmochimica Acta, 75(12): 3318-3334
[8] Liu S A, Teng F Z, Yang W, Wu F Y. 2011. High-temperature inter-mineral magnesium isotope fractionation in mantle xenoliths from the North China craton. Earth and Planetary Science Letters, 308(1): 131-140
[9] Pogge Von Strandmann P A E, Burton K W, James R H, Van Calsteren P, Gislason S R, Sigfússon B. 2008. The influence of weathering processes on riverine magnesium isotopes in a basaltic terrain. Earth and Planetary Science Letters, 276(1): 187-197
[10] Pogge Von Strandmann P A E, Elliott T, Marschall H R, Coath C, Lai Y J, Jeffcoate A B, Ionov D A. 2011. Variations of Li and Mg isotope ratios in bulk chondrites and mantle xenoliths. Geochimica et Cosmochimica Acta, 75(18): 5247-5268
[11] Teng F Z, Wadhwa M, Helz R T. 2007. Investigation of magnesium isotope fractionation during basalt differentiation: Implications for a chondritic composition of the terrestrial mantle. Earth and Planetary Science Letters, 261(1): 84-92
[12] Teng F Z, Yang W. 2014. Comparison of factors affecting the accuracy of high-precision magnesium isotope analysis by multi-collector inductively coupled plasma mass spectrometry. Rapid Communications in Mass Spectrometry, 28(1): 19-24
[13] Tipper E T, Louvat P, Capmas F, Galy A, Gaillardet J. 2008. Accuracy of stable Mg and Ca isotope data obtained by MC-ICP-MS using the standard addition method. Chemical Geology, 257(1): 65-75
[14] Wombacker F, Eisenhauer A, Heuser A, Weyer S. 2009. Combined chemical separation of Mg, Ca and Fe from geological reference samples for stable isotope determination by MC-ICPMS and double-spike TIMS. Geophysical Research Abstracts, 8: 06219
[15] Xiao Y, Teng F Z, Su B X, Hu Y, Zhou M F, Zhu B, Shi R D, Huang Q S, Gong X H, He Y S. 2016. Iron and magnesium isotopic constraints on the origin of chemical heterogeneity in podiform chromitite from the Luobusa ophiolite, Tibet. Geochemistry, Geophysics, Geosystems, 17(3): 940-953. doi: 10.1002/2015GC006223
[16] Young E D, Tonui E, Manning C E, Schauble E, Macris C A. 2009. Spinel-olivine magnesium isotope thermometry in the mantle and implications for the Mg isotopic composition of Earth. Earth and Planetary Science Letters, 288(3): 524-533
[17] Zhou M F, Robinson P T, Su B X, Gao J F, Li J W, Yang J S, Malpas J. 2014. Compositions of chromite, associated minerals, and parental magmas of podiform chromite deposits: The role of slab contamination of asthenospheric melts in suprasubduction zone environments. Gondwana Research, 26(1): 262-283
[18] 何学贤, 朱祥坤, 李世珍, 唐索寒. 2008. 多接收器等离子体质谱(MC-ICP-MS)测定Mg同位素方法研究. 岩石矿物学杂志, 27(5): 441-448