土耳其西部新生代斑岩成矿岩浆氧逸度-挥发份元素组成演化特征*——来自锆石和磷灰石组成约束

赵俊兴 苏本勋 UYSAL Ibrahim AYDIN Faruk 肖燕 SEN Cuneyt回凯旋, 3 秦克章, 3

1. 中国科学院矿产资源研究重点实验室，中国科学院地质与地球物理研究所，北京 100029

2. 中国科学院地球科学研究院，北京 100029

3. 中国科学院大学地球与行星科学学院，北京 100049

4. Department of Geological Engineering, Karadeniz Technical University, Trabzon 61080

5. 岩石圈演化国家重点实验室，中国科学院地质与地球物理研究所，北京 100029

图1 特提斯土耳其段主要斑岩型矿床分布图(据修改)其中土耳其西部主要为四个成矿阶段：1-俯冲阶段；2-碰撞阶段(Tavanl带)；3-晚碰撞-后碰撞阶段(Biga半岛成矿带)；4-后碰撞伸展阶段(Afyon-Konya成矿带);成矿带划分来自Kucu et al. (2019). 本次研究选择2-4区域内典型矿床开展研究Fig.1 Distribution map of main porphyry deposits of Tethys Turkish section (modified after 2009)The western Turkish is mainly composed of four ore-forming stages: 1-subduction stage; 2-collision stage (Tavanl belt); 3-from late-collision to post-collision stage (Biga Peninsula metallogenic belt); 4-post-collision extension stage (Afyon-Konya metallogenic belt); the division of metallogenic belt refers to Kucu et al. (2019). The typical deposits of 2～4 region were selected to research in this study

1 土耳其西部斑岩矿床构造背景

2 土耳其西部新生代典型斑岩成矿系统基本特征

2.1 Tavanl带始新世斑岩成矿系统

2.2 Biga半岛始新世至渐新世斑岩成矿系统

2.3 Afyon-Konya中新世-上新世斑岩金-钼-铜矿及浅成低温成矿系统

3 研究对象和测试方法

锆石和磷灰石在分选后制靶，并拍摄透反射光照片和选择无裂缝和矿物-流体包裹体颗粒区域，以便后续开展矿物主微量元素测试工作。锆石和磷灰石的阴极发光图像主要在中国科学院地质与地球物理研究所使用配置MonoCL4的Nova NanoSEM 450场发射扫描电子显微镜(FSEM)完成。锆石微量元素含量中国科学院地质与地球物理研究所岩石圈演化重点实验室的Agilent 7500a四级杆电感耦合等离子体质谱仪上完成(与193nm ArF准分子激光剥蚀仪联机)。每个时间分辨分析数据包括大约20～30s的空白信号和50s的样品信号。对分析数据的离线处理(包括对样品和空白信号的选择、仪器灵敏度漂移校正、元素含量及U-Th-Pb同位素比值和年龄计算)采用软件ICPMSDataCal(Liuetal., 2008)完成。详细的仪器操作条件和数据处理方法同Liuetal.(2008, 2010)。锆石微量元素选点主要为先前锆石U-Pb定年测试点，激光束斑选择为44～60μm，略大于之前的年龄测定束斑。锆石微量元素含量利用NIST610、USGS参考玻璃BCR-2G与BIR-1G作为多外标、Zr作内标的方法进行定量计算(Liuetal., 2010)。锆石微量元素结果见于表2。

表2 土耳其西部新生代斑岩矿床中岩浆岩锆石LA-ICPMS微量元素数据(×10-6)Table 2 The analytical data of LA-ICPMS zircon trace elements from magmatic rocks in the Cenozoic porphyry ore deposits, western Turkey (×10-6)

续表2Continued Table 2

磷灰石主量元素在中国科学院地质与地球物理研究所(IGGCAS)使用JEOL-JXA8100电子探针完成。本次研究工作条件为加速电压15kV，束流10nA，探针束斑大小为3～5μm，Si、Fe、Mn、Mg、Na、S和Cl的计数时间为20s，F的计数时间为40s，Ca和P的计数时间为10s。选点过程中我们只将感兴趣的磷灰石颗粒在电子束下进行移动，以防止F、Cl、Ca和P成分的迁移。Stormeretal.(1993)指出这些元素在EPMA分析中的Kα-X射线计数率与电子束的照射时间密切相关。所有数据均采用ZAF进行校正，所有结果见表3。

表3 土耳其西部新生代斑岩矿床中岩浆岩磷灰石EPMA主量元素数据(wt%)Table 3 The analytical data of EPMA apatite major elements from magmatic rocks in the Cenozoic porphyry ore deposits, western Turkey (wt%)

