岩石学报  2018, Vol. 34 Issue (9): 2716-2732   PDF    
鄂东南铜绿山矿床石榴子石显微结构及微区成分对成矿过程的指示
纪敏1,2 , 赵新福2 , 曾丽平2 , 范田纬腾2     
1. 中国科学院壳幔物质与环境重点实验室, 中国科学技术大学地球与空间科学学院, 合肥 230026;
2. 中国地质大学地质过程和矿产资源国家重点实验室, 中国地质大学资源学院, 武汉 430074
摘要:鄂东南矿集区铜绿山矿床是典型的矽卡岩型铜铁多金属矿床,矿体产出在铜绿山岩体与三叠系碳酸盐地层的接触带。尽管本矿床的研究程度很高,但对早期成矿流体的成分与演化及矿质富集沉淀等过程的精细制约依然比较欠缺。石榴子石在铜绿山矿床中分布广泛,本文对不同产状的石榴子石利用SEM显微结构、EPMA主量元素和LA-ICPMS微量元素分析,去探讨石榴子石的生长动力学及其对成矿过程的指示。大理岩和镁质外矽卡岩中的石榴子石均以钙铝榴石为主,端元成分比较均一,∑REE含量低(3.01×10-6~14.34×10-6),具有弱的Eu异常、轻微富集LREE的模式。钙质外矽卡岩和内矽卡岩中的石榴子石以钙铁榴石为主,端元成分变化较大,单颗粒从核部到边部Fe含量具有增加的趋势,其∑REE含量较高(22.71×10-6~806.8×10-6),具有明显的正Eu异常、极度富集LREE亏损HREE的模式。大理岩中的钙铝榴石环带不发育,表明在晶体生长过程中,界面反应速率很慢且物质迁移以扩散为主;而镁质外矽卡岩中的钙铝榴石可见振荡环带,表明界面反应速率比以扩散为主要方式的物质迁移快,这与白云石相对于方解石具有较低的反应吉布斯自由能,因此较难被热液消耗有关;钙质外矽卡岩中的钙铁榴石振荡环带清晰,该类石榴子石只有在界面反应速率较快且物质迁移方式以对流为主的情况下才能形成;内矽卡岩中的钙铁榴石的环带杂乱,尽管其形成的动力学模式与钙质外矽卡岩中的类似,但由于其距离岩体更近,热量迁移困难,在结晶生长受阻的同时由于对流提供的充足物质使得原先的环带发生溶解再沉淀进而形成结构混乱的再吸收环带。上述生长动力学模式也能够很好的对应Nb、Ta含量的脱耦和轻重稀土分异的特征。铜绿山不同产状的石榴子石中,外矽卡岩带中钙铁榴石Eu含量最高,而Eu2+与成矿金属元素一样可被Cl-络合的,其含量的高低可能指示了对应的成矿元素在成矿热液中的浓度。石榴子石环带自核部向边部Fe含量具有逐渐增加的趋势,反应干矽卡岩阶段从早到晚,热液中的金属成矿元素含量增加。
关键词: 石榴子石     环带     微量元素     生长动力学     热液演化     铜绿山矽卡岩矿床    
Microtexture and geochemistry of garnets from Tonglushan skarn Cu-Fe deposit in the southeastern Hubei metallogenic province: Implications for ore-forming process
JI Min1,2, ZHAO XinFu2, ZENG LiPing2, FAN TianWeiTeng2     
1. CAS Key Laboratory of Crust-Mantle Materials and Environments, School of Earth and Space Sciences, University of Science and Technology of China, Hefei 230026, China;
2. State Key Laboratory of Geological Processes and Mineral Resources, Faculty of Earth Resources, China University of Geosciences, Wuhan 430074, China
Abstract: The Tonglushan Cu-Fe deposit is one of typical skarn-type deposits in southeastern Hubei metallogenic province. 13 ore bodies occur along the contact zone between the Tonglushan quartz diorite pluton and the Lower Triassic carbonate rocks. Although numerous studies have been made, the evolution of the early ore fluids and mechanism of precipitation of ore minerals is still unclear. Garnet is widespread in Tonglushan skarn deposit, and has a record of the composition, properties and evolution of the hydrothermal fluid. We perform SEM microtexture observation, EPMA major elements analysis and LA-ICPMS trace elements analysis to understand the growth kinetics of garnets and its implication for ore forming process. Garnets from marble and magnesian exoskarn zone are mainly grossulars, which are homogeneous and have low ∑REE (3.01×10-6~14.34×10-6) with LREE slightly enriched patterns with weak Eu anomalies. In contrast, andradites are dominate phases in calcareous exoskarn and endoskarn, which have a large range of variation and have Fe content increasing from core to rim in single grains. Their ∑REE are very high (22.71×10-6~806.8×10-6) and exhibit LREE-rich and HREE-depleted patterns with strong positive Eu anomaly. The grossulars in marble have no zonation, indicating that the rate of interface reaction is slower than mass transfer, which is mainly controlled by diffusion during crystal growth. Whereas the magnesian exoskarn grossulars show oscillatory zonings, suggesting that the rate of interface reaction is faster than diffusion dominated mass transfer. The mass transfer is relatively slow, which may be related to the lower reaction Gibbs free energy of dolomite than calcite. The andradites in calcareous exoskarn develop oscillatory zonings, which is related to the process of high rate of interface reaction and the advection dominated mass transfer. The zoning of endoskarn andradites is very chaotic, although it has growth kinetics similar to the andradites in calcareous exoskarn. However, the difficulties of heat transfer limit the crystal growth. Meanwhile, fluid advection is enhanced by dissolving the previous zoning to form absorbed zoning. This pattern of garnets growth kinetics can also explain the characteristics of Nb/Ta content decoupling and differentiation between LREE and HREE. In regard to the highest concentration of Eu in andradites of exoskarn zone, Eu2+, just like iron ion, could combine with Cl- to form stable complex, so the content of which could be an indicator of the corresponding metal ions in hydrothermal fluid. The content of Fe in garnets zonings increases from core to rim, indicating the content of metal ions increase from early to late prograde skarn stage.
Key words: Garnet     Zoning     Trace elements     Growth kinetics     Fluid evolution     Tonglushan skarn deposit    

石榴子石能够在比较宽泛的物理化学的条件下形成,在各类岩石中分布广泛,它具有较高的Lu/Hf和Sm/Nd比值,可应用于同位素地质年代学、地球化学示踪的工作(Gaspar et al., 2008)。石榴子石的端元类型较多,元素和同位素的环带很发育,存在生长环带、扩散环带、再吸收环带等类型,因而在地质温压计、地质速率计上具有广泛的应用(Kohn, 2014)。以铁铝榴石和镁铝榴石为主要端元成分的石榴子石是超高压变质岩、金伯利岩以及岩浆岩中限定岩石形成的年代、温压条件以及演化历史的重要矿物,研究程度较高(Li et al., 2000; Lackey et al., 2012),但以钙铝榴石和钙铁榴石为主要端元的矽卡岩中石榴子石的研究却相对较少。

