矿物岩石地球化学通报  2017, Vol. 36 Issue (2): 213-221   PDF    
造山带碳酸岩起源与深部碳循环
许成 , 曾亮 , 宋文磊 , 冯梦 , 邓淼 , 韦春婉     
北京大学 地球与空间科学学院, 造山带与地壳演化重点实验室, 北京 100871
摘要: 碳酸岩是揭示地幔地球化学动力学的"探针岩石",迄今有关碳酸岩的研究集中在裂谷环境,而鲜见造山带地区碳酸岩研究的报道。在四川攀西喜马拉雅期造山带、秦岭造山带和华北中央造山带内均有碳酸岩产出,且蕴藏了大型稀土和钼矿床,是研究深俯冲、壳-幔作用和深部碳循环的理想实验室。传统观念认为碳酸岩形成于裂谷环境与地幔柱活动相关,而造山带碳酸岩很可能是陆源富碳沉积物俯冲至地幔低程度熔融的产物。中国造山带内碳酸岩地球化学研究均显示了地壳物质对碳酸岩地幔源区的贡献,这暗示地表碳俯冲至深部地幔,交代地幔发生熔融,这不仅为较还原的地幔源区提供富氧成分,还可以使碳酸岩的母体岩浆更富集稀土,形成稀土矿床。
关键词: 碳酸岩      沉积物俯冲      碳循环      造山带     
Orogenic Carbonatite Petrogenesis and Deep Carbon Recycle
XU Cheng, ZENG Liang, SONG Wen-lei, FENG Meng, DENG Miao, WEI Chun-wan     
Laboratory of Orogenic Belts and Crustal Evolution, School of Earth and Space Sciences, Peking University, Beijing 100871, China
Abstract: Carbonatites play a critical role in understanding mantle geochemistry dynamic processes. Previous studies focus mainly on carbonatite in rifting settings, however, carbonatites emplaced in orogenic belts have rarely been reported. The rocks from Qinling, Panxi, and the North China central orogenic belt in China host large REE and Mo deposits that can serve as a perfect "laboratory" for studies of subduction, crust-mantle interaction, and Earth's deep carbon recycle. It is well known that carbonatites in rifting environments are generally associated with mantle plume activities. In Contrast, orogenic carbonatites are considered to be the products of low-degree partical melting of metasomatized mantle by subducted carbon-rich sediments. Geochemistry of carbonatites in the orogenic belts in China have shown the involvement of crustal materials into the mantle. Subduction of crustal carbon can not only provide oxygen-rich compounds to the mantle, but also result in high REE in carbonatite parental magmas, thus forming REE deposits.
Key words: carbonatites     subduction     carbon recycle     orogenic belt    

碳酸岩是主要由方解石、白云石、铁白云石等碳酸盐矿物和少量硅酸盐矿物组成的岩浆岩,碳酸盐矿物比例大于50%,SiO2含量通常小于20% (Le Maitre et al., 2002) 。碳酸岩的研究历史相对较短,最早由Bose (1884) 提出,Högbohm (1895) 首先对瑞典Alnö碳酸岩进行了研究,但是否存在岩浆成因的碳酸岩,一直存在很大争论,直到20世纪60年代初在非洲坦桑尼亚东部发现了正在喷溢的碱性 (钠质) 碳酸岩火山 (Dawson,1962) ,才开始广泛被接受与认同。随着研究的不断深入,发现碳酸岩起源地幔,而且与稀土、铌等稀有战略矿产资源关系密切,因而越来越受到重视 (宋文磊等,2013) 。碳酸岩产出小且少,稳定性较差,容易被风化、剥蚀,而且很容易与其他类型的碳酸盐岩混淆。据Woolley和Kjarsgaard (2008) 的不完全统计,截至于2008年,全世界报道的碳酸岩仅有527处,其中喷发相型47处,独立产出的岩体 (无共生/伴生碱性硅酸岩) 110余处。尽管碳酸岩的产出较少,但其分布非常广泛,遍布了全世界七大洲 (图 1) ,非洲西北近海地带的Canary Islands、Cape Verde等地还发现了几处洋岛背景产出的碳酸岩 (Hoernle et al., 2002) 。

图 1 全球碳酸岩分布图 (引自Woolley和Kjarsgaard,2008) 和中国造山带碳酸岩产出位置 Figure 1 Distribution of carbonatite of the world (modified from Woolley and Kjarsgaard, 2008)

从形成时代上看,从太古代至今都有碳酸岩的产出,其中最古老的碳酸岩距今已约30亿年,位于格陵兰西部的Tupertalik (Bizzarro et al., 2002) ,而最年轻的当属非洲坦桑尼亚的OlDoinyo Lengai碱质火山碳酸岩,20世纪60年代首次观察到其喷溢活动之后,2007年后又发现了一次新的喷溢记录,规模比较小,但活动时间较长,一直到2012年还有活动迹象 (Bosshard-Stadlin et al., 2014) 。碳酸岩的产出有2个较为明显的高峰和低谷,高峰期主要集中于元古代 (古元古代末-新元古代1800~650 Ma) 和白垩纪-新生代。

目前已知的大多数碳酸岩主要产于裂谷背景的板内环境,如德干高原、加拿大地盾 (东南部) 、东非大裂谷系和西伯利亚地盾等地区,其成因被认为与地幔柱活动相关 (Ernst and Bell, 2010) 。近年来,在一些造山带均发现了碳酸岩,这对研究深部碳循环 (地球表层-地幔系统之间) 具有重要价值。碳在地球深部属于不相容元素,含量非常少,它在地球深部的储存形式以及起源过程是目前关注的焦点。随着环境和气候的日益恶化,人们益发关注地球各圈层之间物质和能量的交换,以及地球内部碳 (主要是以CO2的方式) 的释放对地表环境、气候变化的影响,关于地球表层的碳循环 (气圈-水圈-生物圈-土壤圈) 已有大量的研究成果,但由于条件所限,对地球深部的了解还相对较少,大部分认识主要基于高温高压实验工作。因此,加强对碳在地球内部的分布、存储和循环的研究,不仅对于认识地球内部碳的来源,而且对全球规模的碳循环作用的认识,都具有十分重要的作用和意义 (Hirschmann et al., 2009;Dasgupta and Hirschmann, 2010) ,而造山带碳酸岩的研究,是一个重要的突破口。

1 造山带碳酸岩的分布

目前在芬兰、加拿大、巴基斯坦和中国的造山带内均有碳酸岩产出。芬兰楠塔利 (Naantali) 碳酸岩位于波罗的海地台 (Woodard and Hetherington, 2014) ,该区域分布有1.92~1.88 Ga的岛弧火山岩,1.83~1.81 Ga的高温麻粒岩以及由泥质混合岩、玄武质和安山质钙碱性变质火山岩、英云闪长岩、重熔花岗岩组成的混杂岩,显示波罗的海地台经历了早元古代的俯冲和造山过程。楠塔利碳酸岩呈脉体产出,与钾质煌斑岩共生,它们主要由方解石组成,含少量磷灰石和副矿物褐帘石、榍石和透辉石。碳酸岩的锆石和独居石U-Th-Pb年龄为1.8 Ga,显示它们为造山带的产物。加拿大马尼托巴碳酸岩位于北美地台的元古代哈德森中央造山带 (Trans-Hudson Orogen) (Chakhmouradian et al., 2008) 。该造山带是太古代三大板块 (Hearne,Superior和Sask克拉通) 拼合的产物,它们主要由1.9 Ga的变质火山岩+变质沉积岩组成的绿岩带组成。绿岩带被1.8 Ga的辉长岩和同碰撞花岗岩侵入。马尼托巴碳酸岩呈脉体,侵入正长岩内,它们主要由方解石组成,含少量磷灰石和单斜辉石,副矿物包括褐帘石、榍石、锆石、角闪石、金云母、磁铁矿、重晶石和萤石。马尼托巴碳酸岩的成岩年龄并没有报道,但它们具有低的δ13C组成,Chakhmouradian等 (2008) 提出来源于地壳俯冲的物质。在巴基斯坦的喜马拉雅造山带也有碳酸岩的产出 (Tilton et al., 1998) ,主要分布在Loe Shilman和Sillai Patti地区,它们呈脉体产出,主要由方解石组成,黑云母的Rb-Sr年龄为31 Ma。Loe Shilman碳酸岩与钾质岩和煌斑岩密切共生,而Sillai Patti呈独立脉体产出。

