西太平洋下地幔D”层的地震波速度各向异性研究

摘要

位于地幔底部数百公里的D”层是固体地幔和液态外核的动力学、热学和化学边界，它控制着CMB物质、能量的交换，对地幔对流、板块的运动、磁场的变化有关键性的影响。地震学研究表明D”层具有明显的地震各向异性。地震各向异性是理解地球动力学过程的重要工具，携带了大量的关于固体地球演化与动力学过程的潜在信息。在合适的条件下，地震各向异性可以为采用其它方法(如层析成像)观测不到的地幔对流提供地震学上的指示。因此，对D”层地震各向异性位置、强度和方向的研究有助于深化对核幔边界附近区域物理性质的认识。西太平洋地区是地球上现今构造活动非常活跃、复杂的地区之一，其深地幔中可能存在非常复杂的结构，并且和中国大的地质活动背景有关。因此西太平洋下D”层地震各向异性的研究对于认识该区域的地球动力学过程及地球演化具有十分重要的意义。 首先，本文利用1994-2007年间发生在汤加-斐济及周边地区的146个地震，由IRIS36个台站记录到的数据和2000-2005年发生在东亚北部及北太平洋地区的26个地震，由IRIS的18个台站记录到的数据，使用ScS-S的相对走时分析方法研究了西太平洋下D”层的剪切波速度各向异性。我们得到了512个ScS波径向分量和横向分量的分裂时间，并计算了D”层的地震波速度各向异性强度。发现沿近北．南方向传播的地震波和沿近北西-南东方向传播的地震波产生的剪切波分裂模式差别较大：前者以V<,ScSH>＞V<,ScSV>的分裂为主，后者以V<,ScSH>＜V<,ScSV>(V<,ScSH>和V<,ScSV>分别表示水平和垂直极化的ScS震相的速度)的分裂为主，并且后者的各向异性强度明显要大于前者。整体而言，剪切波分裂呈现出V<,ScSH>＜V<,ScSV>的模式，并且V<,ScSH>＜V<,ScSV>的剪切波分裂主要分布在研究区的中东部，V<,ScSH>＞V<,ScSV>的剪切波分裂主要位于研究区的西南部。SoS的分裂时间值从-4.08s到4.53s，绝对值平均值分别为1.26s；各向异性强度值为从-1.96％到2.55％，绝对值平均值为0.61％；分裂时间的分布和各向异性强度的分布趋势一致。 用S-ScS相对走时分析方法研究D”层各向异性可以非常方便地消除震源一侧和台站一侧各向异性的影响，但是只能在径向和横向分量上测量，从而只能检测具有垂直对称轴的横向各向同性。如果要探测具有更一般形式的各向异性，则需要采用更复杂的算法。为此，本文利用来自发生在汤加-斐济俯冲区及周围地区的15个地震，由IRIS的8个台站记录到的数据，利用旋转相关法获取了对西太平洋下D”层采样的ScS震相的分裂参数-快分量的极化方向和快慢分量的延迟时间，并将其结果与相对走时方法的结果进行了比较。发现快波的极化方向(以正北方向为参考，顺时针方向为正)变化很大，从4°变化到170°(存在180°整数倍的不确定性)；快轴方向与径向的夹角从7.48°变化到73.86°，大部分分裂结果的夹角小于45°，说明研究区的大部分区域有V<,ScSV>＞V<,ScSH>；时间延迟从0.2s变化到3.9s，平均值为2.2s；其各向异性的强度值从0.11％变化到2.01％，平均值为1.06％。旋转相关方法与相对走时方法得出的结果对于V<,ScSV和V<,ScSH>相对大小的判断有75％相符，两者得出的延迟时间与各向异性强度存在一些差别。由于互相关计算中波形窗口的选择和时间平移等会引起误差，并且上地幔各向异性校正存在很大的不确定性，D”层剪切波各向异性的旋转相关方法分析结果可能会存在较大的误差，但是它可以在更大的方位角范围内探测快波的极化方向。地震各向异性是地球内部动力学过程的反映，利用地震各向异性的测量可以推测地幔不同部分的主要变形机制和流动特性。但是，地震各向异性的解释非常复杂，解释结果常常互相矛盾。由于震源与接收台站的地理局限性，目前所有的D”层各向异性研究的方位覆盖采样都是有限的，有限的方位覆盖很难区分D”层不同的各向异性模式，因此由分裂测量到D”层各向异性的解释必须作一些必要的假设。基于对西太平洋下D”层采样的剪切波分裂观测数据，本文详细讨论了西太平洋下D”层各向异性可能的流变场机制：西太平洋下D”层内垂直上升流应占支配地位；水平流动构造也可能同时存在，但尺度相对小一些，并主要分布在研究区的西南面。研究区下可能存在上升流所致的不均匀性物质定向排列成的垂直组构；下地幔物质的晶格优选方位也是一种可能的机制：研究区中可能存在偏离水平面的流动，使各向异性晶体或不均匀性物质定向排列形成方位各向异性，可以解释成具有倾斜对称轴的横向各向同性。 (Mg,Fe)SiO<,3>后钙钛矿可能是D”层的主要矿物，具有明显的弹性各向异性，为解释D”层地震波各向异性等现象提供了一种新的途径。后钙钛矿滑动系和所导致的晶格优选方位对于理解观测到的各向异性非常重要，但是人们对于后钙钛矿主滑移系的性质有不一致的认识。因此，针对复杂的D”层应力环境中晶体特征与应力场和晶轴取向的密切关系，采用第一原理模拟三轴应力场中的MgSiO<,3>后钙钛矿相的弹性性质，利用最小能量原理获得稳定的晶体空间取向类型。D”层的晶体取向和应力差对晶体性质产生明显的影响，其能量、弹性常数和波速随晶体取向和应力差的变化而变化：晶体a轴取向平行于最大压缩方向，以b轴为垂直对称轴的三轴应力场中的剪切波横向各向异性大于静水应力场中的剪切波横向各向异性，随着应力差的增加，具有最小能量的晶体空间取向类型比其他取向类型更为稳定，结果支持以(010)面为主滑移面的观点。在D”层垂直上升流区域，[010]水平取向时具有V<,SH>＜V<,SV>的各向异性：在水平流区域，[001]和[010]垂直取向均可产生V<,SH>＞V<,SV>的各向异性。

