岩石学报  2020, Vol. 36 Issue (10): 3137-3151, doi: 10.18654/1000-0569/2020.10.12   PDF    
新疆塔什库尔干温泉地区花岗岩体侵入与新生代构造变形:对东北帕米尔中新世构造演化的启示
蔡志慧, 何碧竹, 刘若涵     
自然资源部深地动力学重点实验室, 中国地质科学院地质研究所, 北京 100037
摘要: 帕米尔高原是受到印度-亚洲大陆碰撞、持续汇聚影响最显著的地区之一,以强烈地壳增厚和缩短、大量断裂和片麻岩穹窿的形成以及广泛的陆内岩浆活动为特征。以往有关帕米尔陆内岩浆岩的讨论多集中于对其地球化学成分及其所指示的构造背景的研究方面,而对岩浆形成与大型新生代构造之间的联系研究较少。本文通过对东北帕米尔塔什库尔干温泉地区新生代花岗岩及其围岩中锆石U-Pb和Ar-Ar年代学研究,结合该地区新生代构造变形分析,揭示岩浆侵入与区域构造变形之间的关系。U-Pb及Ar-Ar测年结果显示温泉地区花岗岩脉形成时代为中新世(11.8±0.2Ma),其及其围岩在10.8±0.1Ma冷却到300℃左右;中新世花岗岩脉中继承锆石及围岩片岩中碎屑锆石U-Pb的年龄分布特征迥异,它们分别具有类似松潘-甘孜地体以及中帕米尔地体的物质来源。花岗岩内部几乎无变形,围岩片岩变形主要体现为近E-W向的伸展构造,反映其形成于拉张为主的构造应力环境。结合区域构造背景,推测温泉中新世花岗岩脉为公格尔山、慕士塔格峰构造单元岩石部分熔融产物,是在东北帕米尔地壳从挤压增厚向局部伸展转换的过程中形成的,此时,公格尔伸展断裂系可能已经开始发育;此后,直到6~4Ma,公格尔伸展断裂系开始快速运动,与之相伴公格尔-慕士塔格片麻岩穹窿快速折返。
关键词: 东北帕米尔    中新世岩浆岩    构造变形    
Emplacement of granitic pluton and Cenozoic deformation in the Wenquan region, Tashkorgan, Xinjiang: The implications for the Miocene tectonic evolution of the Northeast Pamir
CAI ZhiHui, HE BiZhu, LIU RuoHan     
Key Laboratory of Deep-Earth Dynamics of Ministry of Natural Resource, Institute of Geology, Chinese Academy of Geological Sciences, Beijing 100037, China
Abstract: The Pamir Plateau is one of the most important regions to study effects of the India-Asia continental collision and its continuous convergence. It is characterized by significantly crust thickening and shortening,formation of lots of faults and the gneiss domes as well as volumerous intra-continental magmatism. The previous works on the intra-continental magmatic rocks in Pamir area focused on their geochemical compositions and indicative tectonic settings,but few studies on the relationship between the formation of the magmatism and the regional deformation. This paper presents new structural geology,U-Pb and Ar-Ar chronological data in the Wenquan area,Tashkorgan,Northeast Pamir,reveals the relationship between the emplacement of granitic pluton and the regional deformation. The U-Pb dating results showed that the granitic veins in the Wenquan area were formed at 11.8±0.2Ma. They and their wall rocks cooled to Ar-Ar closure temperature at 10.8±0.1Ma. The inherited zircons in Miocene granitic veins and detrital zircons in wall rocks display quite different U-Pb age characteristics. They were most likely derived from the Songpan-Ganzi and Central Pamir terrains,respectively. The granites show a non-deformed feature. The surrounding schists are manifested as the E-W striking extensional structures,reflecting an extensional stress field. Combined with regional tectonic settings,it indicates that the Miocene granitic veins in the Wenquan region are the partial melting products from the structural terrane of Konger Shan and Muztaghata. The granites were formed during the transformation stage from the crustal thickening to local crustal extension in the Northeast Pamir. It is also possible that the Konger extensional fault system was activated during this process. Until 6~4Ma,rapid slipping along the Kongershan extensional fault system and exhumation of the Konger-Muztaghata gneiss dome commenced.
Key words: Northeast Pamir    Miocene magmatic rocks    Deformation    

青藏-帕米尔造山带是世界上最高、最大的高原,是研究板块动力学、区域大地构造以及地表抬升-剥蚀过程的天然实验室。新生代印度-亚洲大陆碰撞、持续汇聚过程是形成青藏-帕米尔造山带的原始动力(Tapponnier et al., 1981, 2001; Burtman and Molnar, 1993; Clark and Royden, 2000; Yin and Harrison, 2000; DeCelles et al., 2002; Royden et al., 2008)。随着印度陆块向北推进,喜马拉雅及青藏高原腹地不断隆升,大量物质侧向挤出,同时形成青藏高原东、西两侧的造山带(Tapponnier et al., 2001)。帕米尔造山带位于青藏高原的西侧,呈向北凸出的弧形,故又称帕米尔弧形构造带(图 1)。帕米尔弧形构造带是印度/亚洲碰撞效应最强烈的地区之一,其在新生代之前地壳厚度仅为25~30km(Burtman and Molnar, 1993),而深部地球物理资料显示现今帕米尔地壳厚度可达70km(Mechie et al., 2012)。中新世火山岩中的捕掳体表明印度-亚洲大陆碰撞初期帕米尔地区岩浆-沉积岩石组合曾埋深到>50~80km(Ducea et al., 2003)。此外,帕米尔弧形构造带目前的宽度明显小于青藏高原主体,有的学者认为帕米尔地壳增厚的同时,通过逆冲推覆等方式发生了强烈缩短变形(Burtman and Molnar, 1993)。

图 1 帕米尔东部构造简图(据Robinson et al., 2012修改) GMD-公格尔-慕士塔格片麻岩穹窿;KKP-卡日巴生-苦子干岩体;KKF-喀喇昆仑断裂;KES-公格尔伸展断裂系;TAF-塔阿西断裂;NK-公格尔山北地区. ①温泉构造单元; ②慕士塔格峰西构造单元; ③慕士塔格南构造单元; ④塔阿西南部构造单元; ⑤公格尔山及慕士塔格峰构造单元; ⑥马尔洋构造单元. B-B′为东北帕米尔NE向剖面,见图 8 Fig. 1 Sketch tectonic maps of the eastern Pamir region (after Robinson et al., 2012)

