地球物理学报  2020, Vol. 63 Issue (10): 3727-3739   PDF    
造山带下岩石圈地幔小尺度减薄的动力学成因
雷天1, 李忠海1, 刘勉2,1     
1. 中国科学院大学地球与行星科学学院, 计算地球动力学重点实验室, 北京 100049;
2. 密苏里大学哥伦比亚分校, 地质学系, MO 65211, 美国
摘要:地质与地球物理观测数据表明青藏高原、安第斯山、以及帕米尔等典型造山高原之下均有明显的岩石圈地幔小尺度/分段式减薄现象.这些小尺度岩石圈减薄难以用经典的拆沉或对流减薄理论来解释,一方面,拆沉预示大尺度岩石圈地幔的剥离过程,而对流减薄则在黏度相对低的地幔岩石圈中发生,其主要以小尺度的局部增厚触发并仅减薄地幔岩石圈的底部区域.另一方面,拆沉或对流减薄模型都预测造山带尺度的地幔岩石圈拆离,都假设造山带岩石圈横向均一,然而实际的造山带岩石圈往往由多个不同的地块构成,块体之间岩性、物性、流变结构可能大有差别,即横向不均一性.这些造山带岩石圈地幔的横向不均一性,能否有效解释观测到的局部小尺度减薄现象?为此,我们构建了一系列高精度动力学数值模型,系统模拟了碰撞造山过程中岩石圈地幔的形变和不稳定性.结果表明,在塑性屈服强度很低的情况下,横向不均一的造山带岩石圈有发生分段式/小尺度减薄的可能性;其主要机理是由位错蠕变与强塑性作用所导致的应变集中使得地块间及壳幔间耦合弱化,从而使得较弱地块的岩石圈地幔在增厚时由于重力不稳定性而产生局部剥离,进而诱发小尺度软流圈上涌.模拟结果可以良好地解释发生在青藏高原东北缘、安第斯中部高原、以及帕米尔高原之下岩石圈的局部小尺度/分段式减薄现象.
关键词: 小尺度岩石圈减薄      造山带      横向不均一性      位错蠕变      塑性屈服      数值模拟     
Dynamics of small-scale removing of mantle lithosphere under orogens
LEI Tian1, LI ZhongHai1, LIU Mian2,1     
1. Key Laboratory of Computational Geodynamics, College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China;
2. Department of Geological Sciences, University of Missouri, Columbia, MO 65211, USA
Abstract: Geological and geophysical data suggests that small-scale fragmental lithospheric thinning, some all the way to the base of the crust, has occurred under some large orogens, including Tibetan Plateau, central Andes, and Pamir. Such localized small-scale removal of mantle lithosphere under orogens cannot be readily explained by the existing models of delamination and convective thinning (dripping)-the formal would predict large-scale peeling of the mantle lithosphere, and latter could occur in small scale (drips) when the viscosity is low, but in that case it tends to occur near the base of the mantle lithosphere rather than thinning the lithosphere all the way to the Moho. On the other hand, the previous models usually predict orogenic scale removal of mantle lithosphere, and it assume laterally homogeneous rheological structure of the lithosphere, but natural orogens often consists of numerous terranes with different lithology, history, and thermal structure, and hence laterally heterogenous rheological structure. Could such rheological heterogeneity explain the observed small-scale lithosphere thinning? To answer this question, we used a high-resolution thermomechanical model to simulate strain evolution and gravitational instability of collisional orogenic lithosphere. Our results indicate that lateral rheological heterogeneity could cause small-scale/fragmental lithospheric thinning under orogens when plastic yielding strength is rather low. The key factors are the combined weakening effects, by power-law viscous rheology and strong plastic yielding, which promote strain localization and lead to small-scale removal of mantle lithosphere. The numerical results provide useful insight into the cause of the observed small-scale lithosphere thinning under northeastern Tibetan Plateau, central Andes, and the Pamir.
Keywords: Small-scale/fragmental lithospheric removal    Orogen    Heterogeneity    Power-law viscous rheology    Plasticity    numerical simulation    
0 引言

造山带下岩石圈地幔常有显著的减薄,如喜马拉雅-青藏高原、中安纳托利亚以及科罗拉多高原(Owens and Zandt, 1997; 李海鸥等, 2008;Göğüş and Pysklywec, 2008; Levander et al., 2011; Fichtner et al., 2013).通常的解释是拆沉或对流减薄(Bird, 1979; Molnar, 1988; Jiménez-Munt et al., 2008; Göğüş et al., 2017; Li et al., 2016).然而,拆沉或对流减薄模型一般用于造山带尺度的减薄,难以解释中安第斯,帕米尔以及青藏高原东北缘下所见的分段式与小尺度的岩石圈减薄(图 1).这些减薄的岩石圈局部而离散,有些地方岩石圈地幔剥离到莫霍面(Beck and Zandt, 2002; Schurr et al., 2006; Kufner et al., 2016; Chapman et al., 2018; Deng et al., 2018), 并在地表导致小区域分布的岩浆事件(Kay et al., 1994; Schurr et al., 2006; Ducea, 2011; Chapman et al., 2018).这些局部的岩石圈减薄,既不同于拆沉或对流减薄模型所预测的空间尺度,也没有类似于拆沉/对流减薄模型所预测的时空迁移性或对称性分布(Göğüş and Pysklywec, 2008; Lei et al., 2019).

