克拉通岩石圈重力不稳定过程的数值模拟及其动力学响应

1、克拉通岩石圈重力不稳定过程的数值模拟及其动力学响应作为大陆内部最古老的组成部分,克拉通一般认为可以抵抗地幔对流的作用而稳定存在长达几十亿年。克拉通地区通常具有低的表面热流和厚的、冷的地幔岩石圈。相比于下覆软流圈地幔而言,克拉通地幔岩石圈在化学组分上是亏损的,因而具有低的化学组分密度从而是有浮力的。熔融亏损作用以及低温使得克拉通岩石圈地幔具有较高的强度（粘性）。克拉通岩石圈的化学组分浮力以及高粘性有助于维持克拉通的稳定性,因而在克拉通的演化过程中起着重要的作用。然而,并不是所有的克拉通都是一直保持稳定的,人们同样发现有一部分克拉通遭到了一定程度的破坏和改造。因此,研究克拉通的稳定与破坏过程是理解

2、大陆演化历史必不可少的组成部分。中国的华北克拉通是克拉通破坏的一个典型的例子。华北克拉通东部自中生代-新生代以来发育有大规模的岩浆活动,而且地球化学和地球物理学研究表明华北克拉通东部不仅是岩石圈的厚度和热状态发生了改变,其岩石圈地幔的化学组分也发生了改变。尽管目前大量观测证据证实华北克拉通遭到了破坏,但是仍然有很多问题还有待解决,例如克拉通破坏的动力学过程,克拉通破坏的时空尺度等等。为解释不同的地球化学观测证据人们提出了不同的克拉通破坏机制,例如拆沉、热-化学侵蚀、熔体-橄榄岩相互作用以及水化作用等。然而这些破坏机制都缺乏地球动力学模型的支持。本文的主要研究目的就是通过利用动力学数值模拟方法研

3、究克拉通岩石圈前人关于克拉通岩石圈重力不稳定过程的研究主要是实验室条件下基于牛顿流体流变性质的讨论。本文建立二维热-化学对流模型,研究非牛顿流体流变条件下克拉通岩石圈在重力不稳定性作用下的动力学演化过程及其引起的表面响应,包括不稳定过程的演化模式、时间和空间尺度,以及地形和热流响应等等。数值结果显示,非牛顿流变性质可以有效的减小克拉通岩石圈的粘性,因而在克拉通岩石圈的不稳定过程中起着非常重要的作用。在非牛顿流变条件下,克拉通岩石圈的不稳定过程表现为分期、多阶段的特征而且持续时间超过100Myr。华北克拉通东部新生代和中生代时期岩浆活动具有幕式分布特征而且持续时间接近100Myr。因此,我们的数

4、值模拟结果可以很好的解释华北克拉通破坏时期的岩浆活动特征。而且我们发现,克拉通岩石圈的不稳定过程开始于岩石圈的浅部,这使得在破坏早期有大量的热的软流圈物质上涌到达下地壳底部,因此可以为下地壳的榴辉岩化提供足够的温度,从而为地球化学观测到的华北地区下地壳榴辉岩的存在提供动力学支持。同时,一部分失稳的克拉通岩石圈物质会返回到岩石圈层并与上升的软流圈地幔物质混合,这也解释了华北地区古老的克拉通岩石圈地幔和新生海洋型岩石圈地幔并存的地球化学观测。因此,本文的模拟结果所给出的岩石圈不稳定过程为华北克拉通破坏的动力学过程提供了很好的解释。地形和热流是岩石圈不稳定过程在表面最直接的响应,主要受岩石圈热-化学结构变化的控制,因此可以用来约束克拉通破坏的动力学过程。数值模拟结果显示,克拉通岩石圈的重力不稳定过程会抬升失稳区域的地形和热流,这与华北克拉通东高西低的残余地形和热流特征基本一致。我们的数值模拟也为定量化的研究克拉通破坏以后古老的岩石圈地幔物质残留提供了一个可行的办法。根据华北克拉通东西部