帕米尔弧形构造带受印度-亚洲大陆碰撞及持续汇聚的影响发生强烈变形不仅体现在地壳的增厚和缩短,还体现在形成许多大型新生代断裂(Brunel et al., 1994; Murphy et al., 2000; Lacassin et al., 2004; Schwab et al., 2004; Valli et al., 2008; Robinson, 2009; Amidon and Hynek, 2010; Cowgill, 2010; Klemperer et al., 2013)和片麻岩穹窿(Robinson et al., 2004, 2007, 2012; Schmidt et al., 2011; Stearns et al., 2013, 2015; Stübner et al., 2013a, b ; Rutte et al., 2017a, b )。这些新生代构造一方面是碰撞造山的结果,另一方面它们也改造着帕米尔地区的构造格局,切断连续地体,为侧向、深部-浅部物质运移提供途径等。例如,在帕米尔弧形构造带东部与青藏高原主体相连接的位置,各地体的宽度明显变窄,并且被许多新生代断裂切断(图 1),使得帕米尔弧形构造带与青藏高原主体中的地体对应关系变得难以确定;帕米尔弧形构造带的中、高级变质岩折返到地表过程主要是通过大型断裂及穹窿构造实现的(Robinson et al., 2004, 2007; Schmidt et al., 2011; Stearns et al., 2013, 2015; Stübner et al., 2013a, b ; Rutte et al., 2017b; Chapman et al., 2018)。另外,帕米尔弧形构造带受印度-亚洲大陆碰撞和持续汇聚影响还体现在区内普遍存在中新世陆内岩浆活动(Jiang et al., 2002, 2012; Ducea et al., 2003;罗照华等, 2003; Hacker et al., 2005, 2017;柯珊等, 2006; Robinson et al., 2007; Chapman et al., 2018)。许多学者围绕这些岩浆岩的地球化学成分及其所指示的岩浆成因和构造背景等展开讨论,指出这些岩浆岩多为中新世时期的钾质-超钾质岩石,其形成与地壳减压熔融、软流圈上涌有关(Jiang et al., 2002, 2012; Hacker et al., 2005, 2017)。然而,有关帕米尔中新世岩浆岩的形成与前面所述该区大型新生代构造之间的联系则很少有人研究。例如,中新世岩浆的形成、侵位与大型断裂或片麻岩穹窿等的形成演化之间是否存在关系?如果有,其相互影响过程如何?造山带是构造变形、岩浆以及变质过程的综合产物,结合区域构造变形对岩浆侵入过程的研究是全面认识造山作用的一个重要内容。选择帕米尔地区与中新世岩浆岩相邻的新生代大型断裂及片麻岩穹窿等关键位置的构造变形研究,进行构造-岩浆活动分析,可以为探讨帕米尔造山过程提供更多信息。

在对帕米尔东北部进行野外考察的基础上,我们选取新疆塔什库尔干温泉地区新生代岩浆岩及其围岩进行构造变形、U-Pb及Ar-Ar年代学测试分析,其目的是阐述岩浆侵入过程与公格尔伸展断裂系及慕士塔格片麻岩穹窿构造形成之间的关系,进一步探讨其对东北帕米尔新生代构造演化的启示,并为印度-亚洲大陆碰撞后东北帕米尔新生代构造演化提供科学依据。

1 地质背景

根据若干大型断裂及古老蛇绿岩残片的分布,可将帕米尔弧形构造带划为北帕米尔、中帕米尔、南帕米尔及兴都库什等几个地体(Burtman and Molnar, 1993; Robinson et al., 2004, 2007; Schmidt et al., 2011; Stübner et al., 2013a)(图 1)。

新疆塔什库尔干地区位于帕米尔弧形构造带的东北部,区内大型新生代构造包括喀喇昆仑断裂、公格尔伸展断裂以及公格尔-慕士塔格片麻岩穹窿等。喀喇昆仑断裂为区内规模最大的断裂,前人围绕其规模、走滑位移量、形成时间、以及对印度亚洲大陆碰撞后物质侧向挤出贡献等方面进行了广泛讨论(Brunel et al., 1994; Murphy et al., 2000; Lacassin et al., 2004; Schwab et al., 2004; Valli et al., 2008; Robinson, 2009; Amidon and Hynek, 2010; Cowgill, 2010; Klemperer et al., 2013)。除喀喇昆仑断裂之外,帕米尔弧形构造带东部的新生代公格尔伸展断裂系(Brunel et al., 1994; Robinson et al., 2004, 2007, 2010, 2012; Cao et al., 2013)规模也较大,长约250km,并被认为与北帕米尔中部的卡拉库尔地堑(Robinson et al., 2004; Amidon and Hynek, 2010)具有近一致的初始活动时间,二者的形成皆晚于早中新世,标志着东北帕米尔与北帕米尔从地壳增厚到局部伸展转换的时间(Robinson et al., 2004; Amidon and Hynek, 2010)。Robinson et al.(2004)提出公格尔山正滑断裂E-W向大规模伸展开始于中新世,并认为其演化模式为由公格尔山向西扩展或由木吉向东扩展。除了上述新生代大型断裂之外,东北帕米尔塔什库尔干地区另外一种显著的新生代构造为片麻岩穹窿。区内的片麻岩穹窿主要为公格尔-慕士塔格片麻岩穹窿,其形成与公格尔山伸展断裂系的演化密切相关(Robinson et al., 2004, 2007, 2010, 2012)。Cao et al.(2013)也强调这一点,并指出慕士塔格-公格尔片麻岩穹窿快速折返时间集中在6~4Ma和3~1Ma。Thiede et al.(2013)从低温热年代学研究的角度出发,指出东北帕米尔片麻岩穹窿形成于同造山、E-W向伸展作用过程中。Cai et al.(2017)对慕士塔格麻岩穹窿西侧和东侧两条剪切带进行了构造变形及年代学分析,指出二者分别具有向W-SW和向E剪切指向,在12~8Ma时位于中上地壳。

除了上述大型断裂及片麻岩穹窿之外,东北帕米尔塔什库尔干地区还包括大量新生代岩浆岩(Jiang et al., 2002, 2012;罗照华等, 2003;柯珊等, 2006),如卡日巴生-苦子干岩体及昝坎岩体(Jiang et al., 2002)。卡日巴生-苦子干岩体属于Sares-Murgab岩体的一部分,位于塔什库尔干河谷西侧,侵入在二叠系及前寒武纪布伦阔勒群中(河南省地质调查院, 2004)。锆石U-Pb同位素年代学研究表明卡日巴生-苦子干岩体形成时代为10~11.9Ma(柯珊等, 2006; Robinson et al., 2007; Jiang et al., 2012;徐晓尹等, 2017),岩浆具有中新世加厚下地壳(Jiang et al., 2002)或地幔(Jiang et al., 2012)来源,或成分不同的岩浆具有地幔和地壳不同来源(罗照华等, 2003)。

