20201129文献翻译_从660公里的边界地形推断出地球不连续的化学分层

从660公里的边界地形推断出地球不连续的化学分层

作者：Wenbo Wu1,2,3, Sidao Ni1*, Jessica C. E. Irving

摘要：固体地球中某些界面的地形或深度变化可以为我们了解地球内部动态（动力学）提供重要的见解。虽然对地球上地幔和下地幔之间的660公里边界的中期和长期地形变化已经得到了充分的研究，但小尺度的测量则更具挑战性。通过使用离散P’P’地震波，我们在某些区域的660公里边界面发现了惊人数量的短长度尺度的地形（变化）。我们的观察要求在短尺度粗糙度高的地区进行化学分层。相比之下，在上地幔410公里的边界上，我们并没有看到如此小规模的地形。我们的发现支持了上地幔和下地幔之间部分阻塞或不完全循环的概念

  全球观测到的660公里的地震不连续面定义了下地幔的顶部，通常被认为包含了矿物林伍德石到布里奇曼石和铁方镁石的相变转换带。这个界面的详细性质约束了整个地幔的化学和动力学特性。有几条线的证据支持边界仅仅是由相变引起的。如果是这样的话，自然的结论是整个地幔在地质时间尺度上(1,2)对流，尽管上地幔和下地幔之间存在着矿物学上的差异。然而，并非所有的观测结果都支持这种不连续现象。其他地球化学和矿物学证据表明存在一种化学界面，这种界面需要某种形式的主要的分层或分区对流，以保持上地幔和下地幔之间的化学差异(3-6)。地震波可以用来测量与边界物理特性有关的不连续面的许多特征，包括陡度、密度、弹性对比和地形变化。
  地形变化为理解这660公里不连续面的性质提供了必要的线索，地形是动力学过程和密度不均匀分布的结果。地球的自由表面和核幔边界（CMB）以1000公里到几公里(12-14)尺度的地形为特征，这导致了660公里不连续面在许多尺度上也可能是粗糙的预期。边界上的粗糙度尺度提供了对地形变化的动力学过程的洞察。
  我们目前对这660公里的不连续面地形的成像来自于反射体波的前驱波(如PP、SS和P’P’)(7,8,10)或转换震相，如P660s(15,16)和S660P(17,18)。这些方法从高达数十公里(7,10)的热变化(图1A)中揭示了大尺度(约1000公里)地形。中等尺度(约100公里)地形(图1A)已经用接收函数(15,16)或转换相位(17,18)绘制出来，因为它们的菲涅尔带较小。大部分的中等尺度和大尺度地形结果被解释为横向温度变化(19,20)，但一些研究表明，660公里界面的地形可能与更复杂的机制有关，而不仅仅是相变(21-23)。

图1。 660公里界面的PSD（功率谱密度）和P’•660•P’的射线路径。(A)二维(2D)自由表面地形(蓝色虚线)和660km界面地形(红色虚线)的PSD（功率谱密度）关于波数k ($k{m}^{-1}$)的函数。蓝色虚线为P = C2Dk−3,C2D = 0.3 m(33)，红色虚线分别为C2D = 10、100和1000 m。将660 km的间断地形的大、中、小横向尺度范围分别标记为“1”、“2”和“3”。由于地震台站的有限分布，中等尺度地形的采样还不够彻底。(B)P’•660•P’的射线路径。我们在赤道上设置了一个假想的源(红星)和接收台站(蓝色三角形)。与大多数在大圆平面(图中赤道)上传播的常规震相相反，非对称散射允许平面外散射波。( C)P’•660•P’(黑线)和P’ surfp '(灰线)震相在相应界面地形上散射的射线路径卡通图。黑色和灰色的虚线分别显示了P’P’和P’660P’的射线路径，它们分别经历了对称反射。

(from 25-2011BSSA_Earle_Scattered-PP.short-distances)