中国造山带发育,也蕴藏了大量的碳酸岩,包括四川攀西喜马拉雅期造山带、秦岭造山带和华北中央造山带 (图 1) 。四川攀西地区北起牦牛坪,向南经里庄至大陆乡,纵贯冕宁、西昌、德昌3县市,长150 km以上的碳酸岩带。碳酸岩主要呈脉体或透镜体产出,与石英正长岩密切共生。牦牛坪碳酸岩产状总体与北东向断裂一致,主脉两侧有细脉侵入正长岩。脉宽90~200 m,沿倾向延伸400 m未尖灭。岩脉主要由方解石组成,次为黑云母、霓石、霓辉石、钠铁闪石和正长石等,常见副矿物有磷灰石和榍石等。岩脉边缘发育霓长石化和钠铁闪石化。蒲广平 (2001) 获碳酸岩脉中镁钠铁闪石K-Ar为32 Ma,与区内正长岩年代 (K-Ar稀释法年龄为28~40 Ma) 相近。里庄碳酸岩规模较小,呈透镜体产于正长岩内,其矿物组成与牦牛坪相似。碳酸岩中辉钼矿Re-Os模式年龄为28 Ma、黑云母Ar-Ar年龄为31 Ma (田世洪等,2008) ,与正长岩的Ar-Ar年龄 (27 Ma) 相近。大陆乡碳酸岩沿走滑断裂形成的构造裂隙侵入,岩体周围角砾岩筒发育,其矿物组成基本与牦牛坪相似,但出现天青石和石英斑晶,其黑云母Ar-Ar年龄为13.2 Ma (田世洪等,2008) 。四川攀西地区位于扬子克拉通的西缘,从元古代至新生代经历了复杂的构造演化 (Hou et al., 2009):元古代原特提斯洋板块向扬子板块下俯冲;古生代的古特提斯洋板块在早二叠世向西的俯冲,导致了扬子板块西缘转变为被动大陆边缘;中二叠世的地慢柱活动,在扬子克拉通的西缘形成了近南北展布的攀西古裂谷;60~45 Ma印度-亚洲大陆板块的大规模碰撞,扬子克拉通西缘进入喜马拉雅期碰撞造山作用,形成了位于印度-亚洲大陆碰撞带东缘的碰撞造山带,发育了一系列新生代走滑断裂系统,因此攀西喜山期碳酸岩的成因与地壳物质的俯冲循环存在密切联系,它们也蕴藏了中国第二大的轻稀土成矿带,稀土矿物主要以氟碳铈矿为主,与萤石和方解石形成稀土矿脉 (Xu et al., 2008, 2012Hou et al., 2009Liu et al., 2015Xie et al., 2015) 。

秦岭造山带是华北与扬子两大地质构造单元的分界和衔接地,在中国大陆地质构造演化中占有重要地位。前人对秦岭造山带进行了大量深入的地质研究,取得了很多重要的研究成果,确定了扬子板块与秦岭微板块和华北板块之间的碰撞时限以及中新生代陆内造山作用与构造演化等 (张国伟等,2001张本仁等,2002路凤香等,2006) 。在志留-泥盆纪,扬子板块和华北板块之间出现秦岭古洋盆,古秦岭洋向北俯冲碰撞,形成了商丹缝合线。在古秦岭洋向北俯冲碰撞的同时,扬子板块北缘演化出独立的秦岭微板块。三叠纪,扬子板块顺时针转动,与华北板块拼合碰撞,导致秦岭洋向北俯冲闭合以及强烈的陆内造山活动。秦岭造山带从南到北依次划分为南秦岭、北秦岭和小秦岭,其中南北端元均有碳酸岩的产出,包括南秦岭的湖北庙垭和小秦岭的陕西华阳川、黄龙铺等地 (Xu et al., 2007, 2010a) 。湖北庙垭碳酸岩位于南秦岭武当地块西南缘,呈岩株状侵入元古界耀岭河群中。碳酸岩与正长岩在空间上共生,主要由方解石组成,含有少量黑云母、霓石、磷灰石和长石等矿物。Xu等 (2014) 获得碳酸岩内独居石的U-Pb年龄为232 Ma。陕西华阴至洛南一线的华阳川、垣头、黄龙铺等地的碳酸岩产于秦岭造山带北缘的小秦岭,碳酸岩主要呈脉状,侵入太古界太华群和中元古界熊耳群的变质岩中,未见同时代的碱性岩出露。碳酸岩主要属方解石碳酸岩,普遍含有石英、霓石、钠铁闪石、金云母、天青石、磷灰石和长石等。Stein等 (1997) Song等 (2015) 获得黄龙铺和垣头碳酸岩内辉钼矿的年龄分别为221 Ma和226 Ma。尽管南北秦岭造山带的碳酸岩为三叠纪造山带的产物,但其产状和成矿物质明显不同。南秦岭庙垭碳酸岩呈岩株状,富含原生的稀土矿物,包括磷灰石、独居石和氟碳铈矿等,是大型轻稀土矿床。小秦岭碳酸岩内富含辉钼矿,是中国大型钼矿床,也是世界上唯一一处碳酸岩型钼矿床 (Xu et al., 2010a) 。

华北板块是世界上最古老的 ( > 38亿年) 克拉通之一,它经历了复杂的地质变迁,在太古宙末 (约25亿年) 形成东、西板块,在古元古代 (19~18.5亿年) 东、西板块沿着中央造山带碰撞拼合 (Zhao et al., 2005) ,最终完成了克拉通化。然而华北克拉通却表现出与世界其他克拉通十分不同的特征:从18亿年至古生代一直保持相对稳定,并存在巨厚的太古宙岩石圈根 (包括下地壳和岩石圈地幔) ,但构造地质学、地幔包体和岩浆岩的岩石学和地球化学以及地球物理 (包括地震和地热) 研究均证明 (吴福元等,2008) ,华北克拉通自中生代以来遭受了强烈的活化改造,发生了大规模的构造变形和岩浆活动,形成了多种类型的盆地,并伴随产生了大量金属矿产和油气资源,使原有的克拉通结构和性质遭到明显的破坏,被认为是世界上古克拉通被移离的最佳实例。因此,华北克拉通破坏 (或岩石圈减薄) 已成为继青藏高原和大别-苏鲁超高压变质带之后,中国又一在国际地球动力学界引起广泛关注的重大科学问题,已成为国际地球科学前沿研究领域。内蒙古丰镇和河北右所堡碳酸岩位于中央造山带内,两地碳酸岩均与辉石岩密切共生,呈脉状产出,富含磷灰石,已达到小型磷矿床。丰镇地区的围岩主要为石榴斜长片麻岩、石英岩、石榴黑云片麻岩、石墨片麻岩、大理岩和钙硅质片麻岩。镜下观察显示,该碳酸岩主要由方解石组成,次要矿物和副矿物主要包括辉石、磷灰石、金云母、尖晶石、钙钛矿、钛铁矿和稀土矿物。右所堡出露围岩主要为英云闪长质片麻岩、片麻岩、辉长质片麻岩和斜长角闪岩,其碳酸岩矿物组合与丰镇相似。