① 河南省地质调查院. 2004.中华人民共和国1/25万塔什库尔干塔吉克自治县幅(J43C003003)、克克吐鲁克幅(J43C003002)、叶城县幅(J43C003004)区域地质调查报告

2 塔什库尔干温泉地区新生代岩体及其围岩的构造变形特征

塔什库尔干温泉新生代岩体位于温泉剖面西段(图 2图 3a),其为卡日巴生-苦子干岩体东侧一部分(图 1)。温泉剖面的岩体主体为花岗岩,主要成分包括长石、石英及云母等。另外在岩体主体的围岩布伦阔勒群片岩中还侵入许多花岗岩脉,其成分与岩体主体花岗岩成分相似。这些较小的花岗岩脉规模大小不等,若干厘米至若干米宽;形态各异,呈条带状(图 3b),透镜状(图 3ce)或“布丁”状(图 3f)。内部无明显变形,矿物颗粒自形且粒度较大,靠近围岩布伦阔勒群片岩位置颗粒变细,具有冷凝边,反映了花岗岩脉与片岩之间的侵入接触关系。

图 2 东北帕米尔塔什库尔干地区地质简图(a)、温泉-库克西克剖面图(b, 剖面位置见图a中A-A′)及(c)E-W向温泉剖面图 WMSZ-慕士塔格西剪切带; KS-库克西克剪切带; SDF-辛迪断裂.同位素年龄来源:(0)本文; (1)Cai et al., 2017;(2)Robinson et al., 2007;(3)徐晓尹等, 2017;(4)柯珊等, 2006;(5)Jiang et al., 2012;(6)Robinson et al., 2004 Fig. 2 Geological map of Tashkorgan area in Northeast Pamir (a), Wenquan-Kukexilik cross-section (b, see A-A′ in the Fig. 1a) and E-W trending Wenquan corss-section (c)

图 3 温泉地区野外地质现象 Fig. 3 Field photos of outcrops in the Wenquan region

温泉地区布伦阔勒群主要为云母长石石英片岩,变质程度较低,主要体现在形成绢云母和绿泥石等低温变质矿物,局部可见原生沉积层理。布伦阔勒群岩石变形主要分为二期:第一期表现为沿片岩层理方向面理化D1(图 3b-h),面理产状以中等-高角度向NE倾为主(图 3b-h)。该期变形作用早于新生代花岗岩侵入,但根据本文数据无法限定其具体时代。第二期为与花岗岩侵入有关的剪切变形D2。在新生代岩体与布伦阔勒群片岩接触位置岩石发生较明显变形。由长石及石英等矿物定向组成的线理产状向E、SE方向倾伏(图 3d)。片岩中可见细小“σ”形石英脉及不对称褶皱形石英脉等,显示其经历正滑韧性剪切过程,剪切指向为E、SE方向(图 3f-h)。花岗岩脉产状明显与片岩的面理(软弱面)相关,多沿着片岩面理侵入、局部位置小角度斜切面理(图 3c, e, f)。而稍大些的花岗岩脉常常表现为部分斜切片岩早期面理(图 3b),反映了片岩面理形成早于花岗岩年龄。有的透镜状花岗岩脉呈“σ”形,显示逆冲运动指向(图 3e)。与其紧邻的片岩局部的不对称褶皱也显示逆冲的运动特征(图 3d),这反映了花岗岩脉由下向上逆冲侵入的过程。与其非常紧邻的片岩中的不对称褶皱由于岩体向上拖曳作用而呈现逆冲特征。有的花岗岩脉侵入过程受到片岩中断裂或节理的影响而产生位移等明显形状改变(如图 3f中的褶皱状花岗岩脉的形态受到小断裂的影响)。“布丁”状花岗岩脉则显示被“拉断”的形态,且“布丁”中各段岩体非矩形而是呈透镜状(图 3f),反映了花岗岩脉的侵入与水平方向的拉张是同时发生的。

显微镜下,温泉花岗岩脉中长石、石英及云母等矿物自形,不显示定向排列特征(图 4a)。布伦阔勒群片岩中矿物定向排列(图 4b-d),石英变形特征为有的边界不平滑(图 4b, c)、具有带状消光(图 4c);有的石英边部显示细小颗粒化、消光位置明显发生改变(图 4b, c),体现了其以亚晶粒旋转和膨凸为主的塑性变形特征。长石主要以碎裂变形为主,斜长石双晶纹较平直(图 4b, c)。石英带状消光位置与岩石面理之间锐角夹角指示下滑的剪切指向(图 4c),为D2变形产物。面理被花岗岩脉斜切或横切的片岩样品显微镜下石英、长石及云母等定向排列,石英明显拉长,剪切指向不明显(图 4d),此期变形早于花岗岩脉的侵入,为D1。有的片岩样品中局部显示石英近等粒状、边部较平直且常与相邻颗粒呈120°夹角,黑云母颗粒碎小,边界不平整(图 4b, d),说明流体的存在,黑云母被溶蚀,岩石经历静态重结晶过程(Evans et al., 2001)。

图 4 温泉样品薄片显微照片 Bt-黑云母; Fsp-长石; Pl-斜长石; Qz-石英 Fig. 4 Microscopic photographs of petrologic thin sections from the Wenquan region

总而言之,塔什库尔干温泉地区新生代构造变形在花岗岩体内部不明显,而主要表现在围岩布伦阔勒群片岩中,其变形主要分为二期:第一期体现为沿片岩层理方向的面理化D1;第二期为韧性下滑剪切D2。D2整体叠加于D1之上,虽然由于岩体侵位造成应变分布不均使得片岩局部显示逆冲特征,但总体反映了此时温泉地区处于拉张的应力环境。

3 塔什库尔干温泉地区岩石同位素年代学特征 3.1 样品描述

采集塔什库尔干温泉地区布伦阔勒群云母长石石英片岩XK4-1-2与XK4-1-8,以及侵入到布伦阔勒群中的花岗岩脉(约50cm宽)样品(XK4-1-3),采样点见图 2c。片岩XK4-1-2与XK4-1-8样品为浅灰色片状变晶结构、块状构造,中等粒度,主要矿物为石英(50%~60%)、长石(25%~35%)、黑云母(15%~18%)及白云母(5%~8%)等。花岗岩脉XK4-1-3呈中粗粒花岗结构、块状构造(图 4a)。主要由长石(40%~50%)、石英(30%~40%)以及黑云母(5%~10%)组成,副矿物主要为铁的氧化物、磷灰石、榍石及锆石等,其与卡日巴生岩体的矿物组成(柯珊等, 2006)类似,属于花岗岩系列。