图 1 对流减薄、拆沉与分段式/小尺度减薄概念图. (a)由Rayleigh-Taylor所驱动的岩石圈地幔对流减薄模式;(b)岩石圈地幔沿下地壳剥离的拆沉模式;(c)岩石圈地幔分段式减薄 Fig. 1 The conceptual models of convective thinning, delamination and small-scale removing of continental lithosphere under orogens. (a) Convective thinning of mantle lithosphere by Rayleigh-Taylor instability; (b) Delamination of mantle lithosphere peeling off from crust; (c) Small-scale removing of mantle lithosphere

岩石圈地幔的这种分段式/小尺度减薄,与经典拆沉/对流减薄模式也存在动力学上的冲突.理论上讲,拆沉与对流减薄均由局部增厚体的Rayleigh-Taylor重力不稳定性所驱动,虽然二者的时间尺度与流变学背景存在着明显差异(Beall et al., 2017; Lei et al., 2019).造山带岩石圈地幔挤压增厚的尺度越大,其重力不稳定性亦越大,故经典拆沉/对流减薄模式通常预测大尺度的岩石圈地幔拆离减薄.同样道理,小尺度的重力不稳定性很难造成岩石圈地幔发生显著的减薄(Houseman and Molnar, 1997; Conrad and Molnar, 1997; Currie et al., 2008; Ducea, 2011).

岩石圈的弱化是其被破坏和减薄的关键条件.位错蠕变与塑性屈服均可导致一定程度的应变集中,并进而弱化岩石圈地幔,从而可导致造山带之下与汇聚板块交接处岩石圈地幔减薄并进而引发拆沉(Lei et al., 2019).大陆碰撞过程中的多地体拼接理论可以解释某些造山高原的大尺度岩石圈俯冲和拆沉问题(Huangfu et al., 2018; 皇甫鹏鹏等, 2019),然而其并未解释造山带之下岩石圈的局部分段式/小尺度减薄.基于前人研究,大陆岩石圈流变学的弱化机制为小尺度/分段式减薄的发生提供了动力学上的可能性,然而仍然存在着一些不确定性.造山带之下岩石圈地幔流变性质横向不均一性,是否会造成其局部分段式、小尺度减薄?产生显著的小尺度减薄,需要什么样的位错蠕变与塑性屈服呢?为了回答这些问题,我们建立了二维高精度热-动力学数值模型,对小尺度减薄的动力过程与流变学条件进行了系统的数值模拟研究.

1 数值模拟算法与模型设计简介 1.1 控制方程

活动造山带之下岩石圈地幔小尺度/分段式减薄问题相关的动力学数值模型,建立在三组控制方程上:物质守恒、动量守恒、及能量守恒方程.基于Marker-in-cell方法,我们在不规则的欧拉网格上针对不同黏塑性模型建立了一系列高精度二维有限差分模型,对这三组控制方程进行相应离散化并近似求解(Gerya and Yuen, 2003).

(1) 不可压缩流体的物质守恒方程

模型中的物质守恒采用了不可压缩流体假设下的连续性方程(考虑扩展的Boussinesq近似):

(1)

式中xy指示模型水平与垂直坐标方向;vxvy分别对应xy方向的速度.

(2) 动量守恒方程(Stokes方程)

(2)

式中P为压力;g为重力加速度;σ′xxσ′xyσ′yy代表偏应力张量分量,并且它们相应的本构关系如下所示:

其中

(3)

式中ηeff表示依赖于组分,温度,压力以及应变率的有效黏滞系数;则代表偏应变率分量.

对于不同岩石圈类型,其密度ρ是压力(P)与温度(T)的函数:

(4)

式中ρ0代表岩石在P0=0.1 MPa和T0=298 K温压条件下的标准密度.αβ分别代表热膨胀系数和可压缩系数.

(3) 能量守恒方程

(5)

式中T为温度,Cp是等压热容;qxqy是热流值;k是热传导系数,其是物质(C),温度(T)以及压力(P)的函数;Hr, HaHs分别代表放射性生热,绝热变压生热和摩擦生热.

1.2 黏-塑性流变学性质

(1) 黏性流变模型

岩石圈地幔的黏性流变特征由位错蠕变定义:

(6)

式中表示应变率张量第二主不变量;η0表示物质常数;E, Vn则分别表示活化能,活化体积和偏应力幂级(表 1).R代表气体常数.

表 1 数值模型采用的黏性流变参数a Table 1 Viscous flow laws used in the numerical experimentsa

(2) Drucker-Prager塑性屈服准则

广义Drucker-Prager塑性屈服准则(Ranalli, 1995)被应用于我们的数值模型中来模拟岩石圈地幔的黏塑性流变行为

(7)

式中ηplastic代表Drucker-Prager塑性等效黏滞系数,σyield表示屈服应力,P为压力,C0P=0时的岩石剩余强度,φdry表示无水化作用岩石的内摩擦角.λ表示塑性屈服系数,其控制着存在流体/熔体活动作用下岩石的脆性强度:, 由于流体/熔体压力Pfluid/melt并不能从模型中直接计算获得,因此我们采用‘sin(φdry)λ’的范围值来表征岩石的有效摩擦系数(Gerya and Meilick, 2011; Li et al., 2010, 2016; Vogt et al., 2012).值得注意的是,应变弱化效应被考虑在塑性流变模型中,这表明岩石的内聚强度C0与有效摩擦系数sin(φdry)λ会随着应变的增加而降低(表 2).

