拉萨地体巴松措片岩变质条件与时代：对其早侏罗世地壳增厚的约束

Metamorphic P-T Conditions and Age of Schist in Basong Tso， Southern Tibet ： Implications for Early Jurassic Crustal Thickening of the Lhasa Terrane

摘要

南、北拉萨地体于晚三叠世—早侏罗世碰撞拼贴，但其碰撞后的变质作用时限、条件与动力学过程仍缺乏精确约束。位于二者之交的巴松措地区是研究该过程的理想区域，然而已有研究对其碰撞后约200 Ma以来的变质记录缺乏系统认识。本文对巴松措地区石榴子石云母片岩开展了详细的岩相学、全岩地球化学、锆石U-Pb年代学与相平衡模拟研究。结果表明：（1）变质成因锆石边部U-Pb加权平均年龄为（191.84 ± 0.84）Ma（MSWD = 0.25，n = 12），限定了该期变质作用发生于早侏罗世。（2）相平衡模拟获得峰期变质条件为589.5 ℃和12.4 kbar，并揭示了一条具顺时针演化特征的P-T轨迹，指示了增压升温后近等温折返的过程。（3）全岩地球化学特征指示原岩为大陆碰撞背景下的泥质沉积岩。本研究在巴松措地区识别出（191.84 ± 0.84）Ma的高压变质事件，其较高的压力与顺时针P-T轨迹共同揭示了南、北拉萨地体碰撞后地壳持续增厚与热松弛滞后的动力学过程，为完善拉萨地体构造演化模型提供了关键的温压-年代学约束。

Abstract

The collision between the South and North Lhasa terranes occurred from the Late Triassic to the Early Jurassic， yet the timing， conditions， and dynamic processes of post-collisional metamorphism remain poorly constrained. The Basong Tso region， located at the boundary between the two terranes， is an ideal area for studying this process. However， existing studies lack a systematic understanding of the metamorphic record since ca. 200 Ma （i.e.， after the collision）. In this paper， we present detailed petrographic， whole-rock geochemical， zircon U-Pb geochronological， and phase equilibrium modeling analyses on garnet-mica schist from the Basong Tso region. The results show that： （1） The weighted mean U-Pb age of metamorphic zircon rims is （191.84 ± 0.84） Ma （MSWD = 0.25， n = 12）， constraining the timing of this metamorphic event to the Early Jurassic. （2） Phase equilibrium modeling yields peak metamorphic conditions of 589.5 ℃ and 12.4 kbar and reveals a clockwise P-T path， indicating near-isothermal exhumation following an increase in pressure and temperature. （3） Whole-rock geochemical characteristics suggest that the protolith formed in a continental collision setting. This study identifies a （191.84 ± 0.84） Ma high-pressure metamorphic event in the Basong Tso region. The elevated pressure and clockwise P-T path collectively reveal the dynamic process of continuous crustal thickening and delayed thermal relaxation after the collision of the South and North Lhasa terranes， providing critical P-T-t constraints for refining the tectonic evolution model of the Lhasa terrane.

Graphical abstract

关键词

拉萨地体 / 高压变质 / 变质相图 / 地壳增厚 / 相平衡模拟

Key words

Lhasa Terrane / high-pressure metamorphic / metamorphic phase diagram / crustal thickening / phase equilibrium modeling

0 引 言

拉萨地体作为青藏高原的核心组成部分，起源于冈瓦纳大陆北缘，其演化历史经历了多期的裂解、漂移，最终向北拼贴至欧亚大陆南缘^{［1 - 9］}。拉萨地体中部的唐加—松多俯冲增生杂岩带（即洛巴堆—米拉山断裂带）将其进一步划分为北拉萨和南拉萨两个主要构造单元^［10］。现有研究表明，拉萨地体经历了复杂的变质演化历史，主要包括四个阶段：新元古代和晚古生代洋壳俯冲背景下的高压变质作用；早古生代和早中生代陆-陆碰撞带的高压和中压变质作用；晚白垩世与洋中脊俯冲相关的中压变质作用；以及新生代两期与地壳加厚有关的中压变质作用^［6，11］。这些变质事件记录了特提斯洋演化及印度—亚洲碰撞造山的关键信息，对构建青藏高原南部的构造演化模型具有重要意义。

