印度-亚洲大陆碰撞是地质历史时期最重大的构造事件之一，形成“世界屋脊”青藏高原和雄伟的喜马拉雅山脉(图1)，对亚洲乃至全球的海陆分布格局、气候环境和生物多样性都产生了深远的影响。青藏高原作为研究陆-陆碰撞最为理想的天然实验室，对印度-亚洲大陆初始碰撞时间和碰撞动力学过程的研究长期以来一直是固体地球科学领域的前沿和热点。印度-亚洲大陆初始碰撞时间是理解喜马拉雅造山作用的起点，而印度-亚洲大陆碰撞动力学过程则是重建大陆碰撞造山带和青藏高原演化过程的关键。近年来自沉积物源方面的证据表明，约60 Ma亚洲碎屑物质首次沉积到印度北缘被动大陆边缘(DeCelles等, 2014; Wu等, 2014; Hu等, 2015)，标志着印度-亚洲大陆初始碰撞，该初始碰撞时间被国内外广泛接受(丁林等, 2017; 胡修棉等, 2017; Kapp and DeCelles, 2019; van Hinsbergen等, 2019; Parsons等, 2020)。
然而，对印度-亚洲大陆碰撞动力学过程的理解还存在诸多不同看法，学界提出了多个模型，可归结为单阶段碰撞模型和双阶段碰撞模型两大类。其中，单阶段碰撞模型以大印度大陆模型(Greater India model)为代表，又被称为continental Greater India model (Westerweel等, 2019)。双阶段碰撞模型包括大印度洋盆假说(Greater India Basin hypothesis) (van Hinsbergen等, 2012)和弧陆碰撞模型(the island arc-continent collision model) (Aitchison等, 2007)。最近，Kapp and DeCelles (2019)提出了新版的双阶段弧陆碰撞模型(the India-arc collision with the Xigaze backarc ocean basin)，这里称之为日喀则弧后洋盆模型。这些模型各有优点，但都只能或多或少地解释目前已观察到的部分地质事实。
为准确限定印度-亚洲大陆碰撞的时间、位置和动力学过程，在国家自然科学基金和中科院国际合作伙伴计划项目的资助下，中国科学院地质与地球物理研究所博士研究生袁杰、邓成龙研究员、朱日祥院士、郭正堂院士、贺怀宇研究员、李仕虎博士、沈中山博士研究生、秦华峰副研究员，以及首都师范大学杨振宇教授，荷兰乌特勒支大学Wout Krijgsman教授，南京大学胡修棉教授、安慰博士，中国科学院青藏高原研究所丁林院士共同合作，对藏南特提斯喜马拉雅地块江孜地区的上白垩统床得组大洋红层和萨嘎地区的古新统桑单林组深水红色硅质页岩(图1)开展构造古地磁学、岩石磁学、磁性地层学和岩相学(图2，图3)等综合研究，获得了可靠的古地磁数据，揭示了特提斯喜马拉雅地块在白垩纪晚期-古近纪初期快速向北漂移的特征，提出了“北印度海”(North India Sea)假说，并在此基础上构建了两阶段的印度-亚洲大陆碰撞动力学新模型(图4，图5)。
图1 研究区地质与地形。(a)青藏高原及其周边地形图；(b)印度与亚洲大陆碰撞带地质简图，修改自Yin(2006)；(c)江孜地区地质简图，修改自Chen等(2006)；(d)萨嘎地区地质简图
首先，本研究从江孜地区床得组白垩纪大洋红层中获得了通过倒转检验、并由碎屑赤铁矿携带的原生剩磁，应用E/I法(Tauxe and Kent, 2004)和等温剩磁(IRM)各向异性法(Hodych and Buchan, 1994)进行了磁倾角校正，从而将特提斯喜马拉雅地块中东部约75 Ma的古地磁极确定为40.8°N/256.3°E，A95 = 1.8°，古纬度为19.4° ± 1.8°S，而此时亚洲大陆南缘的拉萨地块位于约13.7°N。因此，约75 Ma时，特提斯喜马拉雅地块和拉萨地块仍然被纬向宽度为3600 km左右的新特提斯洋所分隔(图4，图5)。其次，从萨嘎地区桑单林组深水红色硅质页岩中获得了通过倒转检验和褶皱检验、并由碎屑赤铁矿携带的原生剩磁，应用IRM各向异性法(Hodych and Buchan, 1994)进行了磁倾角校正，从而获得特提斯喜马拉雅地块中东部约61 Ma的古地磁极为74.0°N/278.5°E，A95 = 2.5°，古纬度为13.7° ± 2.5°N。
图2 代表性样品岩相学与岩石磁学分析。CL043-1、SE12-1和MB4-8分别来自才狼把剖面、桑单林剖面和木巴拉剖面。(a-c)扫描电子显微镜图片；(d-f)磁滞回线；(g-i)矫顽力谱分析
图3 古地磁结果。(a-i)代表性样品系统热退磁正交投影图、等面积投影图及剩磁衰减图；(j-l)高温剩磁分量(HTC)分别在地理坐标和地层坐标下的等面积投影图
图4 印度-亚洲大陆碰撞动力学演化，分别示意(a)约75 Ma、(b)约61 Ma、(c)约53 Ma和(d)约48 Ma的古地理模式。BNSZ为班公湖-怒江缝合带，WZFZ为Wallaby-Zenith断裂带。黑色星号表示计算大印度延伸尺寸的参考点(29.3°N, 85.3°E)
上述高质量的古地磁结果准确限定了特提斯喜马拉雅地块的位置，即在约75 Ma位于19.4° ± 1.8°S，而在约61 Ma已经向北快速移动到13.7° ± 2.5°N。因此，75～61 Ma期间，特提斯喜马拉雅地块的北向漂移速率为260.1 mm/year。而通过印度大陆的视极移曲线等方法计算得到印度克拉通在80～60 Ma期间的北向漂移速率仅为99.6 mm/year。特提斯喜马拉雅地块与印度克拉通存在明显的漂移速率差异，意味着在约75 Ma后特提斯喜马拉雅地块裂解并离开印度大陆，在两者之间形成一个拉分盆地，本研究称之为“北印度海”(North India Sea) (图4，图5)。本研究进一步探讨了特提斯喜马拉雅地块的裂解机制，认为Reunion地幔柱上涌对印度被动大陆边缘岩石圈减薄和新特提斯洋壳岩石圈长时间俯冲产生的拖拽力等过程的共同作用下，导致特提斯喜马拉雅地块裂解并张开形成“北印度海”。特提斯喜马拉雅地块快速向北漂移，在约61 Ma与亚洲大陆南缘的拉萨地块发生初始碰撞，碰撞位置位于约14°N。随后，印度与特提斯喜马拉雅地块在约53～48 Ma发生第二阶段碰撞，该碰撞自西向东穿时性发生，导致“北印度海”自西向东逐渐关闭(图4)。
沿低喜马拉雅带在巴基斯坦、印度、尼泊尔等地和喜马拉雅东部地区广泛发育了上白垩统-始新统或古新统-始新统浅海相地层，本研究认为是“北印度海”存在的证据。根据区域大地构造分析，本研究推测，主中央逆冲断裂带很可能是第二阶段洋壳俯冲的古地理位置。
图5 晚白垩世以来印度-亚洲大陆碰撞过程的阶段划分(a)及其与印度-亚洲大陆汇聚速率变化历史(b)的关联。图b的数据来自Cande and Stegman(2011)
本研究提出的“北印度海”假说及两阶段的印度-亚洲大陆碰撞动力学新模型与印度-亚洲的汇聚速率变化历史可以完美匹配(图5)。例如，从70 Ma到63 Ma，汇聚速率从约80 mm/year快速增加到约180 mm/year，本研究认为与Reunion地幔柱的推动作用有关；在63－61 Ma期间汇聚速率从约180 mm/year急剧下降到约110 mm/year，特提斯喜马拉雅地块和拉萨地块在61Ma发生的碰撞可很好地解释该汇聚速率变化；在61－53 Ma期间汇聚速率由约110 mm/year缓慢上升到约130 mm/year，对应于“北印度海”的逐步收缩；在53～48 Ma期间汇聚速率显著降低，则与特提斯喜马拉雅地块和印度在此期间发生穿时性碰撞导致北印度海自西向东逐步关闭的过程相符。
上述假说和模型还协调了一系列地质证据。本研究的古地磁结果表明，特提斯喜马拉雅地块与拉萨地块在61 Ma发生第一阶段的大陆碰撞，碰撞位置约为14°N，该碰撞最为直接的地质证据就是亚洲碎屑物质在60 Ma左右首次到达印度北缘被动大陆边缘(DeCelles等, 2014; Wu等, 2014; Hu等, 2015)。
同时，对于第二阶段的印度-亚洲大陆碰撞，即53～48 Ma期间特提斯喜马拉雅地块和印度之间的穿时性碰撞，也得到大量地质证据的支持。例如：(1)亚洲来源的碎屑锆石在始新世早期(约55 Ma)已经沉积到印度大陆小喜马拉雅(sub-Himalaya)巴基斯坦Balakot地区(Ding等, 2016)，而低喜马拉雅带(Lesser Himalaya)中部尼泊尔地区获得特提斯喜马拉雅冲断带碎屑的年代不晚于45 Ma (DeCelles等, 2004; Najman等, 2005)，可支持本次研究提出的53～48 Ma期间发生的第二阶段穿时性大陆碰撞。(2)近年来的研究显示，特提斯喜马拉雅海相地层自西向东穿时性消失，在西部Zanskar地区消失时间为52～50 Ma，在中-东部的定日-岗巴地区为43～41 Ma，在东部的堆纳地区为35 Ma (Rowley, 1996; 胡修棉等, 2017)，这些地质现象也支持53～48 Ma期间发生的第二阶段穿时性大陆碰撞。(3)藏东贡觉盆地在52～48 Ma发生约30°的顺时针旋转(Li等, 2020)，同时贡觉盆地和可可西里盆地共同记录了54～52 Ma期间急剧增加的沉积速率变化特征(Li等, 2020; Jin等, 2018)。(4)孟加拉盆地在始新世中期来自喜马拉雅的碎屑也急剧增加(Alam等, 2003)，也是对53～48 Ma期间发生的第二阶段大陆碰撞的响应。(5)特提斯喜马拉雅带、高喜马拉雅带和低喜马拉雅带都记录了各种独立的证据，例如，特提斯喜马拉雅带在始新世-渐新世记录了地壳缩短、快速折返和岩浆作用，高喜马拉雅带记录了同期构造埋藏作用、高级变质作用和深熔作用以及印度克拉通发育同期前陆盆地(DeCelles等, 2014; 丁林等, 2017)。以上地质证据都表明特提斯喜马拉雅地块和印度大陆在48 Ma已经完成碰撞。
本研究提出的“北印度海”假说及两阶段的印度-亚洲大陆碰撞动力学新模型对印度-亚洲大陆碰撞过程、青藏高原隆升和变形、亚洲古地理和生物多样性格局演化都提供了重要制约，还为喜马拉雅造山作用与全球气候变化关联的构造-气候关系研究及古气候数值模拟提供了关键的边界条件。
研究成果发表于National Science Review。（Yuan J, Yang Z Y*, Deng C L*, Krijgsman W, Hu X M, Li S H, Shen Z S, Qin H F, An W, He H Y, Ding L, Guo Z T, Zhu R X. Rapid drift of the Tethyan Himalaya terrane before two-stage India-Asia collision[J]. National Science Review, 2020: nwaa173. DOI: 10.1093/nsr/nwaa173）（原文链接）