表 2 数值模型采用的物质参数a Table 2 Material properties used in the numerical experimentsa

黏-塑性流变模型的有效黏滞系数被定义为韧性与塑性黏滞系数中的最小值:

(8)

1.3 模型设计与边界条件

我们采用大尺度几何模型(4000 km×400 km)来研究碰撞造山带之下岩石圈地幔的小尺度/分段式减薄问题.模型中的碰撞区域与其他区域分别以1 km×1 km和5 km×1 km分辨率网格所表征.虽然模型长度达到4000 km,然而我们所关注的变形区域却集中在非常窄的尺度内(~250 km),因此地球曲率的影响可以忽略.

模型的初始岩石圈由20 km厚的上地壳,15 km厚的下地壳以及105 km厚的岩石圈地幔组成.模型中岩石圈包括强地体与弱地体(图 2),模型中不同块体所对应的流变参数在表 1表 2中给出.参考模型中,我们假设地幔岩石圈与下伏软流圈具有相同的参考密度ρ0(表 2),因此模型中地幔岩石圈的重力不稳定性主要由增厚过程中的温压变化所导致,如方程(4)所描述.岩石圈中的初始温度结构(图 2a中白线)侧向均一,在深度上从表面0 ℃线性增温至岩石圈地幔底部的1350 ℃(Turcotte and Schubert, 2002).软流圈的初始绝热温度梯度为0.5 ℃/km.为了使得岩石圈顶部可以自由隆升或沉降,我们在模型顶部采用了“sticky air”黏性空气层(Gerya and Yuen, 2003),其具体表现为位于模型地壳之上的一个初始10 km厚的低黏性层.

图 2 初始数值模型与边界条件示例.(a)实际模型(4000×400 km)中的一部分(2000×400 km)物质场.白色线条为间隔为200 ℃的等温线.边界条件用黄色箭头表示,代表向模型左右两侧强地体岩石圈地幔内施加的初始的汇聚速率(Vx=0.5+0.5 cm·a-1~2.5+2.5 cm·a-1). (b)模型中颜色代表不同岩石类型,1, 2:强地体大陆上, 下地壳;3, 4:弱地体大陆上, 下地壳;5:强地体岩石圈地幔;6:弱地体岩石圈地幔;7:软流圈 Fig. 2 Reference model and boundary conditions. (a) Composition field. Enlargement (2000×400 km) of the numerical box (4000×400 km) is shown. White lines are isotherms in ℃ with the interval of 200 ℃. Initial convergence rates (Vx=0.5+0.5 cm·a-1~2.5+2.5 cm·a-1) is imposed to the small internal domains in the left and right parts of the strong-terrane, respectively (yellow arrow). (b) Colors indicate the different rock types, 1, 2: continental upper and lower crust of the strong terrane; 3, 4: continental upper and lower crust of the weak terrane; 5: lithospheric mantle of the strong terrane; 6: lithospheric mantle of the weak terrane; 7: the asthenosphere

模型的顶部及左右边界均采用了自由滑动的速度边界条件,与之不同,模型底部边界则采用了渗透性边界条件,其通过在模型底部近无限深度施加外部自由滑动边界条件来保证计算区域内部的整体物质守恒(Li et al., 2016).类似于一般的自由滑动边界条件,施加于模型底部的外部自由滑动边界条件具体表现为以下形式:ϑvx/ϑy=0和ϑvy/ϑy=-vyyexternal,式中Δyexternal表示模型底部边界到外部自由滑动边界之间的垂直距离(图 2a).此外,相向汇聚速率(Vx)额外施加于岩石圈左右两侧强地体之上, 其用来表征岩石圈所受到的挤压汇聚作用(图 2a).

模型顶部初始采用了0 ℃的温度边界条件,与此同时,模型两侧边界采用绝热条件.进一步地,固定温度(2000 ℃)被施加于距模型底部1000 km的外部边界之上,从而允许模型底部能够随着时间的演化而动态地计算其温度扰动(Li et al., 2016).

1.4 横向不均一性流变结构

造山带之下岩石圈地幔的横向不均一性在模型中用强-弱相间的岩石圈来代表(图 2a).我们分别采用“湿石英”与“斜长石An75”的流变参数来表征模型中强地体所对应的上、下地壳;同时,用“无水橄榄岩”的流变参数来表征其岩石圈地幔的流变强度(表 1表 2).弱地体的上、下地壳均由“湿石英岩”表征;其地幔岩石圈的强度则由“含水橄榄岩”的流变性质决定.

2 数值模拟结果

我们的模拟旨在探讨岩石圈地幔的横向不均一性能否引起小尺度减薄.对比横向均一性模型,研究导致应变集中的位错蠕变和塑性屈服在岩石圈局部减薄中的作用.

2.1 横向不均一性下的应变集中

理论上,位错蠕变(方程(6))与塑性屈服(方程(7))均有助于局部应变集中,并能够弱化相应区域的岩石圈地幔.为了系统探讨由汇聚作用产生的应变集中是否能足够弱化岩石圈地幔并使其发生明显减薄,我们分别考虑了以下两种情况:(1)位错蠕变;(2)位错蠕变与塑性屈服的耦合效应.

2.1.1 位错蠕变的应变集中