巴松措地区位于南、北拉萨地体交界处，完整记录了松多古特提斯洋闭合至南、北拉萨地体碰撞拼合的地质过程，是研究这一关键阶段的理想窗口。Weller等^［12］在巴松措地区识别出了两条性质不同的变质带：北部的扎拉（Zhala）变质带由副片麻岩和花岗岩组成，记录了198 ~ 192 Ma期间Buchan型伸展背景下的低压高温变质作用；南部的巴松措（Basong Tso）变质带则由剪切片岩和正片麻岩组成，记录了204 ~ 201 Ma期间Barrovian型大陆碰撞背景下的中等温压变质作用。然而，目前的年代学与岩石学资料尚不足以完整重建南、北拉萨地体碰撞后的精细热力学过程，特别是对于碰撞后地壳是否经历了持续的增厚作用及其对应的P-T演化轨迹，仍缺乏系统的约束^［13］。

为了探讨上述科学问题，本研究选取巴松措地区对温压条件变化极为敏感的石榴子石云母片岩为研究对象，综合运用岩相学、全岩地球化学、锆石U-Pb年代学以及相平衡模拟方法，精确限定该区域变质作用的时限与P-T轨迹。本研究旨在通过解析变质作用的温压条件与动力学机制，揭示南、北拉萨地体碰撞后地壳热演化的深部过程，为完善拉萨地体碰撞造山模型提供新的关键证据。

1 地质背景

青藏高原由多个自北向南依次拼贴的地体构成，主要包括昆仑地体、松潘—甘孜地体、羌塘地体、拉萨地体和喜马拉雅地体。这些地体之间被一系列主缝合带所分隔，如阿尼玛卿缝合带、金沙江缝合带（JS）、班公湖—怒江缝合带（BN）以及印度河—雅鲁藏布江缝合带（IYZ）（图1（a））。拉萨地体位于青藏高原南部，南北介于班公湖—怒江缝合带与雅鲁藏布江缝合带之间，呈近东西向展布，南北宽100 ~ 300 km，东西延伸超过2000 km^［5，14］（图1（b））。其内部可识别出三条近东西向展布的中高级变质带^［6］：北部为北拉萨地体内的新元古代高压（HP）变质带；中部为位于南、北拉萨地体边界附近的晚二叠世高压（HP）—三叠世中压（MP）变质带；南部为南拉萨地体南缘的晚白垩世高温/中压（HT/MP）—始新世中压（MP）变质带。

巴松措地区（图1（c））主要由两条近东西向展布的晚三叠世—早侏罗世变质带组成，这两条变质带被约5 km宽的变沉积岩（变浊积岩）分隔，其界限可能为断层^［15］。北部带主要由扎拉杂岩组成，该杂岩带主要由副片麻岩和泥质混合岩组成，伴有少量大理岩和角闪岩。该杂岩带记录了早侏罗世（约192 Ma）的变质峰期，前人根据副片麻岩中低Th/U锆石边部的U-Pb分析结果，限定其峰期变质时代为192 Ma^［15］。扎拉杂岩带内发育部分变形的钾长石巨斑S型花岗岩侵入体，表明其构造-热演化历史复杂。南部带由巴松措杂岩组成，该杂岩带主要由含蓝晶石的泥质混合岩和正片麻岩构成，发育强烈的糜棱岩化构造，其内正片麻岩的原岩结晶时代为197 ~ 203 Ma，副片麻岩记录的变质时代为204 ~ 193 Ma。两个高级变质带均经历了复杂的多期变形，区域上以区域动力热变质作用为主。

2 采样与分析方法

2.1 样品与岩相学

本研究样品（LCBR21-05）采自南拉萨地体巴松措杂岩带内（图1（c）），岩性为石榴子石云母片岩，采样点GPS坐标为30°00′27.05″N，94°13′00.13″E。样品具典型的斑状变晶结构和片状构造（图2（a））。变斑晶主要为自形-半自形石榴子石，粒径最大可达3 ~ 5 mm，多数颗粒发育明显的网格状裂纹（图2（b-d））。基质由石英（约40%）、白云母（约20%）、黑云母（约15%）及少量斜长石等矿物组成。该矿物组合是进行相平衡模拟以约束变质温压条件的理想对象^{［18 - 19］}。

2.2 分析方法

2.2.1 矿物主量元素分析

矿物的原位主量元素分析在武汉上谱分析科技有限责任公司完成，采用日本电子JXA-8230型电子探针（EPMA）。分析条件为：加速电压15 kV，束流2 nA，束斑直径1 μm。采用ZAF方法进行数据校正，标样为SPI公司提供的系列矿物及单质标样。