续表3Continued Table 3

续表3Continued Table 3

4 测试结果

4.1 锆石微量元素组成

本次研究共测定333颗锆石微量元素数据，结合全岩数据(赵俊兴等，未发表数据)和Lietal.(2019)的Geo-fO2来计算锆石钛饱和温度和氧逸度(表2)。结果显示土耳其西部新生代斑岩矿床的成矿岩浆具有如下的特征：(1)锆石微量元素表明土耳其西部新生代斑岩矿床成矿岩浆大部分落于高水含量-高氧逸度区间(图2)；(2)岩浆岩锆石具有相似的稀土元素配分曲线(图3)，且暗色包体锆石较同期岩体有较高的稀土总量(图3a, d)；(3)绝大部分锆石的EuN/EuN*位于0.4～0.7之间，斑岩Mo和斑岩Mo-Cu矿相关锆石在结晶时熔体受斜长石结晶影响更大，具有更低EuN/EuN*比值(图4a)；(4)利用Ferry and Watson (2007)的锆石Ti饱和温度计(假设TiO2活度为0.7)估算土耳其西部斑岩成矿及其相关岩体的形成温度在650～900℃之间，锆石Ti温度显示岩体结晶温度表现为：斑岩Au矿>斑岩Cu-Mo=斑岩Mo-Cu矿>斑岩Mo矿，呈现递减趋势 (图4b)；(5)锆石微量元素显示(Yb/Gd和Hf/Y-Th/U，图4c, d)锆石在岩浆房结晶过程中，其微量元素组成大部分受角闪石±榍石±磷灰石分离结晶控制，该分离结晶控制趋势与全岩微量元素指示演化受角闪石结晶相控制一致(图5)；(6)斑岩Au和斑岩Cu-Mo相关岩体锆石形成时氧逸度高于斑岩Mo和斑岩Mo-Cu矿岩体，结合各岩体形成年龄，锆石氧逸度Ce4+/Ce3+比值呈现出先降低(斑岩Cu-Mo到斑岩Mo)后升高(斑岩Mo-Cu到斑岩Au)的趋势(图4d)。

图2 锆石(Ce/Nd)/Y-10000×(Eu/Eu*)/Y图解其中锆石微量元素组成圈定高氧逸度-高含水量成矿岩体(△FMQ>2和H2O<12%)与低氧逸度-低含水量(△FMQ<2和H2O<5%)非成矿岩体范围来自Lu et al. (2016)Fig.2 Zircon (Ce/Nd)/Y vs. 10000×(Eu/Eu*)/Y diagramThe zircon trace element ranges to define high oxygen fugacity-high water content fertile magma (△FMQ>2 and H2O<12%) and low oxygen fugacity-low water content (△FMQ<2 and H2O<5%) infertile suite refer to Lu et al. (2016)

图3 土耳其西部新生代斑岩矿床岩浆岩锆石球粒陨石标准化稀土元素配分图解(标准化值据Sun and McDonough, 1989)Fig.3 Chondrite-normalization REE patterns of western Turkish Neozoic porphyry deposit magmatic rock zircon (normalization values after Sun and McDonough, 1989)

图4 锆石Hf-EuN/EuN*(a)、Hf-锆石Ti饱和温度(b)、Th/U-Yb/Gd(c)和Th/U-Hf/Y(d)之间演化关系图a中氧化条件演化趋势来自Dilles et al. (2015)，图c和d中分离结晶趋势来自Lee et al. (2021). 图例见图2Fig.4 The evolution relationship of zircon Hf vs. EuN/EuN* (a), Hf vs. zircon Ti saturation temperature (b), Th/U vs. Yb/Gd (c) and Th/U vs. Hf/Y (d)The evolution trend of oxidizing condition in Fig.4a refers to Dilles et al. (2015), the fractional crystallization trend in Fig.4c, d refers to Lee et al. (2021). Legends refer to Fig.2

图5 土耳其西部斑岩矿床中锆石年龄和锆石Th/U、Ti温度、EuN/EuN*和Ce4+/Ce3+比值之间演化关系Fig.5 The evolution relationship of zircon age and zircon Th/U, Ti-in-zircon temperature, EuN/EuN* and Ce4+/Ce3+ ratio in the porphyry ore deposits, western Turkey

4.2 磷灰石主量元素组成

图6 斑岩矿床中磷灰石中F含量-F/Cl比值(a)和Cl含量-F/Cl比值(b)之间的关系Fig.6 The relationship of apatite F content vs. F/Cl ratios (a) and Cl content vs. F/Cl ratios (b) in porphyry deposits

5 讨论