（from 24-2012GRL_Wu_asymmetrical-scattering)

  为了确定的660公里和410公里的地震不连续面的小尺度(约10公里)的地形变化，我们使用了短周期波散射。这种策略在上地壳是成功的，其中（上地壳中）P’surfP’波（震相）是由 上地壳中小尺度自由表面地形和/或上地壳非均质体 非对称(平面外：此平面指的是震源、台站、地心三点构成的面)后向散射产生的(24)。我们使用了P’•d•P’震相(25)，其中，（d表示）粗糙界面位于地球表面以下深度d处(图1、B和C)，而PKP在近焦散距离处的双重放大可以显著增强微弱信号。
  我们选择了小角震中距离(大约0°到40°)的地震波形，在这种距离下，非对称散射波P’•d•P’很容易观察到（24-26）。由于这660公里的不连续面比自由表面的电阻抗要低得多(27)，我们仔细地选择了低噪音和记录强震的地震台(表S1)，以增加观测到清晰P’•660•P’的机会。我们还使用了高质量的小孔径地震阵列，如IL (Eielson)、YK (Yellow Knife)和as (Alice Springs)来提高信噪比(SNR)和进行聚束分析。这三个阵列由于其卓越的性能被广泛用于检测各种类型的微弱信号(25,26,28,29)。

(from Supplementary)

(from obspy manual)
  本段说明图2A  原始地震数据通常是有噪声的，所以我们使用短周期滤波(1.5到2.5 Hz，除非另有说明，否则整个过程都使用)来增加信噪比。滤波后，数据中出现了两个强P’•d•P’震相信号(图2A)。到时约2280 s的最强信号主要是P’surfp ‘，部分能量可能来自P3KP(24)震相。另一个清晰的信号比P’surfp ‘的理论到达时间早150到160 s，这个150- 160-s的提前时间与P’•660•P’和P’surfp ‘之间的时间差（理论时间差）是一致的，这意味着前面的明显信号可能是P’•660•P’。如果我们取短周期地震波形的平滑包络线(图2A，底部轨迹)，前面的信号和P’surfp ‘都是纺锤体，这表明它们来源于散射(12,30)，尽管P’•660•P’的信噪比并不高。
  本段说明图2B  为了确认可能的反(背)向散射地震信号的性质，我们分析了鄂霍次克海下不同地震的IL阵列的地震图(图2B)。在叠加结果中，P’surfp ‘和P’•660•P‘震相都清晰可见，P3KP在两者之间到达(图2B)。P’•660•P‘震相包络的到间(~2130秒)和整体形状与LPAZ (La Paz，即图2A的最下面的图)的类似。此外，所有的19个地震波形的平滑包络线的都非常相似，这就要求地震相一致和稳健（表明这是一个明显的震相，而不是随机的噪音，因为不同的地震事件都有类似的波形特征），而不是局部噪声。叠加结果得到可靠的观测结果，信噪比大于2。
  本段说明图2C  我们采集了三个台阵（AS、IL、YK）记录的10次地震数据(表S1)，并重复了叠加过程(如图2c，叠加地震图，图S2，详细数据)。我们检查了国家地震信息中心(NEIC)的地震目录(31)，并排除（去除）了在所关注的时间窗内任何来自余震的污染信号的可能性。在角度震中距(0°)到~40°范围内，P’•660•P’震相几乎恒定在在P’ surfp ‘震相到达之前大约160秒(图2C)。 恒定的旅行时间，与角震中距离无关，是非对称反向散射(脱离大圆平面散射)的特有特征(24)，与其他通常观测到的沿大圆路径行进的地震震相相反(例如，图2C, P3KP)。因此，我们确认了该信号为P’•660•P’。然而，我们没有看到从410km散射的410•410•P震相的信号。