2.2.2 全岩地球化学分析

样品的全岩主量元素分析在武汉上谱分析科技有限责任公司完成。全岩主量元素分析仪器使用日本理学（Rigaku）生产的ZSX PrimusⅡ型波长色散X射线荧光光谱仪（XRF），4.0 kW端窗铑靶X射线光管，测试条件为电压：50 kV，电流：60 mA，主量各元素分析谱线均为Kα，标准曲线使用国家标准物质岩石系列GBW07101-14、土壤系列：GSS07401-08、水系沉积物系列和GBW07302-12建立。数据校正采用理论α系数法，测试相对标准偏差（RSD）< 2%。

全岩微量元素含量在武汉上谱分析科技有限责任公司利用Agilent7700e ICP-MS分析完成，仪器工作参数为：RF功率为1550 W，等离子体气流速15 L/min，采样深度7.5 mm，检测模式为模拟和脉冲两种模式，采集时间100 ms，矩管：屏蔽矩。测试精度为5%（RSD），准确度为10%（相对偏差（RD）），内标采用In。用于ICP-MS分析的样品处理如下：（1）将200目样品置于105 ℃烘箱中烘干12小时；（2）准确称取粉末样品50 mg置于Teflon溶样弹中；（3）先后依次缓慢加入1 ml高纯HNO₃和1 ml高纯HF；（4）将Teflon溶样弹放入钢套，拧紧后置于190 ℃烘箱中加热24小时以上；（5）待溶样弹冷却，开盖后置于140 ℃电热板上蒸干，然后加入1 ml HNO₃并再次蒸干；（6）加入1 ml高纯HNO₃、1 ml MQ水和1ml 内标In（浓度为1 ppm），再次将Teflon溶样弹放入钢套，拧紧后置于190 ℃烘箱中加热12小时以上；（7）将溶液转入聚乙烯料瓶中，并用2% HNO₃稀释至100 g以备ICP-MS测试。

2.2.3 锆石U-Pb年代学与微量元素分析

锆石U-Pb定年及微量元素测试在中国地质大学（北京）矿物激光微区分析实验室（MILMA Lab）利用LA-ICP-MS方法完成。激光剥蚀系统为193 nm ArF准分子激光器，剥蚀斑束直径35 μm。采用国际标准锆石91500和GJ-1作为外部标准进行U-Pb分馏校正，NIST SRM 610玻璃标样用于微量元素含量校准。数据处理采用ICPMSDataCal软件。本次研究共对样品中的40颗锆石进行了分析。

2.2.4 相平衡模拟

相平衡模拟使用GeoPS软件（3.6.5.16878版本），采用Holland和 Powell^［20］的热力学数据库及相应的矿物活度-成分（a-X）模型。模拟所采用的具体固溶体模型包括：石榴子石（Gt；White^［21］）、黑云母（Bio；White等^［21］）、白云母（Ms；White等^［21］）、斜长石（Pl；Holland和Powell^［22］）、钾长石（Kf；Waldbaum&Thompson^［23］）以及绿泥石（Chl；Holland&Powell^［24］）等。

模拟基于样品的全岩成分各氧化物（扣除烧失量归一化至100%后对TFe₂O₃简化处理，对CaO含量依据全岩P₂O₅含量和磷灰石分子式进行重新计算以扣除部分磷灰石中CaO的含量）及独立组分H₂O进行，以约束其稳定的矿物组合及P-T条件。

变质峰期水以流体相或含水矿物结构水的形式存在，是热力学系统的必要组分，且在中-高压变质作用的进变质过程中，含水矿物（如白云母）的脱水反应会持续释放H₂O，且流体不易逃逸，系统中通常存在游离流体相。因此尽管全岩成分已扣除烧失量，但模拟中假设水为饱和状态（过剩纯H₂O流体相）。

样品全岩主量元素分析中总铁以TFe₂O₃形式给出，未分别测试FeO和Fe₂O₃。基于对样品的岩相学观察（图2（a-d））未见磁铁矿、赤铁矿等富Fe³⁺矿物，指示岩石形成于相对还原的环境，且即使存在少量的Fe³⁺，其主要赋存于副矿物（如磁铁矿）中，对硅酸盐矿物平衡的影响可忽略。因此，在相平衡模拟过程中，将TFe₂O₃简化处理，全部换算为FeO（换算系数为0.8998），即假设所有铁均为Fe²⁺。