5.1 土耳其西部新生代斑岩系统岩浆氧逸度演化趋势

通过斑岩成矿带的构造演化历史和岩浆岩微量元素-同位素研究，能够获得俯冲-碰撞过程中地壳属性(Sillitoe, 2018)、地壳厚度演化(如, (La/Yb)N比值, Profetaetal., 2015; Zhuetal., 2017; Y/Sr比值, Dengetal., 2018; 岩浆金属元素含量和俯冲带地壳厚度关系, Chiaradia, 2014)、斑岩成矿带的岩石圈结构(如, 同位素填图, Houetal., 2015b; 地幔或下地壳包体研究, Wangetal., 2016, Xuetal., 2017)、俯冲板块边界信息(如, 配合板块汇聚速度和地壳缩短估算, Wangetal., 2015)等宏观要素。

本次研究得到土耳其西部新生代斑岩系统的岩浆氧逸度自晚古新世-早始新世碰撞阶段到后碰撞环境经历先降低后升高的过程(图5d，Ce4+/Ce3+比值的变化关系)。该特征与西藏冈底斯碰撞环境斑岩矿床(Wangetal., 2014a, b; Yangetal., 2016)和安第斯型俯冲环境斑岩矿床(Caoetal., 2016; Lietal., 2017; Angereretal., 2018)的氧逸度演化趋势不同。西藏冈底斯碰撞环境和安第斯俯冲型成矿岩浆演化趋势表明为成矿岩浆含水量和氧逸度均显著升高，地壳加厚过程可能是控制其氧逸度升高的重要因素(Tangetal., 2020)。因此，估算土耳其西部斑岩矿床成矿岩体及其同期岩浆岩形成时地壳厚度，并结合锆石氧逸度评价地壳厚度演化过程对岩浆氧逸度的控制过程，是解释土耳其西部新生代斑岩系统岩浆氧逸度演化趋势的重要研究内容。

本次研究使用Huetal.(2017)在碰撞造山带岩浆岩研究中提出的地壳厚度估算公式。Huetal.(2017)利用全岩Sr/Y和(La/Yb)N比值以及Crust 1.0模型获得碰撞带地壳厚度(H)计算公式估算地壳厚度：H=0.67×Sr/Y+28.21和H=27.78×ln[0.34×(La/Yb)N]，并指出与前人在俯冲岛弧带岩浆岩中获得的地壳厚度估算模型相比(Chiaradia, 2014; Profetaetal., 2015)，在相同地壳厚度下碰撞环境的中酸性岩浆岩(SiO2含量为55%～72%)具有较低的Sr/Y和较高的 (La/Yb)N比值(Huetal., 2017)。我们结合前人在土耳其西部地区发表的新生代岩浆岩数据(选择SiO2含量在55%～72%)，使用该模型对土耳其西部50Ma以来的斑岩矿床及同期岩浆岩形成时地壳厚度进行估算(图7)。估算结果表明三条成矿带内成矿与非成矿岩体尽管估算值略有不同，但随时间演化趋势基本一致，即Tavanl带和Biga半岛成矿带在始新世到渐新世岩浆岩形成时地壳厚度估算值略微升高，Afyon-Konya成矿带早中新世地壳依然有与Biga半岛成矿带相似的估算厚度。但到中晚中新世Afyon-Konya成矿带的地壳厚度明显增大。而针对本次研究所选择的样品，全岩的(La/Yb)N和估算地壳厚度(图8a, b)显示从始新世到早中新世有着地壳逐渐加厚的趋势。而在渐新世-中新世的Tepeoba斑岩Mo矿、Pnarba斑岩Mo-Cu矿和斑岩Au矿中，从成矿前到成矿期岩浆演化过程明显经历地壳减薄事件(图8a, b，随时间明显降低的(La/Yb)N)。在该过程中Pnarba斑岩Mo-Cu矿和斑岩Au矿的岩浆氧逸度明显上升(图8c)，说明与地壳减薄相关的深部过程可能是控制土耳其西部中新世高氧逸度成矿岩浆形成的重要机制。同时，全岩微量元素(图9a，b)和锆石微量元素(图4d)指示，土耳其西部斑岩矿床成矿岩浆地球化学特征可能主要受到以角闪石为主的分离结晶作用控制。因此，在该地壳减薄过程中可能发生基性岩浆注入过程，使得成矿期岩浆氧逸度升高。