3.2 测试方法

将样品粉碎后用重液和磁力分选方法初步选出锆石, 然后在显微镜下挑选出较完整的、典型的颗粒,将其粘置于环氧树脂中制靶,打磨抛光,使锆石内部充分暴露,接着进行阴极发光照相。锆石的阴极发光照相在中国地质科学院地质研究所完成。LA-MC-ICP-MS锆石U-Pb同位素年代学分析在中国地质科学院矿产资源研究所MC-ICP-MS实验室完成,分析所用仪器为Finnigan Neptune型MC-ICP-MS及与之配套的Newwave UP 213激光剥蚀系统。锆石U-Pb测年激光剥蚀所用斑束直径为25μm,频率为10Hz,能量密度约为2.5J/cm2,以He为载气。信号较小的207Pb、206Pb、204Pb、202Hg用离子计数器(multi-ion-counters)接收,208Pb、232Th和238U信号用法拉第杯接收,实现了所有目标同位素信号的同时接收并且不同质量数的峰基本上都是平坦的,进而获得高精度的数据。详细实验测试过程可参见侯可军等(2009)。采用ICPMSDataCal程序(Liu et al., 2010)进行数据处理。采用ISOPLOT程序(Ludwig, 2003)进行谐和图绘制及加权平均年龄计算等。

用来进行Ar-Ar测年的黑云母分选方法与锆石分选方法相似:利用重液和磁选方法初步选出黑云母,接着在显微镜下进一步挑选,然后仔细清洗选出的黑云母。将清洗后的样品封入石英瓶中送至中国原子能科学研究院进行中子照射。照射总时间为1440分钟,同时接受中子照射的还有用来监控的标样。Ar-Ar同位素测年分析在中国地质科学院地质研究所同位素热年代学实验室完成。使用石墨炉对样品进行阶段升温加热,每阶段加热10分钟,净化30分钟。采用HelixMC多接收稀有气体质谱仪进行质谱分析,每组峰值均采集20个数据。所有数据在回归到时间零点值后再进同位素校正。详细实验流程参考文献(张彦等,2006)。用ISOPLOT程序(Ludwig, 2003)进行坪年龄计算。

3.3 测试结果

花岗岩脉XK4-1-3中的锆石颗粒中等-较大粒度,长轴约50~150μm,浅褐色-浅黄色,晶体自形程度较高,多为柱状-长柱状。锆石阴极发光图像显示大多数锆石具有明显的岩浆结晶振荡环带(图 5a),说明其岩浆成因;有的锆石具有较大、形状不规则的继承锆石核(图 5a)。13个具有振荡环得锆石测点的Th/U皆大于0.1(表 1),其206Pb/238U表面年龄的加权平均值为11.8±0.2Ma(MSWD=1.19)(图 5b)。31颗继承锆石的年龄范围为208~2420Ma(图 5c表 1)。

图 5 温泉地区样品中锆石的阴极发光图像(a)和LA-ICP-MS U-Pb谐和图(b-e) Fig. 5 Cathodoluminescence images (a) and LA-ICP-MS U-Pb concordia diagrams (b-e) of representative zircons from Wenquan region

表 1 东北帕米尔温泉地区样品LA-ICP-MS锆石U-Pb测年结果 Table 1 LA-ICP-MS zircon U-Pb dating results for samples from Wenquan region, Northeast Pamir

温泉地区2个布伦阔勒群片岩样品(XK4-1-2和XK4-1-8)中的锆石呈浅褐色-浅黄色-无色,长轴约80~150μm,晶体具有一定磨圆度,少数可见柱状、短柱状自形(图 5a),表明锆石经历了一定距离搬运,为碎屑锆石。阴极发光图显示多数锆石具有岩浆生长振荡环带,少部分具弱环带或颜色较均匀(图 5a)。2个片岩样品的锆石U-Pb测试显示,绝大多数锆石Th/U>0.1(表 1)。样品XK4-1-2中22颗碎屑锆石年龄分布范围为544~3587Ma (图 5d)。样品XK4-1-8中20颗碎屑锆石年龄分布范围为515~3098Ma (图 5e)。

片岩样品XK4-1-8中的黑云母Ar-Ar同位素测年分析数据列于表 2。840~1400℃连续9个温度阶段获得的加权平均坪年龄为10.8±0.1Ma(MSWD=1.02)(图 6表 2),10.8 Ma代表了冷却到黑云母Ar-Ar同位素体系封闭温度(约300℃)的时代。

表 2 温泉地区片岩样品白云母40Ar-39Ar阶段升温加热测年结果 Table 2 The results of 40Ar-39Ar stepwise heating dating of muscovite for the Wenquan schsit sample

图 6 温泉样品X4-1-8黑云母40Ar-39Ar阶段升温年龄谱 Fig. 6 40Ar-39Ar age spectrum of biotite in the Sample XK4-1-8 from Wenquan region
4 讨论 4.1 塔什库尔干温泉地区花岗岩脉中继承锆石及围岩片岩中碎屑锆石的源区

本文2个温泉布伦阔勒片岩样品XK4-1-2和XK4-1-8中最小碎屑锆石年龄分别为544Ma和512Ma。综合这2个样品测试结果可见:温泉地区布伦阔勒群片岩沉积时代晚于512Ma,具有500~550Ma、750~800Ma及~1000Ma的碎屑锆石年龄峰,另外也有较多年龄为2500~2600Ma的锆石(图 7a)。而花岗岩脉中继承锆石年龄分布特征与其围岩布伦阔勒群片岩碎屑锆石年龄分布特征明显不同,主要表现在:具有200~300Ma、450~550Ma及750Ma的年龄峰(图 7b)。与塔什库尔干其它地区岩石碎屑锆石年龄谱对比可见:温泉地区布伦阔勒群的锆石年龄分布特征(图 7a)与慕士塔格峰西侧(图 1-②)布伦阔勒群副片麻岩(Robinson et al., 2012)以及塔什库尔干河谷东塔阿西断裂南部地区(图 1-④)副片麻岩中锆石年龄分布特征(Zhang et al., 2018)(图 7c)相似,不同点是温泉及塔阿西断裂南部碎屑岩的变质程度较低,并未记录~200 Ma的变质年龄(Yang et al., 2010; Robinson et al., 2012; Zhang et al., 2018;本文)(图 7c)。而温泉地区花岗岩脉中继承锆石年龄分布特征(图 7b)与公格尔山及慕士塔格峰地区(图 1-⑤)及马尔洋地区(图 1-⑥)副片麻岩中的锆石年龄分布特征(Robinson et al., 2012)(图 7c)相似。即温泉地区花岗岩脉并非本地温泉地区沉积岩熔融产物,而是非常可能来源于公格尔山和慕士塔格峰地区或马尔洋地区的沉积岩。