3 分析结果

3.1 矿物化学成分

石榴子石的电子探针成分分析结果列于表1，其端元组分三角投图示于图3。分析结果表明，石榴子石具有明显的成分环带：从核部（Gt-c）到边部（Gt-r），铁铝榴石（Alm）组分由73.9 ~ 75.4%保持主导地位；钙铝榴石（Grs）组分从核部约5%增至边部9.7%，显著升高；而镁铝榴石（Prp）组分则从14.0 ~ 15.2%降至12.1 ~ 14.0%呈现轻微的降低趋势。锰铝榴石（Sps）组分含量极低（ < 1%）。

白云母的电子探针分析结果（数据详见表2）显示其成分相对均一。SiO₂含量为46.2 ~ 48.3%，Al₂O₃为34.3 ~ 36.1%，K₂O为9.2 ~ 10.1%，Na₂O为0.9 ~ 1.2%，TiO₂为0.3 ~ 0.9%，FeO为0.9~1.4%，MgO为0.8 ~ 1.1%。

3.2 全岩地球化学特征

样品LCBR21-05的全岩主、微量元素分析结果见表3和4，其稀土元素（REE）球粒陨石标准化图解和微量元素原始地幔标准化蛛网图示于图4。其主量元素特征为：SiO₂含量为66.91%，高Al₂O₃（15.52%），低CaO（1.06%）、MgO（2.26%），具有高Al₂O₃/TiO₂（20.29）比值和高K₂O/Na₂O（1.90）比值。值得注意的是，MnO含量仅为0.07%，属于泥质岩中Mn含量的低范围值，这一极低的Mn含量意味着Mn在体系中的热力学影响有限。球粒陨石标准化图解（图4（a））显示，样品∑REE较高，La（56.7 × 10^-6）、Ce（111 × 10^-6）等轻稀土元素（LREE）富集，Yb（3.05 × 10^-6）、Lu（0.48 × 10^-6）等重稀土元素（HREE）亏损，（La/Yb）_N比值为16.79，轻、重稀土分异强烈，呈典型右倾分布模式；Eu/Eu*为0.21，呈现显著的Eu负异常（Eu/Eu* < 1），与林彦蒿等^［25］在巴松措附近报道的片麻岩相比，该样品的右倾趋势更明显，Eu负异常强度更大。在其它微量元素特征方面，样品富集Rb（144 × 10^-6）、Ba（617 × 10^-6）等大离子亲石元素（LILE），相对亏损Nb（20.6 × 10^-6）、Ta（1.18 × 10^-6）等高场强元素（HFSE）。在原始地幔标准化图解（图4（b））中显示，与林彦蒿等^［25］报道的片麻岩相比，该片岩LIFE含量更高，Ti亏损更显著，HFSE含量相近；Nb/Ta比值为17.46，接近大陆地壳平均比值17.5，Zr/Hf比值为33.50符合大陆地壳比值范围^［26］，Th/Yb为7.44，Nb/Yb为6.76，与大陆碰撞构造背景的变质岩地球化学特征相符^［27］。

3.3 锆石U-Pb年代学结果

样品中的锆石多呈半自形-自形长柱状，长宽比多为2∶1至3∶1，粒径约100~150 μm。阴极发光（CL）图像显示，大部分锆石具有核-边结构（图5（a））。锆石核部通常发育振荡环带或结构复杂，其锆石U-Pb分析结果如表5所示，显示出约1060 Ma和约1029 Ma两组年龄（图5（b-c）），代表中元古代末期的继承锆石年龄。锆石CL图像边部弱分带/亮白色，Th/U比值普遍低于0.1，并且具有低REE含量，这些证据均指示其为变质成因^［29-33］。对边部变质增生区域的定年结果集中分布在（191.84±0.84） Ma，其谐和图与加权平均年龄高度一致（图5（d-e）），此年龄被解释为本研究样品所经历主期变质作用的时间。

3.4 相平衡模拟和变质作用P-T条件