图7 土耳其西部斑岩矿床岩体及同时期岩浆岩地壳厚度估算与形成时代演化图图中地壳厚度估算公式来自Hu et al. (2017)，并分别利用全岩(La/Yb)N(a)和Sr/Y比值(b)进行估算，其中Tavanl带岩体数据来自Altunkaynak (2007)，Kucu et al. (2019)和Güraslan and Altunkaynak (2019)；Biga半岛成矿带数据来自Kucu et al. (2019)；Afyon-Konya成矿带数据来自 et al. (2015), Deliba et al. (2017)和Kucu et al.(2019)Fig.7 The evolution diagrams of western Turkish porphyry deposits rock and contemporary magmatic rock crustal thickness estimation and formation ageThe computational formula of crustal thickness refers to Hu et al. (2017), whole rock (La/Yb)N (a) and Sr/Y (b) ratio was used to calculate respectively, the Tavanl belt rock values refer to Altunkaynak (2007), Kucu et al. (2019) and Güraslan and Altunkaynak (2019); Biga Peninsula metallogenic belt values refer to Kucu et al. (2019); Afyon-Konya metallogenic belt values refer to et al. (2015), Deliba et al. (2017) and Kucu et al. (2019)

图8 斑岩矿床中各类岩体锆石年龄-(La/Yb)N(a)、锆石年龄-地壳估算厚度(b)和岩浆氧逸度-地壳估算厚度(c)图解△FMQ通过锆石计算获得，以平均值和1StDev表示Fig.8 Magmatic zircon age vs. whole-rock (La/Yb)N (a), zircon age vs. crustal estimation thickness (b), magmatic oxygen fugacity (vs. crustal estimation thickness diagrams△FMQ was obtained by zircon calculation, and showed in terms of average and 1StDev

图9 本次研究的斑岩矿床中各类岩体SiO2-Dy/Yb(a)和SiO2-(La/Yb)N(b)图解其中角闪石和石榴石演化趋势来自Davidson et al. (2007)，全岩数据来自赵俊兴等(未发表数据)Fig.9 Different types of magmatic rock SiO2 vs. Dy/Yb (a) and SiO2 vs. (La/Yb)N (b) diagrams in this studyThe evolution trends of amphibole and garnet refer to Davidson et al. (2007). The whole-rock geochemical compositions are cited from Zhao et al. (unpublished data)

5.2 岩浆中挥发份F-Cl-S演化特征——来自磷灰石组成的约束

斑岩型矿床的金属元素主要以Cl、S络合物的形式搬运，以金属硫化物的形式沉淀。如Cu与Au在熔体中低氧逸度条件下以氯络合物形式迁移，金在高氧逸度条件下以硫氢络合物形式迁移，岩浆系统高氯含量有助于提高铜和金的溶解度(Zajacz and Halter, 2009)。而从斑岩铜矿到斑岩钼矿，岩体酸性程度常与磷灰石F/Cl比值成正比(Mathez and Webster, 2005; 陈雷等, 2017; Zhaoetal., 2020)。因此，研究斑岩矿床成矿物质富集-演化过程，有助于理解挥发分F、Cl、S以及成矿元素本身的地球化学行为(Richards, 2015b)。鉴于岩浆磷灰石和共存熔体处于平衡状态，可使用磷灰石和平衡熔体实验分配系数从磷灰石成分估算在给定P-T条件下熔体的F和Cl含量(Stocketal., 2016)。近期研究表明，在磷灰石和熔体之间存在F和Cl的非能斯特分配，在确定熔体中F、Cl和/或H2O的浓度时，各种实验结果和热力学模型利用磷灰石和熔体之间不同的F和Cl交换系数(KD)进行估算(Pan and Fleet, 2002; Rikeretal., 2018; Li and Costa, 2020)。然而，该方法要求已知熔体中的三种成分之一(F或Cl或H2O含量)，以便使用交换系数确定岩浆的挥发分及水含量。但由于熔体中的F、Cl、H2O含量常需要利用较难找寻的熔体包裹体或火山玻璃进行测定，该类研究对象较难找寻。因此，在本研究中，我们使用卤族元素在磷灰石与熔体间的分配系数值来估算岩浆的F和Cl含量。考虑到岩浆成分(全岩成分熔体的摩尔A/CNK和N/NK)和估计的温度-压力条件，我们使用磷灰石与英安质-流纹质熔体之间的Cl和F分配系数(A/CNK为0.90～1.05，N/NK为0.60～0.77，Websteretal., 2009; Dohertyetal., 2014)，其所代表的压力和温度范围是50～200MPa和900～924℃。长英质硅酸盐熔体的F浓度方程式为(Websteretal., 2009)：CF(%)=[(磷灰石中的XF)-0.12]/3.02，其中XF表示磷灰石中的摩尔F含量。对于熔体中Cl含量的计算，我们使用Baoetal.(2016)的公式，该公式利用了Websteretal.(2009)和Dohertyetal.(2014)数据提供了一个简单的基于权重的能斯特分配系数方法：CCl(熔体，%)=CCl(磷灰石，%)×0.16052。我们利用该公式估算熔体Cl含量。对于岩浆的硫含量，岩浆磷灰石中SO3含量主要受温度、氧化状态和硅酸盐熔体中硫含量的控制(Pengetal., 1997; Parat and Holtz, 2005)。考虑到熔体中硫含量的变化，早期Pengetal.(1997)利用与温度相关的磷灰石-熔体分配系数，开展岩浆相对硫含量的估算，公式为ln(S磷灰石/熔体(%))=21130/T(单位为开氏度)-16.2。Paratetal.(2011)总结前人发表的实验数据，给出了SO3磷灰石(%)=0.157×lnSO3熔体(%)+0.9834(r2=0.68)。我们假设各岩浆磷灰石从硫酸盐饱和熔体结晶，利用Paratetal.(2011)的方法估算岩浆硫含量。虽然两种计算方法会获得不同的岩浆硫含量，但仍可判断其总体演化趋势。

图10 土耳其西部斑岩矿床中磷灰石的F-Cl含量(a)，及锆石年龄和对应岩体磷灰石估算的岩浆F(b)、Cl(c)和S(d)含量的演化关系其中岩浆中氟含量计算来自Webster et al. (2009)，岩浆中氯含量计算来自Webster et al. (2009)和Doherty et al. (2014)，岩浆中硫含量计算来自Parat et al. (2011)，图例见图7Fig.10 The evolution relationship of apatite F-Cl content (a), zircon age and magmatic F (b), Cl (c), S (d) content calculated from apatite in the corresponding rock in western Turkish porphyry depositThe calculation of F content in magmatic refers to Webster et al. (2009), Cl content calculation refers to Webster et al. (2009) and Doherty et al. (2014), S content calculation refers to Parat et al. (2011). Legends refer to Fig.7

斑岩成矿岩浆是否富硫，目前主要利用岩体中的岩浆矿物(岩浆硬石膏和磷灰石)和同期火山岩中熔融包裹体开展研究。如在成矿相关的斑岩中发现岩浆硬石膏，表明其形成时岩浆的氧逸度至少大于NNO+0.5，岩浆中硫含量为0.05%～1%(Carrol and Rutherford, 1987; Paratetal., 2011; Xiaoetal., 2012)。而利用高硫含量的磷灰石和硫在磷灰石和熔体间的分配系数，计算平衡岩浆中的硫含量可达1000×10-6以上(Streck and Dilles, 1998; Liaghat and Tosdal, 2008)，而通过矿床内未发生硫化物饱和的基性端元中熔体包裹体测定硫含量为100×10-6～900×10-6(Seoetal., 2011; Grondahl and Zajacz, 2017; Zhang and Audétat, 2017)，属于正常的岩浆硫含量范围。由于高氧逸度使得岩浆中的硫以+6价形式赋存，并记录在结晶的磷灰石中(Paratetal., 2011)，因此前人在俯冲与碰撞环境斑岩成矿带区域对比研究中发现成矿期磷灰石常具有较高的SO3含量。如中亚造山带哈萨克斯坦地区斑岩铜矿研究(Caoetal., 2016)指出，俯冲作用可以加厚新生下地壳和岩石圈地幔，并使楔形地幔氧逸度不断升高，较高的氧逸度使得成矿期岩石磷灰石较成矿前具有明显较高的SO3含量。西藏班公湖-怒江带俯冲期斑岩铜金矿(Lietal., 2021)和冈底斯后碰撞环境斑岩铜钼矿(Tangetal., 2020)中，岩浆的高氧逸度使得成矿岩体中岩浆磷灰石具有较高的硫含量。上述成矿带的研究与我们在Tavanl带、Biga半岛斑岩成矿带和Afyon-Konya成矿带斑岩矿床中观察到高氧逸度(图5d)条件下的磷灰石高硫组成相一致。表明在俯冲、碰撞-后碰撞和伸展环境的斑岩矿床形成过程中，成矿期高氧逸度岩浆中的磷灰石常具有较高的SO3含量。

5.3 触发土耳其西部斑岩矿床形成的动力学机制——与新特提斯单向汇聚过程有关

根据斑岩矿床的形成时代、成矿岩浆的性质及与区内(或重大)地质事件的时空耦合关系，可以表征斑岩矿床形成时的动力学机制，进而约束板块之间俯冲、碰撞和陆内伸展等一系列构造过程及细节(Tosdal and Richards, 2001)。例如在洋-陆或洋-洋俯冲环境中，若成矿有关的埃达克质岩石，与高镁安山岩/富铌玄武岩、A型花岗岩具有密切的时空关系，则其形成可能与洋脊俯冲有关(Cookeetal., 2005)。依据安第斯位于平板俯冲正上方或缓俯冲向陡俯冲过渡区域的众多超大型斑岩铜矿建立的“平板俯冲模型”(Kay and Mpodozis, 2001; Reichetal., 2003)，斑岩铜矿的形成可能与俯冲带结构和俯冲板片的几何形态有关，如俯冲板片上存在洋底高原(Gutscheretal., 2000; Rosenbaumetal., 2005)，或俯冲板片上的转换断层(Richards and Holm, 2013)均可触发斑岩铜矿的形成。俯冲带极性反转(Solomon, 1990)或者碰撞使得俯冲作用中止的阶段(Sillitoe, 1997)，能够在陆缘弧后背景下形成富金斑岩铜矿。在增生和碰撞造山带内，斑岩矿床形成的时空分布特征可以指示造山带增生及小陆块拼合过程(如中亚造山带西段，秦克章, 2000; Wanetal., 2011, 2017)。通过斑岩矿床的形成时代和区域上成岩-成矿时代朝缝合带方向变年轻的现象，可能说明深部发生俯冲板片回转过程(slab roll-back，如保加利亚地区俯冲相关矿床, Von Quadtetal., 2005; 西藏冈底斯主碰撞阶段斑岩矿床, Zhaoetal., 2014)。在碰撞型斑岩铜矿则集中分布在后碰撞阶段(Houetal., 2004)或挤压向伸展转换阶段(Qinetal., 2005)，可能与俯冲大陆板片断离作用有关(slab break-off, De Boorderetal., 1998)。陆内环境的斑岩铜金矿或者斑岩钼矿，则可根据其形成时代判断与陆内伸展背景或者走滑转换背景相关(Wangetal., 2006; Richards, 2009)。更为重要的是，不同的构造动力学背景之间可以相互叠置衍生出复合模型，用来解释某些大型-超大型斑岩矿床及矿集区形成机制，如Gaoetal.(2018)利用“洋脊俯冲和俯冲改造的基性下地壳”模型解释中亚造山带的部分斑岩矿床形成机制。Dengetal.(2021)提出三江地区新生代大尺度岩石圈构造变形——地壳差异性旋转造成壳幔圈层非耦合变形，活化早期预富集的成矿物质，形成北衙和玉龙巨型斑岩-矽卡岩矿集区。由此可见，当斑岩矿床形成于俯冲-碰撞-后碰撞一系列构造演化过程时，触发斑岩矿床形成的动力学机制往往是较为复杂的。

鉴于特提斯构造域鲜明的地质演化特征是南侧冈瓦纳大陆不断裂解、块体向北漂移并与劳亚大陆持续聚合(吴福元等, 2020)，该特征被形象称为“特提斯号”单向列车(Wanetal., 2019)。而这种板块单向裂解-聚合过程在土耳其-伊朗地区尤为突出(engör and Yilmaz, 1981; Stampfli and Borel, 2004)。数值模拟研究(Menantetal., 2016)显示板块动力学和相关的软流圈影响着该地区弧与弧后区域岩浆活动的形成、运移和排布，而形成的熔体将弱化岩石圈并促进岩石圈的变形及其与软流圈的耦合。该三维板块动力学模型强调俯冲动力学过程与相关的软流圈流动对矿床在空间和时间上分布的影响，以及在伸展的弧后环境中变质穹窿对含矿流体循环-集中的重要性(Menantetal., 2018)。由此可见，这种俯冲体制是控制包括土耳其在内的白垩纪晚期以来地中海东部俯冲带的地球动力学和成矿演化重要动力学因素。在该种向北单向俯冲体制控制下，土耳其西部斑岩型矿床形成时代具有向南逐渐变年轻的特点。这反映其形成与分布受控于碰撞后俯冲的新特提斯洋板片(Vardar洋)后撤-回转过程(slab roll-back)，即地幔楔热结构在板片后撤-板片回转过程中发生改变，热的软流圈地幔回流触发地幔楔部分熔融，形成Tavanl带与Biga半岛斑岩成矿带的斑岩型矿床。这种俯冲的洋壳发生板片回转和断裂过程，在土耳其以西的Apuseni-Banat-Timok-Srednogorie成矿带上形成时代向南迁移岛弧性质的火山岩和花岗岩，von Quadtetal.(2005)认为这些与钙碱性岩浆有关的斑岩型Cu-Au矿和浅成低温矿床可能与碰撞后俯冲板片的回转过程有关。而在特提斯东段西藏冈底斯斑岩矿床中，仅在古新世-始新世岩浆活动有关的斑岩矿床具有相似时代分布特征(Zhaoetal., 2014, 2016)，并被认为是与俯冲的新特提斯洋壳板片后撤的构造过程相关。这一时期，陆陆碰撞造成俯冲的新特提斯洋板片发生回转，引发的软流圈回流改变了地幔楔的热结构，并诱发同时期藏南林子宗组岛弧火山岩向南迁移(Leeetal., 2009)。由此可见，这种新特提斯洋向北俯冲过程中由于陆-陆碰撞形成板片后撤-板片回转过程是特提斯中-东段形成的斑岩矿床的一类重要触发机制。

此后土耳其西部地区经历了地壳隆升过程(图7和图8)，在中新世Afyon-Konya带形成了以PnarbaMo-Cu矿和Au矿为代表的斑岩型矿床。新近研究指出(Rabayrol and Hart, 2021)，在西安纳托利亚中新世成矿作用主要集中在晚中新世(15～9Ma)，之前是成矿较差的早中新世岩浆(21～18Ma)，之后为贫矿的上新世至第四纪碱性岩浆。这种成矿集中期与中-东安纳托利亚地区岩浆岩有明显差别，即中东部早期为成矿岩浆，随后岩浆向含水少、氧化程度低和以地幔组分为主的源区特征演化。这种整体上岩浆活动特征与中新世俯冲的爱琴海与塞浦路斯板块回转-撕裂，及阿拉伯板块断裂过程有关，该过程往往造成贫水的软流圈地幔上涌和回流，进而减薄了安纳托利亚地区深部岩石圈地幔(Rabayrol and Hart, 2021)。此时控制该区域岩浆演化和斑岩矿床形成的动力学条件已由北侧Vardar洋俯冲消减造成的Tauride-Anatolide板块碰撞，转为由土耳其以南爱琴海板片俯冲控制的伸展区域。我们的研究指出，在这一中新世深部过程响应下同时期土耳其西部地区浅部地壳经历了明显减薄过程(图8a, b)，深部富集的岩石圈地幔减压发生部分熔融(Richards, 2009; Bakeretal., 2016)，形成的岩浆向着更高氧化程度演化(图5d与图8c中绿色与黄色)，这类高氧逸度的岩浆条件使得硫元素更容易保存在熔体中(如最高的熔体硫估计量)，为后期成矿提供有利条件。该动力学推论与前人利用地震数据获得土耳其西部地壳厚度(Zhuetal., 2006)与Moho面深度研究结果相一致(Karabulutetal., 2013)，即从Izmir-Ankara-Erzincan缝合带到Menders Massif再到地中海沿岸，地壳厚度从36km到28～30km减薄至25km，Moho面深度从32km到25km到20km。

因此，我们认为触发土耳其西部新生代斑岩矿床形成的动力学机制是：在新特提斯洋向北单向汇聚的背景下，北部始新世-渐新世斑岩矿床受控于碰撞后俯冲的新特提斯洋板片(Vardar洋)后撤-回转-断离过程；南部中新世斑岩矿床的形成受控于爱琴海板片俯冲控制的地壳伸展-减薄过程。

5.4 土耳其西部斑岩型矿床形成的深部过程

当关注到斑岩矿床形成的深部壳幔相互作用时，研究的着力点便集中到两个概念模型：一个是地幔楔部分熔融形成的玄武质岩浆MASH过程(熔融-同化-储存-均一过程，Hildreth and Moorbath, 1988; Richards, 2003; Lietal., 2016)，普遍用来解释俯冲环境斑岩系统中酸性岩浆的形成；另一个是深部地壳热区过程(Deep Crustal Hot Zone过程，简称DCHZ过程，Annenetal., 2006; Kempetal., 2007)，该模型强调岩浆成分取决于玄武质岩浆的侵位深度及与其平衡的地温梯度，其中富水玄武质母岩浆形成于玄武岩分离结晶和地壳围岩部分熔融过程(由冷却玄武岩的热和水造成, Annenetal., 2006)，目前有学者利用该模式解释特提斯中段碰撞环境斑岩矿床(Dengetal., 2018; Wanetal., 2018) 和我国华南晚中生代与花岗岩有关矿床(Wangetal., 2021) 形成的深部过程。但应用这两个概念模型有一定的适用条件：MASH模式需要有持续的俯冲作用(俯冲流体造成地幔楔部分熔融)并提供稳定的热和基性岩浆，以便在壳幔边界附近停留和均一化。而DCHZ模式需要考虑到斑岩矿床的形成过程中发生基性岩浆注入过程(Hattori and Keith, 2001)。 因此，MASH模式更容易解释北侧Izmir-Ankara-Erzincan缝合带附近的Tavanl与Biga半岛斑岩成矿带斑岩矿床的形成，因为白垩纪-古新世的Pontides斑岩铜矿带是由Vardar洋向北持续俯冲形成的，该时期岩浆岩的弧型岩浆特征2009; Kucuetal., 2019)表明至少在古新世该区域深部仍存在持续的俯冲作用。加之板片回转-断裂形成广泛的中酸性岩浆作用和混合作用，指示深部MASH模式可能是控制该区域斑岩矿床形成的重要深部过程。而对于Afyon-Konya带中Pnarba斑岩Mo-Cu矿(18～17Ma)和斑岩Au矿(16～12Ma)，本研究指出形成该时期其高氧逸度成矿岩浆可能直接继承于富集的岩石圈地幔部分熔融。伸展过程中地壳减薄造成富集地幔部分熔融形成玄武质岩浆，并且该时期熔体Cl含量明显高于Tavanl与Biga半岛斑岩矿床中熔体Cl含量(图10c)，表明岩浆磷灰石结晶的熔体环境存在明显差异。由于深部DCHZ带内底侵的玄武质岩浆能够同时为浅部岩浆房提供S和Cl(Zhuetal., 2018)，本文更倾向于利用DCHZ模式解释中新世Afyon-Konya带斑岩矿床形成的深部过程。但在DCHZ过程中底侵玄武岩注入岩浆房的深度和时限则需要下一步更为精细的研究确定。