石榴子石是矽卡岩中代表性的矿物,如果在其生长时没有发生强烈的构造变形,一般都能形成非常好的晶型,常发育显著的振荡环带(Jamtveit et al., 1992),环带的结构与主微量元素、同位素的变化特征受到内部晶体化学与外界物质供给和物理化学条件的影响(Shore and Fowler, 1996; Jamtveit and Hervig, 1994; Hu et al., 2014),因而能够完整的记录矿物结晶时热液的组成、性质和演化(Jamtveit and Hervig, 1994; Smith et al., 2004; Gaspar et al., 2008; Zhai et al., 2014)。Jamtveit等(Jamtveit, 1991; Jamtveit and Andersen, 1992, 1993; Jamtveit et al., 1992, 1993, 1995)对矽卡岩中的石榴子石进行了系统的研究工作,将石榴子石环带结构与地球化学特征同生长速率、热液演化以及构造背景等联系起来。在大多数矽卡岩中,石榴子石的特点是核部更富钙铝榴石、边部更富钙铁榴石端元,早期晶体生长速率与流体的运移速率相对较慢,晚期的晶体生长速度与流体运移速率相对较快(Jamtveit and Hervig, 1994; Meinert et al., 2005)。Somarin (2004)总结了矽卡岩型矿床中石榴子石端元组成与矿种的关系,认为石榴子石的主量元素特征可以作为铜矿化的指示剂。振荡环带发育的石榴子石中存在着O同位素和微量元素突变的层,这些突变层的形成受到生长时的化学不平衡、大尺度的流体流动行为和构造活动的影响(Chernoff and Carlson, 1999; Crowe et al., 2001; D'Errico et al., 2012; Ismail et al., 2014),因而石榴子石的同位素和微量元素指标可以用来反演这些地质过程。多数矽卡岩矿床的蚀变分带非常明显,不同蚀变带中石榴子石的特征差别较大,靠近岩体的内矽卡岩中石榴子石的数量较辉石多且颜色较深,而靠近围岩的外矽卡岩中石榴子石的数量较辉石少且颜色较浅(Meinert et al., 2005)。作为空间分布广泛的早期矽卡岩矿物,石榴子石的主微量元素特征能够用来反演矽卡岩矿床形成过程中早期成矿流体性质和时空演化特征(Dziggel et al., 2009; Zhai et al., 2014; 朱乔乔等, 2014; Peng et al., 2015),但目前对矽卡岩矿床不同蚀变分带中产出的石榴子石的矿物学、地球化学组成和生长动力学特征等方面进行综合对比研究的工作并不多,而这些工作能够对早期成矿热液的时空演化规律和成矿过程进行有效的制约。

鄂东南矿集区铜绿山矿床是我国著名的矽卡岩型铜铁矿床,矿床的形成与~140Ma侵入到三叠系碳酸盐地层中的石英闪长岩相关(Xie et al., 2008, 2011; Li et al., 2009, 2010a, b),辉钼矿Re-Os、热液榍石U-Pb以及云母Ar-Ar定年确定成矿时代为136~140Ma(Xie et al., 2007, 2011; 谢桂青等, 2009; Li et al., 2010a, 2011),与成岩时代相近。但更年轻的岩浆作用(120~130Ma)及热液活动事件(~121Ma)也有报导,表明矿床可能经历过后期局部的热液改造作用(Li et al., 2010a; Xie et al., 2011, 2012)。硫化物和硫酸盐的S同位素数据表明成矿物质主要来自于岩浆,但在形成过程之中有膏岩层的参与(余元昌等, 1985; 王彦博等, 2011; Xie et al., 2015)。Zhao et al. (2012)对不同成矿期和成矿阶段形成的矿物系统的开展了H-O同位素和流体包裹体的研究,发现主要起源于岩浆热液的干矽卡岩阶段的流体具有高温、高盐度的特征;随着持续的降压和降温以及大气降水的加入,湿矽卡岩阶段流体的盐度逐渐降低,导致磁铁矿沉淀成矿;而具有岩浆起源同位素特征的石英硫化物阶段流体的温度和盐度变化范围较大,流体沸腾作用及大气降水加入导致硫化物沉淀;最后碳酸盐阶段流体的温度和盐度都非常低,可能与大气降水的大量加入相关。Xie et al.(2011, 2016)分析硫化物中的Re含量以及He-Ar同位素,认为成矿热液中有地幔组分以及少量大气水的贡献。尽管前人在成矿流体的性质和演化方面做了大量的工作,但是这些工作主要是反映时间跨度较大的成矿阶段之间的演化特征,对于如早期干矽卡岩阶段热液流体的演化特征、蚀变带不同位置流体性质的差异及其对铁质沉淀的影响等成矿过程中的细节问题缺乏制约。本文在对铜绿山矿床中蚀变分带不同位置的石榴子石进行详细的岩相学结构观察的基础上,通过电子探针(EPMA)和激光剥蚀电感耦合等离子质谱(LA-ICPMS)进行原位高精度的主微量元素的分析,揭示石榴子石的环带成分特征,进一步探讨REE加入矿物晶格的机制、石榴子石生长动力学模式、成矿热液性质和时空演化特征及其对成矿过程的指示。

① 余元昌, 李刚, 肖国荃, 杨洪朝, 薛迪康, 刘云龙, 李家珍. 1985.湖北省大冶县铜绿山接触交代铜铁矿床.大冶:湖北省鄂东南地质大队, 1-102

1 地质背景及矿床地质特征

鄂东南矿集区位于长江中下游成矿带的最西端,在大地构造位置上位于扬子板块北缘、大别造山带的西南部(翟裕生等, 1992; 周涛发等, 2012)。处于麻城-团风断裂(NE向)、襄樊-广济断裂(NW向)和高桥-阳新断裂(近EW向)圈并成的三角地带中(图 1)。鄂东南地区地层出露相对齐全,除缺失的中、下泥盆统和上侏罗统外,其他时代的地层发育基本完整、层序齐全(余元昌等, 1985)。区域的构造比较复杂,印支期主要形成区域性的构造规模大的NW-NWW向的复式褶皱、走向断裂,对区域成岩成矿影响较大,而燕山期主要形成规模较小的NE-NNE向的褶皱、断裂(舒全安等, 1992)。区内中生代的岩浆活动频繁,从北到南分布有鄂城、铁山、金山店、阳新、灵乡和殷祖六大岩体以及铜山口花岗闪长岩体、丰山洞花岗斑岩体等小岩体以及金牛盆地中的火山岩(余元昌等, 1985)。鄂东南地区矿产资源丰富,是我国重要的铁铜矿产的产地,此外还伴生有铅、锌、银、金、钨、钼等矿种(舒全安等, 1992)。

图 1 鄂东南矿集区地质简图(据舒全安等, 1992修改) Fig. 1 Simplified geological map of southeastern Hubei metallogenic province (modified after Shu et al., 1992)

铜绿山矽卡岩型铜铁多金属矿床是鄂东南矿集区最大、最富的矿床,主要出产Cu、Fe、Au和Ag,此外还有Co、Ga、In、Se、Te、Mo和Re等矿种(Li et al., 2014)。矿床包含13个矿体(图 2),根据这些矿体在地表展布的方向,划分出北北东向(北偏东22°)、北东东向(北偏东60°)和北北西向(北偏西15°~20°)三个矿带。矿体或呈舒缓的舌状或呈陡倾斜状分布在大理岩之下、穿插大理岩或者超覆于大理岩之上,而大理岩主要是呈捕虏体状分布在石英二长闪长岩体内(余元昌等, 1985)。铜绿山矿床的接触交代变质作用强烈,矽卡岩空间上的分带非常发育,根据原岩的不同可以划分为钙矽卡岩类为主和镁矽卡岩类为主的分带。此外还可以根据矿物组合的差异从岩体到围岩划分为石英闪长岩→内矽卡岩(斜长石岩-石榴子石矽卡岩)→外矽卡岩(透辉石矽卡岩)→矽卡岩化的大理岩(余元昌等, 1985)。

图 2 铜绿山矿床地质简图(a, 据Li et al., 2014)及31号(b)、5号(c)和4号(d)勘探线剖面中采样点分布简图(据魏克涛等, 2010修改) Fig. 2 Geological map of Tonglushan skarn Cu-Fe (Au) deposit (a, after Li et al., 2014), also showing locations of the studied samples in the cross section of exploration line 31 (b), 5 (c) and 4 (d)

① 魏克涛, 胡清乐, 黄智辉等. 2010.湖北省大冶市铜录山铜矿接替资源勘查(深部普查)报告.大冶:湖北省鄂东南地质大队, 1-111

铜绿山矿床中主要的金属矿物有磁铁矿、赤铁矿、黄铜矿、斑铜矿等,共伴生的非金属矿物有石榴子石(主要为钙铝榴石和钙铁榴石端元,如果钙铝端元成分大于钙铁榴石的端元成分则称该石榴子石为钙铝榴石,反之则称为钙铁榴石)、透辉石、阳起石、金云母、绿帘石、方解石等。通过野外观察以及岩相学工作,根据主要矿物组合、结构和构造关系,将铜绿山矿床的成矿过程划分为Ⅰ矽卡岩期(ⅰ干矽卡岩阶段、ⅱ湿矽卡岩阶段和ⅲ氧化物阶段)、Ⅱ石英硫化物期(ⅳ石英硫化物阶段和ⅴ碳酸盐阶段),其中氧化物阶段和石英硫化物阶段分别是铁矿和铜矿的主成矿期,而石榴子石是在干矽卡岩阶段形成的早期矽卡岩矿物之一。

2 样品描述及分析方法

本文系统的采集了铜绿山矿床中不同产状的石榴子石样品,包括矽卡岩化大理岩(14TLS64)、镁质外矽卡岩(14TLS22)、钙质外矽卡岩(14TLS56-1、14TLS78)和内矽卡岩(14TLS75、14TLS89)中的石榴子石(图 3)。其中14TLS64和14TLS56-1采集于ZK3102中735m处矽卡岩化大理岩和728m处透辉石矽卡岩(图 2b),14TLS78和14TLS75采集于TZK402中833m处的透辉石矽卡岩和850m处石榴子石矽卡岩,14TLS89采集于ZK504中726m处的石榴子石矽卡岩,而14TLS22采自铜绿山露天采坑附近的镁矽卡岩带的透辉石矽卡岩。自围岩到岩体,产出石榴子石含量逐渐增加、颗粒颜色逐渐变深,这与Meinert et al. (2005)总结的多数矽卡岩矿床的分带特征一致。矽卡岩化大理岩中主要矿物组合为方解石、石榴子石和少量的透辉石,淡黄色的钙铝榴石呈条带状分布在白色的大理岩中(图 3a),半自形-他形粒状结构,粒径0.6~1.3mm,裂理比较发育;镁质外矽卡岩带主要为透辉石矽卡岩,主要矿物有透辉石、钙铝榴石、尖晶石、绿帘石、绿泥石等(图 3b),黄绿色的半自形-他形的钙铝榴石颗粒与透辉石和镁铁尖晶石共生,大小0.8~3mm,被绿帘石、阳起石、绿泥石等交代。钙质外矽卡岩带中矿化的透辉石矽卡岩中的主要矿物有透辉石、石榴子石、绿帘石、阳起石、磁铁矿、黄铜矿等,钙铁榴石呈浅黄色(图 3c),自形-半自形粒状结构,粒径0.8~5mm,后期的绿帘石、阳起石、磁铁矿、黄铜矿、斑铜矿、方解石等沿着颗粒间隙和裂理生长。内矽卡岩带中矿化的石榴子石矽卡岩中以石榴子石为主,含少量的透辉石(图 3d),褐红色的钙铁榴石呈他形粒状结构,粒径0.5~2.5mm,与透辉石共生,粒间生长着阳起石、磁铁矿、黄铁矿等。

图 3 铜绿山矿床不同产状的石榴子石 (a)矽卡岩化大理岩中的条带状钙铝榴石(14TLS64);(b)镁质外矽卡岩带透辉石矽卡岩中粒状的钙铝榴石(14TLS22);(c)钙质外矽卡岩带透辉石矽卡岩中粒状的钙铁榴石(14TLS56-1);(d)钙质内矽卡岩带石榴子石中块状的钙铁榴石(14TLS89). Di-透辉石;Grt-石榴子石;Spl-尖晶石;Ep-绿帘石;Cc-方解石;Chl-绿泥石;Mt-磁铁矿;Ccp-黄铜矿 Fig. 3 Photographs of representative garnet samples from different zonations of the Tonglushan deposit (a) banded grossular in skarnization marble (14TLS64); (b) granular grossular in diopside skarn of magnesian exoskarn zone (14TLS22); (c) granular andradite in diopside skarn of calcareous exoskarn zone (14TLS56-1); (d) massive andradite in garnet skarn of endoskarn zone (14TLS89). Di-diopside; Grt-garnet; Spl-spinel; Ep-epidote; Cc-calcite; Chl-chlorite; Mt-magnetite; Ccp-chalcopyrite

石榴子石主量元素的电子探针(EPMA)分析在武汉理工大学材料与研究测试中心电子探针室完成,仪器型号为JEOL JXA-8230,配备LaB6电子枪、4道波谱仪。实验条件为:加速电压20kV,电流20nA,束斑直径5μm,所有元素的信号采集时间均为15s,背景时间均为5s,修正方法ZAF,检测限优于100×10-6。电子探针分析直接得到的是氧化物重量百分比,由原始分析结果计算石榴子石晶体化学计量式及端元组成的算法据Locock (2008),由于石榴子石端元成分众多,本文考虑常见的钙铝榴石(Grs)、钙铁榴石(Adr)、铁铝榴石(Alm)、锰铝榴石(Pyp)、镁铝榴石(Sps)这五种类型,而如富集Fe2+的钙钛铁榴石等少见类型未考虑。计算的端元组成即为某一种类型的含量在这五种类型中所占的百分比。

石榴子石微量元素的激光剥蚀电感耦合等离子质谱(LA-ICPMS)分析在中国地质大学(武汉)地质过程与矿产资源国家重点实验室矿床地球化学分室完成。电感耦合等离子体质谱仪由热电公司制造,型号为iCAP Qc。激光剥蚀系统为瑞索公司制造的Resonetics-M50。ArF准分子激光发生器产生193nm深紫外光束,经匀化光路聚焦于矿物表面。激光束斑直径为33μm,频率为10Hz,剥蚀时间55s,以高纯He气为载气,与Ar气和少量N2气混合后进入质谱仪。测试过程中以NIST SRM 612作为信号漂移矫正,以NIST玻璃BCR-2G,BHVO-2G,BIR-1G作为外标,以Si为内标测定主量和微量元素含量。测试数据利用ICPMSDataCal 9.9离线处理(Liu et al., 2010)。在进行EPMA和LA-ICPMS分析时,要求同样的分析点尽量一致,但两种方法的束斑相差较大,两者分析的主量元素大部分在误差范围内一致,但也有少量点存在差别。但考虑到激光剥蚀的面积较大,因此在讨论微量元素和主量元素的相关关系时,采用的是激光分析得出的主量元素数据(数据表略)。

3 分析结果 3.1 岩相学特征

矽卡岩化大理岩中的钙铝榴石蚀变较严重(图 4a图 5a),包含的矿物包裹体较多,环带不发育;镁质外矽卡岩中的钙铝榴石表面很干净(图 4b),具有波状消光的特征(图 4e),局部可见压力侵蚀结构(图 4g),振荡环带比较清楚(图 5b);两个不同接触带中钙质外矽卡岩中的钙铁榴石的均具有聚片双晶的特征,在单偏光下就看到清晰的振荡环带(图 4c),其中颜色较深的环带表现为全消光(图 4f),较浅的环带具有约20°的消光角。有的环带核部富Al(图 5c),有的核部富Fe(图 5f)。BSE图像揭示出可以反应晶体生长速率凹槽状结构(图 5d; Jamtveit et al., 1992),此外BSE图像上还可观察到较暗的层的生长面{211}比较发育,而较亮的层{211}和{110}生长面均比较发育,从核部到边部,{110}面的长度越来越大(图 5c, d);两个不同接触带中内矽卡岩中的石榴子石在单偏光下均存在颜色深浅有差异的区域(图 4d),但是在BSE图像中看到其核部的环带非常混乱,颗粒边缘局部发育小规模的规律环带(图 5e, f)。

图 4 铜绿山矿床中石榴子石的光学和结构特征 (a)大理岩中的钙铝榴石(单偏光);(b)与尖晶石共生的镁质外矽卡岩中的钙铝榴石(单偏光);(c)振荡环带发育的钙质外矽卡岩中的钙铁榴石(单偏光);(d)内矽卡岩中的钙铁榴石颗粒中颜色深浅不一(单偏光);(e)具有波状消光特征的钙铝榴石(正交光);(f)钙铁榴石振荡环带中不同的层消光性质不同,右上方白色方框为左下方对应区域物台旋转20°后全消光的特征(正交光);(g)钙铝榴石中发育的压力侵蚀结构(单偏光);(h)钙铝榴石中较多数量的次生流体包裹体(单偏光);(i)钙铁榴石中含有较多的原生流体包裹体(单偏光). Sp-尖晶石;Ap-磷灰石;Ser-绢云母;Qz-石英; Adr-钙铁榴石; Grs-钙铝榴石; P-原生流体包裹体; S-次生流体包裹体 Fig. 4 Photomicrographs showing the optical and textural characteristics of garnets from the Tonglushan deposit (a) grossular in marble (PPL); (b) magnesian exoskarn grossular co-existing with spinel (PPL); (c) andradite in calcareous exoskarn showing oscillatory zonings under PPL; (d) endoskarn andradite grains with different colors (PPL); (e) wavy extinction developed in grossular (XPL); (f) layers in andradite having contrast extinction property, white box in the right-up showing complete extinction after the substage rotating 20° (XPL); (g) pressure corrosion texture in grossular (PPL); (h, i) in grossular, the type of fluid inclusions is dominated by secondary fluid inclusion, while in andradite is primary fluid inclusion. Sp-spinel; Ap-apatite; Ser-sericite; Qz-quartz; Adr-andradite; Grs-grossularite; P-primary fluid inclusion; S-secondary fluid inclusion

图 5 铜绿山矿床中石榴子石的BSE图像揭示的结构和成分特征 (a)大理岩中的钙铝榴石没有明显环带;(b)镁质外矽卡岩中的钙铝榴石具有振荡环带;(c-f)钙质外矽卡岩中的钙铁榴石振荡环带非常发育,(c、d)中可观察到晶面{110}和{211}不同的结构,(d)中有反应生长速率加快的凹槽状结构;(g、h)为内矽卡岩中的钙铁榴石,核部环带非常的混乱,边部些许发育振荡环带. (b、c、e、f)中的实心和空心的圆形的点代表暗层和亮层的电子探针分析点,方形的代表激光分析点 Fig. 5 BSE images showing the internal textures of garnets from the Tonglushan deposit (a) grossular in marble with no obviously zoning; (b) oscillatory zonings developed in the magnesian exoskarn grossular; (c-f) regular oscillatory zonings in andradite of calcareous exoskarn, the texture difference of growth face {110} and {211} in single grain and hopper-like texture showing in (c, d); (g, h) endoskarn andradite with chaotic zonings in the core and little oscillatory zoning in the rim. The solid and hollow spots in (b, c, e, f) represent the EPMA analysis points of dark and light layers in BSE image, while quadrate spots denote LA-ICPMS analysis points

大理岩和镁质外矽卡岩中的钙铝榴石中孤立分布的原生类型流体包裹体数量较少,主要为一些线状分布的次生类型(图 4h)。而钙质外矽卡岩和内矽卡岩中钙铁榴石中的流体包裹体数量较多(图 4i),且多为原生类型。

3.2 主量元素特征

6件石榴子石样品的电子探针主量元素代表性的分析结果见表 1,计算的端元成分见图 6a。不同产状的石榴子石中的SiO2、CaO、FeOT和Al2O3的含量均较高,SiO2和CaO的含量变化很小,但FeOT和Al2O3的含量差别较大,且FeOT和Al2O3的含量存在负相关关系。同一个颗粒的石榴子石在BSE图像中的亮度存在差异(图 5e),FeOT含量高的层较亮而Al2O3含量较高的层较暗。大理岩和镁质外矽卡岩中石榴子石主要是钙铝榴石端元,而钙质外矽卡岩和内矽卡岩中的主要为钙铁榴石端元。同一个样品的石榴子石的端元组成存在差异,其中钙质外矽卡岩中的差异最大,钙铁榴石端元能相差42.20%,而两个钙铝榴石样品的端元组成比较均一,钙铝榴石端元最大只相差14.15%。同一石榴子石颗粒核边的元素含量也存在差异,大致呈现出从核到边Fe含量逐渐上升,Al含量逐渐下降或保持不变的趋势(图 6b-e)。其中采集于不同接触带的两个钙质外矽卡岩(14TLS56-1和15TLS78)和两个内矽卡岩样品(14TLS75和14TLS89)具有类似的端元组成和端元变化特征(图 6),因而可以进行对比讨论。

表 1 铜绿山矿床中石榴子石代表性的主量元素组成(wt%) Table 1 Representative major element compositions of garnets from Tonglushan deposit (wt%)

图 6 铜绿山矿床中石榴子石的端元组成(a)和环带中核边成分的变化(b-e) (b1、b2)为14TLS-22镁质外矽卡岩中的钙铝榴石; (c1、c2)为14TLS-56-1钙质外矽卡岩中的钙铁榴石; (d、e)为14TLS-78钙质外矽卡岩中的钙铁榴石.箭头指示从核到边的变化. Grs-钙铝榴石;Adr-钙铁榴石;Alm-铁铝榴石;Sps-锰铝榴石;Pyp-镁铝榴石 Fig. 6 Calculated end member compositions of garnets from the Tonglushan deposit (a) The changing trend of garnet compositions from core to rim (b-e), (b1, b2) are granular in magnesian exoskarn zone (14TLS-22), while (c1, c2) are calcareous exoskarn andradite (14TLS-56-1), (d, e) all refer to calcareous exoskarn andradite (14TLS-78). The direction of arrow represents the change of composition from core to rim. Grs-grossular; Adr-andradite; Alm-almandine; Sps-spessartine; Pyp-pyrope
3.3 微量元素特征

石榴子石的LA-ICPMS微量元素代表性的分析结果见表 2图 7,各种产状的石榴子石中大离子亲石元素的含量(如Rb、Ba、Sr)非常低(图 7a),有的甚至低于检出限。相较而言,高场强元素的含量(如Nb、Ta、Zr)要高(图 7b, c),对于地球化学性质很相近的Nb和Ta,钙铁榴石中的Nb和Ta含量很相近,而钙铝榴石中的Nb和Ta的含量出现此高彼低的脱耦现象。U的含量与石榴子石中Fe的含量正相关,钙铁榴石中的U含量很高(图 7k)。

表 2 铜绿山矿床中石榴子石代表性的微量元素组成(×10-6) Table 2 Representative trace element compositions of garnets from Tonglushan deposit (×10-6)

图 7 铜绿山矿床中石榴子石的微量元素组成 (a-c) Sr、Nb、Ta的含量与钙铝榴石端元含量的关系;(d、e) REE球粒陨石标准化模式图(标准化值据Sun and McDonough, 1989);(f-h) ∑REE含量、轻重稀土的分异和Eu异常与钙铝榴石端元组成或含量比值Fe/(Fe+Al)的关系;(i) ∑REE含量与石榴子石晶格中四配位的Al含量的关系;(j、k) U含量与∑REE含量和钙铝榴石端元组成的关系 Fig. 7 Trace elements composition of garnets from Tonglushan deposit (a-c) Plot of the content of Sr, Nb, Ta with Grs; (d, e) Chondrite normalized REE patterns (REE values of Chondrite after Sun and McDonough, 1989); (f-h) plot of ∑REE, LREE/HREE and Eu anomalies with the end member of granular or the content ration of Fe/(Fe+Al); (i) plot of ∑REE with the Al in the four-coordination; (j, k) Plot of U content with ∑REE and Grs

对于稀土元素特征,不同产状的石榴子石相差较大(图 7d, e)。大理岩和镁质外矽卡岩中的钙铝榴石的∑REE(3.01×10-6~14.3×10-6)较低,具有较平坦的REE分配模式,轻微的富集LREE,具有弱的Eu异常。而钙质外矽卡岩和内矽卡岩中的钙铁榴石的∑REE(22.7×10-6~807×10-6)较高,具有陡峭的REE分配模式,多数极度富集LREE,具有明显的正Eu异常,值得注意的是14TLS56-1中有一颗粒中相对富Al的核部∑REE为149.5×10-6,LREE/HREE为3.71,而边部的暗层∑REE为158.6×10-6,LREE/HREE为90.8,这两者∑REE相似,但是轻重稀土元素分异差别较大。∑REE总量、Eu异常和轻重稀土的分异程度均与石榴子石中Fe和Al的含量相关(图 7f-h),其中富Fe的铁铝榴石相对于钙铝榴石具有更高的∑REE、Eu/Eu*和LREE/HREE。微量元素U的含量也与ΣREE近似正相关的关系(图 7j)。

4 讨论 4.1 石榴子石晶格中的稀土元素替代机制及非均质性成因

铜绿山矿床中石榴子石的∑REE与Fe的含量存在正相关的关系,钙铁榴石或者钙铁榴石环带中富Fe的层比钙铝榴石或者钙铁榴石中富Al的层具有更高的∑REE(图 7f)。石榴子石是岛状硅酸盐矿物,其晶体化学通式可以写成X3Y2Z3O12,其中X通常是二价阳离子(如Ca2+、Fe2+、Mg2+、Mn2+以及REE3+、Na+、K+等),Y一般为三价的阳离子(Fe3+、Al3+、Cr3+、V3+等以及Ti4+、Zr4+),Z通常为Si4+(赵珊茸等, 2004; Gaspar et al., 2008)。根据戈尔德施密特定律,REE3+只能通过类质同相的型式替代晶格中与其离子半径相近的八配位的X2+(图 8a)。对于Eu2+,这种替换是等价的,而对于其他的REE,还存在着电价不平衡的问题,为了实现电价平衡,还需其他离子的替代与REE替代耦合。[X2+]-1[REE3+]+1[Si4+]-1[Z3+]+1(YAG式的替换)和[X2+]-2[Na+]+1[REE3+]+1被认为是REE加入石榴子石晶格的两种主要的机制,但也可能存在如[X2+]-3[ ]-1[REE3+]+2少见的缺位替代形式(McIntire, 1963; Smith et al., 2004)。Gaspar et al. (2008)根据Crown Jewel金矿床矽卡岩中石榴子石的∑REE与Al含量呈正相关,而认为YAG式的替换是REE加入石榴子石晶格的机制,而且较多REE的加入导致了钙铝榴石消光的各向异性。但是铜绿山石榴子石与此相反,∑REE与Fe含量大致正相关而与Al含量反相关(图 7i),并且Fe3+不像Al3+能够出现在四次配位的位置上,因此该种替代不太可能发生。而且铜绿山矿床中石榴子石晶格八配位上主要是Ca2+,而且基本上不含有OH,所以Ca2+-Fe2+-Mg2+-Mn2+在立方配位中位置的调整(Griffin et al., 1992)和含水榴石的替代、OH的加入(Rossman and Aines, 1991)并非是导致该区石榴子石不均匀消光的原因,因此六配位上的Al3+-Fe3+有序度(张金民和曹正民, 1988; Griffin et al., 1992)可能影响了铜绿山矿床中石榴子石的非均质性。此外石榴子石电子探针分析得到的Na含量大多数都低于检出限,因而第二种替换形式不太可能形成铜绿山石榴子石的特征。综合看来,铜绿山矿床中石榴子石中的REE可能是通过缺位替代或其他少见的形式而存在的,而这种替代机制与六配位离子的种型存在潜在的联系。

图 8 正三价的稀土元素的在八配位位置上离子半径以及与同在该位置上的Ca2+和Fe2+的半径绝对差值(a, 离子半径值据Shannon, 1976)和铜绿山石英闪长岩体的REE组成以及根据熔体-热液分配系数计算的岩浆脱水脱挥发分形成热液的REE组成(b, 据Flynn and Burnham, 1978) Fig. 8 Plot of ionic radius of trivalent REE in eight-coordination against the absolute percent difference with Ca2+ and Fe2+ in the same co-ordination (a, ionic radius from Shannon, 1976) and REE composition of Tonglushan quartz dirorite and the fluid dehydrated from the magma calculated by the partition coefficient between fluid and melt (b, after Flynn and Burnham, 1978)

不同于高压、超高压变质岩中以铁铝榴石端元为主的石榴子石在球粒陨石标准化后强烈富集HREE的特征(Skublov and Drugova, 2000),铜绿山矿床中的石榴子石主要的端元组成为钙铝榴石和钙铁榴石,均具有相对富集LREE亏损HREE的特点。石榴子石中REE模式受到晶体化学等内部因素和热液组成与形成时的物理化学条件等外部因素的影响(Jamtveit and Hervig, 1994; Smith et al., 2004)。从晶体化学的角度来看,在钙铝榴石和钙铁榴石中,REE3+替代的是八配位上的Ca2+,其中LREE的离子半径与Ca2+更加的接近(图 8a),因而替代更容易发生。而以铁铝榴石端元为主的石榴子石中,八配位的Mn2+、Fe2+、Mg2+的半径要较Ca2+小,HREE的离子半径与Fe2+相差较小,更容易发生类质同象。从热液组成与物理化学条件来看,根据赵海杰等(2010)测定的铜绿山石英闪长岩体的稀土元素组成和Flynn and Burnham (1978)在1.25kbar和800℃的条件下实验岩石学测定的热液与熔体的分配系数,可以估算出铜绿山岩体形成时岩浆作用后期脱水脱挥发分形成热液的REE组成具有富集LREE、亏损HREE的特征(图 8b)。从形成热液到热液与围岩发生接触交代作用形成矽卡岩矿物期间,REE在热液中以络合物的形式(如LaF2+、LaCl2+、LaOH2+等种型)发生迁移,在形成石榴子石时,络合物首先要发生分解释放出REE离子,因而REE络合物种型及其稳定性也影响了石榴子石中REE模式。根据铜绿山石榴子石的流体包裹体的研究(余元昌等, 1985; Zhao et al., 2012),发现其在偏酸性的热液中结晶,形成的温度为350~550℃。Smith et al. (2004)对热液中进行REE络合物的行为进行了理论计算和地球化学模拟,发现在偏酸性条件下LREE络合物的种型是以氯络合物为主(如LaCl2+),而HREE络合物以氟络合物为主(如LuF2+)。在500bar和350℃时(接触交代作用发生的温压条件),氟络合物的稳定性要高于氯络合物,因而在氯络合物发生分解时氟络合物仍然稳定,而氯络合物是LREE重要的载体,因而此时热液中游离的LREE离子浓度要更高,更易被石榴子石生长所利用。

4.2 铜绿山矿床石榴子石的生长模式

矿物从与其化学组成不同的流体中结晶生长是一个多相反应,如石榴子石从热液中结晶至少涉及到两个不同的相,这一反应包括成核、界面反应、物质迁移和热传递这三个步骤(Zhang, 2008)。界面反应是晶核与介质的界面处发生的分子或离子的附着/脱附的过程,界面反应速率可以理解为晶体生长或溶解速率。物质迁移指的是分子或离子等在介质中通过扩散或者对流的方式发生迁移。物质迁移和界面反应是相继进行的,因而速率慢的步骤是决速步,影响着整个多相反应的进程。

铜绿山矿床不同位置产出的石榴子石的岩相学特征存在差异,这种结构的差异与生长速率相关。钙铝榴石中保留有压力侵蚀结构,而钙铁榴石中却没有发现,这种特殊的结构反应了较慢的生长速率(Gaspar et al., 2008; Zhai et al., 2014)。在环带发育的钙铁榴石中,在边部可以观察到凹槽状结构,这种结构Jamtveit et al. (1992)解释为生长速率的逐渐加快,而从核部到边部,钙铁榴石端元组成增加,暗示含Fe高的钙铁榴石的生长速率较快。从核部到边部,{110}面的长度越来越大,根据布拉维法则,面网密度越小的晶面,生长速率越快(赵珊茸等, 2004),说明在结晶过程中{211}面的生长速率要快于{110}面,而高指数{211}面代表了钙铁榴石亚稳态的结构面,指示较快的生长速率(Jamtveit et al., 1993)。钙铁榴石中较多的原生包裹体数目也是其生长速率相对于钙铝榴石较快的一个佐证(Gaspar et al., 2005, 2008)。

在主量元素组成上,钙铝榴石和钙铁榴石也存在着较大的差异:钙铝榴石的端元组成变化很小,而钙铁榴石的端元组成变化要大的多;外矽卡岩中钙铝榴石和钙铁榴石具有成分周期性振荡的环带,而大理岩中的钙铝榴石和内矽卡岩中的钙铁榴石却不发育。振荡环带的形成与热液体系成分周期性的变化相关(Smith et al., 2004),假设在一批次富Fe外部流体注入到下一批次的外部流体注入这一个周期内体系保持封闭,由于富Fe流体注入,使得该体系中的Fe发生过饱和从而结晶出钙铁榴石层,富Fe层的生长使得体系中Fe的过饱和程度降低而Al的过饱和程度增加,从而在钙铁榴石层外结晶出钙铝榴石,Fe、Al的过饱和程度此消彼涨,导致振荡环带的形成。

根据不同位置的石榴子石的结构特征和主量元素差异,建立以下生长动力学模式。值得说明的是尽管我们样品来自三个不同地点的接触带,外矽卡岩及内矽卡岩中相对应的石榴子石(如14TLS56-1与14TLS78)的结构、主量元素及其变化趋势均非常的相似,虽然微量元素的含量存在差别,但是微量元素的比值等指标(LREE/HREE、Eu/Eu*等)却很类似,说明这两个接触带的围岩存在一些差异,但这个差异对于基于结构特征和主量元素差异建立的生长动力学模式影响较小,因此可以将不同产状的石榴子石综合在一起进行讨论。

(1) 矽卡岩化大理岩中的钙铝榴石没有环带特征,表明在其生长的过程中,界面反应的速率要小于物质迁移的速率,否则将会出现从核到边成分逐渐减少的类似于瑞利分馏形成的环带;

(2) 而镁质矽卡岩中的钙铝榴石具有明显的振荡环带,且从核到边Fe的含量上升,在这典型的生长环带的形成过程中,界面反应的速率大于物质迁移的速率,因而热液成分的变化能够记录在环带中,但是其成分变化范围很小,说明物质迁移主要以速率较慢的扩散为主。这与白云岩和灰岩的化学性质相关,ΔrGDol(-2161.3KJ·mol-1) < ΔrGCc(-1128.5KJ·mol-1)(Robie and Hemmingway, 1995),意味着在相同条件下热液更容易与灰岩发生反应,导致发生接触交代作用后物质迁移方式的差别,钙质外矽卡岩以对流为主,而镁质外矽卡岩以扩散为主;

(3) 钙质外矽卡岩中钙铁榴石,其振荡环带非常的发育,而且端元组成变化很大,表明晶体形成时的界面反应速率大于物质迁移速率,而且物质迁移以对流为主;

(4) 内矽卡岩中钙铁榴石端元成分变化较大,该晶体的生长过程中的界面反应速率与物质迁移方式与钙质外矽卡岩中钙铁榴石类似。但出现的是混乱无规律的环带而不是振荡环带,这与热传递过程密切相关。结晶是一个放热的过程,在靠近岩体的位置上,热量不仅迁移不出去反而源源不断的传递过来,因此抑制了结晶的进行,同时热液中的成分以对流的形式发生迁移,物质供应充足,使得先期形成的晶体发生溶解再沉淀形成复杂的再吸收环带。到了岩体侵位晚期,热量耗散,在边部形成少量的振荡环带(图 7g, h)。

铜绿山矿床中不同位置的石榴子石的生长模式如图 9所示。岩浆侵位早期带来的热和少量脱水脱挥发分形成的热液分别使得围岩发生接触热变质和接触交代作用(图 9a),接触热变质使得灰岩变质成大理岩、白云岩变成白云质大理岩,而少量的热液与围岩发生反应,使得体系中的Al发生过饱和,在大理岩中成核并结晶出钙铝榴石(Grs1);随着岩浆的继续上侵、冷却,批次的大量的脱水脱挥发分,随着反应前锋的推进,早期包裹有钙铝榴石的大理岩被消耗完,此时物质迁移的方式由扩散变为对流,同时流体的成分也可能变的更加的富Fe,由于成核相对于界面反应需要更大的过饱和度,因此新生长出的钙铁榴石可能以原先的钙铝榴石为核继续生长(Exo-Adr1)(图 7e),如果Fe的过饱和程度高的话,也可能是新成核进而继续生长(Exo-Adr2)(图 7f);而在白云质大理岩处,由于热液较难与白云岩反应完全,因而会通过扩散的方式形成钙铝榴石,流体成分周期性的改变也会导致振荡环带的形成(MgExo-Grs),由于所处的位置既具有充足的物质又具有合适的热量,此后如果流体成分周期性的变化,则会导致振荡环带的形成;而靠近岩体的位置,尽管具有充足的物质供给,但是由于过多的热量供给限制了结晶过程,而且还会使得原先生长的晶体发生溶解再沉淀,形成复杂的环带(Endo-Adr);在热液向围岩迁移的过程中,存在物质和热量的耗散,到达围岩时可能会像岩浆早期一样,物质通过扩散的方式发生迁移也形成钙铝榴石(Grs2)。

图 9 铜绿山矿床中石榴子石的生长模式图 (a)、(b)分别表示岩体侵位早期、晚期的主要岩石单元以及石榴子石的生长情况 Fig. 9 Growth pattern of garnets from different positions in Tonglushan deposit (a), (b) refer to the distributions of rocks and the conditions of garnets growth in the early and late stage of magma emplacement respectively

此外,上述生长动力学模式也能够很好的反映微量元素指标的差异。钙铁榴石中的化学性质很相似的Nb和Ta的含量几乎一致,而钙铝榴石中出现了Nb和Ta含量脱耦的现象,钙铁榴石形成过程中的界面反应速率很快且大于以对流为主的物质迁移速率,而在物质迁移速率中,对流又较扩散更快。因而关系到钙铝榴石形成的较慢的界面反应速率或者扩散过程,使得性质很相近Nb和Ta元素的活动性差异扩大,进而导致含量的脱耦;钙铁榴石相对于钙铝榴石的LREE/HREE要高,LREE在流体中的活动性要较HREE更强,较快的界面反应速率或物质迁移速率使得原本活动性相差较大的LREE与HREE的差异更大,所以钙铁榴石相对更亏损HREE。因而14TLS56-1钙铁榴石中核部暗层与边部暗层的∑REE相似、LREE/HREE差别较大的特征,反映的是早期核部钙铁榴石生长速率慢于晚期边部。Smith et al. (2004)认为石榴子石中高U的含量反应了较快的界面反应速率,铜绿山钙铁榴石的界面反应速率更快,其U的含量也较高。

本文研究的对象为石英闪长岩与三叠纪的碳酸盐岩发生接触交代作用形成的石榴子石,记录的是早期大规模的岩浆-热液活动事件。没有包括Li et al. (2010)和Xie et al. (2011)等认为铜绿山矿床形成后,后期岩脉的局部热液活动事件的影响。

4.3 对成矿过程的指示

热液中石榴子石的生长要求相应的元素浓度达到高度的过饱和,从而发生成核作用,紧接着通过物质迁移和界面反应实现晶体的生长,因而石榴子石环带的结构和成分能够记录晶体生长时的热液的成分与物理化学条件等信息(Jamtveit and Hervig, 1994)。在岩浆侵位早期,脱水脱挥发分形成的热液与围岩反应后呈现出Al过饱和的状态,因而结晶出大理岩中的钙铝榴石或外矽卡岩中钙铁榴石的核部。随着岩浆批次的脱水脱挥发分,围岩不断被消耗完,反应前锋不断的向前推进,石榴子石不断的结晶,热液在原先Al过饱和的状态和Fe过饱和的状态之间相互转换,在某些热液Fe过饱和程度特别高的位置上,当Al过饱和程度降低时结晶出的富Fe的层其Fe含量不会再发生变化(图 6b),但大多数环带中呈现的是Fe含量不断增加的趋势(图 6c-e),表明在干矽卡岩阶段从早到晚热液中的Fe含量是不断增加的。

从微量元素上来看,Eu2+是弱阳离子,容易与Cl-形成络合物(Allen and Seyfried, 2005),如果热液中Cl的含量逐渐增加,Eu2+的溶解度也随之增加,而成矿金属元素主要是以Cl络合物(氯化亚铁、氯化铁)形式从岩浆运移到热液中,因而热液中Eu2+含量的增加也意味着金属成矿元素含量的增加。钙铁榴石(尤其是外矽卡岩中)具有明显的正Eu异常,表明其形成时的热液具有较高的金属成矿元素浓度。

据此,我们可以推断出铜绿山矽卡岩矿床形成过程中,早矽卡岩阶段从早到晚,热液中的金属成矿元素含量增加;而在早矽卡岩阶段的围岩、外矽卡岩带、内矽卡岩带中,以外矽卡岩带的热液中的成矿元素浓度最高。含矿热液在合适的构造和物理化学条件下,金属成矿元素发生沉淀成矿。

5 结论

(1) 在铜绿山矿床中,接触带不同位置产出不同类型的石榴子石,大理岩中的为不发育环带的钙铝榴石,镁质外矽卡岩中为振荡环带明显的钙铝榴石,钙质外矽卡岩中的是振荡环带非常清晰的钙铁榴石,而内矽卡岩中常为环带混乱的钙铁榴石;

(2) REE在石榴子石晶格中的替代机制还需进一步研究,目前两种主要替代机制在铜绿山矿床中的石榴子石都不适用,可能存在空缺替代或者其他少见并且与六配位上的离子种型相关的替代机制;

(3) 根据结构和主量元素变化,可以建立不同产状的石榴子石的生长动力学模式,大理岩中的钙铝榴石形成过程中的界面反应速率较慢、物质迁移以扩散为主,镁质外矽卡岩中的钙铝榴石形成过程中界面反应速率较快、物质迁移以扩散为主,钙质的外矽卡岩中钙铁榴石的形成过程中的界面反应速率较快、物质迁移以对流为主,而内矽卡岩中的钙铁榴石的形成的动力学条件与钙质外矽卡岩中的相似,但是受到热迁移限制;

(4) 石榴子石能够很好的反演早期成矿热液的组成、性质、时空演化,进而可以对铁质沉淀有所指示,有助于精细的认识矽卡岩型矿床早期成矿过程,对于找矿具有启发性的意义。

致谢      野外地质工作得到了鄂东南地质大队魏克涛和大冶有色的工程师们的帮助;岩石薄片的制作和处理得到陆建培老师的支持;电子探针和LA-ICPMS分析分别得到了杨梅君老师和赵葵东、陈唯老师的支持;论文的撰写过程中还与胡浩、段壮博士进行了讨论;两位审稿专家提出了许多宝贵的意见;在此一并表示衷心感谢。

参考文献
Allen DE and Seyfried Jr WE. 2005. REE controls in ultramafic hosted MOR hydrothermal systems:An experimental study at elevated temperature and pressure. Geochimica et Cosmochimica Acta, 69(3): 675-683. DOI:10.1016/j.gca.2004.07.016
Chernoff CB and Carlson WD. 1999. Trace element zoning as a record of chemical disequilibrium during garnet growth. Geology, 27(6): 555-558. DOI:10.1130/0091-7613(1999)027<0555:TEZAAR>2.3.CO;2
Crowe DE, Riciputi LR, Bezenek S and Ignatiev A. 2001. Oxygen isotope and trace element zoning in hydrothermal garnets:Windows into large-scale fluid-flow behavior. Geology, 29(6): 479-482. DOI:10.1130/0091-7613(2001)029<0479:OIATEZ>2.0.CO;2
D'Errico ME, Lackey JS, Surpless BE, Loewy SL, Wooden JL, Barnes JD, Strickland A and Valley JW. 2012. A detailed record of shallow hydrothermal fluid flow in the Sierra Nevada magmatic arc from low-δ18O skarn garnets. Geology, 40(8): 763-766.
Dziggel A, Wulff K, Kolb J, Meyer FM and Lahaye Y. 2009. Significance of oscillatory and bell-shaped growth zoning in hydrothermal garnet:Evidence from the Navachab gold deposit, Namibia. Chemical Geology, 262(3-4): 262-276. DOI:10.1016/j.chemgeo.2009.01.027
Flynn RT and Burnham CW. 1978. An experimental determination of rare earth partition coefficients between a chloride containing vapor phase and silicate melts. Geochimica et Cosmochimica Acta, 42(6): 685-701.
Gaspar M, Knaack C, Meinert LD and Ottonello G. 2005. REE distribution in grandites from the Crown Jewel gold skarn deposit:A LA-ICP-MS study. Geochimica et Cosmochimica Acta, 69(Suppl.): 78.
Gaspar M, Knaack C, Meinert LD and Moretti R. 2008. REE in skarn systems:A LA-ICP-MS study of garnets from the Crown Jewel gold deposit. Geochimica et Cosmochimica Acta, 72(1): 185-205. DOI:10.1016/j.gca.2007.09.033
Griffin DT, Hatch DM, Phillips WR and Kulaksiz S. 1992. Crystal chemistry and symmetry of a birefringent tetragonal pyralspite 75-grandite 25 garnet. American Mineralogist, 77(3-4): 399-406.
Hu H, Li JW, Lentz D, Ren Z, Zhao XF, Deng XD and Hall D. 2014. Dissolution-reprecipitation process of magnetite from the Chengchao iron deposit:Insights into ore genesis and implication for in-situ chemical analysis of magnetite. Ore Geology Reviews, 57: 393-405. DOI:10.1016/j.oregeorev.2013.07.008
Ismail R, Ciobanu CL, Cook NJ, Teale GS, Giles D, Mumm AS and Wade B. 2014. Rare earths and other trace elements in minerals from skarn assemblages, Hillside iron oxide-copper-gold deposit, Yorke Peninsula, South Australia. Lithos, 184-187: 456-477. DOI:10.1016/j.lithos.2013.07.023
Jamtveit B. 1991. Oscillatory zonation patterns in hydrothermal grossular-andradite garnet:Nonlinear dynamics in regions of immiscibility. American Mineralogist, 76(7-8): 1319-1327.
Jamtveit B and Andersen TB. 1992. Morphological instabilities during rapid growth of metamorphic garnets. Physics and Chemistry of Minerals, 19(3): 176-184.
Jamtveit B, Grorud HF and Bucher-Nurminen K. 1992. Contact metamorphism of layered carbonate-shale sequences in the Oslo Rift. Ⅱ:Migration of isotopic and reaction fronts around cooling plutons. Earth and Planetary Science Letters, 114(1): 131-148. DOI:10.1016/0012-821X(92)90156-P
Jamtveit B and Andersen TB. 1993. Contact metamorphism of layered shale-carbonate sequences in the Oslo Rift; Ⅲ, The nature of skarn-forming fluids. Economic Geology, 88(7): 1830-1849. DOI:10.2113/gsecongeo.88.7.1830
Jamtveit B, Wogelius RA and Fraser DG. 1993. Zonation patterns of skarn garnets:Records of hydrothermal system evolution. Geology, 21(2): 113-116. DOI:10.1130/0091-7613(1993)021<0113:ZPOSGR>2.3.CO;2
Jamtveit B and Hervig RL. 1994. Constraints on transport and kinetics in hydrothermal systems from zoned garnet crystals. Science, 263(5146): 505-508. DOI:10.1126/science.263.5146.505
Jamtveit B, Agnarsdottir KV and Wood BJ. 1995. On the origin of zoned grossular-andradite garnets in hydrothermal systems. European Journal of Mineralogy, 7(6): 1399-1410. DOI:10.1127/ejm/7/6/1399
Kohn MJ. 2014. Geochemical zoning in metamorphic minerals. In: Rudnick RL (ed.). Treatise on Geochemistry. 2nd Edition. Oxford: Elsevier, 4: 249-280
Lackey JS, Romero GA, Bouvier AS and Valley JW. 2012. Dynamic growth of garnet in granitic magmas. Geology, 40(2): 171-174. DOI:10.1130/G32349.1
Li JW, Zhao XF, Zhou MF, Ma CQ, de Souza ZS and Vasconcelos P. 2009. Late Mesozoic magmatism from the Daye region, eastern China:U-Pb ages, petrogenesis, and geodynamic implications. Contributions to Mineralogy and Petrology, 157(3): 383-409. DOI:10.1007/s00410-008-0341-x
Li JW, Deng XD, Zhou MF, Liu YS, Zhao XF and Guo JL. 2010a. Laser ablation ICP-MS titanite U-Th-Pb dating of hydrothermal ore deposits:A case study of the Tonglushan Cu-Fe-Au skarn deposit, SE Hubei Province, China. Chemical Geology, 270(1-4): 56-67. DOI:10.1016/j.chemgeo.2009.11.005
Li JW, Vasconcelos PM, Zhou MF, Deng XD, Cohen B, Bi SJ, Zhao XF and Selby D. 2014. Longevity of magmatic-hydrothermal systems in the Daye Cu-Fe-Au district, eastern China with implications for mineral exploration. Ore Geology Reviews, 57: 375-392. DOI:10.1016/j.oregeorev.2013.08.002
Li SG, Jagoutz E, Chen YZ and Li QL. 2000. Sm-Nd and Rb-Sr isotopic chronology and cooling history of ultrahigh pressure metamorphic rocks and their country rocks at Shuanghe in the Dabie Mountains, Central China. Geochimica et Cosmochimica Acta, 64(6): 1077-1093. DOI:10.1016/S0016-7037(99)00319-1
Li XH, Li WX, Li QL, Wang XC, Liu Y and Yang YH. 2010b. Petrogenesis and tectonic significance of the~850Ma Gangbian alkaline complex in South China:Evidence from in situ zircon U-Pb dating, Hf-O isotopes and whole-rock geochemistry. Lithos, 114(1-2): 1-15. DOI:10.1016/j.lithos.2009.07.011
Liu YS, Hu ZC, Zong KQ, Gao CG, Gao S, Xu J and Chen HH. 2010. Reappraisement and refinement of zircon U-Pb isotope and trace element analyses by LA-ICP-MS. Chinese Science Bulletin, 55(15): 1535-1546. DOI:10.1007/s11434-010-3052-4
Locock AJ. 2008. An Excel spreadsheet to recast analyses of garnet into end-member components, and a synopsis of the crystal chemistry of natural silicate garnets. Computers & Geosciences, 34(12): 1769-1780.
McIntire W. 1963. Trace element partition coefficients:A review of theory and applications to geology. Geochimica et Cosmochimica Acta, 27(12): 1209-1264. DOI:10.1016/0016-7037(63)90049-8
Meinert LD, Dipple GM and Nicolescu S. 2005. World skarn deposits. In: Hedenquist JW, Thompson JFH, Goldfarb RJ and Richards JP (eds.). Economic Geology 100th Anniversary Volume. Littleton, Colorado: Society of Economic Geologists, Inc., 299-336
Peng HJ, Zhang CQ, Mao JW, Santosh M, Zhou YM and Hou L. 2015. Garnets in porphyry-skarn systems:A LA-ICP-MS, fluid inclusion, and stable isotope study of garnets from the Hongniu-Hongshan copper deposit, Zhongdian area, NW Yunnan Province, China. Journal of Asian Earth Sciences, 103: 229-251. DOI:10.1016/j.jseaes.2014.10.020
Robie RA and Hemingway BS. 1995. Thermodynamic properties of minerals and related substances at 298.15K and 1bar (105 pascals) pressure and at higher temperatures. United States Geological Survey Bulletin, No. 2131
Rossman GR and Aines RD. 1991. The hydrous components in garnets:Grossular-hydrogrossular. American Mineralogist, 76(7-8): 1153-1164.
Shannon RD. 1976. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides. Acta Crystallographica Section A, 32(5): 751-767. DOI:10.1107/S0567739476001551
Shore M and Fowler AD. 1996. Oscillatory zoning in minerals:A common phenomenon. The Canadian Mineralogist, 34(6): 1111-1126.
Shu QA, Chen PL and Cheng JR. 1992. The Geology of Iron and Copper Deposits in Eastern Hubei Province. Beijing: Metallurgical Industry Press: 1-532.
Skublov SG and Drugova GM. 2000. REE distribution in metamorphic garnets. Herald DGGGMS RAS, 2(5): 60-61.
Smith MP, Henderson P, Jeffries TER, Long J and Williams CT. 2004. The rare earth elements and uranium in garnets from the Beinn an Dubhaich Aureole, Skye, Scotland, UK:Constraints on processes in a dynamic hydrothermal system. Journal of Petrology, 45(3): 457-484. DOI:10.1093/petrology/egg087
Somarin AK. 2004. Garnet composition as an indicator of Cu mineralization:Evidence from skarn deposits of NW Iran. Journal of Geochemical exploration, 81(1-3): 47-57. DOI:10.1016/S0375-6742(03)00212-7
Sun SS and McDonough W. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. In: Saunders AD and Norry MJ (eds.). Magmatism in the Ocean Basins. Geological Society, London, Special Publications, 42(1): 313-345
Wang YB, Gu XX, Zhang ZB and Dong SY. 2011. Isotopic geochemistry of the Tonglushan Cu-Fe skarn type deposit in Hubei Province. Geoscience, 25(4): 730-739.
Xie GQ, Mao JW, Li RL, Qu WJ, Pirajno F and Du AD. 2007. Re-Os molybdenite and Ar-Ar phlogopite dating of Cu-Fe-Au-Mo (W) deposits in southeastern Hubei, China. Mineralogy and Petrology, 90(3-4): 249-270. DOI:10.1007/s00710-006-0176-y
Xie GQ, Mao JW, Li RL and Bierlein FP. 2008. Geochemistry and Nd-Sr isotopic studies of Late Mesozoic granitoids in the southeastern Hubei Province, Middle-Lower Yangtze River belt, Eastern China:Petrogenesis and tectonic setting. Lithos, 104(1-4): 216-230. DOI:10.1016/j.lithos.2007.12.008
Xie GQ, Zhao HJ, Zhao CS, Li XQ, Hou KJ and Pan HJ. 2009. Re-Os dating of molybdenite from Tonglushan ore district in southeastern Hubei Province, Middle-Lower Yangtze River belt and its geological significance. Mineral Deposits, 28(3): 227-239.
Xie GQ, Mao JW, Zhao HJ, Wei KT, Jin SG, Pan HJ and Ke YF. 2011. Timing of skarn deposit formation of the Tonglushan ore district, southeastern Hubei Province, Middle-Lower Yangtze River Valley metallogenic belt and its implications. Ore Geology Reviews, 43(1): 62-77. DOI:10.1016/j.oregeorev.2011.05.005
Xie GQ, Mao JW, Zhao HJ, Duan C and Yao L. 2012. Zircon U-Pb and phlogopite 40Ar-39Ar age of the Chengchao and Jinshandian skarn Fe deposits, Southeast Hubei Province, Middle-Lower Yangtze River Valley metallogenic belt, China. Mineralium Deposita, 47(6): 633-652. DOI:10.1007/s00126-011-0367-2
Xie GQ, Mao JW, Zhu QQ, Yao L, Li YH, Li W and Zhao HJ. 2015. Geochemical constraints on Cu-Fe and Fe skarn deposits in the Edong district, Middle-Lower Yangtze River metallogenic belt, China. Ore Geology Reviews, 64: 425-444. DOI:10.1016/j.oregeorev.2014.08.005
Xie GQ, Mao JW, Li W, Zhu QQ, Liu HB, Jia GH, Li YH, Li JJ and Zhang J. 2016. Different proportion of mantle-derived noble gases in the Cu-Fe and Fe skarn deposits:He-Ar isotopic constraint in the Edong district, eastern China. Ore Geology Reviews, 72: 343-354. DOI:10.1016/j.oregeorev.2015.08.004
Zhai DG, Liu JJ, Zhang HY, Wang JP, Su L, Yang XA and Wu SH. 2014. Origin of oscillatory zoned garnets from the Xieertala Fe-Zn skarn deposit, northern China:In situ LA-ICP-MS evidence. Lithos, 190-191: 279-291. DOI:10.1016/j.lithos.2013.12.017
Zhai YS, Yao SZ and Lin XD. 1992. Metallogeny of Fe-Cu-(Au) ore Deposits in Mid-Lower Changjiang. Beijing: Geological Publishing House: 1-120.
Zhang JM and Cao ZM. 1988. An investigation of brie fringence of garnets from a skarn deposite. Scientia Geologica Sinica, (3): 229-238.
Zhang YX. 2008. Geochemical Kinetics. New Jersey: Princeton University Press: 325-441.
Zhao HJ, Mao JW, Xiang JF, Zhou ZH, Wei KT and Ke YF. 2010. Mineralogy and Sr-Nd-Pb isotopic compositions of quartz diorite in Tonglushan deposit, Hubei Province. Acta Petrologica Sinica, 26(3): 768-784.
Zhao HJ, Xie GQ, Wei KT and Ke YF. 2012. Mineral compositions and fluid evolution of the Tonglushan skarn Cu-Fe deposit, SE Hubei, east-central China. International Geology Review, 54(7): 737-764. DOI:10.1080/00206814.2011.569418
Zhao SR, Bian QJ and Ling QC. 2004. Crystallography and Mineralogy. Beijing: Higher Education Press: 313-326.
Zhou TF, Fan Y, Yuan F and Zhong GX. 2012. Progress of geological study in the Middle-Lower Yangtze River Valley metallogenic belt. Acta Petrologica Sinica, 28(10): 3051-3066.
Zhu QQ, Xie GQ, Li W, Zhang F, Wang J, Zhang P and Yu BF. 2014. In situ analysis of garnets from the Jingshandian iron skarn deposit, Hubei Province, and its geological implications. Geology in China, 41(6): 1944-1963.
舒全安, 陈培良, 程建荣. 1992. 鄂东铁铜矿产地质. 北京: 冶金工业出版社: 1-532.
王彦博, 顾雪祥, 张宗保, 董树义. 2011. 湖北铜绿山矽卡岩型铜铁矿床同位素地球化学研究. 现代地质, 25(4): 730-739. DOI:10.3969/j.issn.1000-8527.2011.04.014
谢桂青, 赵海杰, 赵财胜, 李向前, 侯可军, 潘怀军. 2009. 鄂东南铜绿山矿田矽卡岩型铜铁金矿床的辉钼矿Re-Os同位素年龄及其地质意义. 矿床地质, 28(3): 227-239. DOI:10.3969/j.issn.0258-7106.2009.03.001
翟裕生, 姚书振, 林新多. 1992. 长江中下游地区铁铜(金)成矿规律. 北京: 地质出版社: 1-120.
张金民, 曹正民. 1988. 夕卡岩中石榴石的双折射性之研究. 地质科学, (3): 229-238.
赵海杰, 毛景文, 向君峰, 周振华, 魏克涛, 柯于富. 2010. 湖北铜绿山矿床石英闪长岩的矿物学及Sr-Nd-Pb同位素特征. 岩石学报, 26(3): 768-784.
赵珊茸, 边秋娟, 凌其聪. 2004. 结晶学及矿物学. 北京: 高等教育出版社: 313-326.
周涛发, 范裕, 袁峰, 钟国雄. 2012. 长江中下游成矿带地质与矿产研究进展. 岩石学报, 28(10): 3051-3066.
朱乔乔, 谢桂青, 李伟, 张帆, 王建, 张平, 于炳飞. 2014. 湖北金山店大型矽卡岩型铁矿石榴子石原位微区分析及其地质意义. 中国地质, 41(6): 1944-1963. DOI:10.3969/j.issn.1000-3657.2014.06.012