2 碳酸岩的成因争论

碳酸岩通常以小的岩脉、岩墙、岩株、岩锥产出,与金伯利岩、煌斑岩、辉长岩、辉石岩、霓霞岩、霞石岩、粗面岩、响岩、黄长岩和正长岩等碱性硅酸岩构成环状杂岩体,也有独立产出的碳酸岩。大量同位素地球化学和高温高压实验证实碳酸岩起源于地幔,其岩浆平衡压力大于3 GPa (Wyllie and Huang, 1975) 。与硅酸盐熔体相比,碳酸岩熔体相对富含挥发分,而固液温度、黏度 (1.5×10-2~5×10-3 Ps) 和密度 (2.2 g/cm3) 均较低 (Dobson et al., 1996) ,极易流动,富含不相容元素,扩散速率很大,在其上升穿越相对较“冷”的岩石圈地幔时,可与地幔二辉橄榄岩等进行反应,促使地幔交代作用发生,在改变地幔矿物学特征和化学成分的同时,释放富含CO2的流体 (Dasgupta et al., 2009) 。上世纪80年代以来,碳酸岩的研究快速发展,特别是在同位素和实验岩石学方面取得了很多重要进展和新认识: (1) 碳酸岩普遍富含Sr和Nd,可以缓冲地壳物质混染的影响,而保留和记录地幔源区的Sr、Nd同位素组成 (Bell,1998) ; (2) 碳酸岩熔体是引起地幔交代作用和地幔地球化学不均一性的重要介质之一,0.1%碳酸岩溶体的交代作用,将能导致地幔源区30%~60%不相容元素的活化迁移 (Dasgupta et al., 2009) ; (3) 碳酸岩的成因与地幔热柱活动有关 (Marty et al., 1998) ; (4) 碳酸岩源区形成于HIMU-EM1地幔端元的混合作用 (Bell and Simonetti, 2010) 。然而碳酸岩与共生碱性岩的成因联系一直是争论的焦点,包括:直接来源于地幔橄榄岩的低程度部分熔融作用,与碱性岩并无成因联系 (Wallace and Green, 1988) ;富CO2的硅酸岩岩浆经过液态不混溶 (Kjarsgaard and Hamilton, 1989) 或分离结晶作用 (Lee and Wyllie, 1998) 形成碳酸岩和共生的碱性岩。

相对裂谷碳酸岩的研究,造山带碳酸岩与碱性岩的成因关系研究较少。Chakhmouradian等 (2008) 发现北美板块马尼托巴碳酸岩的Ba/Mn、Ba/La、Nb/Th值低于共生的正长岩,这与液态不混溶实验比值相反,因此,提出两类岩石并无成因联系。尽管许成等 (Xu et al., 2003) 提出攀西牦牛坪碳酸岩与正长岩有相似的年龄和Sr-Nd同位素组成,两者为液态不混溶的产物。但整个攀西地区碳酸岩都富集REE,但液态不混溶实验显示REE优先进入硅酸岩熔体相内,碳酸岩熔体并不富集REE (Veksler et al., 1998) 。值得注意的是,地表侵入的碳酸岩都是经历了不同程度的分离结晶作用,其化学组成并不能代表母体岩浆 (Xu et al., 2007) ,而实验结果来源于初始熔体。南秦岭庙垭碳酸岩是目前少有的呈岩株状侵入的岩体,是研究造山带碳酸岩与共生碱性岩成因联系的理想实验室。我们利用SHRIMP和LA-ICPMS分别确定了南秦岭湖北庙垭碳酸岩里的独居石U-Pb年龄234 Ma,与碳酸岩共生的正长岩内发现了两类岩浆锆石,其中相对粗晶的锆石U-Pb年龄766 Ma,细粒岩浆锆石年龄U-Pb年龄147 Ma (Xu et al., 2014, 2015) 。在同一岩体内,发现了多种岩浆锆石的年龄,显示锆石来源的多源性,很难代表真实的岩体就位年龄。镜下观察发现,碳酸岩和正长岩内均有许多石英-长石细脉侵入。这也暗示可能外来锆石的混入。而碳酸岩是贫硅岩浆,其是否能结晶出锆石一直存在争论 (Le Bas,2006) ,扫描电镜观察也未在碳酸岩内发现初始的锆石 (Xu et al., 2015) ,而独居石是该岩体内早期结晶矿物,也是定年的理想对象。

液态不混溶和分离结晶是碳酸岩与共生碱性岩的主要成因模式。庙垭碳酸岩富含原生的稀土矿物,它们的Pb含量也与正长岩相近,而液态不混熔实验显示REE,Pb优先配分入硅酸岩岩浆 (Veksler et al., 1998, 2012Martin et al., 2013) ,这与庙垭碳酸岩富含REE和Pb的特征不相符。大量长石的分离结晶,将导致残余熔体具有低Ga/Rb和Cr、V、Ni含量 (Zajacz and Halter, 2007) ,而庙垭碳酸岩的Ga/Rb值 (0.7~8.8) 高于正长岩 (0.2~0.6) ,两者的Cr、V、Ni含量也相近,因此并不能用分离结晶模式来解释两者的成因联系。此外,庙垭碳酸岩的初始 (87Sr/86Sr)i为0.7036~0.7040,εNd为0.6~-1.1,而正长岩的初始 (87Sr/86Sr)i为0.701~0.7034,εNd为9.4~10.1。同位素特征也不支持两者有成因联系。笔者等在碳酸岩和正长岩内均发现有初始的磷灰石 (Xu et al., 2015) 。磷灰石富含不相容元素REE、Th、U、Sr和Ba等,它们在碳酸岩和硅酸岩岩浆体系的配分系数已有报道。因此,根据物质平衡和瑞利分离结晶公式,通过磷灰石的化学组成,可以计算它们的母体岩浆成分:

式中:C0CS分别为初始岩浆和结晶固相的化学成分;DS为固相在熔体中的配分系数;F为分离出的熔体比例。

基于上述公式的模拟结果显示,庙垭碳酸岩质岩浆比正长岩岩浆具有高的Sr和REE,低的Ba、Th、U、Nb和Ta含量 (图 2) 。初始碳酸岩熔体具有强的轻稀土富集,缺少Eu异常,而正长岩母体具有相对低的轻稀土富集程度和明显的Eu负异常。在长石/硅酸岩体系,Sr和Eu是相容元素,而Th和U是不相容元素,大量长石的分离结晶,将导致残余富碳酸质的硅酸岩岩浆亏损Sr和Eu,富集Th和U,从而使最后分离的碳酸岩岩浆继承上述4个元素的组成特征,这与模拟的结果并不相同。模拟的碳酸岩比正长岩岩浆具有更高的稀土组成,这与液态不混溶实验明显不同 (图 2) 。因此,尽管庙垭碳酸岩与正长岩在空间上共生,但两者并无成因联系。磷灰石是碳酸岩和碱性岩内最常见的副矿物,它们富含不相容元素,其微区成分的研究工作对揭示碳酸岩的成因具有重要意义。