图 7 锆石U-Pb直方统计图 Fig. 7 Zircon U-Pb cumulative age probability plots

结合东北帕米尔地区碎屑锆石U-Pb年龄及Sr-Nd同位素数据与青藏高原其他地区碎屑锆石年龄数据,Robinson et al.(2012)将东北帕米尔塔什库尔干地区进行了构造单元划分(图 1图 7):慕士塔格峰西构造单元②与慕士塔格南构造单元③同属于中帕米尔地体;公格尔山和慕士塔格峰构造单元⑤与马尔洋构造单元⑥属于同一地体,具有与麻扎-松潘-甘孜地体相同的前中生代基底;塔阿西断裂南部构造单元④与羌塘地体具有相同的前中生代基底(Robinson et al., 2012)。据此划分温泉构造单元①属于中帕米尔地体(图 1),本文温泉地区布伦阔勒群片岩碎屑锆石年代学特征也支持此观点(图 7a, d)。另外碎屑锆石年代学研究还表明中帕米尔与羌塘地体的碎屑锆石年龄分布相近(Robinson et al., 2012),即二者在其沉积时曾属于同一地体。温泉地区布伦阔勒群片岩与羌塘地体的碎屑锆石年代分布图类似也说明这一点(图 7a, d)。温泉构造单元(属于中帕米尔)与公格尔山及慕士塔格峰以及马尔洋构造单元(属于麻扎-松潘-甘孜地体)之间存在一条中生代逆冲断裂,公格尔山及慕士塔格峰以及马尔洋构造单元位于断裂下盘(Robinson et al., 2012)。东北帕米尔塔什库尔干曾发生印支期古特提斯洋闭合事件(Yang et al., 2010; Robinson et al., 2012; Zhang et al., 2018),二者在此时拼合为一(Robinson et al., 2012)。

4.2 塔什库尔干温泉地区新生代花岗岩侵入与公格尔伸展断裂系的关系

由上可见公格尔伸展断裂系(KES)的E、W两侧的构造单元属于不同地体。然而这些地体主要区别是具有不同的前中生代基底,它们之间的界限为中生代断裂。故此我们推测新生代公格尔伸展断裂系是几乎沿着已存在的印支期断裂基础上开始活动的。又由于本文温泉花岗岩脉中的继承锆石来自断裂东侧,说明公格尔伸展断裂系初始活动时间应不早于花岗岩脉形成时代(11.8Ma)。

卡日巴生-苦子干岩体北部黑云母Ar-Ar坪年龄为11.2Ma,本文温泉地区片岩黑云母Ar-Ar坪年龄为10.8±0.1Ma,说明卡日巴生-苦子干岩体侵位到围岩片岩中、围岩发生下滑剪切之后快速冷却到黑云母封闭温度的时代为10.8~11.2Ma,这之后温泉地区岩石没有发生明显的构造热事件。其他文献中热年代学数据也显示公格尔伸展断裂系主断裂上盘的云母冷却年龄皆大于8.5Ma(Robinson et al., 2007; Cai et al., 2017),而断裂下盘云母冷却年龄明显较低,如慕士塔格峰及公格尔山附近的云母Ar-Ar年龄多小于8Ma,最小年龄为~2Ma(Robinson et al., 2007),即公格尔伸展断裂系E、W两侧的岩石记录的冷却年龄明显不同,可推断公格尔伸展断裂系具有一定规模和位移量。Robinson et al.(2007)根据野外观测及低温热年代学数据提出公格尔断裂系北段的位移量大于南段位移量,北段E-W向水平伸展位移量可达~20km,而南段 < 3km,并认为这种现象可能是由于帕米尔地区弧形弯曲或向外侧逆冲过程导致的(Robinson et al., 2007),而不支持公格尔伸展断裂系为右行喀喇昆仑断裂的北端弯曲伸展转换部分的观点(Ratschbacher et al., 1994; Murphy et al., 2000)。本文数据虽然无法协助判断喀喇昆仑断裂的演化是否是公格尔伸展断裂系形成的主要诱因,但可初步表明在公格尔伸展断裂系形成之前或同时发生地壳部分熔融、形成花岗岩浆,而此部分熔融过程在单纯的走滑剪切环境下难以发生。

4.3 塔什库尔干温泉地区中新世岩浆岩侵入模式

卡日巴生-苦子干岩体整体呈NEN-SWS走向(图 1图 2)、与公格尔断裂系相平行。说明卡日巴生-苦子干岩体的侵入很可能与公格尔断裂系的初始活动位置类似,皆沿着NEN-SWS走向的中生代断裂面发生(Robinson et al., 2012)。虽然有的花岗岩脉侵入引起附近一小部分片岩显示逆冲变形,但岩体“布丁”构造及大部分围岩变形特征表明温泉花岗岩脉是由于该地区拉张伸展过程中形成的。即岩浆侵入可引起局部围岩被动褶皱等变形,在围岩中同时存在伸展和挤压构造,但以伸展构造为主。这种与断裂等构造密切相关的模式称为“被动同构造侵入”(passive syntectonic emplacement),经常在拉张应力环境下发生(Hutton, 1988; Eliasson et al., 2003)。宏观及显微现象表明卡日巴生-苦子干岩体及花岗岩脉内部几乎无变形,与围岩透入性面理构造(D1)没有直接联系,围岩布伦阔勒群片岩面理构造D1表现为岩石压实过程发生矿物重结晶,可能形成于东北帕米尔某个地壳加厚或布伦阔勒群深埋的阶段。花岗岩体内部靠近岩体与围岩边界的位置没有明显的面理,说明岩体的侵入产生微量应变,几乎没有后期构造事件及区域应变对此岩体-围岩边界位置产生作用。花岗岩脉围岩布伦阔勒群片岩黑云母Ar-Ar年龄没有明显年轻于花岗岩脉U-Pb年龄也说明没有后期热事件的改造。

4.4 东北帕米尔中新世构造演化的启示