基于样品全岩成分，通过相平衡模拟构建了其P-T视剖面图（图6（a））。视剖面图中投绘了根据实际矿物成分（白云母和石榴子石的XMg等值线）计算得到的温压计约束。模拟结果显示，薄片中观察到的矿物组合（石榴子石 + 白云母 + 黑云母 + 石英 + 斜长石）稳定在一个相对狭窄的P-T区间如图6（b）红色框中稳定域所示。综合矿物成分等值线交汇区域（图6（b）），可以进一步限定本次变质事件的峰期条件约为589.5 ℃和12.4 kbar。其中，石榴子石的6条XMg等值线由左到右的数值依次为：0.121，0.122，0.123，0.127，0.129和0.144；白云母的6条XMg等值线由下到上的数值依次为：0.51，0.52，0.53，0.59和0.62。模拟揭示顺时针P-T轨迹（图6（b）箭头所示），表现为在压力从约8 kbar显著增加至约12 kbar的过程中，温度仅小幅变化，指示近等温增压过程。

4 讨 论

4.1 原岩性质与构造背景

全岩地球化学分析结果表明，样品LCBR21-05具有高Al₂O₃/TiO₂（20.29）、高K₂O/Na₂O（1.90）比值及高Al₂O₃（15.52%）含量，这些特征与被动大陆边缘或碰撞造山带典型的泥质沉积岩地球化学组成高度一致^［34］。样品强烈富集LREE和LILE（如Rb、Ba），明显亏损HFSE（如Nb、Ta、Ti），且具有显著的Eu负异常（Eu/Eu* = 0.21）（图4）。这些特征共同指示其原岩为成熟度较高的陆源泥质沉积岩，其物源主要来自上地壳物质^［25］。显著的Eu负异常通常归因于源区存在斜长石残留或沉积/变质过程中发生了以斜长石分离为特征的熔体抽取作用^［35］，这与碰撞造山背景下地壳深熔作用的过程相符。此外，样品的Th/Yb（7.44）和Nb/Yb（6.76）比值落在典型大陆碰撞环境范围内^［27］。因此，综合地球化学特征表明，该石榴子石云母片岩的原岩形成于南、北拉萨地体碰撞相关的构造背景，其原岩记录了碰撞造山过程中地壳物质的沉积再循环与变质改造。

4.2 变质时间与P-T演化轨迹

锆石U-Pb定年结果显示，样品中变质成因锆石边部给出了约192 Ma的加权平均年龄（图5（d-e）），这为巴松措地区存在早侏罗世变质事件提供了直接证据。该年龄略晚于Weller等^［12］在南侧巴松措杂岩带报道的峰期变质时代（204 ~ 201 Ma），但与北侧扎拉杂岩带报道的变质时代（198 ~ 192 Ma）在误差范围内完全重叠。这一发现将巴松措杂岩带内与碰撞相关的热事件记录明确延续至约192 Ma，表明南、北拉萨地体碰撞后的构造热活动并非仅限于204 ~ 201 Ma，而是持续到了早侏罗世早期。

相平衡模拟成功限定了本次变质事件的峰期条件为约589.5 ℃和12.4 kbar，并揭示了一条具有顺时针演化特征的P-T轨迹（图6（b））。该轨迹表现为压力显著增加（从约8 kbar至12.4 kbar）的同时，温度仅小幅变化，指示了一个近等温的快速埋藏（增压）过程。石榴子石的成分环带为此提供了独立的矿物学佐证：从核部到边部，钙铝榴石（Grs）组分显著增加，而镁铝榴石（Prp）组分轻微降低（图3，表1）。这种成分变化正是进变质过程中压力迅速升高主导、升温相对滞后的典型响应，符合Barrovian型变质作用中与地壳增厚相关的石榴子石生长模式^{［36 - 38］}。P-T轨迹的顺时针形态及峰期高压特征，共同揭示了变质作用发生在碰撞导致的构造负载与地壳持续增厚的动力学背景之下^［39］。近等温增压（约8→12.4 kbar）与石榴子石边部Grs升高（约5%→9.7%）耦合，反映构造埋藏速率 > 热松弛速率。同时，值得讨论的是在本次相平衡模拟过程中，采用Holland和Powell^［20］的热力学数据库及相应的矿物活度-成分（a-X）模型未专门考虑Mn的影响。本研究样品中锰铝榴石（Sps）组分含量极低（< 1%），且全岩MnO含量仅为0.07%，属于White等^［21］定义的“大多数泥质岩中实测含量”的低范围（0.1 ~ 0.3%）。根据White等^［21］对含Mn泥质岩的系统研究，在本文的Mn含量水平下，Mn的主要作用是扩大石榴子石的稳定域，但对其他常见泥质矿物（如白云母、黑云母）的稳定性影响很小，含Mn体系计算出的相图与不含Mn体系高度相似；即使考虑Mn，石榴子石在大部分P-T范围内的模式含量仍小于1%，对于相平衡计算的整体影响可以忽略。因此，本研究样品极低的Mn含量对于最终P-T条件的约束不产生实质性影响。

