Deep magma evolution in the extensional Youjiang Basin in late Yanshanian period: Evidence from geochemical characteristics of Baiceng ultramafic rock, Guizhou Province
-
摘要:
白层超基性岩产出于扬子地台西南缘,该区域存在多期次的岩浆活动,出露多处超基性岩体。对白层11件超基性岩样品进行了主量和微量元素(包括稀土元素)及铂族元素分析,结果表明,白层超基性岩主量元素变化范围不大,具低钠富钾钙碱性岩石的地球化学特征;该超基性岩富集Ba、Sr、Rb等大离子亲石元素及Nb等高场强元素;稀土元素总量高,球粒陨石标准化曲线呈强烈的右倾趋势,轻、重稀土元素分馏明显,Eu、Ce异常不明显;铂族元素总含量低,且分异不明显。该超基性岩的地球化学特征显示,其形成于燕山晚期右江盆地大规模岩石圈伸展减薄的构造背景下;岩浆来源于低程度部分熔融的超镁铁质地幔,石榴子石及硫化物在熔融过程中残留于源区。岩浆在上升过程中未发生明显的地壳物质混染。岩浆结晶分异程度不高,发生了橄榄石与单斜辉石的结晶分异但未发生斜长石的分离。
Abstract:The ultramafic rocks in Baiceng area are located on the southwest margin of Yangtze block. There exist multi-stage magmatic activities and quite a few ultramafic rocks in this area. Lots of research work on ultramafic rocks has been done in this area, but some important questions remain unsolved. In this paper, the authors analyzed major elements, minor elements (including REE) and PGE for 11 samples of ultramafic rocks from Baiceng area. The results show that their major elements vary in a narrow range, implying the characteristics of high potassium and low sodium. These rocks are also enriched in LILE (Ba, Sr, Rb), HFSE (Nb) and rare earth elements. The enrichment of LREE forms a right inclined curve, which indicates definite fractionation of REE. The anomalies of Eu and Ce are negative. The total content of platinum group element is low, and the differentiation is not significant. Geochemical characteristics of these ultramafic rocks indicate that the rocks were formed under a tectonic setting of large-scale lithospheric extension in Youjiang Basin. The magma was derived from the mantle with a low degree of partial melting, with garnet and sulfide being left over during the melting. Obvious crustal contamination never happened during the process of magma rising. The crystallization and differentiation level of magma was not high. The fractionation of clinopyroxene and olivine occurred during magma intruding while there was no fractionation of plagioclase.
-
北山地区花岗岩类岩石分布广泛,约占该区全部侵入岩的95%[1].虽然前人对本区的花岗岩类进行了大量的研究,积累了丰富的资料[2-19] ,但对北山同一时代花岗岩构造属性的认识却不相同,尤其是华力西期花岗岩,对其构造属性[20-21]和岩浆成因的认识存在分歧.对北山花岗岩地区构造属性、岩浆成因特点和与区域构造演化关系的认识影响了对区域动力学的进一步探讨.另外,花岗岩体的多集中在北山北带,而北山南带很少.目前仅有3 处石炭纪—二叠纪花岗岩的精确定年报道:桥湾北花岗岩303.7±2.4Ma[22]、音凹峡南花岗岩281.7±2.9Ma[23]、石板泉花岗岩280.5±5.5Ma[24].因此,对北山南带的花岗岩体形成时代、岩石成因的研究,对于研究其构造岩浆作用具有重要的意义.近年来,北山南带东段地区发现了一系列与花岗岩类侵入体密切相关的金属矿床,如白山堂铜矿[25]、玉山钨矿[26-28]等,对花岗岩类侵入体的认识在一定程度上制约了该地区的找矿勘探工作.基于此,本文针对北山南带东段梧桐沟马山花岗岩体,进行了详细的地球化学测试和锆石U-Pb 测年,对其形成时代和岩石成因进行约束,并结合区域资料,对其地球动力学意义进行探讨,以期为北山南带构造岩浆研究提供可靠的地质资料.
1. 区域地质背景
研究区梧桐沟—神螺山一带位于北山南带东段,内蒙古西南侧与甘肃省接壤位置.地质背景较复杂,左国朝等[3]、龚全胜等[29]、何世平等[30]认为其属于塔里木板块增生带北缘,并将其进一步划归为磁海-红柳园-白山堂晚古生代陆内裂谷带;刘雪亚等[31]、张新虎[32]、聂凤军等[1]认为其属于哈萨克斯坦板块的南缘.区内以晚古生代石炭纪—二叠纪火山-沉积地层为主,少量前震旦纪基底形成规模较小的推覆体,奥陶系浅变质碎屑岩发育;下石炭统红柳园组可分为上、下2 段,下段为正粒序的碎屑岩组合,底部为砾岩,向上变细,上段为火山熔岩及碎屑岩组合,熔岩包括安山岩和顶部的流纹岩,少量英安岩及流纹质熔结凝灰岩,局部为灰岩,上石炭统芨芨台子组为白云质灰岩;下二叠统双堡塘组为磨拉石组合,由砾岩、含砾砂岩、粗砂岩、杂砂岩组成正粒序层;双堡塘组与红柳园组为平行不整合接触或断层接触,与上覆中二叠统金塔组为整合接触,金塔组以大面积出露玄武岩、枕状玄武岩及少量安山岩、火山-沉积碎屑岩为主,上二叠统方山口组以凝灰质为主,夹少量沉凝灰岩,与金塔组无直接接触关系;差异性沉降地区为中新生界不整合覆盖.侵入岩以晚华力西期为主,主要岩石类型包括英云闪长岩、二长花岗岩、黑云母二长花岗岩、花岗闪长岩,多以岩株或小岩基产出,基性侵入岩十分发育,主要为辉长岩及辉绿岩,前者为小岩株,后者多为岩脉.
2. 岩体地质特征
马山小岩体位于北山南带东段梧桐沟一带(图 1-a),岩体外表似椭球状,长约2.7km, 宽约1km, 出露面积约2.67km2,中心坐标为北纬40°37′45″、东经98°35′30″.侵入下石炭统红柳园组砂岩中,后被下二叠统双堡塘组砾岩、含砾砂岩覆盖.前人粗略地将其归入到华力西中晚期岩浆侵入活动①.
① 甘肃省地质局第二区域地质测量队. 红柳大泉幅1∶200000 区域地质测量报告. 1971.
图 1 马山二长花岗岩体地质简图及大地构造位置图(b 据参考文献[30]修改)1—第四系;2—下二叠统双堡塘组;3—下石炭统红柳园组;4—二叠纪花岗岩;5—二叠纪辉长岩;6—不整合界线;7—锆石U-Pb 采样位置;8—国界;9—板块缝合带;10—早古生代缝合带;11—主干断裂;12—省界;13—隐伏断裂;14—构造单元.Ⅰ—西伯利亚板块;Ⅰ-1—大南湖-雀儿山-狐狸山早古生代活动带陆缘带;Ⅱ—塔里木板块;Ⅱ-1-1 —黄山-红石山-路井晚古生代陆内裂谷带;Ⅱ-1-2— 星星峡-明水-旱山地块;Ⅱ-1-3 —白玉山南-公婆泉-东七一山早古生代晚期活动带陆缘带;Ⅱ-2 —红柳河-洗肠井构造混杂岩带;Ⅱ-3-1 —方山口-营毛沱-鹰咀红山早古生代中期活动陆缘带;Ⅱ-3-2 —花牛山早古生代陆缘裂谷带(裂陷槽);Ⅱ-3-3— 磁海-红柳园-白山堂晚古生代陆内裂谷带;Ⅱ-3-4 —敦煌地块Figure 1. Simplified geological map and tectonic location map of Mashan granite该岩体由花岗闪长岩组成(图 2-a),其内部无明显结构分带和成分分带,为一次岩浆侵入活动的产物.岩石具有典型的花岗结构(图 2-b),矿物成分较简单,主要由斜长石、条纹长石和石英组成.长石0.4~3.8mm, 整体表面褐色粘土化显著,较石英脏,斜长石半自形-自形板状,双晶纹可见;石英呈他形粒状,表面干净透明,波状消光,粒径为0.3~2.7mm, 充填于其他矿物之间.
锆石U-Pb 年龄样品的采集位置为北纬40°37′46″、东经98°35′30″.
3. 分析方法
对采集的样品进行主量、稀土和微量元素及LA-ICP-MS 锆石U-Pb 年龄测定.在进行化学分析前,先在镜下观察薄片以确定样品的适应性.主量、微量元素由中国冶金地质总局西北地质勘查院测试中心测定,采用XRF 荧光光谱样分析方法,分析精度优于5%;微量和稀土元素采用ICP-MS 分析,利用国家一级标准物质进行质量监控(GB/T14506—1993),含量大于10×10-6的元素测试精度为5%,小于10×10-6的测试精度为10%.
用常规方法分选出锆石单矿物,然后在双目镜下根据锆石颜色、自形程度、形态、透明度等进行初步分类,挑选出具有代表性的锆石.将锆石样品用双面胶粘在载玻片上,放上PVC 杯,然后将环氧树脂和固化剂进行充分混合后注入PVC 杯中,待树脂充分固化后将其从载玻片上剥离,打磨和抛光至锆石中心部位暴露,然后拍摄透射光、反射光和阴极发光(CL)图像.最后用体积百分比为3%的HNO3清洗样品制成样品靶备用.
LA-ICP-MS 锆石U-Pb 同位素测定在西安地质矿产研究所实验测试中心完成.测定时根据可见光、CL 和BSE(背散射电子成像)图像选择合适的测点位置,如避免包裹体、裂隙位置等.详细的分析步骤和数据处理方法见参考文献[33-34],用Glit-ter(4.0 版)程序处理原始数据.普通铅扣除时,其组成由Stacey 等[35]的模式给出,年龄加权平均计算及U-Pb 谐和图的绘制采用Isoplot(3.0 版)[36]完成.
4. 锆石U-Pb 年龄
从样品中选取的锆石以浅黄色为主,个别颜色为浅玫瑰色和近于无色,玻璃-金刚光泽.按晶体形态可以分为2 类,一类呈自形双锥柱状,长60~160μm, 宽30~120μm, 长宽比多为2~3,透射光下可见锆石内部发育少量裂纹和包裹体,包裹体以长条形和椭圆形为主,多为小锆石包体,阴极发光图像显示,锆石具有典型的岩浆韵律环带和明暗相间的条带构造等,显示岩浆锆石的特征(图 3);另一类则有溶蚀现象,锆石长20~30μm, 宽100~120μm, 长宽比多为4~5.部分锆石具有残留的核部,为继承核或捕获核,部分锆石颗粒具有窄的浅色边,但核部仍显示出清晰的岩浆环带特征.测得的28 个颗粒的28 个数据中(表 1),21 个206Pb/238U 年龄(不包括1、4、10、12、14、20、28 点)介于270.7±6.0~294.0±8.0Ma 之间,给出的206Pb/238U 年龄加权平均值为281.8 ± 3.2Ma(MSWD=1.15,n=21,2σ)(图 4).U、Th 含量分别为155×10-6~1444×10-6和92×10-6~660×10-6.Th/U 值除2 个样品略低于0.4(0.39、0.31)外,其余的比值均介于0.42~0.99 之间,表现出典型岩浆锆石的特征.该年龄代表了岩体的结晶年龄.此外还获得了3 组年龄:①10、12、20 这3 个点的206Pb/238U 年龄分别为258.4±6.0Ma、258.2±7.6Ma、260.8±6.3Ma, 其年龄加权平均值为259.2±0.1Ma, 锆石具有较规则的外形,CL 图像较亮,具有弱分带特征,为热液变质锆石,可能代表一期热液蚀变活动;② 有2 个点获得了334.4Ma、369.9Ma 的年龄数据.这些锆石略磨圆,锆石环带特征不清楚,疑为继承性锆石;③1、28 这2 个点的年龄数据采用207Pb/206Pb 年龄,分别为2487±35.9Ma、1784.1±40.2Ma, 锆石外形不规则,具磨圆,环带特征不清楚,可能为老基底的捕获锆石.
表 1 马山花岗岩体LA-ICP-MS 锆石(TW01) U-Th-Pb 同位素数据Table 1. LA-ICP-MS zircon U-Th-Pb dating data of granite from Mashan(TW01)分析点 含量/10-6 232Th/238U 207Pb/235U 206Pb/238U 207Pb/235U 206Pb/238U 232Th 238U 比值 1σ 比值 1σ 年龄/Ma 1σ 年龄/Ma 1σ TW01-1 280.96 449.94 0.62 9.3972 0.2158 0.4182 0.0092 2377.6 21.07 2252.3 41.73 TW01-2 523.95 1099.59 0.48 0.3621 0.0089 0.0437 0.0010 313.8 6.64 275.4 5.97 TW01-3 375.06 699.75 0.54 0.3941 0.0100 0.0456 0.0010 337.4 7.29 287.1 6.31 TW01-4 594.83 576.85 1.03 0.8243 0.0212 0.0526 0.0012 610.4 11.79 330.4 7.55 TW01-5 349.29 359.94 0.97 0.3703 0.0110 0.0431 0.0010 319.9 8.13 272 6.38 TW01-6 92.39 154.58 0.60 0.4218 0.0152 0.0456 0.0012 357.3 10.86 287.3 7.47 TW01-7 401.97 1310.93 0.31 0.3223 0.0078 0.0438 0.0010 283.7 6.01 276.4 5.96 TW01-8 191.06 372.5 0.51 0.3463 0.0101 0.0438 0.0010 301.9 7.58 276.3 6.38 TW01-9 292.29 294.14 0.99 0.3421 0.0102 0.0444 0.0011 298.8 7.73 280 6.53 TW01-10 99.79 354.71 0.28 0.3287 0.0093 0.0409 0.0010 288.5 7.14 258.4 5.96 TW01-11 408.08 804.03 0.51 0.3764 0.0095 0.0436 0.0010 324.4 7.02 275 6.06 TW01-12 49.69 106.02 0.47 0.4079 0.0170 0.0409 0.0012 347.4 12.27 258.2 7.6 TW01-13 133.16 156.37 0.85 0.4344 0.0141 0.0433 0.0011 366.3 10 273 6.93 TW01-14 81.57 213.17 0.38 0.5080 0.0149 0.0603 0.0015 417.1 10.03 377.7 8.81 TW01-15 291.58 308.2 0.95 0.3552 0.0104 0.0437 0.0010 308.6 7.79 275.7 6.45 TW01-16 252.32 605.27 0.42 0.3590 0.0094 0.0429 0.0010 311.5 7 270.7 6.04 TW01-17 660.33 1444.13 0.46 0.4373 0.0104 0.0461 0.0010 368.4 7.36 290.5 6.29 TW01-18 216.72 353.49 0.61 0.4315 0.0124 0.0451 0.0011 364.2 8.77 284.6 6.69 TW01-19 477.93 938.69 0.51 0.3581 0.0090 0.0457 0.0010 310.8 6.72 287.9 6.32 TW01-20 212.5 325.28 0.65 0.3001 0.0098 0.0413 0.0010 266.4 7.68 260.8 6.33 TW01-21 108.44 260.92 0.42 0.3980 0.0129 0.0456 0.0012 340.2 9.38 287.2 7.13 TW01-22 158.92 169.8 0.94 0.4177 0.0160 0.0467 0.0013 354.4 11.43 294 8.04 TW01-23 372.22 946.18 0.39 0.3493 0.0087 0.0449 0.0010 304.2 6.53 283.2 6.21 TW01-24 315.24 351.81 0.90 0.4617 0.0139 0.0457 0.0011 385.5 9.65 288.2 7.02 TW01-25 261.36 418.89 0.62 0.3418 0.0098 0.0458 0.0011 298.6 7.41 288.6 6.66 TW01-26 157.57 259.58 0.61 0.3533 0.0108 0.0458 0.0011 307.2 8.13 288.7 6.87 TW01-27 157.3 231.78 0.68 0.3569 0.0112 0.0456 0.0011 309.9 8.4 287.7 6.93 TW01-28 77.52 182.09 0.43 4.3637 0.1037 0.2902 0.0066 1705.5 19.64 1642.3 32.72 5. 岩石地球化学特征
5.1 主量元素
马山花岗岩类的岩石地球化学分析(表 2)表明,SiO2 含量为64.85%~69.75%,平均含量67.35%,属于酸性岩类,在SiO2-(Na2O+K2O)图解(图 5)上,大多数样品点投入花岗闪长岩区域;Na2O 含量为2.66%~4.72%,平均值3.68%,K2O含量为1.90%~2.71%, 平均值为2.23% ,Na2O + K2O 含量为5.13~6.62%,平均值为5.91%, Na2O/K2O 值为0.98~2.48,平均值为1.69;Al2O3含量为13.55%~15.49%,平均值为14.63% ;TiO2 含量为0.41% ~0.76% ,平均值为0.55% ;A/CNK 值为0.96~1.07,平均值为1.0,在A/CNK-A/NK 图解(图 6)中,样品点主要投入到偏铝质岩石范围内,2 个样品投入过铝质范围内;里特曼指数σ为1.15~1.84,为钙碱性系列岩石;在SiO2-K2O 图解(图 7)上,样品点全部投入中钾钙碱性系列;在K2O-Na2O 图解(图 8)上,样品点均投入I 型花岗岩区域.岩石的分异指数(DI)为65.74~76.80,表明原始岩浆的结晶分异程度一般.马山岩体主量元素特征显示,其具有钙碱性序列岩石的属性,属于中等分异程度的I 型花岗岩.
图 5 SiO2-(Na2O +K2O)图解(底图据参考文献[37])Figure 5. SiO2 versus Na2O +K2O diagram图 6 A/CNK-A/NK 图解(底图据参考文献[38])Figure 6. A/CNK versus A/NK diagram图 7 SiO2-K2O 图解(底图据参考文献[39])Figure 7. SiO2 versus K2O diagram图 8 K2O-Na2O 图解(底图据参考文献[40])Figure 8. K2O versus Na2O diagram5.2 微量和稀土元素
(1) 稀土元素
稀土元素分析结果如表 2 所示.样品的稀土元素总量为80.57 × 10-6~125.29 × 10-6,平均值为106.44×10-6,稀土元素总含量较低;在球粒陨石标准化配分模式图(图 9-a)上,具有相对平缓的富集轻稀土元素(LREE)的稀土元素配分模式,重稀土元素(HREE)分馏不明显且相对亏损,但都明显高于球粒陨石丰度10 倍.ΣLREE/ΣHREE 值为1.70~2.43, 平均值为1.99;轻稀土元素内部分异较明显,(La/Sm)N=2.37~3.37,平均值为2.80;La/Yb=5.81~7.23,平均值为6.23,(La/Yb)N=4.17~5.19,平均值为4.47;(Gd/Yb)N=0.85~1.25,δEu=0.65~0.91,平均值为0.80,具有弱负Eu 异常,可能与斜长石在源区残留有关.
表 2 马山花岗岩体主量、微量和稀土元素分析结果Table 2. Analytical data of major, trace, and rare earth elements of granite from Mashan样品号
岩石名称P1/GS1
花岗质碎裂岩P1/GS3
花岗质碎岩P13/GS3
中粗粒花岗闪长岩P13/GS5
中粗粒花岗闪长岩P16/GS1、XT1
黑云母花岗闪长岩P16/GS2、XT2
黑云母花岗闪长岩SiO2 65.48 68.06 69.75 69.21 64.85 66.77 Al2O3 15.00 13.55 14.09 14.33 15.30 15.49 Fe2O3 1.40 1.13 1.95 2.08 3.41 3.01 FeO 3.44 3.02 1.68 1.75 2.22 1.88 TiO2 0.76 0.62 0.41 0.43 0.56 0.53 MgO 2.55 2.05 1.36 1.20 1.78 1.62 CaO 3.92 2.93 3.15 3.57 3.48 3.25 P2O5 0.16 0.16 0.07 0.09 0.12 0.12 K2O 1.90 2.71 2.40 2.29 2.21 1.90 MnO 0.09 0.07 0.07 0.08 0.08 0.07 Na2O 3.23 2.66 3.74 3.67 4.05 4.72 烧失量 1.55 3.72 0.82 0.54 1.60 1.29 总量 99.46 100.69 99.49 99.23 99.66 100.65 A/CNK 1.04 1.07 0.97 0.96 0.99 0.99 NK/A 0.49 0.54 0.62 0.59 0.59 0.63 σ 1.17 1.15 1.41 1.35 1.79 1.84 Rb 88.84 114.35 75.02 170.00 47.34 48.32 Ba 430.78 435.03 370.10 481.29 109.25 283.55 Th 12.55 13.66 8.49 17.14 6.92 4.94 U 2.10 1.30 1.44 3.79 1.68 1.40 Nb 13.36 11.54 5.78 12.36 5.75 4.42 Ta 1.18 1.20 0.62 1.28 0.64 0.49 Sr 255.42 139.70 209.88 195.36 291.65 226.93 Zr 143.82 72.81 127.41 374.10 196.08 185.01 V 101.52 64.53 71.65 23.01 78.67 64.00 Cr 54.31 43.16 22.04 21.90 18.38 22.06 Co 13.38 6.74 6.29 1.34 10.11 7.72 Cu 15.99 16.80 7.21 3.16 22.49 17.54 Zn 76.13 53.91 33.03 41.80 47.13 39.88 Ga 18.83 16.43 15.27 16.56 14.36 13.37 La 18.00 18.59 15.70 12.85 14.47 13.48 Ce 24.17 31.86 31.06 19.91 37.94 34.31 Pr 5.07 4.90 3.65 3.23 4.69 3.94 Nd 19.47 18.51 14.39 12.23 18.34 15.69 Sm 4.14 4.16 3.01 2.83 3.94 3.60 Eu 1.08 0.85 0.86 0.74 0.91 0.90 Gd 3.76 3.84 2.77 2.56 3.31 3.16 Tb 0.74 0.81 0.55 0.48 0.60 0.59 Dy 4.40 4.95 3.57 3.00 3.52 3.57 Ho 0.88 1.00 0.81 0.66 0.75 0.76 Er 2.80 3.34 2.57 2.21 2.08 2.10 Tm 0.40 0.55 0.37 0.34 0.34 0.33 Yb 2.49 3.19 2.70 2.06 2.29 2.28 Lu 0.40 0.45 0.38 0.33 0.36 0.36 Y 23.71 28.29 20.30 17.14 19.75 20.25 REE 111.51 125.29 102.69 80.57 113.29 105.31 LREE 71.93 78.87 68.67 51.79 80.29 71.91 HREE 39.58 46.42 34.02 28.78 33.00 33.39 LREE/HREE 1.82 1.70 2.02 1.80 2.43 2.15 (La/Yb)N 5.19 4.18 4.17 4.47 4.54 4.25 (La/Sm)N 2.81 2.88 3.37 2.93 2.37 2.42 (Gd/Yb)N 1.25 1.00 0.85 1.03 1.20 1.15 δEu 0.84 0.65 0.91 0.84 0.77 0.82 注:元素由中国冶金地质总局西北地质勘查院测试中心测定;δEu=EuN/SQRT(SmN∗GdN);主量元素含量单位为%,微量和稀土元素为10-6 (2) 微量元素
在微量元素原始地幔标准化蛛网图(图 9-b)上,所有样品的稀土元素配分模式和原始地幔标准化图基本平行,显示同源演化的特征;显著亏损Ba、Nb、Ta、P、Ti, Zr、Rb、Th、U、K,富集Zr 和亏损Nb、Ta、Ti, 表明岩浆源区岩石以陆壳组分为主.
岩石具有高Rb、低Sr(均小于400×10-6)、高Yb(大于2 × 10-6)的特征,Sr/Y 值很低(Sr/Y=4.94~14.77, 平均值为10.57),马山花岗岩微量元素Sr、Yb特征与华南低Sr、高Yb 型花岗岩相似.低Sr 、高Yb 型花岗岩为浙闽型花岗岩,与其平衡的是斜长石和角闪石残留相[38],说明形成的压力较低(小于0.8GPa 或1.0GPa),可能是地壳伸展变形减薄达正常地壳厚度(30km 左右)或低于正常地壳厚度的造山后阶段形成的[43].在大洋中脊标准化图解(图 10)上,样品与典型的后碰撞花岗岩显示出类似的趋势,说明马山花岗岩为典型的后碰撞花岗岩.
图 10 马山花岗岩体大洋中脊标准化图解(标准化值据参考文献[42])Figure 10. Ocean ridge granite (ORG) normalized geochemicalpatterns for samples of Mashan granite6. 马山花岗岩体形成的构造背景及其成因
北山地区西邻东天山,东接阿拉善,以阿尔金和星星峡两大走滑断裂为界,位于一个巨大的构造楔形区内,构造过程复杂.其大地构造的归属一直存在争议,一是将北山地区划分为北部的哈萨克斯坦板块和南部的塔里木板块,左国朝等[3]以明水-石板井-小黄山缝合带(早古生代末)为界;龚全胜等[29]、何世平等[30]则以红石山-黑鹰山-六陀山蛇绿混杂岩带(晚古生代)为界;部分学者[31-32]以柳园-大奇山(晚古生代末)和红石山-黑鹰山-六陀山2 条深大断裂为界,将北山从南向北依次划分为塔里木板块、哈萨克斯坦板块和西伯利亚板块;李锦轶等[44]则以星星峡—白玉山—牛圈子—小黄山作为北山地区西伯利亚与中朝2 个古板块之间二叠系的分界线,认为北山南部柳园一带可能代表了二叠纪弧后盆地环境,柳园与白玉山—小黄山之间的区域为早古生代岛弧带.
对区域地质演化,前人有不同的观点,分歧的关键点在于洋盆的最终闭合时间.左国朝等[3]认为,在志留纪末—早泥盆世,北山地区的塔里木-中朝与哈萨克斯坦两大古板块最终拼合,并导致洋盆消失和碰撞造山;刘雪亚等[31]认为,在早二叠世之前,随着北山南带及南天山古洋盆的封闭,敦煌地块北缘的安北-旧寺墩构造带与北山造山带前缘的柳园-大奇山地体碰撞,导致塔里木板块、哈萨克斯坦板块和西伯利亚板块最终拼接;龚全胜等[45]认为,北山地区自泥盆纪开始转化为古亚洲洋构造域演化体系,晚石炭世末实现塔里木板块、哈萨克斯坦板块的最终碰撞对接;何世平等[46]认为,石炭纪末,哈萨克斯坦板块和塔里木板块之间的洋盆最终闭合,形成新的统一大陆,认为在早二叠世研究区洋盆已经闭合,经历了碰撞造山.而朱江[24]则认为,北山南带洋壳俯冲作用可能持续到早二叠世,并在早—中二叠世洋盆最终闭合,认为北山南带石炭纪—二叠纪的岩浆作用形成于俯冲环境,与毛启贵[47]、肖文交[48-49]、郭谦谦[50]等观点一致.
很多学者[3, 51-52]认为,北山地区在二叠纪是一种伸展拉伸的构造背景.左国朝等[3]认为,早二叠世北山南带在拉张背景下形成了峡东-俞井子裂陷槽,但并没有形成真正的大洋;刘明强等[53]指出,音凹峡晚古生代陆内裂谷带为一多旋回裂谷,主要裂谷作用发生于奥陶纪、志留纪、石炭纪和二叠纪,反映了多期次、多旋回性、继承性“开”、“合”演化的特征;龚全胜等[45]认为,北山地区晚石炭世—二叠纪为大陆板块碰撞时期;何世平等[46]认为,早石炭世沿石板山—大奇山—神螺山一带在前震旦纪古老基底上形成北山南部陆内裂谷带,该裂谷带一直发展演化到二叠纪末,并在早二叠世发展到鼎盛时期.
综上所述,结合野外地质特征,下石炭统红柳园组为火山-碎屑岩,火山岩组合为安山岩-英安岩-流纹岩组合;与红柳园组平行不整合接触的下二叠统双堡塘组为一套粗碎屑岩和细碎屑岩组合,中二叠统金塔组以出露海相枕状玄武岩为特征,并伴随大量的辉长岩、辉绿岩等幔源岩浆活动,其与双堡塘组为整合接触关系.多数学者[3, 46, 53]对区内北山南带存在晚古生代洋盆有异议,认为晚古生代为陆内裂谷环境,区内晚石炭世可能存在裂谷回返碰撞,导致石炭系与下二叠统之间的平行不整合接触,早二叠世裂谷进一步拉开,并在中二叠世达到鼎盛.结合地球化学特征,可进一步确定281.8±3.2Ma的马山花岗岩形成于后碰撞伸展环境.
吴泰然[54]将这类花岗岩归为拉张型过渡壳花岗岩(ECG),认为其形成是由于地壳的拉伸减薄,上地幔热物质上涌,使地壳形成一种高温低压的环境,并使地壳发生部分重熔,同时上地幔上涌的热物质沿拉张的裂隙与地壳的热物质发生混染作用,对陆壳的物质进行改造,使之向过渡类型转化,形成拉张型过渡壳的花岗岩.冯继承[22]、张文等[23]研究认为,北山南带音凹峡花岗岩也是该类型的花岗岩体.区域上,前人对北山南带不同地段的石炭纪—二叠纪代表性的中酸性侵入岩体进行了Sr、Nd、Hf 同位素研究,桥湾北花岗岩体[22, 24]的εNd(t)值在-0.40~-0.06 之间,εHf(t)值在-1.2~5.8 之间,是典型的壳幔混合型花岗岩,代表了兴蒙造山带广泛出现的晚古生代—中生代大陆地壳生长现象;音凹峡南花岗岩体也位于北山南带,张文[23]报道的该岩体的εHf(t)值在+4.4~+7.8 之间;石板泉花岗岩体的εNd(t)值在-0.1~-1.6 之间[24].一般认为,具有正εNd(t)和εHf(t)值的花岗质岩石来自于亏损地幔,或由亏损地幔中新增生的年轻地壳物质的部分熔融形成.结合研究区二叠纪大量基性岩浆活动(枕状玄武岩和辉长、辉绿岩的发现),认为北山南带石炭纪—二叠纪(305~280Ma)的花岗岩具有普遍的壳幔混合成因特征.
7. 结论
(1) 马山花岗岩体的LA-ICP-MS 锆石U-Pb年龄为281.8±3.2Ma, 为后碰撞花岗岩.
(2) 马山花岗岩体具有典型的壳幔混合特征,其机制可能是幔源基性岩浆底侵后遭受地壳物质混合或混染.北山南带305~280Ma 之间的花岗岩普遍具有壳幔混合成因的特征.
(3) 马山花岗岩体及北山南带众多壳幔混合成因的花岗岩体说明,北山南带在早二叠世总体为后碰撞伸展环境.
致谢: 感谢审稿专家提出的宝贵意见;本次野外工作、室内数据处理及成文过程中得到峨眉403地质队王疆丽硕士、中国地质调查局成都地质调查中心王秀平的帮助,在此一并表示衷心感谢. -
图 1 贵州黔西南大地构造分区及岩浆岩分布(据参考文献[14]修改)
1—峨眉山玄武岩;2—偏碱性玄武岩;3—辉绿岩;4—超基性岩;Ⅳ-4-1-1—威宁隆起区;Ⅳ-4-1-2—六盘水裂陷槽;Ⅳ-4-2-1—兴义隆起区;Ⅳ-4-2-2—右江裂谷-前陆盆地地区;Ⅳ-4-2-3—黔南坳陷区;Ⅳ-4-1-3(1)—织金穹盆构造变形区;Ⅳ-4-2-2(1)—册亨东西向紧闭褶皱变形区;Ⅳ-4-2-2(2)—望谟北西向褶皱带;Ⅳ-4-2-3(1)—都匀南北向隔槽式褶皱变形区
Figure 1. Geotectonic divisions and distribution of magmatic rocks in southwestern Guizhou Provence
图 4 白层超基性岩微量元素蛛网图(a)和稀土元素配分图(b)
(球粒陨石标准化数据据参考文献[18])
Figure 4. Trace element spider diagram and REE patterns of ultramafic rocks in Baiceng
图 5 白层超基性岩铂族元素原始地幔标准化图
(标准化数据据参考文献[22])
Figure 5. Primitive mantle-normalized PGE patterns for Baiceng ultramafic rocks
图 11 白层超基性岩La/Sm-Sm/Yb图解(底图据参考文献[48])
Figure 11. Diagrams of La/Sm-Sm/Yb of ultramafic rocks in Baiceng
表 1 白层超基性岩主量元素分析结果及特征值
Table 1 Major elements compositions of ultramafic rocks in Baiceng
% 样品号 εk52-1 εk52-2 εk52-3 εk52-4 εk52-5 εk52-6 εk52-8 εk52-9 εk52-10 εk52-11 εk52-12 SiO2 35.44 36.25 35.52 36.31 36.19 34.62 35.47 35.62 36.47 35.43 34.64 Al2O3 11.17 11.91 11.65 11.56 11.89 11.87 11.37 11.56 11.59 11.24 11.21 TFe 7.24 7.52 7.43 7.72 7.73 7.14 7.27 7.88 7.61 7.4 7.38 MgO 7.29 7.78 7.25 7.51 7.26 6.33 7.26 7.78 7.67 7.87 7.21 CaO 16.41 15.59 15.77 14.81 14.64 18.62 16.07 15.12 14.84 15.51 16.53 K2O 4.12 3.65 3.82 4.13 4.18 3.54 4.08 3.67 3.77 4.29 3.88 Na2O 0.56 0.41 0.49 0.57 0.44 0.44 0.55 0.48 0.47 0.5 0.51 TiO2 0.79 0.72 0.86 0.85 0.85 0.7 0.85 0.81 0.87 0.82 0.72 P2O5 2.01 1.89 1.97 1.78 2.07 1.72 1.87 2.04 2.07 2.1 1.87 烧失量 7.87 7.71 8.04 8.24 7.90 8.69 8.32 8.85 8.14 7.72 9.87 总和 99.36 99.27 98.88 98.66 98.40 99.28 99.49 99.80 99.99 98.39 99.93 Mg# 50.11 48.74 48.25 46.75 46.12 47.39 48.29 49.11 49.47 49.25 48.43 Na2O+K2O 4.68 4.06 4.31 4.7 4.62 3.98 4.63 4.15 4.24 4.79 4.39 Na2O/K2O 0.14 0.11 0.13 0.14 0.11 0.12 0.13 0.13 0.12 0.09 0.13 表 2 白层超基性岩微量元素分析结果
Table 2 Trace elements compositions of ultramafic rocks in Baiceng
10-6 样品号 εk52-1 εk52-2 εk52-3 εk52-4 εk52-5 εk52-6 εk52-8 εk52-9 εk52-10 εk52-11 εk52-12 Rb 139.47 127.06 147.02 173.73 147.94 127.21 145.53 155.48 143.46 162.19 172.79 Ba 5621 6062 6469 3699 6140 5946 6425 6338 6148 5903 5111 Th 56.5 60.83 60.64 59.82 65.47 60.85 58.82 63.18 60.15 57.98 62.48 U 7.96 8.53 7.73 8.28 8.46 8.98 8.2 10.06 9.42 7.72 9.03 Sc 24.51 24.72 24.55 22.15 22.65 22.28 23.65 25.19 24.24 24.98 21.12 Ta 5.59 5.74 5.71 5.93 5.79 6.07 5.75 6.36 6.22 5.53 5.79 Co 33.6 34.7 34.8 23.3 34.6 31.6 34.5 36.9 35.3 34.1 32.7 Ni 58.9 60.8 59.6 36.4 55.7 52.2 58.8 60.6 58.7 58.5 52.2 Nb 216.26 248.67 228.66 257.14 259.27 245.91 243.27 285.38 269.35 236.53 245.75 Sr 1420.4 1758.08 1678.65 1680.03 1640.44 1519.68 1617.03 1661.91 1618.57 1629.27 2002.21 Cr 128.06 126.81 129.86 108.68 117.64 110.1 120.17 128.32 125.02 124.78 108.41 Zr 272.96 299.57 276.25 267.79 266.67 250.79 263.64 278.15 279.97 268.38 255.44 Hf 5.18 5.25 5.14 4.75 4.71 4.27 4.72 5.15 4.87 4.83 4.37 La 306 329 324 325 343 307 326 343 332 320 318 Ce 513 549 537 537 553 504 541 572 560 540 517 Pr 50.8 54.1 52.2 52 53.7 49.4 52.7 55 54 53.3 50.5 Nd 170 180 175 170 175 162 174 178 177 176 166 Sm 22 22.9 22.6 21.4 22.1 20.8 22.2 22.9 22.1 22.7 20.4 Eu 5.04 5.44 5.14 5 5.11 4.69 5.15 5.3 5.2 5.21 4.85 Gd 16.4 17.7 16.9 15.9 16.9 16 16.3 17 17.6 17.1 15.7 Tb 1.7 1.78 1.66 1.63 1.7 1.6 1.66 1.75 1.73 1.63 1.55 Dy 6.7 6.93 6.62 6.54 6.61 6.36 6.55 6.87 6.83 6.59 6.02 Ho 1.08 1.15 1.09 1.04 1.06 1.02 1.06 1.14 1.12 1.08 1 Er 2.99 3.19 3.03 2.94 2.99 2.88 2.94 3.23 3.01 2.92 2.85 Tm 0.34 0.37 0.33 0.34 0.36 0.33 0.33 0.38 0.34 0.34 0.33 Yb 2.17 2.22 2.15 2.13 2.13 2.11 2.05 2.29 2.25 2.08 2.06 Lu 0.32 0.34 0.32 0.31 0.31 0.3 0.31 0.34 0.32 0.3 0.29 ΣREE 1098.48 1174.38 1148.08 1141.37 1184.2 1079.45 1152.83 1208.22 1183.5 1149.54 1106.72 LREE 1066.74 1140.71 1115.98 1110.52 1152.15 1048.89 1121.61 1175.21 1150.27 1117.53 1076.87 HREE 31.74 33.67 32.1 30.85 32.05 30.57 31.22 33.01 33.23 32.01 29.85 LREE/HREE 33.61 33.88 34.77 36 35.95 34.32 35.92 35.6 34.62 34.92 36.07 (La/Yb)N 97.77 103.01 104.84 106.12 111.89 101.15 110.47 103.74 102.2 106.9 106.83 (La/Sm)N 8.73 9.03 9.01 9.53 9.74 9.31 9.23 9.39 9.43 8.87 9.77 (Gd/Yb)N 6.25 6.59 6.51 6.2 6.56 6.25 6.59 6.13 6.46 6.78 6.3 δEu 0.78 0.79 0.77 0.79 0.78 0.76 0.79 0.79 0.78 0.78 0.8 δCe 0.92 0.92 0.91 0.91 0.9 0.91 0.92 0.92 0.93 0.92 0.9 表 3 白层超基性岩铂族元素分析结果及特征值
Table 3 PGE compositions and characteristic parameters of Baiceng ultramafic rock
10-9 样品号 Ir Ru Rh Pt Pd ∑PGE PPGE/IPGE Pd/Ir Pd/Ru Pd/Rh Pt/Pd Pt/Pt* εk52-1 0.165 2.69 0.108 0.051 2.62 2.676 0.975 15.941 0.976 24.262 0.020 2.970 εk52-2 0.124 2.34 0.111 0.037 1.01 1.046 0.469 8.142 0.431 9.064 0.036 4.110 εk52-3 0.125 2.36 0.105 0.045 1.78 1.825 0.776 14.272 0.753 16.966 0.025 3.210 εk52-4 0.169 2.31 0.121 0.040 2.28 2.322 0.986 13.494 0.988 18.926 0.018 2.590 εk52-5 0.086 2.19 0.098 0.043 2.38 2.426 1.107 27.710 1.086 24.353 0.018 2.670 εk52-6 0.074 2.74 0.105 0.052 2.46 2.509 0.929 33.393 0.896 23.315 0.021 3.170 εk52-8 0.085 2.46 0.115 0.040 2.26 2.300 0.950 26.663 0.920 19.603 0.018 2.830 εk52-9 0.153 2.13 0.105 0.037 2.47 2.506 1.143 16.108 1.158 23.598 0.015 2.460 εk52-10 0.104 2.44 0.104 0.055 2.45 2.501 1.025 23.486 1.004 23.583 0.022 2.840 εk52-11 0.134 2.45 0.140 0.046 1.73 1.772 0.741 12.863 0.706 12.322 0.026 2.920 εk52-12 0.134 2.18 0.093 0.031 1.52 1.556 0.712 11.342 0.699 16.451 0.021 3.410 -
邓晋福, 苏尚国, 周肃, 等.华北地区燕山期岩石圈减薄的深部过程[J].地学前缘, 2003, 10(3):41-50. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200303004.htm Su W C, Hu R Z, Xia B, et al. Calcite Sm-Nd isochron age of the Shuiyindong Carlin-type gold deposit, Guizhou, China[J]. Chemical Geology, 2009, 258:269-274. doi: 10.1016/j.chemgeo.2008.10.030
Peter S G, Huang J Z, Li Z P, et al. Sedimentary rock-hosted Au deposits of the Dian-Qian-Gui area, Guizhou, and Yunnan Provinces, and Guangxi District, China[J]. Ore Geology Reviews, 2007, 31:170-204. doi: 10.1016/j.oregeorev.2005.03.014
朱赖民, 刘显凡, 金景福, 等.滇-黔-桂微细浸染型金矿床时空分布与成矿流体来源[J].地质科学, 1998, 33(4):463-474. http://www.cnki.com.cn/Article/CJFDTOTAL-DZKX804.008.htm 王砚耕, 索书田, 张明发.黔西南构造与卡林型金矿[M].北京:地质出版社, 1994. 陈懋弘, 毛景文, 吴六灵, 等.滇黔桂矿集区微细浸染型金矿成矿年代学研究[J].桂林工学院学报, 2006, 26(3):334-340. http://www.cnki.com.cn/Article/CJFDTOTAL-GLGX200603004.htm 韩至钧, 王砚耕, 冯济舟, 等.黔西南金矿地质与勘查[M].贵阳:贵州科技出版社, 1999. 陈懋弘, 章伟, 杨宗喜, 等.黔西南白层超基性岩墙锆石SHRIMP U-Pb年龄和Hf同位素组成研究[J].矿床地质, 2009, 28(2):240-250. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ200903003.htm Liu S, Su W C, Hu R Z, et al. Geochronological and geochemical constraints on the petrogenesis of alkaline ultramafic dykes from south-west Guizhou Province, SW China[J]. Lithos, 2010, 114:253-264. doi: 10.1016/j.lithos.2009.08.012
冯光英, 刘燊, 苏文超, 等.黔西南碱性超基性脉岩的铂族元素地球化学[J].矿物学报, 2010, 30(2):207-214. http://www.cnki.com.cn/Article/CJFDTOTAL-KWXB201002007.htm 陈懋弘, 毛景文, 屈文俊, 等.贵州贞丰烂泥沟卡林型金矿床含砷黄铁矿Re-Os同位素测年及地质意义[J].地质论评, 2007, 53(3): 371-382. http://www.cnki.com.cn/Article/CJFDTOTAL-DZLP200703009.htm 胡瑞忠, 彭建堂, 马东升, 等.扬子地块西南缘大面积低温成矿时代[J].矿床地质, 2007, 26(6):583-596. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ200706002.htm 刘建中, 邓一明, 刘川勤, 等.水银洞金矿床包裹体和同位素地球化学研究[J].贵州地质, 2006, 23(1):51-56. http://www.cnki.com.cn/Article/CJFDTOTAL-GZDZ200601011.htm 贵州省地质调查院.贵州省区域地质志[M].北京:地质出版社, 1994. Zhang X C, Spiro B, Halls C, et al. Sediment-hosted disseminated gold deposit in Southwest Guizhou PRC: Their geological setting and orgin in relation to mineralogical, fluid inclusion, and stable-isotope characteristics[J]. International Geology Review, 2003, 45:407-470. doi: 10.2747/0020-6814.45.5.407
梅厚钧.峨眉山玄武岩地球化学特征, IGCAS年度报告 (1980-1981)[M].贵阳:贵州人民出版社, 1980. 路远发. GeoKit:一个用VAB构建的地球化学工具软件包[J].地球化学, 2004, 33(5):459-464. http://www.cnki.com.cn/Article/CJFDTOTAL-DQHX200405003.htm McDonough W F, Sun S S. The composition of the earth[J]. Chemical Geology, 1995, 120: 223-230. doi: 10.1016/0009-2541(94)00140-4
Anders E, Grevese N. Abundances of the elements meteoritic and solar[J]. Geochimica et Cosmochimica Acta, 1989, 53(1):197-214. doi: 10.1016/0016-7037(89)90286-X
张成江, 陈友良. 510-1铀矿床垂直分带规律的发现及其成因意义[J].地质与勘探, 2010, 46(4):434-441. http://www.cnki.com.cn/Article/CJFDTOTAL-DZKT201003009.htm 赵振华.微量元素地球化学原理[M].北京:科学出版社, 1997. Sun S S, McDonough W F. Chemical and isotopic systematics of oceanic basalts:Implications for mantle composition and processes[C]//Saunders A D, Norry M J. Magmatism in The Ocean Basins. Geological Society Special Publication, 1989, 42:313-345.
孙晓明, 熊德信, 王生伟, 等.云南哀牢山金矿带墨江金镍矿床铂族元素 (PGE) 地球化学及其对矿床成因的制约[J].矿床地质, 2006, 25(4):438-446. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ200604007.htm Rehkamper M, HallidayA N, Barfod D, et al. Platinum-Group Element Abundance Patterns in different Mantle Environments[J]. Science, 1997, 278:1595-1598 doi: 10.1126/science.278.5343.1595
张贵山, 温汉捷, 裘愉卓.闽西晚中生代基性岩脉的地球化学研究[J].地球化学, 2004, 3(3):243-253. http://www.cnki.com.cn/Article/CJFDTOTAL-DQHX200403002.htm 张贵山, 温汉捷, 胡瑞忠, 等.闽东南基性岩脉形成的构造应力场地质意义[J].大地构造与成矿学, 2006, 30(2):142-148. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGKD200504001144.htm 李立兴, 李厚明, 崔艳合, 等.河北高寺台含铬超基性岩杂岩体成岩成矿时代及岩石成因[J].岩石学报, 2012, 28(11):3757-3771. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201211026.htm 冯光英, 刘燊, 冯彩霞, 等.吉林红旗岭超基性岩体的锆石U-Pb年龄、Sr-Nd-Hf同位素特征及岩石成因[J].岩石学报, 2011, 27(6): 1594-1606. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201106003.htm 吴华, 李华芹, 莫新华, 等.新疆哈密白石泉铜镍矿区基性-超基性岩的形成时代及其地质意义[J].地质学报, 2005, 79(4):498-502. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200504009.htm Pearce J A, Norry M J. Petrogenetic implications of Ti, Zr, and Nb variations in volcanic rocks[J]. Contributions to Mineralogy and Petrology, 1979, 69:33-47. doi: 10.1007/BF00375192
Meschede M. A method of discriminat between different types of mid-ocean ridge basalts and continental tholeiites with the Nb-Zr-Y diagram[J]. Chemical Geology, 1986, 56:207-218. doi: 10.1016/0009-2541(86)90004-5
杜远生, 黄宏伟, 黄志强, 等.右江盆地晚古生代-三叠纪盆地转换及其构造意义[J].地质科技情报, 2009, 28(6):10-15. http://www.cnki.com.cn/Article/CJFDTOTAL-DZKQ200906002.htm 曾允孚, 刘文均, 陈洪德, 等.华南右江复合盆地的沉积构造演化[J].地质学报, 1995, 69(2):113-124. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE199502001.htm 曾允孚, 刘文均, 陈洪德, 等.右江复合盆地的沉积特征及其构造演化[J].广西地质, 1992, 5(4):1-14. http://www.cnki.com.cn/Article/CJFDTOTAL-GXDZ199204000.htm 陈洪德, 曾允孚.右江沉积盆地的性质及演化讨论[J].岩相古地理, 1990, 1:28-37 http://www.cnki.com.cn/Article/CJFDTOTAL-TTSD199001003.htm 毛景文, 谢桂青, 李晓峰, 等.华南地区中生代大规模成矿作用与岩石圈多阶段伸展[J].地学前缘, 2004, 11(1):45-55. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200401002.htm Chen M H, Mao J W, Li C, et al. Re-Os isochron ages for arsenopyrite from Carlin-like gold deposits in the Yunnan-Guizhou-Guangxi "golden triangle", southwestern China[J]. Ore Geology Reviews, 2015, 64:316-327. doi: 10.1016/j.oregeorev.2014.07.019
Su W C, Hu R Z, Xia B, et al. Calcite Sm-Nd isochron age of the Shuiyindong Carlin-type Au deposit, Guizhou, China[J]. Chemical Geology, 2009, 258:269-274. doi: 10.1016/j.chemgeo.2008.10.030
Taylor S, McLennan S. The Continental Crust: Its Composition and Evolution[M]. New York:Oxford Press, 1985.
Hofmann A, Jochum K, Seufert M, et al. Nb and Pb in oceanic basalts: New constraints on mantle evolution[J]. Earth and Planetary Science Letters, 1986, 79(1/2):33-45. https://www.researchgate.net/publication/222453665_Hofmann_A_W_Jochum_K_P_Seufert_M_White_W_M_Nb_and_Pb_in_oceanic_basalts_new_constraints_on_mantle_evolution_Earth_Planet_Sci_Lett_79_33-45
Munker C, Pfander J A, Weyer S, et al. Evolution of planetary cores and the earth-moon system from Nb/Ta systematics[J]. Science, 2003, 301(5629):84-87. doi: 10.1126/science.1084662
Rehkamper M, Halliday A N, Barfod D, et al. Platinum-Group Element Abundance Patterns in different Mantle Environments[J]. Science, 1997, 278:1595-1598. doi: 10.1126/science.278.5343.1595
Munker C. Nb/Ta fractionation in a Cambrian arc/back system, New Zealand:source constraints and application of refined ICPMS techniques[J]. Chemical Geology, 1998, 56:207-218. https://www.researchgate.net/publication/248360164_NbTa_fractionation_in_a_Cambrian_arcback_arc_system_New_Zealand_Source_constraints_and_application_of_refined_ICPMS_techniques
Naldrett A J, Hoffmann E L, Green A H, et al. The composition of Ni-sulfide ores with particular reference to their content of PGE and Au[J]. Canada Mineral, 1979, 17:403-415. http://canmin.geoscienceworld.org/content/17/2/403
储雪蕾, 李晓林, 徐久华, 等.汉诺坝玄武岩及其地幔橄榄岩、麻粒岩捕掳体的PGE分布特征[J].科学通报, 1999, 44(8):859-863. http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB199908016.htm 赵正, 漆亮, 黄智龙, 等.攀西裂谷南段鸡街碱性超基性岩微量元素和Sr-Nd同位素地球化学及其成因探讨[J].岩石学报, 2012, 28(6):1915-1927. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201206016.htm Barnes S J, Naldrett A J, Gorton M P. The orgin of the Fractionation of plartinum-group elements in terrestrial magmas[J]. Chemical Geology, 1985, 53:303-323. doi: 10.1016/0009-2541(85)90076-2
Genc S C, Tuysuz O. Tectonic setting of the Jurassic bimodal magmatism in the Sakarya zone (central and western Pontides), northern Turkey:A geochemical and isotopic approach[J]. Lithos, 2010, 118(1/2):95-111. https://www.researchgate.net/publication/222136664_Tectonic_setting_of_the_Jurassic_bimodal_magmatism_in_the_Sakarya_Zone_Central_and_Western_Pontides_Northern_Turkey_A_geochemical_and_isotopic_approach
Eggins S, Rudnick R, McDonough W. The composition of peridotites and their minerals:A laser-ablation ICP-MS study[J]. Earth and Planetary Science Letters, 1998, 154(1/4):53-71. http://www.academia.edu/2806189/The_composition_of_peridotites_and_their_minerals_a_laser-ablation_ICP-MS_study
McDonough W F, Stosch H, Ware N. Distribution of the titanium and the rare earth elements between peridotitic minerals[J]. Contributions to Mineralogy and Petrology, 1992, 110(2):321-328. doi: 10.1007/BF00310747
Xu Y G. Distribution of trace element in spinel and garnet peridotites[J]. Science in China (Series D), 2000, 43(2):166-175. doi: 10.1007/BF02878146
Zack T, Foley S, Jenner G. A consistent partition coefficient set for clinopyroxene, amphibole and garnet from laser ablation microprobe analysis of garnet pyroxenites from Kakanui, New Zealand[J]. Neues Jahrbuch Fur Mineralogie Abhandlungen, 1998, 172:23-41. http://www.academia.edu/2837201/A_consistent_partition_coefficient_set_for_clinopyroxene_amphibole_and_garnet_from_laser_ablation_microprobe_analysis_of_garnet_pyroxenites_from_Kakanui_New_Zealand
Bodinier J, Merlet C, Bedini R, et al. Distribution of niobium, tantalum, and other highly incompatible trace elements in the lithospheric mantle:The spinel paradox[J]. Geochimica et Cosmochimica Acta, 1996, 60(3):545-550. doi: 10.1016/0016-7037(95)00431-9
Keays R R. The role of komatiitic magmatism and S-saturation in the formation of ore deposits[J]. Lithos, 1995, 34:1-18. doi: 10.1016/0024-4937(95)90003-9