新生代印度-亚洲陆块碰撞以及持续汇聚使得青藏高原及帕米尔高原地壳明显增厚(Rowley, 1996; Yin and Harrison, 2000; Bell and Sapkota, 2012)。地壳不断加厚并且程度不均,导致重力不均衡,从而引起地壳局部伸展、减薄(Stübner et al., 2013a, b ; Rutte et al., 2017a, b ),同时通过部分熔融作用,形成大量花岗岩浆(如卡日巴生-苦子干岩体及本文温泉片岩中花岗岩脉)。地壳局部伸展减薄的重要渠道是伸展断裂系的活动。伸展断裂可以表现在深层次的韧性剪切,通过应变量的增大导致地壳水平减薄。根据目前资料,我们初步认为在12~11Ma时,公格尔伸展断裂系已经开始发生韧性变形(如本文温泉地区的韧性下滑剪切记录)。此时,在公格尔伸展断裂以及东侧库克西力克断裂的作用下慕士塔格穹窿物质也开始缓慢折返(Sobel et al., 2011; Cai et al., 2017)(图 8a)。6~4Ma时,伴随印度亚洲陆块继续汇聚、塔里木地块的南向下插,公格尔伸展断裂开始快速活动,并在近地表产生明显位移量,慕士塔格片麻岩穹窿也在此时快速折返(图 8b)。塔什库尔干及邻区低温热年代学数据存在明显6~4Ma峰值也可以说明此过程的存在(Cao et al., 2013)。

图 8 东北帕米尔中新世构造演化模式图 剖面位置见图 1中B-B′ Fig. 8 Model for the Miocene tectonic evolution of the Northeast Pamir
5 结论

本研究展示了东北帕米尔塔什库尔干温泉地区岩浆岩及围岩接触带的构造变形现象、U-Pb和Ar-Ar年代学测试结果,得到的主要结论如下:

(1) 东北帕米尔塔什库尔干温泉地区侵入到布伦阔勒群片岩中的花岗岩脉锆石U-Pb年龄为11.8±0.2Ma,与其相邻的卡日巴生-苦子干岩体年龄近一致。

(2) 花岗岩侵入体内部几乎无变形,围岩片岩主要体现为近E-W向的伸展构造,局部少量挤压构造,反映了拉张为主的应力环境。

(3) 温泉中新世花岗岩脉中继承锆石与邻区公格尔山、慕士塔格峰构造单元岩石碎屑锆石的U-Pb年龄分布特征类似,说明其具有类似松潘-甘孜地体的物质来源;围岩布伦阔勒片岩中的碎屑锆石与慕士塔格西地区碎屑锆石的U-Pb年龄分布特征相似,说明其属于中帕米尔地体。

(4) 温泉布伦阔勒片岩黑云母Ar-Ar测年结果为10.8±0.1Ma,说明温泉地区岩石在此时已冷却到大约300℃,并且后期没有发生过较强烈的构造热事件。

(5) 结合东北帕米尔区域构造背景,我们认为温泉地区岩浆岩形成在东北帕米尔地壳挤压增厚向局部伸展转换的过程中。公格尔伸展断裂系可能在此过程开始部分启动。直到6~4Ma,公格尔伸展断裂系开始快速运动、公格尔-慕士塔格片麻岩穹窿快速折返。

致谢      感谢张建新研究员与李广伟教授的意见与建议; 感谢马绪宣博士的建议和帮助。

