7月28日，中国科学院青藏高原研究所丁林院士带领的大陆碰撞与高原隆升团队受《自然综述：地球与环境》(Nature Reviews Earth & Environment)主编邀请，以“青藏高原隆升时间和机制 (Timing and Mechanisms of Tibetan Plateau uplift)”为题，发表综述文章，系统阐述了青藏高原的差异性隆升过程和深部动力学机制。
大陆碰撞-俯冲等深部圈层作用驱动的青藏高原隆升是新生代全球最重要的地质事件之一。高原隆升显著影响了地表圈层—大气圈、水圈/冰冻圈、生物圈和人类圈的耦合作用过程，深刻地影响了亚洲气候动力学、生物多样性、碳循环、现代水资源和大江大河演化（图1），是21世纪地球系统科学研究的前沿阵地。
然而，在大陆碰撞过程中，青藏高原大陆岩石圈变形和地表高程时空变化的机制仍不清楚。近些年，随着青藏高原定量古高度数据的加速产生，学者们逐渐认识到高原具有差异性隆升的特征，部分地区的隆升时间比以前的推测或早或晚，已有的动力学模式均不能完整体现高原隆升过程（图1）。
图1. 青藏高原地貌特征、西风季风系统和主要构造简图
白垩纪构造事件和青藏高原山脉雏形
“青藏高原的完整演化模式必须考虑亚洲在印度-欧亚大陆碰撞之前的构造事件中继承下来的古地貌和岩石圈的不均一性，这对认识高原差异性隆升至关重要。”丁林说。
通过详细分析青藏高原白垩纪海陆转换、构造变形、岩浆和低温热年代学证据，研究团队提出，拉萨-羌塘地体的碰撞以及随后的拉萨岩石圈向北俯冲导致了分水岭山脉的初步生长；南部新特提斯洋的持续俯冲，在约9500万年前将冈底斯地区隆升至海平面之上，从而形成与现今的安第斯山相似的发展过程，称之为安第斯型冈底斯山，并在藏南地区形成显著的降水效应。此时的青藏高原仅有两条狭窄的山脉，即分水岭山脉和冈底斯山脉，但地表隆起的幅度仍有待量化。
印度-欧亚大陆初始碰撞时间和模式
印度-欧亚板块碰撞时间和方式对于限定印度北缘范围和新生代陆内缩短变形量至关重要，而它们又是约束高原地表隆升幅度和深部动力学机制的关键。
目前，关于新特提斯洋闭合历史的假说包括：大印度洋盆模型（图2a）、洋内俯冲模型（图2c-d）和单阶段俯冲碰撞模型（图2d）。这些假说对大印度（已俯冲消失到青藏高原下的印度部分）的规模以及印度-欧亚大陆初始碰撞时间做出了截然不同的预测。然而，以上模型都是基于丁林院士团队等发现的关键证据—即印度-欧亚大陆碰撞形成的前陆盆地，该盆地在6500万年至5900万年之间开始接受来自冈底斯岛弧区的物源，表明此时印度-欧亚大陆已经开始开始碰撞，大大早于此前国际公认的印度与亚洲大陆在5000万年才初始碰撞的时间。
同时，团队近20年在西藏、巴基斯坦和印度工作还获得了关键认识：（1）目前发育在藏南的冈底斯岩浆弧、弧前盆地、蛇绿岩和海沟形成于拉萨地体南缘的洋-陆俯冲系统，而非形成于远离大陆的洋内俯冲系统；（2）巴基斯坦北部的证据表明约5200万年亚洲物质可以到达印度次大陆，不支持印度-欧亚大陆碰撞前存在洋盆。因此，该文指出，单阶段俯冲碰撞模型是解释印度-欧亚大陆碰撞最简单，也是得到地质证据支持的模型。
令人遗憾的是，古地磁数据指出，如果印度-亚洲大陆在约6000万年发生碰撞，印度和亚洲大陆岩石圈必须在完全的大陆环境中吸收约4000公里缩短量。然而目前地质证据表明亚洲（<1000公里）和喜马拉雅（<1000公里）地壳缩短量不到2000公里。为了协调地壳缩短和古地磁汇聚量之间的不匹配，国际上许多学者提出了大印度洋盆模型和洋内俯冲模型（图2a-c），但是上述模型目前都缺乏地质证据的支持。另外，如果印度大陆在向北漂移过程中发生过约90度°的逆时针旋转，也可简单的吸收2000公里的纬向缩短量。
图2. 印度-欧亚大陆初始碰撞模型. a. 大印度洋盆模型; b. 洋内俯冲模型; c. 弧后盆地模型；d. 单阶段俯冲碰撞模型
青藏高原新生代差异隆升及其动力学机制
解决高原隆升历史的需求极大地促进了古高度计的发展，广泛使用的古高度定量重建技术包括氢/氧同位素、动植物化石、团簇同位素等。这些古高度定量重建技术为大陆变形和高原生长提供了关键信息，可以更加清晰地认识高原差异隆升过程和动力学机制。
结合已有定量古高度结果和深部动力学证据，研究团队恢复了青藏高原自约6000万年以来不同地体地表隆升历史（图3）和岩石圈演化过程（图4），提出青藏高原不同造山带具有差异的隆升历史。5500万年至45 00万年前，由于新特提斯洋俯冲板片的断离，冈底斯造山带隆升到4500米高海拔；4500万年至4000万年，新特提斯板块断离之后，更具浮力的印度岩石圈向北水平楔入，激活羌塘地体南北部缝合带发生陆内俯冲，使得分水岭山脉隆升到5000米的高海拔；此时位于冈底斯造山带和分水岭造山带之间的中央谷地、高原最南部的喜马拉雅造山带以及高原北部还处于小于2000米的低海拔，高原整体形成“两山夹一盆”的地貌特征。4000万年至3000万年，拉萨岩石圈在中央谷地下方拆沉，上地壳缩短、岩浆底垫和软流圈上涌等多种深部地球动力学过程耦合作用，使得中央谷地抬升4500米的目前高度。标志着青藏高原由造山带正式转变为统一高原。250万年至1500万年，由于印度大陆的持续俯冲，喜马拉雅山脉下方俯冲的印度大陆岩石圈及藏北可可西里-昆仑山下方俯冲的欧亚大陆岩石圈先后发生拆沉，喜马拉雅山与昆仑山先后隆升到现代高度，现代意义上的高原形成（图3 & 4）。但是北部地区的隆升历史还存在较大的不确定性，需要更多的定量古高度数据来约束。
图3. 青藏高原不同地体差异隆升历史
图4. 青藏高原不同地体岩石圈演化过程、地表隆升历史和环境演变
现代印度大陆的俯冲
地球物理探测揭示了现今印度和欧亚大陆岩石圈发生了从水平楔入到陡峭俯冲、板块撕裂和断离以及拆沉等各种地球动力学行为（图5），从而指示在整个新生代印度-亚洲大陆碰撞过程中，类似的过程不断发生，最终导致了青藏高原构造变形、岩浆作用和地表隆升的时空差异。
图5. 现今印度-欧亚大陆岩石圈俯冲几何形态
针对青藏高原隆升的时间和机制问题，该文指出，今后需要重点发展的方向包括：（1）解决印度-欧亚大陆汇聚量和地壳缩短之间的不一致性；（2）需要大量高分辨率古高度数据精确限制高原隆升历史；（3）结合数值模拟和地质数据对高原地球系统演化历史进行准确重建；（4）结合地球物理成像技术和地球动力学模拟，阐明大陆岩石圈的循环过程和分布范围，分析探索大陆碰撞如何影响邻近板块边界的构造以及全球规模的地幔对流。
该研究成果获得第二次青藏高原综合科学考察研究（2019QZKK0708）、基金委青藏高原地球系统基础科学中心（41988101）、中国科学院“丝路环境”先导专项（XDA20070301）等项目的联合资助。