值得注意的是，本研究所获峰期压力（12.4 kbar）高于Zhang等^［6］在该区域报道的早期（约204 ~ 201 Ma）Barrovian型变质作用的峰期压力（约9 kbar），而温度（589.5 ℃）则相对更低。这一“更高压、更低温”的特征，可能指示在约192 Ma时期，研究区地壳经历了进一步的垂向加厚（埋深增大），但地温梯度尚未因深部热松弛而同步升高，从而记录了碰撞造山过程持续进行、热松弛相对滞后的瞬态热状态，即本研究的192 Ma事件压力更高（12.4 kbar）但温度相当（约590 ℃），指示同期地壳垂向增厚而非热抬升。

综上，锆石年代学（约192 Ma）与相平衡模拟（顺时针）P-T轨迹共同约束了碰撞后持续增厚的时间-温压框架。

4.3 区域构造意义

本研究在巴松措杂岩带内识别出的约192 Ma高压变质事件，对厘清南、北拉萨地体碰撞后的构造演化具有重要启示。关于巴松措地区两条早中生代变质带的关系，前人存在不同认识：一种观点认为它们分别代表了碰撞后不同构造层次、不同时期的响应，即巴松措杂岩带记录较早的M1事件（204 ~ 201 Ma），扎拉杂岩带记录较晚的M2事件（198 ~ 192 Ma）^［12］；另一种观点则认为二者可能属于同一持续变质热事件的产物^［40］。本研究获得的证据支持后一种观点。首先，192 Ma的变质年龄不仅存在于扎拉杂岩带，也明确存在于巴松措杂岩带内部，表明该时期的热事件具有区域性。其次，该期变质作用表现出典型的与挤压增厚相关的Barrovian型P-T轨迹特征，与早期（204 ~ 201 Ma）事件的动力学背景一致，均为碰撞造山过程的延续。结合陈智佳^［40］在扎拉杂岩带识别出的201 ~ 190 Ma变质年龄，可以认为巴松措地区在约204 ~ 190 Ma期间普遍经历了一期持续的、与南、北拉萨地体碰撞后地壳持续增厚相关的区域变质作用^［41］。

因此，本研究揭示的约192 Ma高压变质事件表明，南、北拉萨地体的碰撞造山作用并非在初始拼贴（204 ~ 201 Ma）后即刻停止，其导致的构造缩短与地壳增厚过程至少持续了十余个百万年，直至早侏罗世早期。这一持续增厚过程是拉萨地体早期构造演化乃至青藏高原原特提斯构造体制终结阶段的重要组成部分，本研究为其提供了关键的温压-年代学约束。

5 结 论

（1）变质时代：锆石U-Pb定年结果显示，样品中变质成因锆石的年龄为约192 Ma，精确厘定了巴松措杂岩带内存在一期早侏罗世变质事件。该年龄将南、北拉萨地体碰撞后该区域与碰撞相关的构造热事件记录明确延续至192 Ma。

（2）变质条件与过程：相平衡模拟定量限定本次变质作用的峰期温压条件为 589.5 ℃/12.4 kbar。模拟揭示的顺时针P-T轨迹，结合石榴子石成分环带特征，共同指示了一个以近等温增压为特征的快速构造埋藏过程，具有典型的Barrovian型变质作用属性，反映了碰撞造山过程中地壳增厚与热松弛滞后的耦合动力学响应。

（3）构造意义：本研究峰期变质压力（12.4 kbar）高于已报道的约204 ~ 201 Ma的同类变质作用压力（约9 kbar），而温度相对更低。指示南、北拉萨地体碰撞后地壳持续垂向增厚（而非热松弛抬升）。支持巴松措地区在204 ~ 190 Ma期间经历持续地壳增厚变质事件，表明碰撞造山作用持续 > 10 Ma，为完善拉萨地体构造演化模型提供了新的约束。

参考文献

[1] ALLÈGRE C J. Cold Gondwana， warm Tethys and the Tibetan Lhasa Block［J］. Nature， 1984， 310： 165 - 166.