目录

文摘

英文文摘

声明

第一章引言

1.1概述

1.2 D”层的地震学观测

1.2.1 D”层的大尺度速度不均匀

1.2.2 D”层间断

1.2.3超低速区

1.2.4大型低速体

1.2.5 D”层散射体

1.2.6 D”层各向异性

1.3 D”层的本质及模型

1.4 D”层各向异性研究存在的问题

1.5本文的主要内容

第二章地震各向异性研究的理论与方法

2.1概述

2.2 Pn波各向异性

2.3面波各向异性

2.3 Ps波转换

2.5剪切波分裂

2.5.1测量剪切波分裂所用的震相

2.5.2各向异性对称系

2.5.3分裂参数的测量

2.6小结

第三章西太平洋下D”层各向异性的相对走时研究

3.1数据

3.2相对走时分析

3.3剪切波分裂的观测结果

3.3.1分裂时间

3.3.2各向异性强度

3.4结论与讨论

第四章西太平洋下D”层各向异性的旋转相关法研究

4.1旋转相关法

4.2数据资料及处理

4.3观测结果

4.4讨论

4.5小结

第五章D”层各向异性的流变场解释

5.1各向异性的描述

5.1.1推测各向异性对称系

5.1.2推测不均匀的各向异性

5.2各向异性的形成机制

5.2.1概述

5.2.2地球不同圈层的各向异性机制

5.3西太平洋下D”层的各向异性机制

5.3 小结

第六章D”层各向异性的后钙钛矿相变解释

6.1后钙钛矿相变

6.2后钙钛矿相变对地震学观测的解释

6.3三轴应力场中MgSiO3后钙钛矿的第一原理模拟

6.3.1计算方法

6.3.2计算流程

6.3.3计算结果

6.3.4结论与讨论

6.4 小结

第七章结论与展望

7.1结论与讨论

7.2展望

参考文献

附录

致谢

攻读博士学位期间发表的文章