庙垭碳酸岩和正长岩初始母体组成由早期结晶的磷灰石模拟计算 (Xu等,2015) 。液态不混溶实验引自Veksler等 (1998) 图 2秦岭造山带庙垭碳酸岩和正长岩初始熔体化学组成比值与液态不混溶实验结果对比 图 2 秦岭造山带庙垭碳酸岩和正长岩初始熔体化学组成比值与液态不混溶实验结果对比 Figure 2 Comparison of initial compositions of the Miaoya carbonatite and syenite with experimental results of liquid immiscibility
3 碳酸岩母体岩浆与深部碳循环

近年来,人们意识到地球深部碳循环长期影响着地表碳的演化和分布。地球内部碳的分布、储存和循环对于认识全球碳循环具有关键的作用。俯冲的板片是地球深部碳循环的主要机制。Dasgupta等 (2013) 通过岩石学研究估算,提出每年有5.4×1013~8.8×1013 g地表碳以沉积作用进入新生洋壳。沉积物和大洋岩石圈中的含碳矿物由俯冲作用 (Hirschmann et al., 2009;Dasgupta and Hirschmann, 2010; Dasgupta et al., 2013) 带入地球深部,随着俯冲的深入,板片被加压加热,板片内的挥发分元素通过脱水和脱碳作用转化为流体相 (主要为H2O和CO2) ,从板片分离,向上迁移并交代上覆地幔楔,导致其熔点降低发生部分熔融形成弧岩浆,弧岩浆在结晶和喷发过程中将其所含的来自板片的碳以CO2的形式返还到大气中 (Burton et al., 2013) 。研究表明最多有40%的碳通过弧岩浆作用返回到地表 (Zhang and Zindler, 1993Dasgupta and Hirschmann, 20062010) ,这意味着有相当一部分碳酸盐矿物在在弧下并不分解而是随板片被带入到地幔的更深处,并通过碳酸岩岩浆作用再返回地表。

大量对全球碳酸岩的稳定和放射性同位素示踪 (Nelson et al., 1988Xu et al., 2014) 结果表明,俯冲循环的地壳物质对碳酸岩的成因起着非常重要的作用。碳酸盐化榴辉岩体系的实验岩石学研究也证实碳酸岩可由俯冲板片直接熔融形成 (Dasgupta et al., 2004Yaxley and Brey, 2004Litasov and Dasgupta et al., 2010) 。另外,近年来对碳酸盐化泥质岩的部分熔融实验也生成了碳酸岩熔体 (Thomsen and Schmidt, 2008Grassi and Schmidt, 2011) 。相对于碳酸盐化橄榄岩和榴辉岩,碳酸盐化的泥质岩形成的碳酸岩质熔体更加富集CO2、碱质 (主要为K) 及不相容元素,并在所有主要的俯冲岩石类型中具有最低的熔点,也解释了一些含钾质幔源岩浆的成因。以上实验研究表明,随着压力的升高,板片中稳定的碳酸盐矿物依次为方解石-白云石-文石-菱镁矿,在大于5 GPa时,白云石转化为菱镁矿,在更高的压力下,菱镁矿始终是熔融残余物中唯一稳定的碳酸盐矿物。实验发现菱镁矿甚至可以稳定存在于下地幔的超高压力和温度条件下 (100 GPa和2700℃) ,在压力约为15 GPa时菱镁矿转变为高压多晶相和菱镁矿Ⅱ (Isshiki et al., 2004) 。然而,最新的实验研究表明 (Thomson et al., 2016) ,在地幔转换带处 (700 km) 板片会发生熔融形成碳酸岩熔体,丢失几乎所有的碳酸盐矿物,使其不能循环到更深的地幔去。

值得注意的是,碳酸盐化橄榄岩体系的熔融均生成了富镁的熔体 (Dalton and Presnall, 1998Gudfinnsson and Presnall, 2005) ,并不能熔出方解石质的碳酸岩熔体,这与地表难测到的碳酸岩绝大多数以方解石碳酸岩为主并不一致。而碳酸盐化榴辉岩体系熔出的碳酸岩熔体含有更多的CaCO3 (50%~80%) ,这与其实验初始物质组成有关,且实验熔体中Mg的含量随着压力的升高而升高 (Dasgupta et al., 2005) 。但在14~20 GPa,碳酸岩熔体又倾向于方解石质,很可能是在这个压力下超硅石榴子石的稳定存在,其相对于Ca会更富集Mg,使与之平衡的熔体富Ca而Mg含量降低 (Keshav et al., 2011) 。这些实验结果表明钙质的碳酸盐化榴辉岩可能是许多地表方解石质碳酸岩的源区 (Hammouda et al., 2015) 。然而是否碳酸岩的成因与俯冲物质有成因联系尚无定论。Bell和Simonetti (2010) 根据全球显生宙以来碳酸岩的同位素组成特征认为,碳酸岩的地幔源区主要是富集地幔EM1和高U-Pb的HIMU端元混合而成,缺少亏损地幔DMM端元;一些碳酸岩显示地壳特征的Sr-Nd-Pb同位素组成,来源于地幔源区长时间的演化 ( > 3 Ga) ,形成高的Rb/Sr、Sm/Nd、U/Pb比值,因此演化成高的Sr-Nd-Pb同位素比值。他们认为缺少造山带碳酸岩的出露是排除地壳物质俯冲与碳酸岩起源的关键因素。攀西地区喜马拉雅期造山带、秦岭造山带和华北中央造山带内碳酸的发现对研究地壳碳酸盐物质的深俯冲与碳酸岩的成因具要重要价值。我们确定了南秦岭庙垭碳酸岩的成岩年龄为晚三叠世的产物,它们的侵位靠近勉略缝合线,在空间和时间上与勉略洋壳的俯冲密切相联,并且庙垭碳酸岩的C同位素高出正常地幔值,它们具有球粒陨石库的Nd同位素组成,其CHUR模式年龄 (190~300 Ma) 与其成岩年龄相近,但Pb同位素变化较大,部位样品显示富集的特征,这些均指示庙垭碳酸岩是俯冲洋壳低程度熔融的产物 (Xu et al., 2014) 。