图2. 从单个台站和台阵观测到的P’.660.P’震相。 (A) LPAZ台站观测到的2003年6月20日南美7.0Mw地震事件波形的P’·660·P’震相。顶部为宽带速度波形数据。中间和底部分别是高频(1.5到2.5 Hz)滤波数据及其平滑包络线。 (B) IL阵列观测到2008年11月24日鄂霍次克海Mw7.3地震事件波形的平滑包络线。黑线为19个组成台站的地震事件的平滑包络线，红线为叠加道。 （C) AS、IL和YK台阵 记录到的11个地震事件叠加平滑包络线。每一道都代表了某一台阵记录到的一个地震事件的叠加数据波形，相当于(B)中的红线。这三条虚线分别是P’•660•P’(红色)、P3KP(橙色)和P’ surfp ‘(蓝色)的预测到达。圆点状的绿色lin表示P’•410•P’的到来。
  台阵数据的一个优点是能够以波慢度的形式分辨地震信号的传播方向信息(28)，这对识别非对称反向散射很有价值(24,25)。利用大孔径(~ 20km) YK和AS台阵研究了P’•660•P’的慢度。较小孔径的台阵(如10公里的IL阵列)可能分辨率不够。（P‘.660.P’震相）信号到达时方位角约为偏离大圆弧±120度，慢度>为1.9 s/度(图S3)。由于PKP波的非对称背散射，这些结果与在震中距约30度距离(24,25)处P’•660•P’震相的慢度一致。

  这种类型的非对称背向散射可以由660公里边界附近的非均质体或660公里不连续面地形产生(图S4)。观测到的P’•660•P’到达时间将非均质体限制在660公里不连续点附近的一个薄层中（如果是非均质体产生的，那么非均质体只能为一个薄层）。较深或较浅的散射体都将导致其传播时间的实质性提前或延迟，这与我们的观察结果相反(图S5)。先前的研究(32)没有观察到P’P’的强前驱波或尾波(前驱波或尾波能量强的话就意味这可能是非均质体的原因产生非对称背向散射的)，因此粗糙的地形似乎是更有可能的候选者。



(from 32-2003JGR_Xu_Precursors.P’P’)

  660公里不连续面的地形变化受P’•660•P’振幅的限制。我们使用幂律方程$P=C_{2D}{k}^{-3}$（来表示这种地形变化模型，也就是分形模型），其中$C_{2D}$表示二维情况下表示660km界面功率谱密度的起伏（粗糙程度），以及射线理论来模拟合成非对称散射（震相）(33)。为了减少不确定性，我们使用P’surfp’作为参考相位来研究P’•660•P’，因为它们的射线路径非常接近。
  作为有记录以来的第二大深震源事件(34)，1994年6月9日的矩级(Mw) 8.2玻利维亚深震源地震提供了一个机会，通过模拟P’•660•P’震相与P’ surfp ‘震相的振幅比值来估计660公里的不连续面的地形(图3)。玻利维亚和秘鲁的16个宽带台站记录了强烈的P’•660•P’震相和更强的P’ surfp ‘震相，其角度震中距从2°到10°(图3C，黑线)。为了得到P’ 660•P’震相与P’ surfp ‘震相的理论比值，使用初步参考地球模型(PREM)(27)合成了P’ 660•P’(图S4A)、P’ 410•P’和P’surfp ‘（三个震相）的包络函数。将远震观测台站SJG (San Juan)的记录到的前50秒直达P波的平滑和平方包络作为经验源时间函数(图S6)，并与合成的平方包络函数进行卷积。利用2000 - 2100 s时窗内平方包络的平均值(图3C，黑色虚线)，来数据中估算出背景噪声的能量级别，然后计算合成包络线中。