在第一个实验中,我们考察位错蠕变对于弱地幔块体的应变集中效应.图 3展示了物质常数η0=5.01×1020Pa4·s,且汇聚速率νx=2.5+2.5 cm·a-1情况下的模型演化过程.结果显示,若造山带地幔的流变性质仅由湿橄榄岩的位错蠕变控制(方程(6)), 则无法产生局部的减薄.此时,汇聚会导致俯冲且产生明显的地表抬升(图 3).由汇聚所产生的应变集中在单一剪切带,贯穿整个弱地幔(图 4a),这主要由造山带之下岩石圈地幔的横向黏滞系数差异所导致(图 4b),在宏观上表现为强地体向弱地体之下的俯冲.

图 3 位错蠕变控制的地体变形特征.如表 1所示,模型中弱地体岩石圈地幔采用含水橄榄岩(n=4)流变特征,其中物质常数η0=5.01×1020Pa4·s.图中颜色代表岩石类型(图 2).白色线条代表间隔为200 ℃的温度场.演化时间如图(a)-(c)所示 Fig. 3 Model evolution with dislocation creep. As in Table 1, the rheology of wet olivine (n=4) is applied for the mantle lithosphere of weak terrane, with η0=5.01×1020Pa4·s. Colors indicate the different rock types (Fig. 2). White lines are isotherms with 200 ℃ intervals. Time (Ma) of evolution is shown in each panel (a)-(c)
图 4 应变率第二主不变量场(a)与有效黏滞系数场(b)演化特征.图中对应的物质场如图 3a所示 Fig. 4 The second invariants of the strain rate (a) and effective viscosity (b) of the model. Corresponding composition field is as shown in Fig. 3a
2.1.2 位错蠕变与塑性屈服的耦合作用导致的应变集中

在这个实验中,我们考察塑性屈服与位错蠕变的耦合作用对于弱地幔块体应变集中的影响.在数值模型中我们以塑性弱化系数λ的变化范围(方程(7))来体现弱地体岩石圈地幔的塑性屈服性质.图 5展示了汇聚速率νx=2.5+2.5 cm·a-1时,物质常数η0=5.01×1020Pa4·s,且λ=0.001时的模型演化过程.与前一实验不同,塑性屈服与位错蠕变的耦合效应允许小尺度的减薄出现.在汇聚作用的早期阶段,位于造山带之下左右两弱块体地幔均发生了局部增厚(图 5a),随着汇聚的进行,小尺度的重力不稳定最终减薄了岩石圈地幔并产生地表抬升(图 5b-5c).之后,同样的减薄过程也出现在位于造山带中部的弱地幔中(图 5d-5e).

图 5 塑性屈服和位错蠕变共同控制的地体变形特征.如表 1所示,模型中弱地体采用含水橄榄岩(n=4),其中物质常数与塑性屈服强度分别为η0=5.01×1020Pa4·s和λ=0.001,图中颜色代表岩石类型(图 2).白色线条代表间隔为200 ℃的温度场.演化时间如图(a-e)所示 Fig. 5 Model evolution with integrated dislocation creep and plasticity. As in Table 1, the rheology of wet olivine (n=4) is applied in the mantle lithosphere of weak terrane, with η0=5.01×1020Pa4·s and λ=0.001. Colors indicate the different rock types (Fig. 2). White lines are isotherms with 200 ℃ intervals. Time (Ma) of evolution is shown in each panel (a-e)