[2] AUDLEY-CHARLES M G. Reconstruction of eastern Gondwanaland［J］. Nature， 1983， 306： 48 - 50.

[3] SENGOR A C. Tectonics of the tethysides： Orogenic collage development in a collisional setting［J］. Annual Review of Earth and Planetary Sciences， 1987， 15： 213 - 244.

[4] KAPP P， AN Y， MANNING C E， et al. Blueschist-bearing metamorphic core complexes in the Qiangtang block reveal deep crustal structure of northern Tibet［J］. Geology， 2000， 28（1）： 19 - 22.

[5] ZHU D C， ZHAO Z D， NIU Y L， et al. The Lhasa Terrane： Record of a microcontinent and its histories of drift and growth［J］. Earth and Planetary Science Letters， 2011， 301（1/2）： 241 - 255.

[6] ZHANG Z M， DONG X， SANTOSH M， et al. Metamorphism and tectonic evolution of the Lhasa terrane， Central Tibet［J］. Gondwana Research， 2014， 25（1）： 170 - 189.

[7] 刘德民， 王杰， 姜淮， 等. 青藏高原形成演化动力机制及其远程效应［J］. 地学前缘， 2024， 31（1）： 154 - 169.

[8] 李子颖， 郭强， 秦明宽， 等. 青藏高原碰撞隆升与中国西北及中亚邻区砂岩铀成矿作用［J］. 矿床地质， 2024， 43（6）： 1336 - 1359.

[9] 余长胜. 拉萨地体的起源及演化——来自前寒武纪基底的制约［D］. 长春： 吉林大学， 2025.

[10] YANG J S， XU Z Q， LI Z L， et al. Discovery of an eclogite belt in the Lhasa block， Tibet： A new border for Paleo-Tethys？［J］. Journal of Asian Earth Sciences， 2009， 34（1）： 76 - 89.

[11] 王向伟， 张保涛， 杨浩强， 等. 青海省海晏县团宝山一带变质岩年代学、地球化学及其地质意义［J］. 现代地质， 2023， 37（3）： 586 - 598.

[12] WELLER O M， STONGE M R， SEARLE M P， et al. Quantifying the P-T-t conditions of north-south Lhasa terrane accretion： New insight into the pre-Himalayan architecture of the Tibetan Plateau［J］. Journal of Metamorphic Geology， 2015， 33（1）： 91 - 113.

[13] 毕先梅， 莫宣学， 刘艳宾. 西藏南部极低级变质岩及其地质与资源意义［J］. 地学前缘， 2024， 31（1）： 201 - 210.

[14] YIN A， HARRISON T M. Geologic evolution of the Himalayan-Tibetan Orogen［J］. Annual Review of Earth and Planetary Sciences， 2000， 28： 211 - 280.

[15] LIN Y H， ZHANG Z M， DONG X， et al. Early Mesozoic metamorphism and tectonic significance of the eastern segment of the Lhasa terrane， south Tibet［J］. Journal of Asian Earth Sciences， 2013， 78： 160 - 183.

[16] GUYNN J， KAPP P， GEHRELS G E， et al. U-Pb geochronology of basement rocks in central Tibet and paleogeographic implications［J］. Journal of Asian Earth Sciences， 2012， 43（1）： 23 - 50.

[17] SHI R D， YANG J S， XU Z Q， et al. The Bangong Lake ophiolite （NW Tibet） and its bearing on the tectonic evolution of the Bangong-Nujiang suture zone［J］. Journal of Asian Earth Sciences， 2008， 32（5/6）： 438 - 457.

[18] 李方兰， 刘学龙， 周云满， 等. 滇西保山地块陡崖铁铜多金属矿床石榴子石年代学及其地球化学特征［J］. 地学前缘， 2024， 31（3）： 113 - 132.

[19] 王晓彤， 刘军， 杨艳， 等. 新疆北部喇嘛苏铜矿床成矿时代与构造背景： 来自石榴子石原位LA-ICP-MS U-Pb测年的证据［J］. 现代地质， 2023， 37（3）： 645 - 661.

[20] HOLLAND T J B， POWELL R. An improved and extended internally consistent thermodynamic dataset for phases of petrological interest， involving a new equation of state for solids［J］. Journal of Metamorphic Geology， 2011， 29（3）： 333 - 383.