图3. 地图、合成和观察到的Mw 8.2玻利维亚地震在近podal距离的P’.660.P’震相。(A) ( C)图中所使用的地震(红星)和地震台站(绿色三角形)位置图。 (B) 这次地震(红星表示，Mw 8.2 Bolivia earthquake.)P’•660•P’震相的采样区域。(灰色区域)。灰色区域从141.4°的PKP焦散距离在660 km的不连续面上141.4度开始，随着距离的增加，由于PKP的振幅衰减，灰色区域逐渐消失。( C) 观察和合成的高频(1.5至2.5 Hz) P’ 660•P’、P’ 410•P’和P’surfp’震相的平滑包络。用归一化振幅绘制速度地震观测图(对应黑线)和合成地震图(彩色线)的平滑包络线。黑色虚线表示地震图的零基线。对于660公里的不连续地形模型，所有合成观测的$C_{2D}$为100米。P 'SurfP由在地壳顶部10公里指数相关函数(自相关长度L = 13公里和扰动均方根= 5%，指数型模型)表示的非均质体(24)和自由表面地形(幂律PSD $C_{2D}$ = 0.3米，分形模型)生成。以P ‘SurfP作为参考震相,$C_{2D}$= 100合成P•660.P’似乎与观察很匹配，而P’•410•P’震相在观察波形中是看不见的或很弱的（合成波形中有而实际观测中没有，就说明是没有410地形变化的）。

  $C_{2D}$ = 100 m合成的P’•660•P’震相(图3C，红线)的起始点（起跳点）和形状都很好地匹配(观测)数据，我们将该值作为对该采样区域的660km界面地形（起伏）的合理估计值。（1.如果660km用分形模型表示）我们只是试图估计$C_{2D}$的数量级，而不是精确的测量(图S7)，因为对地球结构的约束很差，涉及到大量的不确定性(33,35)。这个$C_{2D}=100m$的数量级比自由表面的全球平均$C_{2}D$ ($C_{2D}=0.3m$)大得多(图1A)。（2.如果660km用指数型模型表示）我们还利用非均质体的替代模型（指数型模型）计算出了P’•660•P’，（结果表明）如果存在非均质层，那么（非均质层）必须很薄(<250 km)(图S8)。观察结果中 P’ surfp’ 的持续时间比合成波形中的（P’ surfp）持续时间长。这可能是由于观察波形中存在P3KP震相，但我们没有在合成波形中考虑到这一点（没有考虑可能存在P3KP震相）。
（660km界面用分形模型表示）

（660km界面用指数型模型表示）

  全球性的410公里不连续面橄榄石-瓦兹利石相变转换带有关。与P’·660·P’震相类似，如果410公里界面尖锐且在小的横向尺度上是粗糙的，那么我们预计会看到非对称散射波P’ 410·P’震相。然而，我们未能从单个的地震图中分辨出P’•410•P’震相。为了增强其可观测性，我们叠加所有的平滑包络，但仍然无法识别P’•410•P’。这说明410公里的小尺度地形比660公里的小尺度地形要弱。对这一观测结果的一种解释是，在410公里界面的陡坡面上，小尺度地形非常弱(图3C)。这一解释表明，410公里的不连续面是一个纯粹的相变边界，只有主要是由平稳变化的热结构引起的大尺度和中间尺度的地形变化，。另一种解释是，410公里不连续面(22、36、37)较宽的过渡宽度可以大大降低短周期地震波的反射系数，且当过渡宽度接近地震波的波长时，会减小P’•410•P’震相的振幅。然而，由于其他报告支持一个尖锐的而不是宽的410公里的不连续面(8,38)，410公里的边界似乎主要是由于相变产生。
  P’surfp ‘在上地幔和地壳中移动的距离是P’ 660•P’的两倍，因此上地幔衰减模型中的不确定性可能会使我们对(660km界面)地形的估计产生偏差。此外，因为对地壳中小尺度的不均匀性量化得很差(24)，（导致）我们合成的P’surfp’震相有很大的不确定性。上地幔和地壳结构影响我们对小尺度的660公里不连续面地形的估计。因此，我们使用另一个震相$PKiKP$，它从内核边界(ICB)反射，(将$PKiKP$)作为一个参考震相来估计660公里的不连续面的地形。$PKiKP$常在高频下观测到，其射线路径与$P^{\prime }•660•P^{\prime }$相似。尽管$P^{\prime }•660•P^{\prime }$在下地幔中旅行了4次，而$PKiKP$只这样(下地幔中旅行)做了两次，但下地幔的衰减要比上地幔低得多(27,39)。利用$PKiKP$作为参考震相，用PREM来描述$P^{\prime }•660•P^{\prime }$和$PKiKP$的衰减，通过比较合成和实际观测波形的包络，估算出660公里不连续面的地形。采用高频(3 ~ 4 Hz)滤波后，IL台阵记录到 鄂霍次克海地震（Sea of Okhotsk earthquake） 事件的清晰的$PKiKP$震相，且$P^{\prime }•660•P^{\prime }$的信噪比大于4(图4A和表S1)。与玻利维亚地震(Bolivia earthquake)(图3C)不同，P’ surf P’甚至比P’•660•P’还要弱，这可能是由于上地幔的衰减所致。我们比较 测量（实际观测波形） 和 合成（合成波形） 的$P^{\prime }•660•P^{\prime }$与$PKiKP$震相振幅比(图4B)，并估计$C_{2D}$超过5000m，这比全球自由表面平均地形要大几个数量级(自由表面$C_{2D}$=0.3m)。ICB反射系数或 (CMB的?)地形 的不确定性等因素可能会导致$PKiKP$震相(40)的幅值波动，从而对660km界面地形变化的估计产生偏差。然而，IL台阵记录到的 其他三个地震事件(图4B) 测量的$P^{\prime }•660•P^{\prime }/PKiKP$比值表明有类似的结果$C_{2D}>1000m$。这个$C_{2D}$（结果）比根据 玻利维亚地震Bolivia earthquake (Fig. 3C) 估计的$C_{2D}=100m$高得多(图3C)。这种巨大的差异可能是由于 $C_{2D}$估算的巨大不确定性 和/或660公里界面地形小尺度的的地理差异(图S9)。尽管如此，我们在这两种情况下都观察到了660公里不连续面的强烈地形。