除碳酸岩外,金刚石也是了解地球深部碳循环和碳酸岩起源的有力工具。对全球发现的金刚石研究发现,它们通常含有两类子矿物:超硅石榴子石、透辉石和Ca-Si钙钛矿、Mg-Si钙钛矿、铁方镁石,暗示金刚石主要来源于玄武质和橄榄质母体 (Harte,2010) ,前者起源深度为300~500 km,超基性母体的金刚石起源深度为600~800 km。许多金刚石内超硅石榴子石的稀土配分模式出现了Eu的负异常 (Tappert et al., 2005) ,并发现了富铝的子矿物,显示地壳俯冲物质与地球深部金刚石成因的联系。地震学研究证实俯冲洋壳可以进入下地幔深度 (Fukao et al., 2009) ,Walter等 (2011) 发现巴西金伯利岩筒中的金刚石具有异常负的碳同位素组成,提出了地表碳可随俯冲洋壳进入下地幔深度 ( > 900 km) 。而地表沉积碳酸盐俯冲进入地幔的形式是碳酸岩形成机理的关键。高温高压实验研究发现辉石和石榴子石在7 GPa以上时可以容纳大量的Fe3+,因此上地幔大于250 km时Fe、Ni等金属可以稳定存在 (Rohrbach et al., 2007) ,深部地幔相较还原,Stagno和Frost (2010) 的高温高压实验表明碳酸盐矿物在100~150 km时就开始被还原成金刚石或石墨。Rohrbach和Schmidt (2011) 的实验也认为俯冲的碳酸盐进入受Fe、Ni金属控制氧逸度的上地幔时 (250 km) 会被大量还原成金刚石。然而,Walter等 (2008) 在金刚石内也发现了碳酸质流体包裹体,提出碳酸岩岩浆可以在地幔过渡层甚至下地幔顶部形成。因此,金刚石在地幔的氧化熔融过程被看作碳酸岩岩浆形成的关键,而这一过程需要地幔柱的上升动力,使深部地幔含Fe3+硅酸盐矿物释放出Fe3+,氧化单质碳为碳酸盐,从而高温熔融。Rohrbach和Schmidt (2011) 提出了形成碳酸岩熔体的2个过幔转化带,如果地幔柱起源于下地幔,它们携带金刚石和富Fe3+的MgSi-钙钛矿上升进入660 km的转化带时,将发生相变,释放出Fe3+,氧化金刚石为碳酸岩熔体;如果地幔柱起源于地幔过渡层,金刚石和富Fe3+的超硅石榴子石上升进入250 km的软流圈,将发生相变转化为普通石榴子石,释放出Fe3+,进而氧化金刚石。值得注意的是,这种模式仅适用于解释板内裂谷型碳酸岩的成因,而造山带碳酸岩熔体能否在150 km以下形成还不清楚。

4 俯冲沉积物与稀土富集成矿

当今世界每6项新技术的发明,就有一项离不开稀土元素。它是现代工业的“维生素”,由于具有优良的光电磁等物理特性,能与其他材料组成性能各异、品种繁多的新型材料,大幅度提高产品的质量和性能,稀土已成为当今重要的战略资源。但稀土矿床在世界上的分布极不均匀,主要集中在中国、美国、俄罗斯、澳大利亚、巴西、印度、马来西亚、加拿大、南非、埃及等少数国家。中国是一个名副其实的稀土资源大国,约占世界储量的36% (Humphries,2010) 。由于稀土矿床在全球的分布非常有限,目前已知资源量已经远远不能满足未来全球范围内高新技术产业发展对稀土的需求 (Hatch,2012) ,因此,许多国家把REE资源勘探开发列为重点项目。碳酸岩是所有火成岩中REE含量最高的,但仅少数碳酸岩蕴藏具有经济价值的稀土矿床。据Woolley和Kjarsgaard (2008) 不完全统计,全球已知具有经济价值的碳酸岩型内生稀土矿床少于10处,而磷矿床32处,铌矿床33处,而攀西喜马拉雅期造山带碳酸岩稀土成矿带具有世界上第3大稀土储量,南秦岭庙垭碳酸岩是大型稀土矿床。

尽管碳酸岩稀土含量高,但大多稀土富集在碳酸盐矿物内 (Ingrid,1998) 而缺少稀土矿物。许多实验也证实碳质橄榄岩熔融形成的碳酸质熔体稀土含量低,岩石学证据也显示地幔包体内的碳酸盐矿物含有低的稀土元素 (Ionov and Harmer, 2002Foley et al., 2009) ,因此大多数碳酸岩的稀土来源于地壳内的分离结晶。而造山带碳酸岩的成因与俯冲的沉积物质密切相关。研究发现东印度洋和太平洋3500~6000 m深海沉积物中含有大量的稀土 (Kato et al., 2011Yasukawa et al., 2014) 。Hou等 (2006) 提出在30 Ma以来,攀西地区出现由挤压向转换拉张的过渡阶段,深俯冲至地幔的洋壳沉积物部分熔融作用形成该区域特有的碳酸岩及其大规模稀土成矿。Ling等 (2013) 也提出内蒙古白云鄂博超大型稀土矿床主要来源于古亚洲洋俯冲分异的富稀土的碳酸质流体的交代作用。Xu等 (2010a, 2015) 也发现南秦岭庙垭碳酸岩也具有高的REE组成和LREE富集,但岩体内的方解石具有低的REE含量和平坦的REE配分模式,明显不同于大多数碳酸岩内的碳酸盐矿物。镜下观察发现内岩体内富含初始岩浆成因的富REE的磷灰石和独居石,它们早于方解石结晶。这说明庙垭碳酸岩的初始母体富含稀土,优先结晶稀土矿物,导致主要矿物方解石相对贫稀土,而如果母体内稀土含量相对较少,稀土将优先进入碳酸盐矿物内,并不发育稀土矿物,这也得到了实验证实 (Wyllie et al., 1996) 。

5 结论

碳酸岩是地表出露相对较少的幔源岩石之一,即可形成于裂谷环境,也可以产生于碰撞造山过程中派生的拉张岩石圈断裂带,是研究地幔构造背景和地幔交代作用的“探针岩石”。中国造山带内分布了大量的碳酸岩,也是重要的稀土成矿带,是研究地壳碳酸盐岩俯冲、循环以及与稀土富集联系的天然实验室。尽管碳酸岩在同位素地球化学和高温高压实验方面取得了重要进展,但不相容元素,特别是稀土元素的富集机制还有待深入。多数碳酸岩主要由碳酸盐矿物组成,大量的镁铁质矿物在碳酸岩侵入地表前已经分离结晶,围岩接触带特有的霓长岩化特征也指示挥发分相的分离,因此碳酸岩的全岩化学成分很难代表其母体岩浆的组成。碳酸岩的稀土元素富集过程与起源条件,源区组成以及俯冲的沉积物质是否相关,还不清楚。地幔柱已成功用于解释裂谷碳酸岩的成因,但造山带内碳酸岩的成因并不与地幔柱相关,可能直接与俯冲的碳酸盐矿物在地幔源区发生熔融相关,而俯冲的碳酸岩盐矿物不可避免的带来高氧化的成分,将有可能抑制碳酸盐在深部地幔向单质碳的转化,进而直接发生熔融,在造山带形成碳酸岩岩浆。高氧逸度下的高温高压实验和碳酸岩地幔包体的工作将可能是解释造山带碳酸岩成因,以及深部碳循环和存在形式的关键。