[21] WHITE R W， POWELL R， JOHNSON T E. The effect of Mn on mineral stability in metapelites revisited： New a-x relations for manganese-bearing minerals［J］. Journal of Metamorphic Geology， 2014， 32（8）： 809 - 828.

[22] HOLLAND T， POWELL R. Activity-composition relations for phases in petrological calculations： An asymmetric multicomponent formulation［J］. Contributions to Mineralogy and Petrology， 2003， 145（4）： 492 - 501.

[23] WALDBAUM D R， THOMPSON J B Jr. Mixing properties of sanidine crystalline solutions： Ⅱ. Calculations based on volume data1［J］. American Mineralogist， 1968， 53（11/12）： 2000 - 2017.

[24] HOLLAND T J B， POWELL R. An internally consistent thermodynamic data set for phases of petrological interest［J］. Journal of Metamorphic Geology， 1998， 16（3）： 309 - 343.

[25] 林彦蒿， 张泽明， 董昕. 拉萨地体东南缘始新世早期变质作用及其构造意义［J］. 岩石学报， 2013， 29（6）： 1962 - 1976.

[26] MCDONOUGH W F， SUN S S. The composition of the earth［J］. Chemical Geology， 1995， 120（3/4）： 223 - 253.

[27] PEARCE J A. Geochemical fingerprinting of oceanic basalts with applications to ophiolite classification and the search for Archean oceanic crust［J］. Lithos， 2008， 100（1/2/3/4）： 14 - 48.

[28] SUN S S， MCDONOUGH W F. Chemical and isotopic systematics of oceanic basalts： Implications for mantle composition and processes［J］. Geological Society， London， Special Publications， 1989， 42（1）： 313 - 345.

[29] 杨伟， 邱亮， 木红旭， 等. 西准噶尔谢米斯台辉绿岩锆石U-Pb年代学及其构造意义［J］. 现代地质， 2025， 39（6）： 1461 - 1482.

[30] 张乐民， 张聚全， 吴晓贺， 等. 太行山南段洪山碱性火山-侵入杂岩成因与构造意义： 锆石U-Pb年代学和地球化学的约束［J］. 现代地质， 2025， 39（3）： 752 - 773.

[31] 胡子奇， 张德贤， 刘磊. 束斑直径和能量密度对锆石U-Pb定年准确度的影响研究［J］. 现代地质， 2023， 37（3）： 722 - 732.

[32] 黄宇， 钟世华， 李三忠， 等. 副矿物包裹体和信号采集时间对锆石U-Pb年龄和微量元素分析结果的影响［J］. 地学前缘， 2025， 32（1）： 388 - 400.

[33] 刘晓惠， 刘一珉， 丁林， 等. 中拉萨地体当惹雍措锆石微量元素特征及其对地壳厚度演化的指示［J］. 地学前缘， 2025， 32（1）： 343 - 366.

[34] TAYLOR S， MCLENNAN S. The Continental Crust： Its Composition and Evolution： An Examination of the Geochemical Record Preserved in Sedimentary Rocks［M］. Oxford： Blackwell Scientific Publications， 1985.

[35] RUDNICK R L， GAO S. 3.01-Composition of the Continental Crust［M/OL］ // Holland H D， Turekian K K， eds. Treatise on Geochemistry. Oxford： Pergamon， 2003， 3： 1-64 ［2025-05-18］.

[36] SPEAR F. Metamorphic phase equilibria and pressure-temperature-time paths［M］. Washington， D.C.： Mineralogical Society of America， 1993.

[37] 夏琼霞， 郑永飞. 高压-超高压变质岩石中石榴石的环带和成因［J］. 岩石学报， 2011， 27（2）： 433 - 450.

[38] 杨建军. 榴辉岩石榴子石反应边的震击成因初探［J］. 地学前缘， 2025， 32（6）： 276 - 285.

[39] ENGLAND P C， THOMPSON A B. Pressure： Temperature： Time paths of regional metamorphism I. heat transfer during the evolution of regions of thickened continental crust［J］. Journal of Petrology， 1984， 25（4）： 894 - 928.

[40] 陈智佳. 青藏高原拉萨地体东部早中生代变质作用与构造意义［D］. 北京： 中国地质大学（北京）， 2017.

[41] 张艳斌， 翟明国， 周艳艳， 等. 大陆下地壳［J］. 地学前缘， 2024， 31（1）： 28 - 45.

基金资助

中央高校基本科研业务费专项资金项目(2652023001)