参考文献
Amidon WH and Hynek SA. 2010. Exhumational history of the north central Pamir. Tectonics, 29(5): TC5017
Bell TH and Sapkota J. 2012. Episodic gravitational collapse and migration of the mountain chain during orogenic roll-on in the Himalayas. Journal of Metamorphic Geology, 30(7): 651-666 DOI:10.1111/j.1525-1314.2012.00992.x
Brunel M, Arnaud N, Tapponnier P, Pan Y and Wang Y. 1994. Kongur Shan normal fault: Type example of mountain building assisted by extension (Karakoram fault, eastern Pamir). Geology, 22(8): 707-710 DOI:10.1130/0091-7613(1994)022<0707:KSNFTE>2.3.CO;2
Burtman VS and Molnar P. 1993. Geological and geophysical evidence for deep subduction of continental crust beneath the Pamir. Special Paper of the Geological Society of America, 281(2): 248-251
Cai ZH, Xu ZQ, Cao H, Robinson AC, Li GW and Xu XY. 2017. Miocene exhumation of Northeast Pamir: Deformation and geo/thermochronological evidence from western Muztaghata shear zone and Kuke ductile shear zone. Journal of Structural Geology, 102: 130-146 DOI:10.1016/j.jsg.2017.07.010
Cao K, Bernet M, Wang GC, van der Beek P, Wang A, Zhang KX and Enkelmann E. 2013. Focused Pliocene-Quaternary exhumation of the Eastern Pamir domes, western China. Earth and Planetary Science Letters, 363: 16-26 DOI:10.1016/j.epsl.2012.12.023
Chapman JB, Scoggin SH, Kapp P, Carrapa B, Ducea MN, Worthington J, Oimahmadov I and Gadoev M. 2018. Mesozoic to Cenozoic magmatic history of the Pamir. Earth and Planetary Science Letters, 482: 181-192 DOI:10.1016/j.epsl.2017.10.041
Clark MK and Royden LH. 2000. Topographic ooze: Building the eastern margin of Tibet by lower crustal flow. Geology, 28(8): 703-706 DOI:10.1130/0091-7613(2000)28<703:TOBTEM>2.0.CO;2
Cowgill E. 2010. Cenozoic right-slip faulting along the eastern margin of the Pamir salient, northwestern China. GSA Bulletin, 122(1-2): 145-161 DOI:10.1130/B26520.1
DeCelles PG, Robinson DM and Zandt G. 2002. Implications of shortening in the Himalayan fold-thrust belt for uplift of the Tibetan Plateau. Tectonics, 21(6): 12-1
Ducea MN, Lutkov V, Minaev VT, Hacker B, Ratschbacher L, Luffi P, Schwab M, Gehrels GE, McWilliams M, Vervoort J and Metcalf J. 2003. Building the Pamirs: The view from the underside. Geology, 31(10): 849-852 DOI:10.1130/G19707.1
Eliasson T, Ahlin S and Petersson J. 2003. Emplacement mechanism and thermobarometry of the Sveconorwegian Bohus granite, SW Sweden. GFF, 125(3): 113-130 DOI:10.1080/11035890301253113
Evans B, Renner J and Hirth G. 2001. A few remarks on the kinetics of static grain growth in rocks. International Journal of Earth Sciences, 90(1): 88-103
Gehrels G, Kapp P, DeCelles P, Pullen A, Blakey R, Weislogel A, Ding L, Guynn J, Martin A, McQuarrie N and Yin A. 2011. Detrital zircon geochronology of pre-Tertiary strata in the Tibetan-Himalayan orogen. Tectonics, 30(5): TC5016 DOI:10.1029/2011TC002868
Hacker B, Luffi P, Lutkov V, Minaev V, Ratschbacher L, Plank T, Ducea M, Patiño-Douce A, McWilliams M and Metcalf J. 2005. Near-ultrahigh pressure processing of continental crust: Miocene crustal xenoliths from the Pamir. Journal of Petrology, 46(8): 1661-1687 DOI:10.1093/petrology/egi030
Hacker BR, Ratschbacher L, Rutte D, Stearns MA, Malz N, Stübner K, Kylander-Clark ARC, Pfänder JA and Everson A. 2017. Building the Pamir-Tibet Plateau-Crustal stacking, extensional collapse, and lateral extrusion in the Pamir: 3. Thermobarometry and petrochronology of deep Asian crust. Tectonics, 36: 1743-1766
Hou KJ, Li YH and Tian YR. 2009. In situ U-Pb zircon dating using laser ablation-multi ion counting-ICP-MS. Mineral Deposits, 28(4): 481-492 (in Chinese with English abstract)
Hutton DHW. 1988. Granite emplacement mechanisms and tectonic controls: Inferences from deformation studies. Earth and Environmental Science Transactions of the Royal Society of Edinburgh, 79(2-3): 245-255 DOI:10.1017/S0263593300014255
Jiang YH, Jiang SY, Ling HF, Zhou XR, Rui XJ and Yang WZ. 2002. Petrology and geochemistry of shoshonitic plutons from the western Kunlun orogenic belt, Xinjiang, northwestern China: Implications for granitoid geneses. Lithos, 63(3-4): 165-187
Jiang YH, Liu Z, Jia RY, Liao SY, Zhou Q and Zhao P. 2012. Miocene potassic granite-syenite association in western Tibetan Plateau: Implications for shoshonitic and high Ba-Sr granite genesis. Lithos, 134-135: 146-162 DOI:10.1016/j.lithos.2011.12.012
Ke S, Mo XX, Luo ZH, Zhan HM, Liang T, Li L and Li WT. 2006. Petrogenesis and geochemistry of Cenozoic Taxkorgan alkalic complex and its geological significance. Acta Petrologica Sinica, 22(4): 905-915 (in Chinese with English abstract)
Klemperer SL, Kennedy BM, Sastry SR, Makovsky Y, Harinarayana T and Leech ML. 2013. Mantle fluids in the Karakoram fault: Helium isotope evidence. Earth and Planetary Science Letters, 366: 59-70 DOI:10.1016/j.epsl.2013.01.013
Lacassin R, Valli F, Arnaud N, Leloup PH, Paquette JL, Li HB, Tapponnier P, Chevalier ML, Guillot S, Maheo G and Xu ZQ. 2004. Large-scale geometry, offset and kinematic evolution of the Karakorum fault, Tibet. Earth and Planetary Science Letters, 219(3-4): 255-269 DOI:10.1016/S0012-821X(04)00006-8
Liu YS, Gao S, Hu ZC, Gao CG, Zong KQ and Wang DB. 2010. Continental and oceanic crust recycling-induced melt-peridotite interactions in the Trans-North China Orogen: U-Pb dating, Hf isotopes and trace elements in zircons from mantle xenoliths. Journal of Petrology, 51(1-2): 537-571 DOI:10.1093/petrology/egp082
Ludwig KR. 2003. Isoplot/Ex Version 3.0: A Geochronological Toolkit for Microsoft Excel. Berkeley, CA: Berkeley Geochronology Center
Luo ZH, Mo XX and Ke S. 2003. Ages of Taxkorgan alkaline intrusive complex and their geological implications. Xinjiang Geology, 21(1): 46-50 (in Chinese with English abstract)
Mechie J, Yuan X, Schurr B, Schneider F, Sippl C, Ratschbacher L, Minaev V, Gadoev M, Oimahmadov I, Abdybachaev U, Moldobekov B, Orunbaev S and Negmatullaev S. 2012. Crustal and uppermost mantle velocity structure along a profile across the Pamir and southern Tien Shan as derived from project TIPAGE wide-angle seismic data. Geophysical Journal International, 188(2): 385-407
Murphy MA, Yin A, Kapp P, Harrison TM, Ding L and Guo JH. 2000. Southward propagation of the Karakoram fault system, Southwest Tibet: Timing and magnitude of slip. Geology, 28(5): 451-454 DOI:10.1130/0091-7613(2000)28<451:SPOTKF>2.0.CO;2
Ratschbacher L, Frisch W, Liu GH and Chen CS. 1994. Distributed deformation in southern and western Tibet during and after the India-Asian collision. Journal of Geophysical Research, 99(B10): 19917-19945 DOI:10.1029/94JB00932
Robinson AC, Yin A, Manning CE, Harrison TM, Zhang SH and Wang XF. 2004. Tectonic evolution of the northeastern Pamir: Constraints from the northern portion of the Cenozoic Kongur Shan extensional system, western China. GSA Bulletin, 116(7-8): 953-973
Robinson AC, Yin A, Manning CE, Harrison TM, Zhang SH and Wang XF. 2007. Cenozoic evolution of the eastern Pamir: Implications for strain-accommodation mechanisms at the western end of the Himalayan-Tibetan orogen. GSA Bulletin, 119(7-8): 882-896
Robinson AC. 2009. Evidence against Quaternary slip on the northern Karakorum Fault suggests kinematic reorganization at the western end of the Himalayan-Tibetan orogen. Earth and Planetary Science Letters, 286(1-2): 158-170 DOI:10.1016/j.epsl.2009.06.025
Robinson AC, Yin A and Lovera OM. 2010. The role of footwall deformation and denudation in controlling cooling age patterns of detachment systems: An application to the Kongur Shan extensional system in the Eastern Pamir, China. Tectonophysics, 496: 28-43 DOI:10.1016/j.tecto.2010.10.003
Robinson AC, Ducea M and Lapen TJ. 2012. Detrital zircon and isotopic constraints on the crustal architecture and tectonic evolution of the northeastern Pamir. Tectonics, 31(2): TC2016 DOI:10.1029/2011TC003013
Rowley DB. 1996. Age of initiation of collision between India and Asia: A review of stratigraphic data. Earth and Planetary Science Letters, 145(1-4): 1-13 DOI:10.1016/S0012-821X(96)00201-4
Royden LH, Burchfiel BC and van der Hilst RD. 2008. The geological evolution of the Tibetan Plateau. Science, 321(5892): 1054-1058 DOI:10.1126/science.1155371
Rutte D, Ratschbacher L, Khan J, Stübner K, Hacker BR, Stearns MA, Enkelmann E, Jonckheere R, Pfänder JA, Sperner B and Tichomirowa M. 2017a. Building the Pamir-Tibetan Plateau-Crustal stacking, extensional collapse, and lateral extrusion in the Central Pamir: 2. Timing and rates. Tectonics, 36(3): 385-419
Rutte D, Ratschbacher L, Schneider S, Stübner K, Stearns MA, Gulzar MA and Hacker BR. 2017b. Building the Pamir-Tibetan Plateau-Crustal stacking, extensional collapse, and lateral extrusion in the Central Pamir: 1. Geometry and kinematics. Tectonics, 36(3): 342-384
Schmidt J, Hacker BR, Ratschbacher L, Stübner K, Stearns M, Kylander-Clark A, Cottle JM, Alexander A, Webb G, Gehrels G and Minaev V. 2011. Cenozoic deep crust in the Pamir. Earth and Planetary Science Letters, 312(3-4): 411-421 DOI:10.1016/j.epsl.2011.10.034
Schwab M, Ratschbacher L, Siebel W, McWilliams M, Minaev V, Lutkov V, Chen FK, Stanek K, Nelson B, Frisch W and Wooden JL. 2004. Assembly of the Pamirs: Age and origin of magmatic belts from the southern Tien Shan to the southern Pamirs and their relation to Tibet. Tectonics, 23(4): TC4002
Sobel ER, Schoenbohm LM, Chen J, Thiede R, Stockli DF, Sudo M and Atrecker MR. 2011. Late Miocene-Pliocene deceleration of dextral slip between Pamir and Tarim: Implications for Pamir orogenesis. Earth and Planetary Science Letters, 304(3-4): 369-378 DOI:10.1016/j.epsl.2011.02.012
Stearns MA, Hacker BR, Ratschbacher L, Lee J, Cottle JM and Kylander-Clark A. 2013. Synchronous Oligocene-Miocene metamorphism of the Pamir and the North Himalaya driven by plate-scale dynamics. Geology, 41(10): 1071-1074 DOI:10.1130/G34451.1
Stearns MA, Hacker BR, Ratschbacher L, Rutte D and Kylander-Clark ARC. 2015. Titanite petrochronology of the Pamir gneiss domes: Implications for middle to deep crust exhumation and titanite closure to Pb and Zr diffusion. Tectonics, 34(4): 784-802
Stübner K, Ratschbacher L, Rutte D, Stanek K, Minaev V, Wiesinger M, Gloaguen R and Project TIPAGE Members. 2013a. The giant Shakhdara migmatitic gneiss dome, Pamir, India-Asia collision zone: 1. Geometry and kinematics. Tectonics, 32(4): 948-979
Stübner K, Ratschbacher L, Weise C, Chow J, Hofmann J, Khan J, Rutte D, Sperner B, Pfänder JA, Hacker BR, Dunkl I, Tichomirowa M, Stearns MA and Project TIPAGE Members. 2013b. The giant Shakhdara migmatitic gneiss dome, Pamir, India-Asia collision zone: 2. Timing of dome formation. Tectonics, 32(5): 1404-1431
Tapponnier P, Mattauer M, Proust F and Cassaigneau C. 1981. Mesozoic ophiolites, sutures, and arge-scale tectonic movements in Afghanistan. Earth and Planetary Science Letters, 52(2): 355-371 DOI:10.1016/0012-821X(81)90189-8
Tapponnier P, Xu ZQ, Roger F, Meyer B, Arnaud N, Wittlinger G and Yang JS. 2001. Oblique stepwise rise and growth of the Tibet Plateau. Science, 294(5547): 1671-1677 DOI:10.1126/science.105978
Thiede RC, Sobel ER, Chen J, Schoenbohm LM, Stockli DF, Sudo M and Strecker MR. 2013. Late Cenozoic extension and crustal doming in the India-Eurasia collision zone: New thermochronologic constraints from the NE Chinese Pamir. Tectonics, 32(3): 763-779 DOI:10.1002/tect.20050
Valli F, Leloup PH, Raquette JL, Arnaud N, Li HB, Tapponnier P, Lacassin R, Guillot S, Liu DY, Deloule E, Xu ZQ and Mahéo G. 2008. New U-Th/Pb constraints on timing of shearing and long-term slip-rate on the Karakorum fault. Tectonics, 27(5): TC5007
Xu XY, Cai ZH, Xu ZQ and Cao H. 2017. The Miocene high-K alkaline rocks in Tashkorgan region, northeastern Pamir: Formation mechanism and tectonic implications. Geological Review, 63(3): 616-629 (in Chinese with English abstract)
Yang WQ, Liu L, Cao YT, Wang C, He SP, Li RS and Zhu XH. 2010. Geochronological evidence of Indosinian (high-pressure) metamorphic event and its tectonic significance in Taxkorgan area of the Western Kunlun Mountains, NW China. Science China (Earth Sciences), 53(10): 1445-1459 DOI:10.1007/s11430-010-4081-1
Yin A and Harrison TM. 2000. Geologic evolution of the Himalayan-Tibetan orogen. Annual Review of Earth and Planetary Sciences, 28: 211-280 DOI:10.1146/annurev.earth.28.1.211
Zhang CL, Zou HB, Ye XT and Chen XY. 2018. Tectonic evolution of the NE section of the Pamir Plateau: New evidence from field observations and zircon U-Pb geochronology. Tectonophysics, 723: 27-40 DOI:10.1016/j.tecto.2017.11.036
Zhang Y, Chen W, Chen KL and Liu XY. 2006. Study on the Ar-Ar age spectrum of diagenetic I/S and the mechanism of 39Ar recoil loss: Examples from the clay minerals of P-T boundary in Changxing, Zhejiang province. Geological Review, 52(4): 556-561 (in Chinese with English abstract)
侯可军, 李延河, 田有荣. 2009. LA-MC-ICP-MS锆石微区原位U-Pb定年技术. 矿床地质, 28(4): 481-492.
柯珊, 莫宣学, 罗照华, 詹华明, 梁涛, 李莉, 李文韬. 2006. 塔什库尔干新生代碱性杂岩的地球化学特征及岩石成因. 岩石学报, 22(4): 905-915.
罗照华, 莫宣学, 柯珊. 2003. 塔什库尔干碱性杂岩体形成时代及其地质意义. 新疆地质, 21(1): 46-50.
徐晓尹, 蔡志慧, 许志琴, 曹汇. 2017. 东北帕米尔塔什库尔干中新世高钾碱性岩的成因机制与大地构造意义. 地质论评, 63(3): 616-629.
张彦, 陈文, 陈克龙, 刘新宇. 2006. 成岩混层(I/S)Ar-Ar年龄谱型及39Ar核反冲丢失机理研究——以浙江长兴地区P-T界线粘土岩为例. 地质论评, 52(4): 556-561.