Cenozoic uplift process in Gangdise, Tibet:Evidence from thermal history modeling of apa-tite fission track
-
摘要:
西藏冈底斯南缘中酸性侵入岩的磷灰石裂变径迹年龄在37~25Ma之间,热史模拟过程反映冈底斯经历了3个阶段的抬升演化。40~26Ma的快速冷却抬升阶段:受控于印度-欧亚大陆完全碰撞拼合的影响,并在37~26Ma抬升至现今海拔高度;26~8Ma的剥蚀阶段:受夷平和大型逆冲推覆活动的影响,出现剥蚀和抬升交替过程;8~0Ma的缓慢冷却阶段:受南北向裂谷作用影响,出现内部差异抬升。此外,北部墨竹工卡地区和南部泽当、桑耶地区,西部桑耶地区和东部泽当地区,均具有相似的抬升过程和历史,没有明显差异,暗示冈底斯经历了整体性、较均一的阶段性抬升过程。
Abstract:The apatite fission track age obtained from the intermediate-acid intrusive rocks on the southern margin of the Gangdise is 37~25Ma. Thermal history modeling process shows that the Gangdise belt has experienced three stages of uplift evolution. The First stage (40~26Ma) was a rapid cooling uplifting stage, which was controlled by the India-Eurasia complete collision, and the Gangdise was raised to the elevation height in 37~26Ma. The second stage (26~8Ma) was an erosion stage, affected by the pedimentation and large thrust system activation, with the existence of alternate erosion and uplift process in this area. The third stage (8~0Ma) was a low cooling stage, during which activation of the NS-striking normal fault systems and rifts led to the formation of basin-mountain system and differential uplifting. In addition, there existed no obvious difference between the north and the south and between the west and the east of Gangdise, with all of these places having similar uplift processes and histories. Thermal history modeling process of apatite fission track age shows that the Gangdise belt has experienced a whole, uniform and phase uplifting.
-
Keywords:
- apatite /
- fission track age /
- thermal history modeling /
- Gangdise /
- uplift
-
钦杭结合带是扬子和华夏2个古陆块于新元古代碰撞拼接所形成的板块结合带[1-2],其南西起自广西钦州湾,经湘、赣往北东延伸至浙江杭州湾,总体呈北东向反“S”状弧形展布。钦杭成矿带是沿钦杭结合带两侧扩展而新圈定的一个重要成矿带,成矿地质条件优越。开展地质调查工作以来,相继在该成矿带发现了一批大中型矿床,如锡田锡矿、思委银矿、园珠顶铜钼矿等,找矿突破成果显著[3]。对钦杭成矿带内一些典型矿床进行系统的研究,对了解区域成矿地质条件、矿床成因、矿床分布规律等具有重要意义。
幕阜山地区位于钦杭成矿带中部,是钦杭成矿带西段部署的重点找矿地区之一,已发现的主要矿种有钨、金、铀、铍、铌、钽、铜、铅、锌等,湖南三墩铜铅锌矿位于幕阜山岩体南部[2-3]①。幕阜山岩体呈岩基状产出,出露面积达2530km2,主体部分位于湖南省境内。幕阜山岩体为多期次侵入的复式岩体②③,区域地质调查对其进行了划分,其中燕山早期岩浆活动持续时期较长(189~145Ma),在长达44Ma的地质时期内,除了侵入以粗中粒斑状二长花岗岩为主的岩体外,在岩浆后阶段还有大量花岗伟晶岩脉生成,燕山晚期岩浆活动基本上紧接着燕山早期岩浆活动发生,其活动时期在136~115Ma之间,持续时间21Ma, 燕山晚期岩体划分出第一、第二、第三次侵入体,还伴有部分花岗伟晶岩脉的生成。虽然前人对幕阜山岩体进行了成岩期次划分和同位素年龄分析,但受当时测试条件和方法本身的制约,同位素年龄存在精度不高的问题,且缺乏系统的岩石成因研究。为了获得高精度的成岩年龄数据,本文采用LA-ICP-MS锆石U-Pb定年,并结合Hf同位素和岩石地球化学分析,深入探讨岩石成因,为研究成岩与成矿的关系提供重要的科学依据。
1. 地质背景
三墩铜铅锌矿区位于下扬子陆块江南古岛弧的北缘,幕阜山岩体南部与冷家溪群侵入接触部位、北东向天府山-幕阜山断裂与枫林-浆市断裂复合处北西侧②④。矿区出露地层较简单(图 1),主要为冷家溪群变质岩,位于矿区南部,岩性以黑云母片岩、二云母片岩、石英片岩为主,其中有大量伟晶岩岩脉及石英脉穿插。区内断裂构造发育,尤其以北北西向和近南北向断裂最发育,延伸稳定、规模较大,断裂带内构造透镜体、石英脉、角砾岩等发育。区内岩浆岩十分发育,主要出露幕阜山岩体燕山晚期第一次侵入体,在矿区东北角还出露燕山早期侵入体和燕山晚期第三次侵入体,燕山早期侵入体(Mγ52Pb)岩性为片麻状粗中粒斑状黑云母二长花岗岩,燕山晚期第一次侵入体(γ53-1)岩性为中细粒二云母二长花岗岩,燕山晚期第三次侵入体(γδ53-3)岩性为细粒花岗闪长岩。
矿区内热液活动强烈,硅化显著,尤其在断裂带内硅化现象最明显,往往被石英脉充填,石英脉除少数呈透镜状产出外,大部分呈整脉状充填在断裂带内。区内已探明19条含矿石英脉,其中具有工业矿体的有5条,皆呈脉状、透镜状产出,其走向为北北西向,倾向为北东东向。围岩蚀变除硅化外,还有绢云母化、绿泥石化、萤石化,偶见碳酸盐化、重晶石化等,萤石化与矿化关系密切。一般在断层构造带的中部、断裂构造转折处、两组断裂相交处,硅化、绿泥石化强烈发育处和北北西向断裂(节理)带中,形成较富的矿体②③。
2. 岩石学特征
三墩铜铅锌矿区花岗岩岩石地球化学和测年样品(东经113°45.626′、北纬28°52.622′)均采于幕阜山岩体燕山晚期第一次侵入体的地表露头。通过观察手标本和显微镜下特征(图 2),鉴定三墩铜铅锌矿区花岗岩岩性为二云母二长花岗岩,灰白色,中细粒花岗结构,块状构造,矿物成分为石英(25% ~30%)、钾长石(30% ~35%)、斜长石(25% ~30%)、黑云母(5%~7%)、白云母(5%~8%)。
3. 分析测试方法
主量、微量元素测试在中国地质调查局武汉地质调查中心中南矿产资源监督检查中心进行(前者用X-射线荧光光谱法XRF,后者用ICP-MS法)。地球化学参数计算与图解生成使用GeoKit软件[4]。
花岗岩样品破碎和锆石挑选由廊坊市宇能岩石矿物分选技术服务有限公司完成。锆石阴极发光(CL)图像在中国地质大学(武汉)地质过程与矿产资源国家重点实验室拍摄,并结合透射光和反射光图像观察锆石内部结构。
锆石U-Pb定年分析在中国地质大学(武汉)地质过程与矿产资源国家重点实验室(GPMR)LAICP-MS仪器上完成,ICP-MS型号为Agilent7500a。激光剥蚀系统为GeoLas 2005,所用斑束直径为32μm。对分析数据的离线处理(包括对样品和空白信号的选择、仪器灵敏度漂移校正、元素含量及U-Th-Pb同位素比值和年龄计算)采用软件ICPMSDataCal [5-6]完成。详细的仪器操作条件和数据处理方法同Liu等[5-6]。锆石样品的U-Pb年龄谐和图绘制和年龄加权平均值计算均采用Isoplot 3.0完成[7]。
锆石Hf同位素原位测试在西北大学大陆动力学国家重点实验室的多接受电感耦合等离子质谱仪MC-ICP-MS(Nu Plasma)上完成。分析时采用氦作为剥蚀物质载气,激光束斑直径为42μm, 分析步骤和流程参见徐平等[8]和Yuan等[9]。采用标准锆石91500、MON-1和GJ-1作为外部标样,分析精度和误差用标样来校准,误差为2σ,本次实验测试中91500的分析结果176Hf/177Hf=0.282307 ± 0.000016,与参考值一致。
4. 分析结果
4.1 锆石U-Pb分析结果
本次研究采集的测年样品编号为SD2-2。样品CL图像显示(图 3),锆石形态大部分呈短柱状,晶形比较完整,裂纹不发育,振荡环带发育。从制备好的测年样品中选取锆石进行测试,每颗锆石一个测点,多数测点位于锆石柱体两端,少数测点在柱体中部。锆石Th含量为192×10-6~963×10-6,U含量为1581×10-6~10530×10-6,Th/U值为0.034~0.39,虽然锆石的Th/U值较低,但其具有较高的Th、U含量和存在振荡环带等特征,表明其为岩浆锆石[10]。
测年样品的同位素分析结果见表 1。共测试20个测点,除3号点数据有明显错误删除外,有效测点为19个,其中15号点年龄明显较老,14号点年龄明显较小,其余17个测试点的206Pb/238U年龄值集中分布于129.0~134.9Ma之间,投影点均落在谐和线上(图 4),206Pb/238U年龄加权平均值为131.9±1.1Ma(95%置信度,MSWD=2.3),可代表三墩铜铅锌矿区花岗岩的成岩年龄。结合阴极发光图像分析,15号测点应为继承锆石核,其206Pb/238U年龄为749.5Ma。
表 1 花岗岩LA-ICP-MS锆石U-Th-Pb同位素测试结果Table 1. LA-ICP-MS U-Th-Pb isotopic compositions of zircons for the granite分析
点含量/10-6 Th/U 同位素比值 年龄/Ma Pb Th U 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 1 231 323 9563 0.034 0.04826 0.00120 0.14259 0.00348 0.02115 0.00021 122.3 57.4 135.3 3.1 134.9 1.3 2 155 377 6425 0.059 0.04627 0.00127 0.13697 0.00381 0.02114 0.00023 13.1 63.0 130.3 3.4 134.8 1.5 3 103 243 4266 0.057 0.04430 0.00125 0.13036 0.00375 0.02102 0.00024 error error 124.4 3.4 134.1 1.5 4 174 358 7255 0.049 0.04723 0.00137 0.13570 0.00381 0.02056 0.00022 61.2 66.7 129.2 3.4 131.2 1.4 5 101 240 4245 0.057 0.04742 0.00182 0.13820 0.00545 0.02082 0.00030 77.9 88.9 131.4 4.9 132.8 1.9 6 109 279 4578 0.061 0.05025 0.00160 0.14233 0.00446 0.02024 0.00022 205.6 74.1 135.1 4.0 129.2 1.4 7 84 224 3476 0.064 0.04772 0.00160 0.13724 0.00451 0.02058 0.00023 87.1 -117.6 130.6 4.0 131.3 1.5 8 110 259 4530 0.057 0.04652 0.00146 0.13614 0.00420 0.02097 0.00021 33.4 64.8 129.6 3.8 133.8 1.3 9 81 211 3513 0.060 0.04946 0.00318 0.14136 0.00903 0.02063 0.00034 168.6 154.6 134.3 8.0 131.7 2.1 10 233 376 10530 0.036 0.04686 0.00134 0.13658 0.00453 0.02102 0.00042 42.7 66.7 130.0 4.0 134.1 2.6 11 120 375 5185 0.072 0.04865 0.00152 0.13863 0.00437 0.02049 0.00026 131.6 78.7 131.8 3.9 130.7 1.6 12 79 192 3258 0.059 0.04918 0.00166 0.14283 0.00465 0.02097 0.00021 166.8 79.6 135.6 4.1 133.8 1.3 13 148 399 6305 0.063 0.04854 0.00110 0.13748 0.00307 0.02040 0.00017 124.2 55.6 130.8 2.7 130.2 1.1 14 210 561 9618 0.058 0.04795 0.00103 0.12706 0.00283 0.01910 0.00017 98.2 45.4 121.5 2.5 122.0 1.1 15 392 963 2472 0.390 0.06553 0.00086 1.11877 0.01526 0.12330 0.00078 790.7 27.8 762.4 7.3 749.5 4.5 16 38 286 1581 0.181 0.04736 0.00244 0.13302 0.00626 0.02055 0.00029 77.9 109.3 126.8 5.6 131.1 1.8 17 141 371 6038 0.061 0.04796 0.00111 0.13634 0.00319 0.02040 0.00016 98.2 53.7 129.8 2.9 130.2 1.0 18 138 326 5740 0.057 0.04908 0.00174 0.14255 0.00407 0.02112 0.00021 150.1 83.3 135.3 3.6 134.7 1.3 19 145 358 6175 0.058 0.04931 0.00142 0.14325 0.00423 0.02089 0.00025 161.2 66.7 135.9 3.8 133.2 1.6 20 154 398 6752 0.059 0.04837 0.00138 0.13621 0.00389 0.02022 0.00017 116.8 63.9 129.7 3.5 129.0 1.1 4.2 地球化学分析结果
岩石地球化学分析结果见表 2。花岗岩SiO2含量为72.56%~74.27%(>70%),Al2O3含量为14.05%~14.46%,铝饱和指数(A/CNK)值为1.13~1.28,属强过铝质(>1.1)。碱含量(Na2O+K2O)变化范围为7.71% ~8.81%,其中K2O含量为4.22% ~5.61%,在TAS岩石分类图解(图 5-a)上[11-12],样品点落在亚碱性花岗岩系列。在岩石系列SiO2-K2O图解(图 5-b)上[13],样品点落在高钾钙碱性系列。
表 2 花岗岩主量、微量和稀土元素含量Table 2. Major, trace and rare earth element values of the granite样号
岩性SD1-4
二云母
花岗岩SD1-5
二云母
花岗岩SD2-2
二云母
花岗岩SD3-1
二云母
花岗岩SD4-1
二云母
花岗岩SiO2 72.70 72.56 74.27 73.49 73.49 Al2O3 14.45 14.46 14.18 14.35 14.05 Fe2O3 0.501 0.383 0.250 0.228 0.237 FeO 1.88 2.04 1.73 1.77 2.34 CaO 0.553 1.02 0.475 0.732 0.739 MgO 0.325 0.346 0.184 0.231 0.320 K2O 5.61 4.50 4.77 5.00 4.22 Na2O 3.20 3.66 3.08 3.12 3.49 TiO2 0.125 0.107 0.106 0.127 0.165 P2O5 0.201 0.181 0.267 0.215 0.111 MnO 0.034 0.035 0.050 0.039 0.064 灼失量 0.399 0.277 0.419 0.459 0.460 Cu 13.5 15.0 19.6 17.8 12.7 Pb 65.0 90.3 28.2 43.2 43.5 Zn 108 59.4 47.0 57.9 146 Cr 14.5 3.86 17.9 7.27 9.28 Ni 2.25 2.29 2.65 19.4 5.07 Co 1.70 1.54 1.20 1.42 2.12 W 2.19 2.14 3.12 3.27 2.29 Mo 0.85 0.55 0.89 Bi 1.04 1.25 1.45 4.82 0.22 Sr 42.4 63.3 18.1 33.3 38.8 Ba 251 199 78.9 163 226 Nb 17.3 14.8 18.6 Ta 2.86 2.57 4.89 Zr 48.0 59.9 51.2 Hf 2.26 2.52 1.97 U 17.2 8.02 8.00 Th 6.55 9.92 10.7 La 22.8 15.7 12.7 17.4 20.7 Ce 36.1 27.2 26.6 35.8 40 Pr 5.5 4.02 3.16 4.37 4.73 Nd 18.9 14 10.9 15.1 16.2 Sm 4.49 3.39 2.67 3.81 3.78 Eu 0.47 0.4 0.25 0.46 0.54 Gd 3.83 2.93 2.26 3.32 3.39 Tb 0.58 0.46 0.45 0.56 0.57 Dy 2.56 2.05 2.5 2.54 2.64 Ho 0.36 0.29 0.44 0.35 0.42 Er 0.87 0.73 1.24 0.81 1.09 Tm 0.12 0.093 0.24 0.11 0.2 Yb 0.74 0.58 1.51 0.74 1.39 Lu 0.094 0.071 0.19 0.091 0.19 Y 9.27 7.49 12.7 9.43 11.8 A/CNK 1.17 1.13 1.28 1.21 1.21 ∑REE 106.68 79.40 77.81 94.89 107.64 δEu 0.35 0.39 0.31 0.40 0.46 δCe 0.79 0.84 1.03 1.01 0.99 注:主量元素含量单位为%,微量和稀土元素含量为10-6 花岗岩稀土元素总量(∑REE)为77.81×10-6~107.64×10-6,平均为80.56×10-6,δEu=0.31~0.46,δCe=0.79~1.03。在微量元素蛛网图(图 6-a)上,富集U、Ta、Pb等元素,亏损Ba、Nb、Sr、Zr、Ti等元素。Nb/Ta值为3.80~6.05,平均为5.20,低于地壳Nb/Ta值,指示源区具有地壳性质。Sr亏损指示斜长石的分离结晶;Ti亏损指示钛铁矿的分离结晶,暗示岩浆物质来源于地壳[14-16]。在稀土元素球粒陨石标准化配分图(图 6-b)上,稀土元素为右倾配分模式,轻、重稀土元素分异强烈,富集轻稀土元素,重稀土元素平坦分布,具有弱负Eu异常特征。
4.3 Hf同位素分析结果
对测年样品锆石19个有效测点中的17个测点(包含15号继承锆石核,除去5号和14号测点)进行了原位Hf同位素分析,除继承锆石用测点的年龄计算外,其余锆石Hf同位素计算所用的年龄为该样品的206Pb/238U年龄加权平均值,分析结果见表 3。测年样品的16个燕山晚期花岗岩锆石的微区原位Hf同位素组成相对均匀,初始176Hf/177Hf值较一致,分布在0.282526~0.282624之间,平均值为0.282589,εHf(t)值集中分布在-5.9~-2.4之间,平均值为-3.7,二阶段模式年龄(T2DM)在1558~1338Ma之间,平均值为1417Ma。15号继承锆石核初始176Hf/177Hf值为0.282453,εHf(t)值为4.8,一阶段模式年龄(T1DM)为1130Ma, 二阶段模式年龄(T2DM)为1355Ma。
表 3 锆石Lu-Hf同位素分析结果Table 3. Lu-Hf isotope data of zircons分析点 176Hf/177Hf
比值2σ 176Lu/177Hf
比值2σ 176Yb/177Hf
比值2σ t/Ma εHf(t) T1DM/Ma T2DM/Ma fLu/Hf 1 0.282554 0.000018 0.000975 0.000003 0.026140 0.000080 131.9 -4.9 988 1496 -0.97 2 0.282570 0.000017 0.000971 0.000008 0.025686 0.000242 131.9 -4.4 966 1460 -0.97 4 0.282526 0.000021 0.000974 0.000004 0.025628 0.000105 131.9 -5.9 1028 1558 -0.97 6 0.282584 0.000016 0.000732 0.000009 0.019666 0.000251 131.9 -3.8 939 1426 -0.98 7 0.282605 0.000019 0.000830 0.000020 0.021906 0.000564 131.9 -3.1 912 1380 -0.98 8 0.282600 0.000014 0.001043 0.000011 0.027627 0.000318 131.9 -3.3 925 1392 -0.97 9 0.282597 0.000019 0.000897 0.000002 0.024262 0.000049 131.9 -3.4 926 1399 -0.97 10 0.282624 0.000016 0.001042 0.000014 0.027703 0.000383 131.9 -2.4 890 1338 -0.97 11 0.282614 0.000028 0.001068 0.000004 0.029513 0.000140 131.9 -2.8 906 1361 -0.97 12 0.282588 0.000017 0.000782 0.000004 0.020700 0.000127 131.9 -3.7 935 1418 -0.98 13 0.282616 0.000016 0.001003 0.000014 0.026896 0.000377 131.9 -2.7 901 1355 -0.97 15 0.282453 0.000018 0.000985 0.000004 0.023484 0.000117 749.5 4.8 1130 1355 -0.97 16 0.282575 0.000015 0.000511 0.000002 0.013413 0.000053 131.9 -4.1 946 1445 -0.98 17 0.282595 0.000019 0.001079 0.000005 0.028965 0.000139 131.9 -3.4 932 1403 -0.97 18 0.282541 0.000021 0.000971 0.000006 0.025835 0.000161 131.9 -5.4 1006 1524 -0.97 19 0.282619 0.000021 0.000940 0.000002 0.024897 0.000063 131.9 -2.6 896 1349 -0.97 20 0.282613 0.000016 0.001165 0.000002 0.031296 0.000053 131.9 -2.8 910 1364 -0.96 5. 讨论
5.1 成岩年龄
幕阜山复式岩体最早的同位素年龄测定是在20世纪80年代区域地质调查时完成的,其中报道有南江桥附近的铷锶(全岩)同位素年龄样品等时线年龄为189±29.7Ma, 黄龙山岩体的白云母K-Ar年龄为139Ma②③,由于当时条件制约,存在同位素年龄精度不高的问题。锆石中Pb的扩散封闭温度高达900℃[17],是目前测定岩浆结晶年龄最理想的矿物之一,锆石结晶年龄代表了岩石成岩年龄[10]。本次选取与成矿关系密切的燕山晚期第一次侵入体开展锆石LA-ICP-MS定年工作,获得了燕山晚期第一次侵入体花岗岩的年龄为131.9±1.1Ma, 为准确限定幕阜山岩体燕山晚期第一次侵入体的年龄提供了重要的年代学依据。
本次工作发现了新元古代继承锆石核的存在,其206Pb/238U年龄为749.5Ma, 与前人研究一致,表明湖南地区的中酸性岩体中存在新元古代—新太古代继承锆石核[18-20]。
5.2 物质来源
近年来的研究发现,锆石原位Hf同位素分析是揭示地壳演化和示踪岩浆源区的重要手段[21]。三墩铜铅锌矿区花岗岩燕山晚期锆石的176Hf/177Hf值为0.282526~0.282624,εHf(t)值为-5.9~-2.4(图 7-a),Hf同位素二阶段模式年龄在1558~1338Ma之间。在t-εHf(t)图(图 7-b)上,花岗岩样品燕山晚期锆石测点均落在球粒陨石和下地壳演化线之间,表明其物质来源为中元古代的古老地壳岩石部分熔融。
继承锆石核的176Hf/177Hf值为0.282453,对应的εHf(t)值为4.8,二阶段模式年龄(T2DM)为1355Ma, 投影点落在亏损地幔和球粒陨石之间,表明其岩浆源区为直接源于亏损地幔分异的新生地壳的迅速重熔,新生地壳年龄为中元古代。
综合证据表明,花岗岩主要物质来源为中元古代地壳岩石的部分熔融,岩浆源区或上升通道可能有新元古代幔源物质加入。燕山晚期锆石和继承锆石的二阶段模式年龄都为中元古代,指示中元古代是该区重要的地壳生长阶段。
5.3 构造背景
钦杭成矿带是一条巨型的构造岩浆活动带,学者们普遍认为其在新元古代完成了拼接,但在中侏罗世又复活[22-23]。中侏罗世是一个重要的构造转折期,华南乃至中国东部发生了特提斯体制向太平洋体制的转变,随着太平洋板块持续俯冲,华南板块受到强烈挤压,地壳缩短加厚,产生一系列北东向断裂,原先不活动的断裂又重新活动;大约在早白垩世太平洋板块运动方向由原来的斜俯冲转向几乎平行大陆边缘运动,中国大陆乃至东亚大陆边缘转换为伸展阶段,区域上对应出现了基性岩浆活动,代表性的有蕉溪岭煌斑岩,其成岩年龄在136Ma左右[24]。沿钦杭结合带,是构造相对薄弱的古板块拼接带,由于强烈挤压活化,壳幔相互作用更加活跃,致使钦杭成矿带成为华南中生代大规模岩浆活动和成矿作用最集中的地区[2,22]。
前人研究表明,岩石圈减薄和软流圈上涌是华南中生代岩浆作用形成的主要机制[25-28],笔者认为,太平洋板块的运动方式制约了华南中生代大规模的岩浆活动和成矿作用,其中中晚侏罗世受太平洋板块挤压加厚的中下地壳发生部分熔融,局部出现岩石圈垮塌,形成大量火山岩浆活动;早白垩世太平洋板块运动方向发生转向,发生了强烈挤压向伸展的构造转换,出现大规模的岩石圈减薄和软流圈上涌事件,加速了中下地壳的熔融,爆发了大规模岩浆活动和成矿作用,幕阜山岩体燕山晚期第一次侵入体就是其成岩事件爆发的响应。三墩铜铅锌矿区花岗岩岩石地球化学特征具有壳源特征,Hf同位素特征表明,其主要物质来源为中元古代地壳岩石的部分熔融,在岩浆熔融的过程中有少量幔源物质的加入。三墩铜铅锌矿区花岗岩可能是由于中下地壳的部分熔融岩浆形成后,混入少量幔源物质上侵形成的。
6. 结论
(1)通过对LA-ICP-MS锆石U-Pb定年,获得了三墩铜铅锌矿区幕阜山岩体燕山晚期第一次侵入体的年龄加权平均值为131.9±1.1Ma, 成岩年龄为燕山晚期,继承锆石核的206Pb/238U年龄为749.5Ma, 暗示区内可能存在新元古代岩浆活动。
(2)通过岩石地球化学研究,确定三墩铜铅锌矿区花岗岩为一套强过铝质高钾钙碱性系列花岗岩,富U、Ta、Pb等元素,贫Ba、Nb、Sr、Zr、Ti。具有向右倾斜的稀土元素配分模式,富集轻稀土元素,重稀土元素分布平坦,具弱负Eu异常。
(3)结合Hf同位素分析,确定三墩铜铅锌矿区花岗岩物质来源主要为中元古代地壳物质重熔,并混入少量的幔源物质。其成岩可能是在大范围的岩石圈伸展背景下,发生了岩石圈拆离、软流圈上涌,中下地壳加热,部分熔融形成花岗岩浆上侵。
-
图 1 研究区地质简图(据参考文献①修改)
YZS—雅鲁藏布江缝合带;BNS—班公湖-怒江缝合带;XJS—金沙江缝合带;MMT—米拉山逆冲断裂;SDT—当雄南部逆冲断裂;GT—冈底斯逆冲断裂;RZT—仁布-泽当逆冲断裂
Figure 1. Simplified geological map of the study area
表 1 研究区磷灰石裂变径迹年龄
Table 1 Apatite fission track age of the study area
样品 岩性 现今海
拔/m粒数 ρs/(105·cm-2)
(Ns)ρi/(105·cm-2)
(Ni)ρd/(105·cm-2)
(Nd)P(χ2)
/%径迹长度/m
(N)中值年龄/a
(±1σ)池年龄/Ma
(±1σ)泽当Ⅰ区 B075-2 花巧闪长岩 3823 28 3.001
(620)21.014
(4342)10.435
(7312)9.8 12.7±2.1
(115)31±3 31±2 B117-1 花巧闪长岩 4019 30 2.165
(364)17.05
(2866)10.45
(7312)71 12.8±2.2
(101)27±2 28±2 B121-1 二长闪长岩 3567 28 1.529
(248)10.031
(1627)10.457
(7312)75 12.5±2.3
(49)33±3 33±3 泽当Ⅱ区 B118-1 石英二长岩 3564 29 2.738
(446)22.55
(3673)10.453
(7312)54.8 13.2±1.8
(61)26±2 26±2 B104-3 石英二长岩 3557 28 1.998
(774)17.106
(6625)10.447
(7312)40.4 13.2±2.1
(97)25±2 25±2 桑耶Ⅲ区 B122-1 二长闪长岩 4329 26 1.36
(156)8.769
(1006)10.433
(7312)89.6 / 34±4 34±4 B122-7 石英二长闪长岩 4309 30 1.421
(125)8.241
(725)10.43
(7312)99.7 12.6±2.4
(34)37±4 37±4 B123-1 石英二长闪长岩 4040 29 1.043
(234)6.682
(1499)10
(7312)10.4 12.6±2.3
(89)32±4 32±3 墨竹工卡Ⅳ区 JM22-2 花岗斑岩 3961 12 0.637
(43)6.501
(439)12.666
(7312)24.5 / 26±5 26±4 JM24-1 花岗斑岩 5022 30 1.095
(150)11.027
(1511)12.593
(7312)84.4 12.3±2.0
(68)26±3 26±3 JM27-1 石英二长闪长玢岩 3965 28 0.460
(171)3.910
(1452)12.519
(7312)68.9 12.3±1.9
(104)31±3 31±3 注:ρs、ρi和ρd分别为自发径迹密度、诱发径迹密度和标准径迹密度;Ns、Ni和Nd分别对应自发径迹数,诱发径迹数和标准径迹数;P(χ2)为χ2检验值 表 2 磷灰石裂变径迹年龄古地表海拔
Table 2 Paleo-altitude of the apatite fission track age sample
样品编号 总抬升/m 古埋深/m 现今样品海拔/m 古地表海拔/m 泽当Ⅰ区 B075-2 1821 -3140 3823 5142 B 117-1 2022 -3140 4019 5137 B 121-1 1582 -3140 3567 5125 泽当Ⅱ区 B 104-3 1882 -3140 3557 4815 B 118-1 1582 -3140 3564 5122 桑耶Ⅲ区 B 122-7 1882 -3140 4239 5497 B 123-1 1702 -3140 4040 5478 墨竹工卡Ⅳ区 JM24-1 1731 -3140 5022 6431 JM27-1 1883 -3140 3956 5213 -
Dewey J F, Shackleton R M, Chang C F, et al. The tectonic evolution of the Tibetan Plateau[J]. Philosophical Transactions of the Royal Society of London A:Mathematical, Physical and Engineering Sciences, 1988, 327(1594):379-413. doi: 10.1098/rsta.1988.0135
Harrison T M, Copeland P, Kidd W S F, et al. Raising Tibet[J]. Science, 1992, 263:1663-1670.
Yin A, Harrison T M, Ryerson F J, et al. Tertiary structural evolution of the Gangdese thrustsystem, southeastern Tibet[J]. Journal of Geophysical Research, 1994, 99(B9):18175-18201. doi: 10.1029/94JB00504
Copeland P, Harrison T M, Kidd W S F, et al. Rapid early Miocene acceleration of uplift in the Gangdese Belt Xizang (southern Tibet), and its bearing on accommodation mechanisms of the India-Asia collision[J]. Earth and Planetary Science Letters, 1987, 86:240-252. doi: 10.1016/0012-821X(87)90224-X
Spicer R A, Harris N B W, Widdowson M W, et al. Constant elevation of Southern Tibet over the past15 million years[J].Nature, 2003, 421:622-624. doi: 10.1038/nature01356
Rowley D B, Currie B S. Palaeo-altimetry of the late Eocene to-Miocene Lunpola basin, central Tibet[J]. Nature, 2006, 439(9):677-681. http://www.ncbi.nlm.nih.gov/pubmed/16467830
崔之久, 高全洲, 刘耕年, 等.夷平面、古岩溶与青藏高原隆升[J].中国科学(D辑), 1996, 26(4):378-386. http://youxian.cnki.com.cn/yxdetail.aspx?filename=YNDL201703004&dbname=CJFDPREP 李吉均, 文世宣, 张青松, 等.青藏高原隆起的时代、幅度和形式的探讨[J].中国科学(A辑), 1979, 6:608-616. http://www.cnki.com.cn/Article/CJFDTOTAL-JAXK197906008.htm 李廷栋.青藏高原隆升的过程和机制[J].地球学报, 1995, 1:1-9. http://youxian.cnki.com.cn/yxdetail.aspx?filename=YNDL201703004&dbname=CJFDPREP 钟大赉, 丁林.青藏高原的隆起过程及其机制探讨[J].中国科学(D辑), 1996, 26(4):289-295. http://youxian.cnki.com.cn/yxdetail.aspx?filename=YNDL201703004&dbname=CJFDPREP 刘顺生, 张峰.西藏南部地区的裂变径迹年龄和上升速度的研究[J].中国科学(B辑), 1987, 9:1000-1010. http://www.cnki.com.cn/Article/CJFDTOTAL-JBXK198709011.htm 江万, 莫宣学, 赵崇贺, 等.矿物裂变径迹年龄与青藏高原隆升速率研究[J].地质力学学报, 1998, 4(1):13-18. http://www.cnki.com.cn/Article/CJFDTOTAL-DZLX801.001.htm 袁万明, 王世成, 李胜荣, 等.西藏冈底斯带构造活动的裂变径迹证据[J].科学通报, 2001, 45(20):1739-1742. doi: 10.3321/j.issn:0023-074X.2001.20.017 Gleadow A J W, Duddy I R. A natural long-term track annealing experiment for apatite[J]. Nuclear Tracks, 1981, 5:169-174. doi: 10.1016/0191-278X(81)90039-1
Nacser C W. The fading of fission tracks in the geologic environment-data from deep drill holes[J]. Nuclear Tracks, 1981, 5:248-250. doi: 10.1016/0191-278X(81)90055-X
Copeland P, Harrison M, Yun P, et al. Thermal evolution of the Gangdese batholith, southern Tibet:A history of episodic unroofing[J]. Tectonics, 1995, 14(2):223-236. doi: 10.1029/94TC01676
陈文寄, 李齐, 郝杰, 等.冈底斯岩带热演化史的MDD模式新证据Ⅱ[J].科学通报, 1998, 43(21):2332-2336. doi: 10.3321/j.issn:0023-074X.1998.21.022 陈文寄, 李齐, 郝杰, 等.冈底斯岩带结晶后的热演化史及其构造含义[J].中国科学(D辑), 1999, 29(1):9-15. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK199901001.htm Harrison T M, Yin A, Grove M. The ZedongWindow:A record of superposedTertiaryconvergencein southeasternTibet[J]. Journal of Geophysical Research, 2000, 105(B8):19211-19230. doi: 10.1029/2000JB900078
袁万明, 候增歉, 李胜荣, 等.雅鲁藏布江逆冲带活动的裂变径迹定年证据[J].科学通报, 2002, 47(2):147-150. http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB200202016.htm 向树元, 张士贞, 胡敬仁, 等.西藏米拉山断裂活动的磷灰石裂变径迹证据[J].地球科学-中国地质大学学报, 2012, 37, 增刊2:39-46. http://www.cnki.com.cn/Article/CJFDTOTAL-DQKX2012S2011.htm 吴珍汉, 江万, 周继荣, 等.青藏高原腹地典型岩体热历史与构造-地貌演化过程的热年代学分析[J].地质学报, 2001, 75(4):468-476. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200104007.htm 丁林, 钟大赉, 潘裕生, 等.东喜马拉雅构造结上新世以来快速抬升的裂变径迹证据[J].科学通报, 1995, 40(16):1497-1500. doi: 10.3321/j.issn:0023-074X.1995.16.018 赵志丹, 莫宣学, 郭铁鹰, 等.西藏南部岩体裂变径迹年龄与高原隆升[J].自然科学进展, 2003, 13(8):877-880. http://www.cnki.com.cn/Article/CJFDTOTAL-ZKJZ200308019.htm 王立成, 魏玉帅.西藏羌塘盆地白垩纪中期构造事件的磷灰石裂变径迹证据[J].岩石学报, 2013, 29(3):1039-1047. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201303025.htm Lu L, Zhao Z, Wu Z H, et al. Fission track thermochronology evidence for the Cretaceous and Paleogene tectonic event of Nyainrong Microcontinent, Tibet[J]. Acta Geologica Sinica, 2015, 89(1):133-144. doi: 10.1111/1755-6724.12400
Yin A, Harrison T M, Murphy M A, et al. Tertiary deformation history of southeastern and southwestern Tibet during the Indo-Asian collision[J]. Geological Society of America Bulletin, 1999, 111(11):1644-1664. doi: 10.1130/0016-7606(1999)111<1644:TDHOSA>2.3.CO;2
Wu Z H, Hu D G, Ye P S. Thrusting of the North Lhasa Block in the Tibetan Plateau[J]. Acta Geologica Sinica, 2004, 78(1):246-259. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzxw200401031&dbname=CJFD&dbcode=CJFQ
叶培盛. 拉萨地块中部蛇绿岩与逆冲推覆构造[D]. 中国地质科学院博士学位论文, 2004. http://cdmd.cnki.com.cn/Article/CDMD-82501-2007213421.htm Yuan W M, Zhang X T, Dong J Q, et al. A new vision of the intracontinental evolution of the eastern Kunlun Mountains, Northern Qinghai-Tibet Plateau, China[J]. Radiation Measurements, 2003, 36:357-362. doi: 10.1016/S1350-4487(03)00151-3
Bellemans F, Decorte F, Denhaute P V. Composition of SRM and CN U-doped glasses:significance for their use as thermal neutron fluence monitors in fission track dating[J]. Radiation Measurements, 1995, 24(2):153-160. doi: 10.1016/1350-4487(94)00100-F
Yuan W M, Dong J Q, Carter A, et al. Mesozoic-Tertiary exhumation history of the Altai Mountains, northern Xinjiang, China:constraints from apatite fission track data[J]. Tectonophysics, 2006, 412:183-193. doi: 10.1016/j.tecto.2005.09.007
Hurford A J, Green P F. The zeta age calibration of fission-track dating[J]. Chemical Geology, 1983, 41:285-317. doi: 10.1016/S0009-2541(83)80026-6
Green P F, Duddy I R, Gleadow A J W, et al. Thermal annealing of fission tracks in apatite:1. A qualitative description[J]. Chemical Geology:Isotope Geoscience section, 1986, 59:237-253. doi: 10.1016/0168-9622(86)90074-6
Green P F. A New Look at Statistics in Fission Track Dating[J]. Nucl Tracks, 1981, 5:77-86. doi: 10.1016/0191-278X(81)90029-9
Galbraith R F, Laslett G M. Statistical models for mixed fissiontrack ages[J]. Nuclear Tracks and Radiation Measurements, 1993, 2(4):459-470.
Brandon M. Decomposition of mixed grain age distributions using BINOMFIT[J].On Track, 2002, 24:13-18.
赵珍, 吴珍汉, 胡道功, 等.青藏高原冈底斯带南缘泽当多金属矿田多期岩浆活动及年代意义[J].地球学报, 2014, 35(6):703-712. doi: 10.3975/cagsb.2014.06.05 Armstrong P A.Thermochronometers in sedimentary basins[J]. Reviews in Mineralogy and Geochemistry, 2005, 58(1):499-525. doi: 10.2138/rmg.2005.58.19
Ketcham R A, Donelick R A, Carlson W D. Variability of apatite fission-track annealing kinetics Ⅲ:Extrapolation to geological time scales[J]. American Mineralogist, 1999, 84:1235-1255. doi: 10.2138/am-1999-0903
Ketcham R A, Donelick R A, Donelick M B. AFT Solve:A program for multi-kinetic modeling of apatite fission-track data[J]. Geological Materials Research, 2000, 2(1):1-32. http://ammin.geoscienceworld.org/content/88/5-6/929
陈文寄, 李齐, 周新华, 等.西藏高原南部两次快速冷却事件的构造含义[J].地震地质, 1996, 18(2):109-115. http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ602.002.htm Quidelleur X, Grove M, Lovera O. Thermal evolution and slip history of the Renbu Zedong Thrust, southeastern Tibet[J]. Journal of Geophysical Research, 1997, 102(B2):2659-2679. doi: 10.1029/96JB02483
Guillot S, Hodges K, Lefort P, et al. New constraints on the age of the Manasluleucogranite:Evidence for episodic tectonic denudation in the Central Himalayas[J]. Geology, 1994, 22(6):559-562. doi: 10.1130/0091-7613(1994)022<0559:NCOTAO>2.3.CO;2
Coleman M, Hodges K. Evidence for Tibetan plateau uplift before 14Myr ago from a new minimum age for east-west extension[J]. Nature, 1995, 374(6517):49-52. doi: 10.1038/374049a0
Murphy M A, Harrison T M. Relationship between leucogranites and the Qomolangma detachment in the Rongbuk valley, South Tibet[J]. Geology, 1999, 27(9):831-834. doi: 10.1130/0091-7613(1999)027<0831:RBLATQ>2.3.CO;2
丁林, 岳雅慧, 蔡福龙, 等.西藏拉萨地块高镁超钾质火山岩及对南北向裂谷形成时间和切割深度的制约[J].地质学报, 2006, 80(9):1252-1261. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200609003.htm Harrison T M, Copeland P, Kidd W S F, et al. Activation of the Nyainqentanghla shear zone:implications for uplift of the southern Tibetan Plateau[J]. Tectonics, 1995, 14(3):658-676. doi: 10.1029/95TC00608
吴珍汉, 江万, 周继荣, 等.青藏高原腹地典型岩体热历史与构造-地貌演化过程的热年代学分析[J].地质学报, 2001, 75(4):468-476. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200104007.htm 吴珍汉, 江万, 吴中海, 等.青藏高原腹地典型盆-山构造形成时代[J].地球学报, 2002, 23(4):289-294. http://www.cnki.com.cn/Article/CJFDTOTAL-DQXB200204000.htm 袁万明, 董金泉, 保增宽, 等.西藏冈底斯地块尼木地区新第三纪构造热史的磷灰石裂变径迹约束[J].原子能科学技术, 2008, 42(6):570-573. http://www.cnki.com.cn/Article/CJFDTOTAL-YZJS200806023.htm 祝嵩. 雅鲁藏布江河谷地貌与地质环境演化[D]. 中国地质科学院博士学位论文, 2012. http://cdmd.cnki.com.cn/Article/CDMD-82501-1012371268.htm Yin A, Harrison T M. Geologic evolution of the Himalayan-Tibetan orogen[J]. Annu. Rev. Earth Planet Sci., 2000, 28:211-280. doi: 10.1146/annurev.earth.28.1.211
吴珍汉, 叶培盛, 胡道功, 等.拉萨地块北部逆冲推覆构造系统[J].地质论评, 2003, 49(1):74-81. http://www.cnki.com.cn/Article/CJFDTOTAL-DZLP200301010.htm Currie B S, Rowley D, Tabor N. Middle Miocene paleoaltimetry of southern Tibet:Implications for the role of mantle thickening and delamination in the Himalayan orogeny[J].Geology, 2005, 33(3):181-184. doi: 10.1130/G21170.1
Copeland P, Harrison T M, Yun P. et al. Thermal evolution of the Gangdese batholith, southern Tibet:a history of episodic unroofing[J]. Tectonics, 1995, 14(2):223-236. doi: 10.1029/94TC01676
Wang Q, Wyman D A, Xu J F, et al. Eocene melting of subducting continental crust and early uplifting of central Tibet:Evidence from central-western Qiangtang high-K calc-alkaline andesites, dacitesandrhyolites[J]. Earth and Planetary Science Letters, 2008, 272(1/2):158-171. http://www.sciencedirect.com/science/article/pii/S0012821X0800277X
Wang Q, Wyman D A, Li Z X, et al. Eocene north-south trending dikes in central Tibet:New constraints on the timing of eastwest extension with implications for early plateau uplift?[J]. Earth and Planetary Science Letters, 2010, 298:205-216. doi: 10.1016/j.epsl.2010.07.046
吴珍汉, 吴中海, 胡道功, 等.青藏高原古大湖与夷平面的关系及高原面形成演化过程[J].现代地质, 2009, 23(6):993-1002. http://www.cnki.com.cn/Article/CJFDTOTAL-XDDZ200906001.htm 李亚林, 王成善, 伊海生, 等.西藏北部新生代大型逆冲推覆构造与唐古拉山的隆起[J].地质学报, 2006, 80(8):1118-1130. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200608009.htm 中国地质科学院成都地质矿产研究所. 青藏高原及邻区地质图(1: 1500000). 1988.