6 结论

(1)锆石微量元素表明土耳其西部新生代斑岩矿床成矿岩浆大部分落于高水含量-高氧逸度区间；锆石Ti温度计估算土耳其西部斑岩成矿及其相关岩体的形成温度在650～900℃之间，结晶温度表现为斑岩Au矿>斑岩Cu-Mo=斑岩Mo-Cu矿>斑岩Mo矿，呈现递减趋势；锆石微量元素显示(Yb/Gd和Hf/Y-Th/U)其组成大部分受角闪石±榍石±磷灰石分离结晶控制；从始新世到中新世，锆石氧逸度Ce4+/Ce3+比值呈现出先降低(斑岩Cu-Mo到斑岩Mo)后升高(斑岩Mo-Cu到斑岩Au)的趋势，且斑岩Au和斑岩Cu-Mo相关岩体锆石形成时氧逸度高于斑岩Mo和斑岩Mo-Cu矿岩体。

(2)始新世-渐新世斑岩型Cu-Mo与Mo矿床中岩浆相对具有高F和低Cl组成，中新世伸展环境形成的斑岩Mo-Cu矿和斑岩Au矿中岩浆Cl含量普遍较高，其中与成矿岩体同期的暗色包体或基性岩脉中磷灰石计算获得的岩浆硫含量均大于侵入体对应岩浆的硫含量，表明基性岩浆注入可能为岩浆房提供硫。

