利用双差层析成像方法反演阿拉斯加地区岩石圈速度结构

周少贤, 薛梅

周少贤,薛梅. 2022. 利用双差层析成像方法反演阿拉斯加地区岩石圈速度结构. 地震学报,44(3):374−387. DOI: 10.11939/jass.20210122
引用本文: 周少贤,薛梅. 2022. 利用双差层析成像方法反演阿拉斯加地区岩石圈速度结构. 地震学报,44(3):374−387. DOI: 10.11939/jass.20210122
Zhou S X,Xue M. 2022. Lithospheric velocity structure of Alaska revealed by double difference tomography. Acta Seismologica Sinica44(3):374−387. DOI: 10.11939/jass.20210122
Citation: Zhou S X,Xue M. 2022. Lithospheric velocity structure of Alaska revealed by double difference tomography. Acta Seismologica Sinica44(3):374−387. DOI: 10.11939/jass.20210122

利用双差层析成像方法反演阿拉斯加地区岩石圈速度结构

基金项目: 国家自然科学基金面上项目(42076064)和上海佘山地球物理国家野外科学观测研究站开放基金(2020K03)共同资助
详细信息
    作者简介:

    周少贤,在读硕士研究生,主要从事双差层析成像和横波分裂方面的研究,e-mail:zhoushaoxian@tongji.edu.cn

    通讯作者:

    薛梅,博士,副教授,主要从事西太平洋边缘海内部构造和地震背景噪声机制研究,e-mail:meixue@tongji.edu.cn

  • 中图分类号: P315.2