参考文献
[] Bell K. 1998. Radiogenic isotope constraints on relationships between carbonatites and associated silicate rocks-a brief review. Journal of Petrology, 39(11-12): 1987–1996. DOI:10.1093/petroj/39.11-12.1987
[] Bell K, Simonetti A. 2010. Source of parental melts to carbonatites-critical isotopic constraints. Mineralogy and Petrology, 98: 77–89. DOI:10.1007/s00710-009-0059-0
[] Bizzarro M, Simonetti A, Stevenson R K, David J. 2002. Hf isotope evidence for a hidden mantle reservoir. Geology, 30(9): 771–774. DOI:10.1130/0091-7613(2002)030<0771:HIEFAH>2.0.CO;2
[] Bose P N. 1884. Geology of the lower Narbada valley between nimáwar and káwant. Geological Survey of India, 21: 1–72.
[] Bosshard-Stadlin S A, Mattsson H B, Keller J R. 2014. Magma mixing and forced exsolution of CO2 during the explosive 2007-2008 eruption of Oldoinyo Lengai (Tanzania) . Journal of Volcanology and Geothermal Research, 285: 229–246. DOI:10.1016/j.jvolgeores.2014.08.017
[] Burton M R, Sawyer G M, Granieri D. 2013. Deep carbon emissions from volcanoes. Reviews in Mineralogy and Geochemistry, 75: 323–354. DOI:10.2138/rmg.2013.75.11
[] Chakhmouradian A R, Mumin A H, Demény A, Elliott B. 2008. Postorogenic carbonatites at Eden Lake, Trans-Hudson Orogen (northern Manitoba, Canada) . Geological setting, mineralogy and geochemistry. Lithos, 103 (3-4): 503-526
[] Dalton J A, Presnall D C. 1998. The continuum of primary carbonatitic-kimberlitic melt compositions in equilibrium with lherzolite:Data from the system CaO-MgO-Al2O3-SiO2-CO2 at 6 GPa. Journal of Petrology, 39(11-12): 1953–1964. DOI:10.1093/petroj/39.11-12.1953
[] Dasgupta R, Hirschmann M M, Dellas N. 2005. The effect of bulk composition on the solidus of carbonated eclogite from partial melting experiments at 3 GPa. Contributions to Mineralogy and Petrology, 149(3): 288–305. DOI:10.1007/s00410-004-0649-0
[] Dasgupta R, Hirschmann M M, Withers A C. 2004. Deep global cycling of carbon constrained by the solidus of anhydrous, carbonated eclogite under upper mantle conditions. Earth and Planetary Science Letters, 227(1-2): 73–85. DOI:10.1016/j.epsl.2004.08.004
[] Dasgupta R, Hirschmann M M. 2006. Melting in the Earth's deep upper mantle caused by carbon dioxide. Nature, 440(7084): 659–662. DOI:10.1038/nature04612
[] Dasgupta R, Hirschmann M M, McDonough W F, Spiegelman M, Withers A C. 2009. Trace element partitioning between garnet lherzolite and carbonatite at 6.6 and 8.6 GPa with applications to the geochemistry of the mantle and of mantle-derived melts. Chemical Geology, 262: 57–77. DOI:10.1016/j.chemgeo.2009.02.004
[] Dasgupta R, Hirschmann M M. 2010. The deep carbon cycle and melting in Earth's interior. Earth and Planetary Science Letters, 298(1-2): 1–13. DOI:10.1016/j.epsl.2010.06.039
[] Dasgupta R, Mallik A, Tsuno K, Withers A C, Hirth G, Hirschmann M M. 2013. Carbon-dioxide-rich silicate melt in the Earth's upper mantle. Nature, 493(7431): 211–215. DOI:10.1038/nature11731
[] Dawson J B. 1962. Sodium carbonate lavas from Oldoinyo Lengai, Tanganyika. Nature, 195(4846): 1075–1076. DOI:10.1038/1951075a0
[] Dobson D P, Jones A P, Rabe R, Sekine T, Kurita K, Taniguchi T, Kondo T, Kato T, Shimomura O, Urakawa S. 1996. In-situ measurement of viscosity and density of carbonate melts at high pressure. Earth and Planetary Science Letters, 143(1-4): 207–215. DOI:10.1016/0012-821X(96)00139-2
[] Ernst R E, Bell K. 2010. Large igneous provinces (LIPs) and carbonatites. Mineralogy and Petrology, 98(1-4): 55–76. DOI:10.1007/s00710-009-0074-1
[] Foley S F, Yaxley G M, Rosenthal A, Buhre S, Kiseeva E S, Rapp R P, Jacob D E. 2009. The composition of near-solidus melts of peridotite in the presence of CO2 and H2O between 40 and 60 kbar. Lithos, 112(S1): 274–283.
[] Fukao Y, Obayashi M, Nakakuki T, the Deep Slab Project Group. 2009. Stagnant slab:A review. Annual Review of Earth and Planetary Sciences, 37(1): 19–46. DOI:10.1146/annurev.earth.36.031207.124224
[] Grassi D, Schmidt M W. 2011. The melting of carbonated pelites from 70 to 700 km depth. Journal of Petrology, 52(4): 765–789. DOI:10.1093/petrology/egr002
[] Gudfinnsson G H, Presnall D C. 2005. Continuous gradations among primary carbonatitic, kimberlitic, melilititic, basaltic, picritic, and komatiitic melts in equilibrium with garnet lherzolite at 3~8 GPa. Journal of Petrology, 46(8): 1645–1659. DOI:10.1093/petrology/egi029
[] Hammouda T, Keshav S. 2015. Melting in the mantle in the presence of carbon:Review of experiments and discussion on the origin of carbonaties. Chemical Geology, 418: 171–188. DOI:10.1016/j.chemgeo.2015.05.018
[] Harte B. 2010. Diamond formation in the deep mantle:The record of mineral inclusions and their distribution in relation to mantle dehydration zones. Mineralogical Magazine, 74(2): 189–215. DOI:10.1180/minmag.2010.074.2.189
[] Hatch G P. 2012. Dynamics in the global market for rare earths. Elements, 8(5): 341–346. DOI:10.2113/gselements.8.5.341
[] Hoernle K A, Tilton G R, Le Bas M J, Duggen S, Garbe-Schönberg D. 2002. Geochemistry of oceanic carbonatites compared with continental carbonatites:Mantle recycling of oceanic crustal carbonate. Contributions to Mineralogy and Petrology, 142(5): 520–542. DOI:10.1007/s004100100308
[] Högbohm A G. 1895. Über das Nephelinsyenitgebiet auf der Insel Alnö. Geology Fören Stockholm Förh, 17: 100–160. DOI:10.1080/11035899509453263
[] Hou Z Q, Tian S H, Yuan Z X, Xie Y L, Yin S P, Yi L S, Fei H C, Yang Z M. 2006. The Himalayan collision zone carbonatites in western Sichuan, SW China:Petrogenesis, mantle source and tectonic implication. Earth and Planetary Science Letters, 244(1-2): 234–250. DOI:10.1016/j.epsl.2006.01.052
[] Hou Z Q, Tian S H, Xie Y L, Yang Z S, Yuan Z X, Yin S P, Yi L S, Fei H C, Zou T R, Bai G, Li X Y. 2009. The Himalayan Mianning-Dechang REE belt associated with carbonatite-alkaline complexes, eastern Indo-Asian collision zone, SW China. Ore Geology Reviews, 36(1-3): 65–89. DOI:10.1016/j.oregeorev.2009.03.001
[] Humphries M. 2010. Rare Earth Elements:The Global Supply Chain. Congressional Research Service.
[] Ingrid H K. 1998. Rare earth elements in sövitic carbonatites and their mineral phases. Journal of Petrology, 39(11-12): 2105–2121. DOI:10.1093/petroj/39.11-12.2105
[] Ionov D, Harmer R E. 2002. Trace element distribution in calcite-dolomite carbonatites from Spitskop:Inferences for differentiation of carbonatite magmas and the origin of carbonates in mantle xenoliths. Earth and Planetary Science Letters, 198(3-4): 495–510. DOI:10.1016/S0012-821X(02)00532-0
[] Isshiki M, Irifune T, Hirose K, Ono S, Ohishi Y, Watanuki T, Nishibori E, Takata M, Sakata M. 2004. Stability of magnesite and its high-pressure form in the lowermost mantle. Nature, 427(6969): 60–63. DOI:10.1038/nature02181
[] Kato Y, Fujinaga K, Nakamura K, Takaya Y, Kitamura K, Ohta J, Toda R, Nakashima T, Iwamori H. 2011. Deep-sea mud in the Pacific Ocean as a potential resource for rare-earth elements. Nature Geoscience, 4(8): 535–539. DOI:10.1038/ngeo1185
[] Keshav S, Gudfinnsson G H, Presnall D C. 2011. Melting phase relations of simplified carbonated peridotite at 12-26 GPa in the systems CaO-MgO-SiO2 and CaO-MgO-Al2O3-SiO2-CO2:Highly calcic magmas in the Transition Zone of the Earth. Journal of Petrology, 52: 2265–2291. DOI:10.1093/petrology/egr048
[] Kjarsgaard B A, Hamilton D L. 1989. Carbonatite origin and diversity. Nature, 324(6216): 547–548.
[] Le Maitre R W, Streckeisen A, Zanettin B, Le Bas M J, Bonin B, Bateman P. 2002. Igneous Rocks: A Classification and Glossary of Terms: Recommendations of the International Union of Geological Sciences Subcommission on the Systematics of Igneous Rocks. Cambridge, U. K.: Cambridge University Press.
[] Le Bas M J. 2006. Re-interpretation of zircon date in a carbonatite dyke at the Bayan Obo giant REE-Fe-Nb deposit, China. Acta Petrologica Sinica, 22(2): 517–518.
[] Lee W J, Wyllie P J. 1998. Processes of crustal carbonatite formation by liquid immiscibility and differentiation, elucidated by model systems. Journal of Petrology, 39(11-12): 2005–2013. DOI:10.1093/petroj/39.11-12.2005
[] Ling M X, Liu Y L, Williams I S, Teng F Z, Yang X Y, Ding X, Wei G J, Xie L H, Deng W F, Sun W D. 2013. Formation of the world's largest REE deposit through protracted fluxing of carbonatite by subduction-derived fluids. Scientific Reports, 3: 1776.
[] Litasov K, Ohtani E. 2010. The solidus of carbonated eclogite in the system CaO-Al2O3-MgO-SiO2-Na2O-CO2 to 32 GPa and carbonatite liquid in the deep mantle. Earth and Planetary Science Letters, 295(1-2): 115–126. DOI:10.1016/j.epsl.2010.03.030
[] Liu Y, Hou Z Q, Tian S H, Zhang Q C, Zhu Z M, Liu J H. 2015. Zircon U-Pb ages of the Mianning-Dechang syenites, Sichuan Province, southwestern China:Constraints on the giant REE mineralization belt and its regional geological setting. Ore Geology Reviews, 64: 554–568. DOI:10.1016/j.oregeorev.2014.03.017
[] Martin L H J, Schmidt M W, Mattsson H B, Guenther D. 2013. Element partitioning between immiscible carbonatite and silicate melts for dry and H2O-bearing systems at 1~3 GPa. Journal of Petrology, 54(11): 2301–2338. DOI:10.1093/petrology/egt048
[] Marty B, Tolstikhin I, Kamensky I G, Nivin V, Balaganskaya E, Zimmermann J L. 1998. Plume-derived rare gases in 380 Ma carbonatites from the Kola region (Russia) and the argon isotopic composition in the deep mantle. Earth and Planetary Science Letters, 164(1-2): 179–192. DOI:10.1016/S0012-821X(98)00202-7
[] Nelson D R, Chivas A R, Chappell B W, McCulloch M T. 1988. Geochemcal and kotopic systematics in cartonatites and implications for the evdution of ocean-island sources. Geochimica et Cosmochimica Acta, 52: 1–17. DOI:10.1016/0016-7037(88)90051-8
[] Rohrbach A, Ballhaus C, Golla-Schindler U, Ulmer P, Kamenetsky V S, Kuzmin D V. 2007. Metal saturation in the upper mantle. Nature, 449(7161): 456–458. DOI:10.1038/nature06183
[] Rohrbach A, Schmidt M W. 2011. Redox freezing and melting in the Earth's deep mantle resulting from carbon-iron redox coupling. Nature, 472(7342): 209–212. DOI:10.1038/nature09899
[] Song W L, Xu C, Qi L, Zhou L, Wang L J, Kynicky J. 2015. Genesis of Si-rich carbonatites in Huanglongpu Mo deposit, Lesser Qinling orogen, China and significance for Mo mineralization. Ore Geology Reviews, 64: 756–765. DOI:10.1016/j.oregeorev.2014.04.003
[] Stagno V, Frost D J. 2010. Carbon speciation in the asthenosphere:Experimental measurements of the redox conditions at which carbonate-bearing melts coexist with graphite or diamond in peridotite assemblages. Earth and Planetary Science Letters, 300: 72–84. DOI:10.1016/j.epsl.2010.09.038
[] Stein H J, Markey R J, Morgan M J, Du A, Sun Y. 1997. Highly precise and accurate Re-Os ages for molybdenite from the East Qinling-Dabie molybdenum belt, Shaanxi province, China. Economic Geology, 92(7-8): 827–835. DOI:10.2113/gsecongeo.92.7-8.827
[] Tappert R, Stachel T, Harris J W, Muehlenbachs K, Ludwig T, Brey G P. 2005. Diamonds from Jagersfontein (South Africa): Messengers from the sublithospheric mantle. Contributions to Mineralogy and Petrology, 150(5): 505–522. DOI:10.1007/s00410-005-0035-6
[] Thomsen T B, Schmidt M W. 2008. Melting of carbonated pelites at 2.5-5.0 GPa, silicate-carbonatite liquid immiscibility, and potassium-carbon metasomatism of the mantle. Earth and Planetary Science Letters, 267(1-2): 17–31. DOI:10.1016/j.epsl.2007.11.027
[] Thomson A R, Walter M J, Kohn C S, Brooker A R. 2016. Slab melting as a barrier to deep carbon subduction. Nature, 529(7584): 76–79. DOI:10.1038/nature16174
[] Tilton G R, Bryce J G, Mateen A. 1998. Pb-Sr-Nd isotope data from 30 and 300 Ma collision zone carbonatites in Northwest Pakistan. Journal of Petrology, 39(11-12): 1865–1874. DOI:10.1093/petroj/39.11-12.1865
[] Veksler I V, Petibon C, Jenner G A, Dorfman A M, Dingwell D B. 1998. Trace element partitioning in immiscible silicate-carbonate liquid systems:An initial experimental study using a centrifuge autoclave. Journal of Petrology, 39(11-12): 2095–2104. DOI:10.1093/petroj/39.11-12.2095
[] Veksler I V, Dorfman A M, Dulski P, Kamenetsky V S, Danyushevsky L V, Jeffries T, Dingwell D B. 2012. Partitioning of elements between silicate melt and immiscible fluoride, chloride, carbonate, phosphate and sulfate melts, with implications to the origin of natrocarbonatite. Geochimica et Cosmochimica Acta, 79: 20–40. DOI:10.1016/j.gca.2011.11.035
[] Wallace M E, Green D H. 1988. An experimental determination of primary carbonatite magma composition. Nature, 335(6188): 343–346. DOI:10.1038/335343a0
[] Walter M J, Bulanova G P, Armstrong L S. 2008. Primary carbonatite melt from deeply subducted oceanic crust. Nature, 454: 622–626. DOI:10.1038/nature07132
[] Walter M J, Kohn S C, Araujo D, Bulanova G P, Smith C B, Gaillou E, Wang J, Steele A, Shirey S B. 2011. Deep mantle cycling of oceanic crust:Evidence from diamonds and their mineral inclusions. Science, 334(6052): 54–57. DOI:10.1126/science.1209300
[] Woodard J, Hetherington C J. 2014. Carbonatite in a post-collisional tectonic setting:Geochronology and emplacement conditions at Naantali, SW Finland. Precambrian Research, 240: 94–107. DOI:10.1016/j.precamres.2013.10.017
[] Woolley A R, Kjarsgaard B A. 2008. Paragenetic types of carbonatite as indicated by the diversity and relative abundances of associated silicate rocks:Evidence from a global database. Canadian Mineralogist, 46(4): 741–752. DOI:10.3749/canmin.46.4.741
[] Wyllie P J, Huang W L. 1975. Peridotite, kimberlite, and carbonatite explained in the system CaO-MgO-SiO2-CO2. Geology, 3: 621–624. DOI:10.1130/0091-7613(1975)3<621:PKACEI>2.0.CO;2
[] Wyllie P J, Jones A P, Deng J. 1996. Rare earth elements in carbonate-rich melts from mantle to crust. In: Jones A P, Wall F, Williams C T, eds. Rare Earth Minerals: Chemistry, Origin and Ore Deposits. Netherlands: Springer
[] Xie Y L, Li Y X, Hou Z Q, Cooke D R, Danyushevsky L, Dominy S C, Yin S P. 2015. A model for carbonatite hosted REE mineralisation-the Mianning-Dechang REE belt, western Sichuan Province, China. Ore Geology Reviews, 70: 595–612. DOI:10.1016/j.oregeorev.2014.10.027
[] Xu C, Huang Z L, Liu C Q, Qi L, Li W B, Guan T. 2003. Geochemistry of carbonatites in Maoniuping REE deposit, Sichuan Province, China. Science in China (D) , 46: 246–256. DOI:10.1360/03yd9023
[] Xu C, Campbell I H, Allen C M, Huang Z L, Qi L, Zhang H Zhang G S. 2007. Flat rare earth element patterns as an indicator of cumulate processes in the Lesser Qinling carbonatites, China. Lithos, 95(3-4): 267–278. DOI:10.1016/j.lithos.2006.07.016
[] Xu C, Campbell IH, Kynicky J, Allen C M, Chen Y J, Huang Z L, Qi L. 2008. Comparison of the Daluxiang and Maoniuping carbonatitic REE deposits with Bayan Obo REE deposits, China. Lithos, 106(1-2): 12–24. DOI:10.1016/j.lithos.2008.06.005
[] Xu C, Kynicky J, Chakhmouradian A N, Qi L, Song W L. 2010a. A unique Mo deposit associated with carbonatites in the Qinling orogenic belt, central China. Lithos, 118(1-2): 50–60. DOI:10.1016/j.lithos.2010.03.013
[] Xu C, Kynicky J, Chakhmouradian A N, Campbell I H, Allen C M. 2010b. Trace-element modeling of the magmatic evolution of rare-earth-rich carbonatite from the Miaoya deposit, central China. Lithos, 118(1-2): 145–155. DOI:10.1016/j.lithos.2010.04.003
[] Xu C, Taylor R X, Li W B, Kynicky J, Chakhmouradian A R, Song W L. 2012. Comparison of fluorite geochemistry from REE deposits in the Panxi region and Bayan Obo, China. Journal of Asian Earth Sciences, 57: 76–89. DOI:10.1016/j.jseaes.2012.06.007
[] Xu C, Chakhmouradian A N, Taylor R X, Kynicky, J, Li W B, Song W L, Fletcher I R. 2014. Origin of carbonatites in the South Qinling orogen:Implications for crustal recycling and timing of collision between the South and North China Blocks. Geochimica et Cosmochimica Acta, 143: 189–206. DOI:10.1016/j.gca.2014.03.041
[] Xu C, Kynicky J, Chakhmouradian A R, Li X H, Song W L. 2015. A case example of the importance of multi-analytical approach in deciphering carbonatite petrogenesis in South Qinling orogen:Miaoya rate-metal deposit, central China. Lithos, 227: 107–121. DOI:10.1016/j.lithos.2015.03.024
[] Yasukawa K, Liu H J, Fujinaga K, Machida S, Haraguchi S, Ishii T, Nakamura K, Kato Y. 2014. Geochemistry and mineralogy of REY-rich mud in the eastern Indian Ocean. Journal of Asian Earth Sciences, 93: 25–36. DOI:10.1016/j.jseaes.2014.07.005
[] Yaxley G M, Brey G P. 2004. Phase relations of carbonate-bearing eclogite assemblages from 2.5 to 5.5 GPa:Implications for petrogenesis of carbonatites. Contributions to Mineralogy and Petrology, 146(5): 606–619. DOI:10.1007/s00410-003-0517-3
[] Zajacz Z, Halter W. 2007. LA-ICPMS analyses of silicate melt inclusions in co-precipitated minerals:Quantification, data analysis and mineral/melt partitioning. Geochimica et Cosmochimica Acta, 71(4): 1021–1040. DOI:10.1016/j.gca.2006.11.001
[] Zhang Y X, Zindler A. 1993. Distribution and evolution of carbon and nitrogen in Earth. Earth and Planetary Science Letters, 117(3-4): 331–345. DOI:10.1016/0012-821X(93)90088-Q
[] Zhao G C, Sun M, Wilde S A, Li S Z. 2005. Late Archean to Paleoproterozoic evolution of the North China Craton:Key issues revisited. Precambrian Research, 136(2): 177–202. DOI:10.1016/j.precamres.2004.10.002
[] 路凤香, 张本仁, 韩吟文, 钟增球, 凌文黎, 张宏飞, 郑建平, 侯青叶. 2006. 秦岭-大别—苏鲁地区岩石圈三维化学结构特征. 北京: 地质出版社.
[] 蒲广平. 2001. 攀西地区稀土成矿历史演化与喜山期成矿基本特征. 见: 陈毓川, 王登红. 喜马拉雅期内生成矿作用研究. 北京: 地震出版社: 104-116
[] 宋文磊, 许成, 王林均, 吴敏, 曾亮, 王丽泽, 冯梦. 2013. 与碳酸岩-碱性杂岩体相关的内生稀土矿床的矿化特征及稀土富集机制评述. 北京大学学报 (自然科学版) , 49: 725–740.
[] 田世洪, 侯增谦, 杨竹森, 陈文, 杨志明, 袁忠信, 谢玉玲, 费红彩, 尹淑苹, 刘英超, 李政, 李小渝. 2008. 川西冕宁-德昌REE成矿带成矿年代学研究:热液系统维系时限和构造控矿模型约束. 矿床地质, 27(2): 177–187.
[] 吴福元, 徐义刚, 高山, 郑建平. 2008. 华北岩石圈减薄与克拉通破坏研究的主要学术争论. 岩石学报, 24(6): 1145–1174.
[] 张本仁, 高山, 张宏飞, 韩吟文. 2002. 秦岭造山带地球化学. 北京: 科学出版社.
[] 张国伟, 张本仁, 袁学诚, 肖庆辉. 2001. 秦岭造山带与大陆动力学. 北京: 科学出版社.