摘要：地壳缩短导致青藏隆升造山是普遍的认识.然而,在青藏东部,越来越多的观测数据和研究支持了中下地壳流与隆升造山的关系.目前,地壳缩短造山机制和中下地壳流造山机制仍然处于争论之中.本文建立了二维黏弹塑性有限元模型,模拟了龙门山断层带的多个地震循环的应变与变形,探讨了无与有中下地壳流情况下,地壳地表的位移、速度与变形的分布和演化;以及有中下地壳流情况下,不同流动范围、速度与黏度对模型结果的影响;并结合地形变观测数据的约束,推测了青藏东缘中下地壳流的流动状态.模拟结果显示,通过对比有和无中下地壳流的模拟结果,发现青藏东部震间的地表垂向速度在变形样式及数值上存在较大差异,即存在地壳流的地表垂向抬升速率显著大于无地壳流;震间在龙门山断层西侧附近产生的垂向凸状隆起随中下地壳流的速度、黏度及通道长度的变化而变化.此外,本文研究结果对青藏其他地区可能存在的地壳流的研究也具有一定的参考意义.

关键词：青藏高原东缘 地壳缩短 中下地壳流 黏弹塑性有限元模拟 

单位：中国科学院计算地球动力学重点实验室; 中国科学院大学地球与行星科学学院; 北京100049