图 6展示了模型的应变率第二主不变量场(图 6a-6b)与有效黏滞系数场(图 6c-d)的演化.与前一实验不同,塑性屈服使得由汇聚导致的应变集中更加显著,在模型中呈现为多条局部共轭剪切的应变集中带(图 6a-6b).这足以弱化其所在区域的岩石圈地幔(图 6c-6d),并产生小尺度减薄.这一应变集中效应反映了位错蠕变与塑性屈服耦合作用下的应变率与有效黏滞系数ηeff的正反馈效应(方程(6)-(8).

图 6 应变率第二主不变量场(a-b)与有效黏滞系数场(c-d)的演化特征.图中对应的物质场如图 5b, d所示 Fig. 6 The second invariants of the strain rate (a-b) and effective viscosity (c-d) of the model.Corresponding composition field is as shown in Fig. 5b, d
2.1.3 不同塑性屈服强度的效应

我们进一步考察不同塑性屈服强度的作用.图 7展示了汇聚速率在νx=2.5+2.5 cm·a-1时,物质常数η0=5.01×1020Pa4·s,且λ=0.001~∞情况下的模型演化特征.这里λ=∞表明此时模型中的含水橄榄岩地幔未经历塑性屈服,其流变特征仅由位错蠕变机制所控制,即图 3的情况.结果表明,塑性屈服强度对于局部地幔的减薄至关重要.只有经历了强烈塑性屈服的地幔才会发生明显的小尺度减薄(图 7).若造山带地幔未经历塑性屈服,即λ=∞,俯冲将主导此处岩石圈地幔的失稳(图 7a).另一方面,若地幔的塑性屈服较弱,即λ=0.1,则会以岩石圈增厚的形式发生失稳,并不会发生减薄(图 7b).而当λ=0.01时,强塑性屈服会在弱地幔底部形成局部增厚,并最终与剩余地幔拆离.值得注意的是,这一过程并不会造成显著的地幔减薄(图 7c).只有当塑性屈服强度降低至λ=0.001时,才会形成岩石圈小尺度减薄(图 7d).

图 7 不同塑性屈服条件下的地体变形特征(λ=0.001~∞). λ=∞表明岩石圈地幔无塑性屈服作用 Fig. 7 Model evolution of terranes with different plastic yielding conditions (λ=0.001~∞).λ=∞ indicates no plastic yielding applied in the mantle lithosphere
2.2 小尺度地幔减薄的流变学条件

我们用不同流变学条件进行了系统的数值模拟研究.若造山带之下岩石圈地幔存在横向不均一性,则基于位错蠕变和不同塑性屈服强度的失稳结果由图 8表示.结果表明局部小尺度/分段式的减薄过程发生在塑性屈服强度很低(λ=0.001)的情况下,与汇聚速率无关.当塑性屈服强度λ=0.01~∞时,局部小尺度/分段式的减薄过程均不容易发生.

图 8 岩石圈地幔小尺度减薄过程的流变学条件.模型中弱地体地幔均采用基于位错蠕变机制的含水橄榄岩(n=4)的流变性质(表 1).图中λ=∞表明岩石圈地幔无塑性屈服作用,其流变性质仅由含水橄榄岩本身位错蠕变机制控制 Fig. 8 Rheological conditions of small-scale removing of mantle lithosphere.Power-law viscous rheology of wet olivine (Table 1) is applied for the weak terranes. In the diagrams, λ=∞ indicates no plastic yielding is applied in the mantle lithosphere of weak terrane
2.3 对比模型:横向均一性下的应变集中

此实验中,我们采用了横向均一性模型设置(Lei et al., 2019)及位错蠕变与塑性屈服耦合的流变性质.图 9展示了汇聚速率νx=2.5+2.5 cm·a-1时,参考黏滞系数η0=5.01×1020Pa4·s,且λ=0.001情况下横向均一的弱地幔减薄的特点.结果显示塑性屈服与位错蠕变的耦合效应使得弱地幔块体发生快速而完全的拆沉,其归因于位错蠕变与塑性屈服耦合作用所产生的应变集中效应(方程(6)-(8)),最终导致整个岩石圈地幔发生拆沉,与前文的小尺度、分段式的岩石圈减薄(图 5)显著不同.

图 9 塑性屈服和位错蠕变共同控制的横向均一地体变形特征.如表 1所示,模型中横向均一弱地体采用含水橄榄岩(n=4),其中物质常数与塑性屈服强度分别为η0=5.01×1020Pa4·s和λ=0.001,图中颜色代表岩石类型(图 2).白色线条代表间隔为200 ℃的温度场.演化时间如图(a-d)所示 Fig. 9 Model evolution of homogeneous weak terrane with integrated dislocation creep and plasticity. As in Table 1, the rheology of wet olivine (n=4) is applied for the mantle lithosphere of homogeneous weak terrane, with η0=5.01×1020Pa4·s and λ=0.001. Colors indicate the different rock types (Fig. 2). White lines are isotherms with 200 ℃ intervals. Time (Ma) of evolution is shown in each panel (a-d)
3 讨论 3.1 小尺度岩石圈减薄的控制因素

通过数值模拟,我们系统探讨了造山带下岩石圈地幔小尺度/分段式减薄的流变学条件.这里,我们总结不同流变学条件与其所导致的重力失稳过程间的动力学联系.

首先,造山带岩石圈地幔流变性质的横向不均一性为其分段式/小尺度减薄过程的发生提供了基本的动力学条件.这样,由汇聚所导致的应变集中会优先出现在较弱岩石圈地幔一侧(含水橄榄岩).此效应会进一步弱化岩石圈地幔并最终导致局部的减薄(图 5图 6).

而位错蠕变与塑性屈服作用本身对地幔块体的重力失稳类型至关重要(Lei et al., 2019).它们决定局部扰动是否会有效地减薄岩石圈地幔,位错蠕变本身并不足以导致局部岩石圈减薄, 其应变集中效应只导致强地体地幔向弱地体地幔之下的俯冲(图 3图 4).但塑性屈服的加入则允许局部扰动自发地形成于造山带之下,强烈弱化岩石圈地幔,使其发生失稳.这个扰动是否足以减薄弱化的岩石圈地幔,取决于塑性屈服强度.我们的模型结果表明,明显的减薄只可能发生在位错蠕变与强塑性屈服作用下的弱地幔岩石圈块体中.

3.2 地质启示

鉴别碰撞造山带之下岩石圈地幔减薄方式的地质观测通常是由地表隆升历史与岩浆活动分布两方面证据提供(Göğüş and Pysklywec, 2008; Li et al., 2016; Huangfu et al., 2018Lei et al., 2019).然而,与分段式/小尺度减薄相关的地表隆升历史与岩浆演化分布均无法体现出类似于拆沉/对流减薄所呈现的明显迁移性/对称性特征.

另一方面,由局部增厚所导致的软流圈物质的绝热上涌一直是用于判断造山带之下岩石圈地幔减薄与否的关键指标(Kay and Kay, 1993; Ducea and Saleeby, 1998),然而触发软流圈发生明显上涌所需要的局部重力扰动尺度往往较大,相比之下,更小尺度的重力扰动一般需要大约数十个百万年才能造成类似的软流圈上涌规模并进而引发相应的岩浆活动(Ducea, 2011),这似乎难以解释某些造山高原所出现的局部的快速(约2~3 Ma内)岩浆事件(Kay and Kay, 1993; Kay et al., 1994; Chapman et al., 2018).我们的数值模型结果表明,造山带下岩石圈地幔的流变横向差异,加以强烈塑性作用,允许局部重力扰动快速发展并在约2~3 Ma内产生明显的软流圈上涌(图 5图 7).这种塑性作用可来自于与俯冲/碰撞过程相关的一系列流体/熔体活动.

3.2.1 青藏高原东北缘

分段式/小尺度减薄过程可能已经/正在青藏高原东北缘之下的岩石圈地幔中进行(Deng et al., 2018).此区域由南向北由不同地块拼接而成(图 10),根据Deng等(2018)的地震波联合反演结果显示,位于这些造山带之下的地壳与地幔岩石圈内的流变强度存在明显的横向不均一性, 并呈现出显著的分段移除特征(图 10a),这些特征均能较好地与低汇聚速率下(vx=0.5+0.5 cm·a-1)考虑强塑性屈服条件(λ=0.001)的数值模型结果相吻合(图 10b).松潘-甘孜与南昆仑-柴达木地体岩石圈地幔已发生完整的移除, 而位于北昆仑-柴达木地体以及中祁连造山带下仍保留着部分残余岩石圈地幔.相比之下,北祁连造山带下地幔岩石圈可能随着部分下地壳一起发生了完整的移除并导致了局部莫霍界面深度的突然降低,之后软流圈物质在此处发生绝热上涌并进一步导致了 < 6 Ma时间尺度下的岩浆活动事件(Xia et al., 2011).与此同时,流变强度最强的华北克拉通岩石圈地幔则较为完整,未发生显著的地幔减薄.

图 10 青藏高原东北缘分段式、小尺度减薄的机制分析. (a)青藏高原东北缘岩石圈地幔变形示意图(修改自Deng et al., 2018),其中SGT:松潘-甘孜地体;KQT:昆仑-柴达木地体;CQL:中祁连;NQL:北祁连;NCC:华北克拉通.竖直虚线代表区分不同地块之间的主要断层, KF:昆仑断层;SQS:南祁连断层;HF:海原断裂;GF:古浪断裂. ML:地幔岩石圈;RML:残余地幔岩石圈. (b)基于含水橄榄岩位错蠕变(n=4)与强塑性屈服作用(λ=0.001)下模型结果.汇聚率Vx=0.5+0.5 cm·a-1 Fig. 10 Fragmental and small-scale lithospheric removal under the northeastern Tibetan Plateau. (a) Deformation of mantle lithosphere under Northeastern Tibetan Plateau (modified from Deng et al., 2018), SGT: Songpan-Garzê Terrain; KQT: Kunlun-Qaidam Terrain; CQL: Central Qilian; NQL: North Qilian; NCC: North China Craton. Vertical dashed lines represent the main fault motions, KF: Kunlun fault; SQS: South Qilian suture; HF: Haiyuan fault. ML: Mantle lithosphere; RML: Residual mantle lithosphere. (b) Model result based on power-law viscous rheology of wet olivine (n=4) and strongly plastic yielding (λ=0.001), with convergence rates of Vx=0.5+0.5 cm·a-1
3.2.2 其他的造山带

类似的局部减薄过程可能普遍发生在造山高原之下.Altiplano高原与东科迪勒拉造山带之间的岩石圈地幔内,地震波数据表明位于此处的岩石圈地幔已发生明显的局部移除,而与其相邻的阿尔蒂普拉诺之下岩石圈地幔却保存完整.进一步地,基于岩石地球化学方面的证据也表明位于此区域以及普纳高原之下岩石圈地幔已经/正在经历与小尺度移除(~数十公里)过程相关的熔融事件(Ducea et al., 2013).这种岩石圈地幔部分移除的方式可能与纳斯卡板片低角度俯冲所导致的流体/熔体活动以及巴西克拉通之下岩石圈地幔的西向汇聚相关(Beck and Zandt, 2002),此两种过程促使局部应变集中在阿尔蒂普拉诺与东科迪勒拉造山带之间的岩石圈地幔区域并最终导致其发生移除.

另一个例子是帕米尔高原,有地震学数据表明位于帕米尔与辛都库什之间交界地区之下岩石圈地幔发生了较为明显地小尺度移除现象(Kufner et al., 2016),可能与印-亚大陆碰撞所导致的一系列流体/熔体弱化作用相关.在印度大陆岩石圈碰撞汇聚之前,冈瓦纳地体(当今由帕米尔, 喀喇昆仑以及青藏岩石圈组成)岩石圈就已经历了与古特提斯洋闭合事件相关的一系列俯冲,拼合与岛弧生长等构造活动(Kumar et al., 2001),这些过程均可以弱化其岩石圈地幔并使其发生塑性屈服,导致局部岩石圈地幔移除.帕米尔高原的岩浆岩区域也有局部的分布,其中一个例子便是分布于中帕米尔高原的Vanj岩浆岩.相应地球化学证据表明Vanj岩浆事件曾发生于至今42~36 Ma的始新世时期,极有可能由发生于岩石圈地幔的小尺度熔融过程导致(Chen et al., 2013; Chapman et al., 2018).

4 结论

利用二维高精度数值模型,我们系统研究了碰撞造山带下岩石圈地幔发生小尺度/分段式减薄过程的动力学机制,得到以下三个结论:

(1) 造山带岩石圈地幔流变性质的横向不均一性能够使得由汇聚所导致的局部应变集中在流变特征显著差异处.这为岩石圈地幔分段式/小尺度减薄提供了必要的动力学条件.

(2) 地幔岩石圈在有水、部分熔融时会发生强烈的塑性屈服作用,位错蠕变机制与强塑性屈服作用(λ=0.001时)共同促使发生于此处的应变集中以多条局部共扼剪切带的形式呈现.此种流变弱化机制允许局部重力失稳快速发育,并于此处最终发生减薄.

(3) 小尺度、分段式岩石圈减薄模型可用于解释青藏高原东北缘及其他造山带的深部观测数据.

致谢  感谢皇甫鹏鹏博士的有益讨论和帮助.感谢两位审稿人的建设性意见.
References
Beall A P, Moresi L, Stern T. 2017. Dripping or delamination? A range of mechanisms for removing the lower crust or lithosphere. Geophysical Journal International, 210(2): 671-692. DOI:10.1093/gji/ggx202
Beck S L, Zandt G. 2002. The nature of orogenic crust in the central Andes. Journal of Geophysical Research:Solid Earth, 107(B10): ESE 7-1-ESE 7-16. DOI:10.1029/2000JB000124
Bird P. 1979. Continental delamination and the Colorado Plateau. Journal of Geophysical Research:Solid Earth, 84(B13): 7561-7571. DOI:10.1029/JB084iB13p07561
Bittner D, Schmeling H. 1995. Numerical modeling of melting processes and induced diapirism in the lower crust. Geophysical Journal International, 123(1): 59-70. DOI:10.1111/j.1365-246X.1995.tb06661.x
Chapman J B, Scoggin S H, Kapp P, et al. 2018. Mesozoic to Cenozoic magmatic history of the Pamir. Earth and Planetary Science Letters, 482: 181-192. DOI:10.1016/j.epsl.2017.10.041
Chen J L, Wu J B, Xu J F, et al. 2013. Geochemistry of Eocene high-Mg# adakitic rocks in the northern Qiangtang terrane, central Tibet:implications for early uplift of the plateau. GSA Bulletin, 125(11-12): 1800-1819. DOI:10.1130/B30755.1
Clauser C, Huenges E. 1995. Thermal conductivity of rocks and minerals.//Ahrens T J ed. Rock Physics & Phase Relations: A Handbook of Physical Constants. Washington D. C.: AGU Reference Shelf, 105-126.
Conrad C P, Molnar P. 1997. The growth of Rayleigh-Taylor-type instabilities in the lithosphere for various rheological and density structures. Geophysical Journal International, 129(1): 95-112. DOI:10.1111/j.1365-246X.1997.tb00939.x
Currie C A, Huismans R S, Beaumont C. 2008. Thinning of continental backarc lithosphere by flow-induced gravitational instability. Earth and Planetary Science Letters, 269(3-4): 436-447. DOI:10.1016/j.epsl.2008.02.037
Deng Y F, Li J T, Song X D, et al. 2018. Joint Inversion for Lithospheric Structures:Implications for the Growth and Deformation in Northeastern Tibetan Plateau. Geophysical Research Letters, 45(9): 3951-3958. DOI:10.1029/2018GL077486
Ducea M, Saleeby J. 1998. A case for delamination of the deep batholithic crust beneath the Sierra Nevada, California. International Geology Review, 40(1): 78-93. DOI:10.1080/00206819809465199
Ducea M N. 2011. Fingerprinting orogrnic delamination. Geology, 39(2): 191-192.
Ducea M N, Seclaman A C, Murray K E, et al. 2013. Mantle-drip magmatism beneath the Altiplano-Puna plateau, central Andes. Geology, 41(8): 915-918. DOI:10.1130/G34509.1
Fichtner A, Saygin E, Taymaz T, et al. 2013. The deep structure of the North Anatolian fault zone. Earth and Planetary Science Letters, 373: 109-117. DOI:10.1016/j.epsl.2013.04.027
Gerya T V, Yuen D A. 2003. Characteristics-based marker-in-cell method with conservative finite-differences schemes for modeling geological flows with strongly variable transport properties. Physics of the Earth and Planetary Interiors, 140(4): 293-318. DOI:10.1016/j.pepi.2003.09.006
Gerya T V, Meilick F I. 2011. Geodynamic regimes of subduction under an active margin:Effects of rheological weakening by fluids and melts. Journal of Metamorphic Geology, 29(1): 7-31. DOI:10.1111/j.1525-1314.2010.00904.x
Göğüş O H, Pysklywec R N. 2008. Near-surface diagnostics of dripping or delaminating lithosphere. Journal of Geophysics Research:Solid Earth, 113(B11): B11404. DOI:10.1029/2007JB005123
Göğüş O H, Pysklywec R N, Şengör A M C, et al. 2017. Drip tectonics and the enigmatic uplift of the Central Anatolian Plateau. Nature Communications, 8: 1538. DOI:10.1038/s41467-017-01611-3
Houseman G A, Molnar P. 1997. Gravitational (Rayleigh-Taylor) instability of a layer with non-linear viscosity and convective thinning of continental lithosphere. Geophysical Journal of the Royal Astronomical Society, 128(1): 125-150. DOI:10.1111/j.1365-246X.1997.tb04075.x
Huangfu P P, Li Z H, Gerya T, et al. 2018. Multi-terrane structure controls the contrasting lithospheric evolution beneath the western and central-eastern Tibetan plateau. Nature Communications, 9: 3780. DOI:10.1038/s41467-018-06233-x
Huangfu P P, Li Z H, Fan W M. 2019. Modes and dynamics of multi-terrane deformation during continental collision. Chinese Journal of Geophysics (in Chinese), 62(4): 1230-1243.
Ji S C, Zhao P L. 1993. Flow laws of multiphase rocks calculated from experimental data on the constituent phases. Earth and Planetary Science Letters, 117(1-2): 181-187. DOI:10.1016/0012-821X(93)90125-S
Jiménez-Munt I, Fernàndez M, Vergés J, et al. 2008. Lithosphere structure underneath the Tibetan Plateau inferred from elevation, gravity and geoid anomalies. Earth and Planetary Science Letters, 267(1-2): 276-289. DOI:10.1016/j.epsl.2007.11.045
Kay R W, Kay S M. 1993. Delamination and delamination magmatism. Tectonophysics, 219(1-3): 177-189. DOI:10.1016/0040-1951(93)90295-U
Kay S M, Coira B, Viramonte J. 1994. Young mafic back arc volcanic rocks as indicators of continental lithospheric delamination beneath the Argentine Puna plateau, central Andes. Journal of Geophysical Research:Solid Earth, 99(B12): 24323-24339. DOI:10.1029/94JB00896
Kirby S H. 1983. Rheology of the lithosphere. Reviews of Geophysics, 21(6): 1458-1487. DOI:10.1029/RG021i006p01458
Kirby S H, Kronenberg A K. 1987. Rheology of the lithosphere:Selected topics. Reviews of Geophysics, 25(6): 1219-1244. DOI:10.1029/RG025i006p01219
Kufner S K, Schurr B, Sippl C, et al. 2016. Deep India meets deep Asia:Lithospheric indentation, delamination and break-off under Pamir and Hindu Kush (Central Asia). Earth and Planetary Science Letters, 435: 171-184. DOI:10.1016/j.epsl.2015.11.046
Kumar M R, Saul J, Sarkar D, et al. 2001. Crustal structure of the Indian shield:new constraints from teleseismic receiver functions. Geophysical Research Letters, 28(7): 1339-1342. DOI:10.1029/2000GL012310
Lei T, Li Z H, Liu M. 2019. Removing mantle lithosphere under orogens:delamination versus convective thinning. Geophysical Journal International, 219(2): 877-896. DOI:10.1093/gji/ggz329
Levander A, Schmandt B, Miller M S, et al. 2011. Continuing Colorado plateau uplift by delamination-style convective lithospheric downwelling. Nature, 472(7344): 461-465. DOI:10.1038/nature10001
Li H O, Xu X W, Jiang M. 2008. Deep dynamical processes in the central-southern Qinghai-Tibet Plateau-Receiver functions and travel-time residuals analysis of north Hi-climb. Science in China Series D:Earth Sciences, 51(9): 1297-1305. DOI:10.1007/s11430-008-0096-2
Li Z H, Gerya T V, Burg J P. 2010. Influence of tectonic overpressure on P-T paths of HP-UHP rocks in continental collision zones:thermomechanical modelling. Journal of Metamorphic Geology, 28(3): 227-247. DOI:10.1111/j.1525-1314.2009.00864.x
Li Z H, Liu M, Gerya T V. 2016. Lithosphere delamination in continental collisional orogens:A systematic numerical study. Journal of Geophysical Research:Solid Earth, 121(7): 5186-5211. DOI:10.1002/2016JB013106
Molnar P. 1988. A review of geophysical constraints on the deep structure of the Tibetan Plateau, the Himalaya and the Karakoram, and their tectonic implications. Philosophical Transactions of the Royal Society A Mathematical Physical and Engineering Sciences, 326(1589): 33-88.
Owens T J, Zandt G. 1997. Implications of crustal property variations for models of Tibetan plateau evolution. Nature, 387(6628): 37-43. DOI:10.1038/387037a0
Ranalli G, Murphy D C. 1987. Rheological stratification of the lithosphere. Tectonophysics, 132(4): 281-295. DOI:10.1016/0040-1951(87)90348-9
Ranalli G. 1995. Rheology of the Earth. 2nd ed. London: Chapman Hall.
Schurr B, Rietbrock A, Asch G, et al. 2006. Evidence for lithospheric detachment in the central Andes from local earthquake tomography. Tectonophysics, 415(1-4): 203-223. DOI:10.1016/j.tecto.2005.12.007
Turcotte D L, Schubert G. 2002. Geodynamics. 2nd ed. New York: Cambridge University Press: 456.
Vogt K, Gerya T V, Castro A. 2012. Crustal growth at active continental margins:Numerical modeling. Physics of the Earth and Planetary Interiors, 192-193: 1-20. DOI:10.1016/j.pepi.2011.12.003
Xia L Q, Li X M, Ma Z P, et al. 2011. Cenozoic volcanism and tectonic evolution of the Tibetan Plateau. Gondwana Research, 19(4): 850-866. DOI:10.1016/j.gr.2010.09.005
皇甫鹏鹏, 李忠海, 范蔚茗. 2019. 大陆碰撞过程中的多地体变形模式. 地球物理学报, 62(4): 1230-1243. DOI:10.6038/cjg2019M0607
李海欧, 徐锡伟, 姜枚. 2008. 青藏高原中南部的深部地球动力学过程-Hi-climb剖面北段接收函数和走时残差分析. 中国科学D辑:地球科学, 38(5): 622-629.