图4. $P^{\prime }•660•P^{\prime }$、$PKiKP$及它们的振幅比。(A) IL台阵记录的2012年8月14日鄂霍次克海地震的高频(3至4 Hz) PKiKP(上)和P’•660•P’(下)速度地震图的平滑包络。红色的道表示叠加的结果。(B)四次事件的$P^{\prime }•660•P^{\prime }/PKiKP$振幅比。利用射线理论计算了$P^{\prime }•660•P^{\prime }$(图S4)和$PKiKP$的包络函数。红线代表根据四个不同的$C_{2D}$系数(10、100、1000和5000)地形模型预测$P^{\prime }•660•P^{\prime }/PKiKP$振幅比。实测$P^{\prime }•660•P^{\prime }/PKiKP$振幅比(黑色圆圈)高于红线$C_{2D}$= 1000m,分布在$C_{2D}$= 5000附近。

  本研究中$P^{\prime }•660•P^{\prime }$的采样区域不仅包括目前存在的（正在发生的？）深俯冲带(如图3B，日本)，还包括未发现重大俯冲的区域(图S10和S11)。在前者环境中(正在发生俯冲的区域)，停滞在660-km界面处的板片(6)会积累大量的化学非均质（体），从而导致660-km的小尺度地形和体积非均质性。在后一种情况下(未发现大的俯冲)，来自远古板块的累积海洋地壳可能在660公里界面上浮起，形成一个充满化学异质性的石榴子岩层(41)。


  660公里边界的小尺度地形或不太可能的薄层的非均匀体，最好可以用化学原因来解释。由于当前（正在进行的）俯冲和/或古代板块形成的海洋地壳积累，在660公里的界面上方形成一个重力稳定的石榴子岩层是可能的。这一假设由Ringwood(42)提出，Irifune和Tsuchiya(41)进一步讨论。一些地区缺乏660公里界面的小尺度地形，这意味着一个全球不连续的化学层。我们的观察支持了将俯冲板块描述为过渡区的瞬态特征的模拟，而过渡区最终会渗透到下地幔中。它们还支持了部分被阻断的上地幔到下地幔环流(43)，它有效地在层状和全地幔对流之间交替。这种类型的模型比层状或全地幔端部的情况更为复杂，但在地球历史上增加物质交换的时间依赖程度，可能有助于统一地幔的地球化学、地球动力学、地震学和岩石学观测。