(3)通过区域动力学分析和地壳厚度估算，本文认为触发土耳其西部新生代斑岩矿床形成的动力学机制是：在新特提斯洋向北单向汇聚的背景下，北部始新世-渐新世斑岩矿床受控于碰撞后俯冲的新特提斯洋板片(Vardar洋)后撤-回转-断离过程；南部中新世斑岩矿床的形成受控于爱琴海板片俯冲控制的地壳伸展-减薄过程。北侧Izmir-Ankara-Erzincan缝合带附近的Tavanl与Biga半岛斑岩成矿带始新世-渐新世斑岩矿床的形成与熔融-同化-储存-均一过程(MASH)有关；深部地壳热区过程(DCHZ)与中新世Afyon-Konya带斑岩矿床有关。

致谢野外工作得到中国科学院地质与地球物理研究所李继磊副研究员、褚杨副研究员、冯振天博士、刘霞博士、景揭俊博士的帮助；测试工作得到中国科学院地质与地球物理研究所余可龙博士和贾立辉工程师的协助；成文过程受益于“特提斯地球动力系统重大研究计划”研究团队的学术交流，并得到了中国科学院地质与地球物理研究所万博研究员、曹明坚研究员、胡方泱副研究员和邹心宇博士，中国科学院青藏高原研究所李金祥研究员的帮助、指导和启发；审稿人的宝贵修改意见和建议，以及期刊编辑的精心修改，使文章得以完善。在此笔者及其团队一并致以诚挚的谢意。

在笔者学生时期对斑岩系统构造背景研究过程中，深得李继亮老师的帮助和教导，谨以此文纪念李继亮老师。