Lithospheric velocity structure of Alaska revealed by double difference tomography

  • 摘要: 阿拉斯加地区由不同地质时期的地体向北增生而成,经历了漫长的构造演化,地质构造复杂。ANF (Array Network Facility)网站新近提供的来自USArray地震台网记录的地震台站观测数据填补了阿拉斯加地区西部和北部的观测空白,本文选取该区域中345个台站记录的5 638个地震事件的P波、S波到时数据,采用区域双差地震层析成像方法反演得到了该地区的岩石圈三维P波速度模型和地震重定位结果。研究结果显示:阿拉斯加西部太平洋板块的俯冲倾角较大,深部地幔楔表现为P波低速异常,推测由俯冲板块顶部脱水产生的流体释放到地幔楔并触发部分熔融所致,这些熔融物质上升到达地表形成阿留申火山岛链;中部亚库塔特(Yakutat)地体与太平洋板块发生耦合,俯冲倾角减小,一方面使地壳压应力增加,引起地壳增厚和楚加奇(Chugach)山脉隆升,另一方面导致该处地幔楔降温从而使产生的熔体减少,并随着地壳压应力的增加部分地壳裂隙闭合,阻断了熔体上升至地表,从而形成迪那利(Denali)火山空区;亚库塔特地体与东部兰格尔(Wrangell)火山区之间存在较明显的分界,兰格尔火山区下方的低速区(与岩浆活动对应)集中于西北侧,火山区的岩浆来源可能与环形地幔流沿太平洋—亚库塔特板块边缘的上升流相关。这些结果表明,阿拉斯加地区深部复杂的地球动力学过程导致了其地表复杂的地质构造。
    Abstract: Alaska region is formed by the northward accretions of terranes from different geological periods and has experienced extensive internal deformation and metamorphism, and the geological structure is complex. The Array Network Facility (ANF) website recently provides new observational data from the seismic network of USAarray, filling the observation gap in the west and north of Alaska. Based on P and S wave arrivals of 5 638 events recorded by 345 stations from ANF, this study relocates earthquakes and images the 3D lithospheric P-wave velocity structure beneath Alaska simultaneously by regional double difference tomography. The results reveal larger dip angle of subducting Pacific Plate and low-vP anomalies in the mantle wedge beneath western Alaska. These observations reflect the subduction process, during which the dehydration of the subduting Pacific Plate releases fluids into the mantle wedge, triggers partial melting, and generates melts, which was then transported to the surface by the upwelling flow so as to form Aleutian volcanic arc. In central Alaska, the coupling between the Yakutat terrane and the Pacific Plate reduces the subduction dip. On the one hand, the shallow subduction of the Yakutat terrane increases the compressive stress of the crust, causing the crustal thickening and uplifting of Chugach mountains. On the other hand, it cools the mantle wedge reducing magma generation, which are then combined with the closure of crust fractures resulted from the increase of the crustal stress, blocking the supply of melt to the surface, and finally leading to the formation of Denali volcanic gap. In addition, there is a clear boundary between Yakutat terrane and Wrangell volcanic field in the east, and the low velocity zone corresponding to magmatic activity in the region is concentrated in the northwest. The magma source may be related to upwelling of the toroidal mantle flow around the Pacific-Yakutat slab edge. These results suggest that the complex geodynamic processes in deep Alaska lead to the complex geological structure on the surface.
  • 地震是由沿断层的应力迅速转化为非弹性应变(岩石破坏)而引起的,断层的破坏可以表现为小地震或大地震、抗震滑移或火山地震。实验室研究表明,由于应力会引起裂缝性质的变化,故地震对应力场很敏感。获取地下介质的动态变化一直是热门的话题,越来越多的研究表明,与地震和火山爆发等构造事件相关的地震波其速度会随时间发生变化(Patanè et al,2006Lei et al,2011Obermann et al,2014Koulakov et al,2016)。目前,研究地震波速随时间变化的方法较多,且应用广泛。其中,地震层析成像是可以提供高精度空间分辨率的方法之一,该方法通常仅分析一组地震波到时数据并建立一个速度模型,假设速度结构不随时间变化,使用此方法可需独立地反演不同时间段的数据集,并假定得到的模型差异代表了真实的速度变化。然而,Julian等(2008)Julian和Foulger (2010)的研究表明此方法得出的变化可能不真实,由于地震事件和台站分布的不确定性导致了不同时段的地震射线分布对速度模型采样不同,为消除震源位置不确定性和仪器误差所产生的非真实变化,可以使用复杂的迭代方法,将所得到的一个时段的速度模型当作初始模型,然后同时反演两个时段的数据,依此可以确定真正的速度变化。在此基础上,Qian等(2018)基于双差层析成像方法(Zhang,Thurber,2003)提出了一种新的时移层析成像方法,通过最小化模型的差异和不同时段的到时残差来同时反演多个数据集,使用不同时段的事件对构建的到时差来克服由部分观测误差和不同时段地震射线分布差异造成的影响,因此该方法可以确定高精度空间分辨率的地震波速随时间的变化。Qian等(2018)将该方法应用于美国西部地下煤矿开采诱发微震的数据并得到了高精度的与开采过程相关的速度变化。通常在中强地震发生前后,震源区及其附近的地震分布和地震定位所使用的台站差异很大,所以使用此方法还可以消除以上差异所造成的非真实变化并得到可靠的结果。

    云南盈江地处云南省西部的中缅交界地区,位于以顺时针旋转的喜马拉雅东构造的东南端,其东邻腾冲火山区。该地区地质构造复杂(图1),地震活动频繁(肖卓,高原,2016),特别是2008年汶川MS8.0地震以来区域内地震活动更加频繁,这些地震活动可能与印度板块深俯冲所形成的“大地幔楔”结构相关(Lei et al, 20092019;Lei,Zhao ,2016)。该地区台站分布稀疏,台站间空隙角较大,且仅有一侧有台站分布,因此对其地震数据进行分析存在一定局限。此外,该地区的历史地震记录不够完备,虽然云南省地震台网的建立使该地区的地震监控能力得以提升,但台站对地震的包围性较差的问题仍然存在,致使该地区地下介质特征的相关研究较少。

    图  1  研究所使用的台站分布
    Figure  1.  The distribution of stations used in the study

    2014年5月24日4时49分在云南省德宏傣族景颇族自治州盈江县发生MS5.6地震,震中位置(25.0ºN,97.8ºE),震源深度为12 km,5月30日9时20分再次发生MS6.1地震,震中位置(25.0ºN,97.8ºE),震源深度为12 km。两次地震震中距约为5 km,其中MS6.1地震震中更靠近卡场—大竹寨断裂,且该地震的极震区烈度为Ⅷ度,共造成45人受伤,多处房屋倒塌倾斜,经济损失十分严重(张彦琪等,2017)。此次双震产生的大量余震数据,为该地区的相关研究提供了基础,许力生等(2014)运用经验格林函数技术提取了两次地震的震源时间函数,获得了破裂过程,上述研究仍缺乏针对盈江地区地下介质特性的研究成果;徐甫坤等(2015)杨婷等(2016)对余震序列进行了重定位;肖卓和高原(2016)利用尾波干涉技术分析了地震序列持续期间地壳介质的物性变化;孙楠等(2017)对近场流动台站记录的大量波形数据进行横波分裂研究;秦敏等(2018)分析了余震序列应力降的时空演变过程。鉴于此,本文拟基于云南数字地震台网2009—2019年间的地震到时记录,利用基于双差层析成像的时移层析成像方法,拟得到2014年云南盈江两次地震前后震源区地下P波速度的时空变化,并结合已有研究结果探讨引起该变化的原因,以期了解此次双震的孕育和发生过程。

    根据杨婷等(2016)和地震应急现场科考结果可知,2014年盈江两次地震的发生与卡场—大竹寨断裂相关,但震源区更靠近苏典断裂。而2008年8月21日和2008年9月3日苏典—盈江断裂上分别发生了MS6.0和MS5.0地震,为了消除这两次中强地震及其余震对本研究的影响,将这两次中强震的余震结束时间作为研究的起始时间来进行数据搜集,最终选取了2009—2019年云南地震台网记录的发生在卡场—大竹寨断裂及苏典断裂附近区域的地震观测报告数据(图2)。盈江MS5.6地震后,云南省地震局在震中附近架设了3个流动台(卡场台、勐弄台和昔马台),盈江MS6.1地震后又新增了2个流动台(新城台和盏西台),震源区的台站密度和地震定位的精确度得到了提升,然而盈江地震前震中附近的台站分布不够密集且尚未对地震形成包围,对地震定位有一定的影响。最终,本文选取了22个台站(图1)记录的6 124次地震事件(图2),为了保证反演结果的可靠性,每个事件需被震中距小于200 km的5个以上台站记录到。同时使用初始一维速度模型(表1)绘制了理论走时曲线和P波震相走时曲线,如图3所示。可以看出在所选取的震中距范围内P波走时呈现出明确的线性趋势,震相的离散度小,表明原始震相观测数据具有较高的可靠性。

    表  1  初始一维速度模型
    Table  1.  The initial 1D velocity model
    深度/km速度/(km·s−1
    05.000
    25.367
    55.803
    76.033
    106.126
    146.400
    186.600
    256.800
    307.200
    388.000
    下载: 导出CSV 
    | 显示表格
    图  2  研究区内地震分布及网格划分
    Figure  2.  Seismic stations used in this study and grid nodes division in the studied area
    图  3  P波走时曲线
    Figure  3.  P-wave travel times

    基于双差层析成像的时移层析成像方法可以确定真正的速度变化,即同时反演多个数据集,最小化不同时段速度模型之间的差异。在双差层析成像方法中,用来构建事件对的事件均来自同一时段,为了求解两个时段间的速度变化,需选取两个不同时段的事件来构建事件对,将双差层析成像方法(Zhang,Thurber,20032006)进行转变。对于同一观测台站k,分属于两个时间段的事件i和事件j的到时差之差为

    $$ r_k^i {\text{-}} r_k^j {\text{=}} \sum\limits_{l {\text{=}} 1}^3 {\frac{{{{{\text{∂}}}} T_k^i}}{{{\text{∂}} x_l^i}}} \Delta x_l^i {\text{+}} \Delta {\tau ^i} {\text{+}} \int\nolimits_i^k {\delta {u_1}} {\rm{d}}s {\text{-}} \Bigg(\sum\limits_{l {\text{=}} 1}^3 {\frac{{{\text{∂}} T_k^j}}{{{\text{∂}} x_l^i}}} \Delta x_l^i {\text{+}} \Delta {\tau ^j} {\text{+}} \int\nolimits_j^k {\delta {u_2}{\rm{d}}s} \Bigg){\text{,}} $$ (1)

    式中,δu1δu2为时段1和时段2的慢度参数;ij为时段1和时段2中构建事件对的事件;Tk为发震时刻;∆τ为事件的发震时刻扰动量;xl为三分量的位置扰动,其中l=1,2,3。

    首先,利用双差层析成像方法得到时段1的三维速度模型;然后,将其作为初始模型解决采样问题;最后,使用基于式(1)的时移层析成像方法得到时段1和时段2之间的速度变化。在此过程中仅使用由分属于两个时段的事件所构建的事件对的到时差数据,事件对的构建要求与双差层析成像方法一致,需由一定距离范围内的事件所构成,并至少被一定数量的相同台站记录到。保持时段1的慢度u1固定,到时差数据用来升级时段2的慢度u2。通过同时反演两个时段的数据,可以确定真实的变化,并能提供高空间分辨率的结果。

    为了获得两次盈江地震发生前后震源区的P波速度变化,得到两次中强震的发生所造成的影响,将收集到的数据划分为4个时段(表2)进行研究,其中t1为前震期,t2MS5.6地震发生后至MS6.1地震发生前,t3MS6.1地震发生后两个月内,t4MS6.1地震后三个月至2019年12月31日。盈江地震发生前研究区监控能力较弱,地震记录数据较少,并且各时段的时间长度不同,致使四个时段的地震次数不同。虽然4个时间段的时间窗尺度不一,但本文所使用的基于双差层析成像的时移层析成像方法可以消除由不同时间段的射线分布差异带来的误差。

    表  2  盈江两次中强震四个时段的数据分布
    Table  2.  Data of the four periods for two moderate strong earthquakes in Yingjiang
    时间段时间窗台站个数地震次数绝对到时/条相对到时/条
    t12009−01−01—2014−05−2312679431064208
    t22014−05−24—2014−05−2914947605467628
    t32014−05−30—2014−07−3116398235139160560
    t42014−08−01—2019−12−31211007851778155
    下载: 导出CSV 
    | 显示表格

    首先选取盈江MS6.1主震的震中位置为(25.0ºN,98.7ºE)作为坐标原点来建立坐标系(图2),并综合考虑地震和台站的分布来划分网格。由于地震和台站分布不均匀,在划分网格时应考虑此种不均匀性,同时测试不同的网格间隔分布,最后得出横向网格间隔在震源区为2 km×2 km,外围区域为4 km×4 km较为合理(图2),垂向网格分别位于0,2,5,7,10,14,18,25,30,38 km深度处。

    初始一维速度模型参考了多个地壳模型(房立华等,2011肖卓,高原,2016杨婷等,2016陈飞,2017)后建立,为了获得更加稳定的结果,三维初始速度模型相对于一维速度模型的精度更高,因此本文首先使用初始一维速度模型(表1)且对所有的地震数据进行双差层析成像反演得到三维速度结构和地震重定位结果,然后根据基于双差层析成像的时移层析成像方法的技术路线,将所得到的三维速度模型作为初始速度模型分别对4个时段的地震数据(表2)进行反演,从而得到这4个时段的较为稳定的P波速度结构。

    首先将所有地震事件进行配对,地震对之间的最大距离为30 km,每次地震最多可与20次地震组成地震对来构建相对到时数据。为了保证反演结果的稳定性,对不同的光滑权重和阻尼因子进行了权衡分析,最后确定光滑权重为30,阻尼因子为150。经过20次迭代后,到时残差的均方根由0.387降低到0.037,最终获得6028次地震的震源参数,重定位后的MS5.6和MS6.1地震的震源深度分别为10.22 km和8.6 km,同时得到了研究区的P波速度结构。并对反演结果进行了棋盘测试,结果如图4所示。可以看出,地震分布较为集中的盈江地震震源区的分辨率恢复得较好,由于地震监控能力不同,震源区以西即边境线以外分辨率较差。在5 km和10 km深度处余震分布较多,主要分布于低速异常区内,而在14 km深度处余震分布较少,主要分布在高速异常内,可能与不同深度处的岩石类型有关,在10 km和14 km深度处的震源区及其附近存在相邻的条带状的高速异常区和低速异常区,而盈江地区周围的主要断层(F1F2F3)无明显的P波速度分布特征。由于盈江地区位置偏远,地震监控能力不足,对该地区地下介质性质的研究很少,已有的研究结果也不够全面。黄璟(2012)和杨婷等(2014)分别使用不同的近震走时层析成像方法对滇西地区的三维速度结构做了研究,虽然两文中均未专门对盈江地区做讨论,但根据其结果,在10 km深度处,盈江地区主要表现为低速区,与本文的结果相符合。李飞等(2011)利用Pn波数据反演了川滇地区的Pn波速度结构和各向异性,结果表明盈江地下的Pn波速度偏低,基本也与本文的结果相一致。由已有的滇西地区电性结构、速度结构研究(孙洁等,1989Wang et al,2003)也可知,滇西地区地壳内普遍存在低阻、低速区。

    图  4  不同深度处的P波速度结构(a)和棋盘测试结果(b)
    Figure  4.  Distribution of P wave velocity structure (a) and checkerboard resolution test (b) at different depths

    采用2.1节所得到的三维速度结构作为初始速度模型,重定位结果作为初始地震定位,分别对4个时段进行双差层析成像得到相应的P波速度结构。首先根据4个时段的地震分布来进行地震匹配,其中t1t4时段的地震较为分散,采用相同的匹配参数,地震对之间的最大距离为30 km,每次地震最多可与20次地震组成地震对。t2t3时段的地震集中分布在震源区,采用另一组匹配参数,即地震对之间的最大距离为20 km,每次地震最多可与10次地震组成地震对,最终构建了4个时间段的相对到时数据(表2)。由图5a中可看出盈江地震震源区4个时间段的射线密度均较密集,基本能覆盖震源区。反演过程中,为了得到更加稳定的结果,对4个时段均进行了光滑权重及阻尼参数权衡分析,并建立解的方差与数据方差(图5b)和模型方差与数据方差(图5c)均衡曲线,分别得到最优阻尼因子和最优光滑权重,4个时段的阻尼因子和光滑权重的选取均不相同。从图5c中可看出,由于数据集不同,4个时段的曲线变化趋势不尽相同,但都存在拐点,可得到相应的最优光滑权重。经过18次迭代后使结果收敛,反演前后4个时间段的到时残差均方根均显著降低(表3),说明最终模型的稳定性有了很大的提高。通过对比初始速度模型及最后得到的三维速度模型的均方根残差变化(图5d)可以看出,最后的三维速度模型相对于初始速度模型更加集中且对称地分布于均方根拟合差零值的两侧,表明最终的三维速度模型能够更好地对数据进行拟合。

    图  5  四时段的参数信息
    (a) 二维射线分布;(b) 不同阻尼参数的解的方差和数据方差的均衡曲线;(c) 不同平滑权重参数的模型方差和数据方差的均衡曲线;(d) 不同模型下的观测走时与理论走时的均方根拟合差变化
    Figure  5.  The parameter information of four periods
    (a) Distribution of 2-D P wave ray paths;(b) The trade-off curves of solution variance and data variance for different damping parameters;(c) The trade-off curve of slowness model variance and data variance for a set of smoothing weight parameters;(d) The root-mean-square misfit improvement between observation travel times and theoretical ones based on 3D models
    表  3  4个时段反演前后到时差的均方根残差变化
    Table  3.  The RMS residuals between observed and predicted differential travel times after and before tomography for the four periods
    时段到时差的均方根残差/s下降百分比时段到时差的均方根残差/s下降百分比
    初始三维速度模型最终三维速度模型初始三维速度模型最终三维速度模型
    t1 0.519 0.203 60.9% t3 0.318 0.052 83.6%
    t2 0.317 0.048 85.8% t4 0.404 0.117 71%
    下载: 导出CSV 
    | 显示表格

    本文同时对4个时段的反演结果分辨率进行了棋盘测试,分别选取了5,10,14 km深度处的分辨率测试结果进行展示。由图6可见,对于所有时间段,地震分布比较集中的盈江地震震源区其P波速度恢复结果明显好于其它区域,因为5 km和10 km深度处的地震分布较多,震源区的P波速度恢复结果好于14 km深度处。并且地震集中分布的t2t3t4这3个时段其分辨率恢复较地震分布分散的t1时段好,在t1时段盈江地震震源区也能恢复出部分分辨率。图7显示了整个研究区在不同时段内的不同深度层面的P波速度分布,为了更加直观的展示可信度较高的区域而加入了衍生加权总和(derivative weight sum,缩写为DWS)值分布,该值表示给定网格节点附近所有射线长度(乘以数据权重)的总和。DWS=150 (图7中的红色等值线)基本可以包围棋盘测试恢复较好的盈江地震震源区,等值线分布基本能与棋盘测试结果相对应,表明本文所得到的盈江地震震源区的P波速度结构具有一定的可信度。由于此结果将作为下一步时移层析成像的初始模型,所以此处不再赘述。

    图  6  四时段t1-t4 (a-d)不同深度处的P波速度棋盘测试结果
    Figure  6.  Checkerboard resolution test for the four periods t1-t4 (a-d)
    图  7  四个时间段t1t4 (a−d)不同深度处的P波速度分布
    Figure  7.  The distribution of P wave velocity for the four periods t1t4 (a−d)

    根据时移层析成像方法的技术路线,将2.2节中得到的两个相邻时段中的前一个时段的P波速度结构作为初始模型,两个时段的地震数据进行配对得到的相对到时数据进行升级后,得到一个时段相对于前一个时段的P波速度变化。地震配对时,要求两个地震事件分属于两个时段,由于地震配对会受地震和台站分布的影响,为了构建尽可能多但又不影响结果稳定性的相对到时数据,要求地震对之间的最大距离为50 km,一个时段内的每次事件最多可以与30次事件相连接,最后得到的相对到时数据如表4所示。t1t2时段所使用的台站和地震分布差异较大,配对后所得到的相对到时数据较少,将影响最后结果精度。3个时移层析成像经过30次迭代后均能很好地收敛。经过对比初始三维速度模型和最后的三维速度变化的均方根残差变化(表5)可以看出,最终得到的P波速度变化与所使用的地震数据较为适宜,精度显著提高。同双差层析成像方法一样,使用棋盘测试来分析结果的分辨率,具体如图8所示。对于t1t2时段,射线分布差异大,虽然已选取尽量多的地震对,其分辨率仍然恢复不佳,再加上盈江地区地震监控能力不均,均会影响反演结果的精度;而其余时段分辨率恢复相对较好,受地震数据数量及台站分布的影响,使得棋盘测试结果总体表现不佳。但在一定程度上依然恢复出了盈江地震震源区的分辨率,就此本文将重点讨论盈江地震震源区的P波速度变化特征。

    表  4  时移层析成像所用到的数据
    Table  4.  The data of time-lapse tomography
    时间段地震次数台站/个P波到时差/个
    t2t116261410364
    t3t249291641193
    t4t349892143656
    下载: 导出CSV 
    | 显示表格
    表  5  相邻两个时间段反演前后到时差的均方根残差变化
    Table  5.  The RMS residuals between observed and predicted differential travel times after and before tomography varies for two adjacent time periods
    时间段到时差的均方根残差/s下降百分比
    初始三维模型三维速度变化
    t2t10.2790.09765.2%
    t3t20.1400.01688.6%
    t4t30.1800.01591.7%
    下载: 导出CSV 
    | 显示表格
    图  8  相邻两时间段的P波速度变化棋盘在不同深度的恢复图
    Figure  8.  Recovered checkerboard result of P wave velocity for the two adjacent time periods at different depths (a) t2t1;(b) t3t2;(c) t4t3

    图9主要显示了5,10和14 km深度处的P波速度变化,为突出显示,在图中将DWS≤500的区域阴影化。虽然在一些深度处DWS>500的区域与棋盘测试结果并无精确对应,但基本均包含了盈江地震震源区。从图9a可看出,t2时段相对于t1时段,每个深度处震源区的P波速度均呈下降趋势,且下降幅度较小,其中5 km深度处MS5.6地震震中区下降了约0.5%,10 km和14 km深度处下降的区域延伸至该震中区以南区域,分别下降了0.7%和0.5%左右。图9b显示,t3时段相对于t2时段,5 km深度处的余震分布较多区域的P波速度下降幅度增大,约为1%,而10 km和14 km深度处变化不大,分别为0.6%和0.3%左右。从图9c可以看出,震后至2019年底整个震源区在愈合,无论是浅层还是深层都是P波速度上升,但是上升幅度均未达到之前下降的幅度,其中浅层5 km的P波速度上升幅度更明显,约为0.4%,深层10 km和14 km处均为轻微上升。

    图  9  两个相邻时间段的P波速度变化在不同深度的分布图
    Figure  9.  Distribution of temporal P-wave velocity changes for the two adjacent time periods at different depths
    (a) t2t1;(b) t3t2;(c) t4t3

    综上所述,在MS5.6地震发生后,震源区的P波速度下降,下降幅度较小,区域较小。而在MS6.1地震发生后,震源区的P波速度下降至最大,且对浅层深度的影响较大。在MS6.1地震发生后至2019年底震源区愈合,但愈合程度未达到地震前的水平。

    沿着余震分布优势方向,选取两条剖面位置见图2,其中AA′剖面与卡场—大竹寨断裂平行,BB′剖面与AA′剖面相交且穿过卡场—大竹寨断裂。图10显示了沿着这两条剖面的P波速度变化分布,同时将两个相邻时间段的后一个时间段的地震事件投影到剖面上。从图10a可以看出,在剖面上投影的是MS5.6地震的余震,t2时段相对于t1时段,MS5.6地震上方的浅层深度的P波速度轻微下降,余震主要分布在AA′剖面上MS5.6地震以北区域,靠近A′端的、与MS5.6地震大体相同深度处P波速度也是下降。BB′剖面的卡场—大竹寨断裂以东区域(B′端),不仅余震分布较多,浅层深度处的P波速度也轻微下降。图10b中两个剖面上投影的是MS6.1地震发生后两个月内的余震,在MS6.1地震发生后,整个AA′剖面基本都有余震分布,BB′剖面的卡场—大竹寨断裂以东区域(B′端)的余震分布也很密集,且余震分布深度加深。在两个剖面上,t3时间段相对于t2时间段的P波速度在余震主要分布的深度范围内(<15 km)基本表现为下降,下降幅度大于图10a,其中在AA′剖面上可观察到,在MS6.1地震发生后,原来没有余震分布的MS5.6地震以北区域(A端)不仅余震增加,P波速度下降更为显著,而余震分布密集的BB′剖面的B′端的P波速度下降也很明显。以上说明P波速度变化与余震的发展有关。图10c中剖面上投影的是MS6.1地震发生后至2019年底的地震事件,基本可看到整个区域内所有深度处的P波速度均为上升,说明整个区域都在愈合,只是愈合的程度不及之前地震发生所产生的影响。

    图  10  两个相邻时间段的P波速度变化沿垂直剖面AA′(左)和BB′(右)的分布图
    Figure  10.  Temporal distribution P-wave velocity changes along the cross sections AA′ (left), BB′ (right) for the two adjacent time periods
    (a) t2t1;(b) t3t2;(c) t4t3

    对于P波速度变化的原因,需结合其它的地球物理研究结果进行分析讨论。本研究表明震源区P波速度的下降与余震的发生相关,在5月24日MS5.6地震发生后至5月30日MS6.1地震发生前这6天内,余震主要分布于卡场—大竹寨断裂以东以及MS5.6地震震中以南区域,接近地表的余震数量较少,余震大多分布在5-10 km深度内。本文的重定位结果显示MS5.6地震的震源深度为10.22 km,表明大多数余震分布比MS5.6地震震中浅,MS5.6地震对浅层区域的影响较大,所以P波速度在盈江地震震源区的浅层区域呈轻微下降。5月30日MS6.1地震震中更靠近卡场—大竹寨断裂,重定位结果显示震源深度为8.6 km,比MS5.6地震震源深度要浅。在MS6.1地震发生后两个月内,余震仍主要分布于卡场—大竹寨断裂以东,但该断裂以西区域的余震增多,且延伸的距离更远,余震区范围相较于MS5.6地震更加向北延伸。深度分布方面,虽然MS6.1地震的余震主要分布于2—10 km,但是主震下方仍然分布很多余震。综上,MS6.1地震的发生扩大了横向和纵向上的影响范围,从而P波速度下降的范围也随之扩大,但是在余震的主要分布深度范围内其下降较为显著,并且向北延伸的余震分布区域的P波速度下降幅度更显著,可能与主震的震级有关,MS6.1地震相较于MS5.6地震所产生的能量更大,因此影响范围和程度更大。分析影响地震波速变化的原因,可能是由于地震速度变化与余震分布的时空变化同步相关,余震的动态和静态应力造成震源区介质物理属性的改变而导致地震波速的变化(刘志坤,黄金莉,2010)。

    MS6.1地震发生后,P波速度下降至最低,下降的幅度约不超过1%,这个幅度与肖卓和高原(2016)得出的波速变化约为3‰基本符合,虽然波速变化大小不太一致,但数量级相一致,并且也同样监测出大地震发生后明显的波速下降,表明本文结果可靠。同时肖卓和高原(2016)也监测出MS6.1地震后走时延迟的减小,即出现地震波速的上升现象,这也与本文的观测结果相符合。MS6.1地震发生后至2019年底,整个震源区的P波速度上升,在这个过程中岩石裂纹闭合,受损岩石刚度恢复,但上升的幅度不及之前下降的幅度,然而这也不能说明该区域还在愈合,因为本文的观测结果对应于时间的平均值,震后至2019年底的时间跨度较长,可能还包含其它的物理过程,还需要进一步研究才能得出更加确切的结论。

    盈江地区位于地质构造和应力分布情况复杂的中缅交界地区,一直是地震高发区,地震监控能力分布不均匀,对于监测该地区地下介质特征存在一定的桎梏。2014年盈江两次地震只间隔了6天,短时间内两次中强地震的接连发生,杨婷等(2014)认为是因为盈江地区处于火山区的边缘,地下介质比较破碎,强度较低引起的。但本文并未得出表明介质强度特征的证据,因此还未找到短时间内发生两次地震的原因。未来随着地震监控能力的不断增强,地震数据精度的不断提高以及观测方法的不断升级,在以后的工作中对发震原因的研究将成为重点。

    云南省地震局王光明助理研究员协助绘制了本文提供了画图件,孙楠助理研究员提供了断层数据,审稿专家对本文提出了宝贵的意见,作者在此一并表示感谢。

  • 图  7   不同位置的P波速度垂直剖面图(剖面位置见图6i,白色实线以上代表DWS>1 000的区域)

    Figure  7.   vP models along ten cross sections at different positions (The profile location is shown in Fig. 6i,and the region above the white solid line is the region with DWS>1 000)

    (a) AA′; (b) BB′; (c) CC′; (d) DD′; (e) EE′; (f) FF′; (g) GGG′′; (h) HHH′′; (i) II′; (j) JJ

    图  1   阿拉斯加构造背景图

    两个箭头分别给出了太平洋板块向北美板块俯冲的速度(DeMets,Dixon,1999)和亚库塔特地体向北美板块俯冲的速度(Sauber et al,1997)。F1:廷蒂纳断层;F2:迪那利断层;F3:塔基特纳断层;F4:接触断层;F5:边界山脉断层,下同

    Figure  1.   Tectonic settings of Alaska

    The two arrows represent subduction velocity of Pacific Plate to North America Plate (DeMets,Dixon,1999) and Yakutat Plate to North America Plate (Sauber et al,1997). F1: Tintina fault;F2:Denali fault;F3:Talkeentna fault;F4:Contact fault;F5:Border Range fault,the same below

    图  2   本文所用的地震和台站分布图

    Figure  2.   Distribution of events and stations used in this study

    图  3   筛选前(a)、后(b)的P波和S波时距曲线

    黑线代表拟合的时距曲线,绿线代表观测与拟合的到时差为±10 s,反演中仅保留绿线以内的震相数据

    Figure  3.   Time-distance curves of P and S waves before (a) and after (b) selection

    The black line represents the fitting time-distance curve,and the green line represents the difference of arrival time between observation and fitting is ±10 s,only the seismic phase data within the green lines are retained to ensure the data quality

    图  4   阻尼参数及平滑因子权重均衡曲线

    Figure  4.   Trade-off curves of damping and smoothing weight parameters

    图  5   研究区不同深度水平剖面的vP棋盘测试结果(白色实线内部为DWS>1 000的区域)

    Figure  5.   Checkerboard test results of vP at different depths in the studied area where the white solid line delineates the area with DWS>1000

    图  6   阿拉斯加地区不同深度处P波速度结构

    Ⅰ . 育空盆地;Ⅱ . 塔纳纳盆地;Ⅲ . 铜河盆地;Ⅳ . 库克湾盆地

    Figure  6.   P-wave velocity structure in Alaska at different depths

    Ⅰ . Yukon Flat basin;Ⅱ . Middle-Tanana basin;Ⅲ . Copper River basin;Ⅳ . Cook Inlet basin(a) 10 km;(b) 25 km;(c) 45 km;(d) 60 km;(e) 80 km;(f) 100 km;(g) 120 km;(h) 140 km;(i) 180 km

    表  1   本文所采用的初始一维速度模型(Eberhart-Phillips et al,2006

    Table  1   Initial ID velocity model used in this study (Eberhart-Phillips et al,2006

    深度/kmvP/(km·s−1深度/kmvP/(km·s−1
    0 2.90 48.0 7.87
    2.0 5.10 65.0 8.12
    6.0 6.05 85.0 8.20
    15.0 6.35 110.0 8.24
    24.0 6.60 140.0 8.30
    33.0 7.15 190.0 8.44
    下载: 导出CSV
  • 刘伟,吴庆举,张风雪. 2019. 利用双差层析成像方法反演青藏高原东南缘地壳速度结构[J]. 地震学报,41(2):155–168. doi: 10.11939/jass.20180083

    Liu W,Wu Q J,Zhang F X. 2019. Crustal structure of southeastern Tibetan Plateau inferred from double-difference tomography[J]. Acta Seismologica Sinica,41(2):155–168 (in Chinese).

    莘海亮,曾宪伟,康敏,高级. 2020. 海原弧形构造区地壳三维精细速度结构成像[J]. 地球物理学报,63(3):897–914. doi: 10.6038/cjg2020N0067

    Xin H L,Zeng X W,Kang M,Gao J. 2020. Crustal fine velocity structure of the Haiyuan arcuate tectonic zone from double-difference tomography[J]. Chinese Journal of Geophysics,63(3):897–914 (in Chinese).

    左可桢,陈继锋. 2018. 门源地区地壳三维体波速度结构及地震重定位研究[J]. 地球物理学报,61(7):2788–2801. doi: 10.6038/cjg2018L0537

    Zuo K Z,Chen J F. 2018. 3D body-wave velocity structure of crust and relocation of earthquakes in the Menyuan area[J]. Chinese Journal of Geophysics,61(7):2788–2801 (in Chinese).

    Allam A A,Schulte-Pelkum V,Ben-Zion Y,Tape C,Ruppert N,Ross Z E. 2017. Ten kilometer vertical Moho offset and shallow velocity contrast along the Denali fault zone from double-difference tomography,receiver functions,and fault zone head waves[J]. Tectonophysics,721:56–69. doi: 10.1016/j.tecto.2017.09.003

    Christensen D H,Abers G A. 2010. Seismic anisotropy under central Alaska from SKS splitting observations[J]. J Geophys Res:Solid Earth,115(B4):B04315.

    Christeson G L,Gulick S P S,van Avendonk H J A,Worthington L L,Reece R S,Pavlis T L. 2010. The Yakutat terrane:Dramatic change in crustal thickness across the transition fault,Alaska[J]. Geology,38(10):895–898. doi: 10.1130/G31170.1

    Chuang L,Bostock M,Wech A,Plourde A. 2017. Plateau subduction,intraslab seismicity,and the Denali (Alaska) volcanic gap[J]. Geology,45(7):647–650. doi: 10.1130/G38867.1

    DeMets C,Dixon T H. 1999. New kinematic models for Pacific-North America motion from 3 Ma to present,I:Evidence for steady motion and biases in the NUVEL-1A model[J]. Geophys Res Lett,26(13):1921–1924. doi: 10.1029/1999GL900405

    Eberhart-Phillips D,Christensen D H,Brocher T M,Hansen R,Ruppert N A,Haeussler P J,Abers G A. 2006. Imaging the transition from Aleutian subduction to Yakutat collision in central Alaska,with local earthquakes and active source data[J]. J Geophys Res:Solid Earth,111(B11):B11303.

    Enkelmann E,Zeitler P K,Pavlis T L,Garver J I,Ridgway K D. 2009. Intense localized rock uplift and erosion in the St Elias orogen of Alaska[J]. Nat Geosci,2(5):360–363. doi: 10.1038/ngeo502

    Feng L L,Ritzwoller M H. 2019. A 3-D shear velocity model of the crust and uppermost mantle beneath Alaska including apparent radial anisotropy[J]. J Geophys Res:Solid Earth,124(10):10468–10497. doi: 10.1029/2019JB018122

    Ferris A,Abers G A,Christensen D H,Veenstra E. 2003. High resolution image of the subducted Pacific (?) plate beneath central Alaska,50−150 km depth[J]. Earth Planet Sci Lett,214(3/4):575–588.

    Gou T,Zhao D P,Huang Z C,Wang L S. 2019. Aseismic deep slab and mantle flow beneath Alaska:Insight from anisotropic tomography[J]. J Geophys Res:Solid Earth,124(2):1700–1724. doi: 10.1029/2018JB016639

    Hacker B R,Peacock S M,Abers G A,Holloway S D. 2003. Subduction factory 2. Are intermediate-depth earthquakes in subducting slabs linked to metamorphic dehydration reactions?[J]. J Geophys Res:Solid Earth,108(B1):2030.

    Hansen P C,O’Leary D P. 1993. The use of the L-curve in the regularization of discrete ill-posed problems[J]. SIAM J Sci Comput,14(6):1487–1503. doi: 10.1137/0914086

    Hayes G P,Wald D J,Johnson R L. 2012. Slab1.0:A three-dimensional model of global subduction zone geometries[J]. J Geophys Res:Solid Earth,117(B1):B01302.

    Jadamec M A,Billen M I. 2010. Reconciling surface plate motions with rapid three-dimensional mantle flow around a slab edge[J]. Nature,465(7296):338–341. doi: 10.1038/nature09053

    Jadamec M A,Billen M I. 2012. The role of rheology and slab shape on rapid mantle flow:Three-dimensional numerical models of the Alaska slab edge[J]. J Geophys Res:Solid Earth,117(B2):B02304.

    Jadamec M A,Billen M I,Roeske S M. 2013. Three-dimensional numerical models of flat slab subduction and the Denali fault driving deformation in south-central Alaska[J]. Earth Planet Sci Lett,376:29–42. doi: 10.1016/j.jpgl.2013.06.009

    Jiang C X,Schmandt B,Ward K M,Lin F C,Worthington L L. 2018. Upper mantle seismic structure of Alaska from Rayleigh and S wave tomography[J]. Geophys Res Lett,45(19):10350–10359.

    Kissling E,Lahr J C. 1991. Tomographic image of the Pacific slab under southern Alaska[J]. Eclogae Geol Helv,84(2):297–315.

    Martin-Short R,Allen R,Bastow I D,Porritt R W,Miller M S. 2018. Seismic imaging of the Alaska subduction zone:Implications for slab geometry and volcanism[J]. Geochem Geophys Geosyst,19(11):4541–4560. doi: 10.1029/2018GC007962

    McNamara D E,Pasyanos M E. 2002. Seismological evidence for a sub-volcanic arc mantle wedge beneath the Denali volcanic gap,Alaska[J]. Geophys Res Lett,29(16):1814.

    Miller M S,Moresi L. 2018. Mapping the Alaskan Moho[J]. Seismol Res Lett,89(6):2430–2436. doi: 10.1785/0220180222

    Page R A, Stephens C D, Lahr J C.1989. Seismicity of the Wrangell and Aleutian Wadati-Benioff zones and the North American Plate along the Trans-Alaska crustal transect, Chugach mountains and Copper River basin, southern Alaska[J]. J Geophys Res: Solid Earth, 94(B11): 16059-16082.

    Qi C,Zhao D P,Chen Y. 2007. Search for deep slab segments under Alaska[J]. Phys Earth Planet Inter,165(1/2):68–82.

    Richter D H,Smith J G,Lanphere M A,Dalrymple G B,Reed B L,Shew N. 1990. Age and progression of volcanism,Wrangell volcanic field,Alaska[J]. Bull Volcanol,53(1):29–44. doi: 10.1007/BF00680318

    Rondenay S,Montési L G J,Abers G A. 2010. New geophysical insight into the origin of the Denali volcanic gap[J]. Geophys J Int,182(2):613–630. doi: 10.1111/j.1365-246X.2010.04659.x

    Sauber J,McClusky S,King R. 1997. Relation of ongoing deformation rates to the subduction zone process in southern Alaska[J]. Geophys Res Lett,24(22):2853–2856. doi: 10.1029/97GL52979

    Schultz R A,Aydin A. 1990. Formation of interior basins associated with curved faults in Alaska[J]. Tectonics,9(6):1387–1407. doi: 10.1029/TC009i006p01387

    Thurber C,Eberhart-Phillips D. 1999. Local earthquake tomography with flexible gridding[J]. Comput Geosci,25(7):809–818. doi: 10.1016/S0098-3004(99)00007-2

    Thurber C H,Brocher T M,Zhang H J,Langenheim V E. 2007. Three-dimensional P wave velocity model for the San Francisco Bay region,California[J]. J Geophys Res:Solid Earth,112(B7):B07313.

    van Stiphout T,Kissling E,Wiemer S,Ruppert N. 2009. Magmatic processes in the Alaska subduction zone by combined 3-D b value imaging and targeted seismic tomography[J]. J Geophys Res:Solid Earth,114(B11):B11302.

    Venereau C M A,Martin-Short R,Bastow I D,Allen R M,Kounoudis R. 2019. The role of variable slab dip in driving mantle flow at the eastern edge of the Alaskan subduction margin:Insights from shear-wave splitting[J]. Geochem Geophys Geosyst,20(5):2433–2448.

    Waldhauser F,Ellsworth W L. 2000. A double-difference earthquake location algorithm:Method and application to the northern Hayward fault,California[J]. Bull Seismol Soc Am,90(6):1353–1368. doi: 10.1785/0120000006

    Wang Y,Tape C. 2014. Seismic velocity structure and anisotropy of the Alaska subduction zone based on surface wave tomography[J]. J Geophys Res:Solid Earth,119(12):8845–8865. doi: 10.1002/2014JB011438

    Ward K M. 2015. Ambient noise tomography across the southern Alaskan Cordillera[J]. Geophys Res Lett,42(9):3218–3227. doi: 10.1002/2015GL063613

    Ward K M,Lin F C. 2018. Lithospheric structure across the Alaskan Cordillera from the joint inversion of surface waves and receiver functions[J]. J Geophys Res:Solid Earth,123(10):8780–8797. doi: 10.1029/2018JB015967

    Xin H L,Zhang H J,Kang M,He R Z,Gao L,Gao J. 2019. High-resolution lithospheric velocity structure of continental China by double-difference seismic travel-time tomography[J]. Seismol Res Lett,90(1):229–241. doi: 10.1785/0220180209

    You T,Zhao D P. 2012. Seismic anisotropy and heterogeneity in the Alaska subduction zone[J]. Geophys J Int,190(1):629–649. doi: 10.1111/j.1365-246X.2012.05512.x

    Zhang H,Thurber C H. 2007. Estimating the model resolution matrix for large seismic tomography problems based on Lanczos bidiagonalization with partial reorthogonalization[J]. Geophys J Int,170(1):337–345. doi: 10.1111/j.1365-246X.2007.03418.x

    Zhang H J,Thurber C H. 2003. Double-difference tomography:The method and its application to the Hayward fault,California[J]. Bull Seismol Soc Am,93(5):1875–1889. doi: 10.1785/0120020190

    Zhang Y,Li A B,Hu H. 2019. Crustal structure in Alaska from receiver function analysis[J]. Geophys Res Lett,46(3):1284–1292. doi: 10.1029/2018GL081011

    Zhao D P,Christensen D,Pulpan H. 1995. Tomographic imaging of the Alaska subduction zone[J]. J Geophys Res:Solid Earth,100(B4):6487–6504. doi: 10.1029/95JB00046

图(7)  /  表(1)
计量
  • 文章访问数:  754
  • HTML全文浏览量:  201
  • PDF下载量:  112
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-07-04
  • 修回日期:  2021-09-28
  • 网络出版日期:  2022-04-28
  • 发布日期:  2022-06-26

目录

/

返回文章
返回