Towards global seismic monitoring of underground nuclear explosions using waveform cross correlation. Part I: Grand master events

We propose to use waveform cross correlation techniques for seismic monitoring of the Comprehensive Nuclear-Test-Ban Treaty (CTBT) by proving the ability to build the Global Grid (GG) of master events using high-quality waveform templates at stations of the International Monitoring System. In seismically active regions, the choice of the highest quality and most representative signals from earthquakes as templates for master events is straightforward and one can reduce the current amplitude detection threshold at the International Data Centre (IDC) by a factor of 2 to 3. This gain practically doubles the number of built events, and thus, is crucial for seismic monitoring under the CTBT. However, the coverage by real masters is confined to areas with natural seismicity. Here, we investigate the possibility to populate the GG with real and synthetic master events as a two part study. In Part I, we estimate the effectiveness of cross correlation based on a regular grid filled with grand-master events with replicated high-quality waveform templates.In Part II, we develop and estimate efficiency of synthetic waveform templates for aseismic zones. Both approaches are quantitatively tested using the aftershock sequence of the April 11, 2012 Sumatera earthquake. The cross correlation bulletin previously built using sixteen real master events is a natural benchmark to evaluate the performance of the replicated masters and synthetic waveforms. In Part I, we prove that the replicated grand masters demonstrate the performance at the level of real masters. The Primary IMS Seismic Network includes many array stations, which provide a higher resolution monitoring based on cross correlation.     

2012年6月30日新疆新源-和静MS6.6 地震发震构造初步研究

2012年6月30日新疆维吾尔自治区新源-和静县交界发生M_S6.6地震,该地震是2010年青海玉树7.1级地震和2013年4月20日四川芦山7.0级地震之间中国大陆发生的最大的地震.本文基于新疆数字地震台网记录的此次地震序列震相资料,分别用绝对和相对定位方法联合对其进行重新定位,重新定位后余震展布为NW向,主震位置为43.429°N,84.755°E,深度为21.8 km.基于新疆地震台网记录6.6级地震波形数据,本文用CAP方法反演了震源机制解和震源深度.结果显示：M_S6.6地震震源机制解：节面Ⅰ走向39°,倾角46°,滑动角12°,节面Ⅱ走向301°,倾角81°,滑动角135°;震源深度为21 km,与利用地震震相到时确定的主震震源深度基本一致.主震震源机制解的节面Ⅱ与伊犁盆地北缘断裂走向和倾角基本一致,综合精确定位余震展布和伊犁盆地北缘断裂性质分析认为,新源-和静M_S6.6地震发震构造是伊犁盆地北缘断裂,震源深度为21 km左右,是一个高角的内陆倾滑地震.     

2010年河南太康MS4.6地震序列震源参数的精确确定

精确确定地震的震源机制解、震源深度和序列的相对位置是地震危险性分析的重要基础.2010年10月24日和2011年3月8日,河南太康地区分别发生M_S4.6和M_S4.1显著地震,为分析该弱震、少震区的地震危险性,本文利用CAP方法反演了震源机制解和震源深度,并结合深度震相波形记录进一步确认了深度的可靠性.结果显示,两次地震的震源机制解较为接近,均以走滑为主,深度也均为13km左右.此外,以主震作为参考事件,采用主事件法对余震的水平位置进行了相对定位,定位结果显示余震空间上分布的走向分别与主震震源机制解两个节面走向大致相同.本文结果为研究当地地震危险性提供了一定依据.     

2018年11月26日台湾海峡M_S6.2地震发震构造研究

本文利用福建省地震台网、广东省地震台网和台湾"中央"气象局17个台的宽频带记录,使用CAP方法反演了2018年11月26日台湾海峡M_S6.2地震震源机制解,得到节面1走向/倾角/滑动角为89°/82°/-173°,节面2走向/倾角/滑动角为358°/84°/-7°,最佳拟合深度14km,矩震级5.8.使用双差定位获取了94个M_L2.0以上地震的精定位结果,结果显示,主震位于北纬23.36°,东经118.62°,震源深度10.43km.根据小震分布和构造应力场反演得到余震断层面走向和倾角分别为88°和60°.研究认为,台湾海峡6.2级地震发震构造为近EW向的台湾浅滩断裂,受南海板块张裂拉伸发育而成,孕震过程中有东山隆起东缘断裂的参与,推测在菲律宾板块对欧亚板块NW-SE向挤压碰撞背景下,近EW向的台湾浅滩断裂与近NS向的东山隆起东缘断裂交接部位属于强度薄弱区,最终产生高倾角右旋走滑错动而引发地震,余震主要沿台湾浅滩断裂分布.     

结合波形互相关的双差定位方法在2008年汶川地震余震序列中的应用

选取了汶川地震主震后的2008年5月12日——2009年8月31日, 震级为3.0le;MSle;5.0的余震4240次．利用波形互相关方法得到其P波到时,用双差定位方法对其进行定位,最终得到了2441次重新定位的结果．统计定位误差（两倍标准偏差）在E-W方向为0.4 km,N-S方向为0.4 km,垂直方向为0.7 km．定位结果表明,汶川地震的余震深度集中在10——20 km,震中分布与龙门山中央断裂带的走向关系密切．沿龙门山断裂的地震分布具有明显的分段性,西南段呈水平带状分布,东北段接近垂直分布,且在北川附近存在深度突变．这与龙门山断裂的地震在西南段多表现为逆冲,东北段多表现为走滑的现象相吻合．在深度剖面上地震的空间分布存在分立的特征,通过对比前人在此地区浅层的地震剖面资料, 发现地震空间分布与已探知的浅部断层有较好的对应关系．      

The 20 April 2013 Lushan,Sichuan, mainshock,and its aftershock sequence: tectonic implications

Using the double-difference relocation algorithm, we relocated the 20 April 2013 Lushan, Sichuan, earthquake (M _S 7.0), and its 4,567 aftershocks recorded during the period between 20 April and May 3, 2013. Our results showed that most aftershocks are relocated between 10 and 20 km depths, but some large aftershocks were relocated around 30 km depth and small events extended upward near the surface. Vertical cross sections illustrate a shovel-shaped fault plane with a variable dip angle from the southwest to northeast along the fault. Furthermore, the dip angle of the fault plane is smaller around the mainshock than that in the surrounding areas along the fault. These results suggest that it may be easy to generate the strong earthquake in the place having a small dip angle of the fault, which is somewhat similar to the genesis of the 2008 Wenchuan earthquake. The Lushan mainshock is underlain by the seismically anomalous layers with low-V_P, low-V_S, and high-Poisson’s ratio anomalies, possibly suggesting that the fluid-filled fractured rock matrices might significantly reduce the effective normal stress on the fault plane to bring the brittle failure. The seismic gap between Lushan and Wenchuan aftershocks is suspected to be vulnerable to future seismic risks at greater depths, if any.     

玉树地震震源区速度结构与余震分布的关系

利用玉树震区21个应急流动地震台站和青海省地震台网固定地震台站的观测数据,采用双差层析成像方法,对2010年4月14日至6月15期间发生的地震进行了重定位,并反演得到了玉树地震震源区的三维速度结构.重定位结果揭示余震主要沿NW向成窄带状分布在断层的两侧,表明脆性破裂应力释放主要集中于一个狭窄的区域内.在西北端,余震偏离玉树—甘孜断裂分布,在SW向也有分布,推测可能与南西向次级断裂有关.双差层析成像得到的速度结构在浅部与地表地质构造相一致,中上地壳的速度结构显示巴颜喀拉地块为高速异常,羌塘地块为低速异常.玉树地震余震分布与特定的速度结构存在相关性：主震发生在高低速过渡带偏高速体的一侧,余震主要分布在高速体外围,高速体内部几乎没有余震分布.一般说来,中上地壳的高速体通常具有较高的强度,可以积累较强的孕震能量.主震发生后,高速体内积累的弹性能量向周边释放,可能是导致高速体周边余震发生的主要原因.     

2015尼泊尔M_S8.1地震中等余震震源机制研究

本研究利用西藏台网记录的波形数据,采用gCAP方法反演了2015年4月25日尼泊尔MS8.1大震5次中等余震(5.0≤MS≤6.5)及西藏定日MS5.9地震震源机制解.结果显示,6次地震包含2个正断、2个走滑及2个逆冲型地震.其中2个正断型地震位于主震的东北方向,即发震断层的上盘,表明该区域受到主震同震位移的影响,表现出应力拉张的变化特征;2个走滑型地震在主震破裂的东南方向上,说明随着破裂往东南方向延伸,余震的走滑分量增强;另外2个逆冲型地震位于5月12日MS7.5强余震区域,与MS7.5地震的滑移状态一致,可能与主震同震位移引起该区域处于应力挤压状态密切相关.这些结果表明,尼泊尔MS8.1主震发生后,由于同震位移的影响,不同区域处于不同的应力状态,从而使中等余震表现出不同的震源类型.     

2003年青海德令哈6.7级地震序列的震源机制解及其构造含义

利用8个流动数字地震台和国家数字地震台站的地震波形记录,测量了2003年4月17日青海德令哈6.7级地震及其主要余震的直达P波、SV波、SH波的初动方向和振幅比,应用Snoke(2003)的测定震源机制解的格点尝试法,测定出德令哈地震序列的48个2.4级以上地震的震源机制解.搜集分析了美国哈佛大学测定的德令哈6.7级主震和2004年二期地震活动中的7个地震的震源机制解.基于余震空间分布特征和对震源机制解特征的分析,讨论了德令哈地震序列的可能断层活动方式和地震的构造含义.结果表明,主震和大部分余震都是沿NWW-SEE走向的逆断层错动,北边的上盘可能沿低角度向北倾的断层面向南仰冲;个别正断层余震可能是震源区挤压变形弧顶区附近发生的局部张性破裂;在二期地震活动中,逆断层和走滑断层都有,走滑断层地震主要发生在震源区东侧.德令哈地震活动是青藏高原东北缘NWW-SEE向延伸的挤压带继续处于隆升活动中的表现,这一继承性新构造运动是德令哈地震序列的可能发震原因.     

Perspectives of Cross-Correlation in Seismic Monitoring at the International Data Centre

We demonstrate that several techniques based on waveform cross-correlation are able to significantly reduce the detection threshold of seismic sources worldwide and to improve the reliability of arrivals by a more accurate estimation of their defining parameters. A master event and the events it can find using waveform cross-correlation at array stations of the International Monitoring System (IMS) have to be close. For the purposes of the International Data Centre (IDC), one can use the spatial closeness of the master and slave events in order to construct a new automatic processing pipeline: all qualified arrivals detected using cross-correlation are associated with events matching the current IDC event definition criteria (EDC) in a local association procedure. Considering the repeating character of global seismicity, more than 90 % of events in the reviewed event bulletin (REB) can be built in this automatic processing. Due to the reduced detection threshold, waveform cross-correlation may increase the number of valid REB events by a factor of 1.5–2.0. Therefore, the new pipeline may produce a more comprehensive bulletin than the current pipeline—the goal of seismic monitoring. The analysts’ experience with the cross correlation event list (XSEL) shows that the workload of interactive processing might be reduced by a factor of two or even more. Since cross-correlation produces a comprehensive list of detections for a given master event, no additional arrivals from primary stations are expected to be associated with the XSEL events. The number of false alarms, relative to the number of events rejected from the standard event list 3 (SEL3) in the current interactive processing—can also be reduced by the use of several powerful filters. The principal filter is the difference between the arrival times of the master and newly built events at three or more primary stations, which should lie in a narrow range of a few seconds. In this study, one event at a distance of about 2,000 km from the main shock was formed by three stations, with the stations and both events on the same great circle. Such spurious events are rejected by checking consistency between detections at stations at different back azimuths from the source region. Two additional effective pre-filters are f–k analysis and F _prob based on correlation traces instead of original waveforms. Overall, waveform cross-correlation is able to improve the REB completeness, to reduce the workload related to IDC interactive analysis, and to provide a precise tool for quality check for both arrivals and events. Some major improvements in automatic and interactive processing achieved by cross-correlation are illustrated using an aftershock sequence from a large continental earthquake. Exploring this sequence, we describe schematically the next steps for the development of a processing pipeline parallel to the existing IDC one in order to improve the quality of the REB together with the reduction of the magnitude threshold.     

The Lemnos 8 January 2013 (M w = 5.7) earthquake: fault slip,aftershock properties and static stress transfer modeling in the north Aegean Sea

We investigate mainshock slip distribution and aftershock activity of the 8 January 2013 M _w?=?5.7 Lemnos earthquake, north Aegean Sea. We analyse the seismic waveforms to better understand the spatio-temporal characteristics of earthquake rupture within the seismogenic layer of the crust. Peak slip values range from 50 to 64 cm and mean slip values range from 10 to 12 cm. The slip patches of the event extend over an area of dimensions 16?×?16 km². We also relocate aftershock catalog locations to image seismic fault dimensions and test earthquake transfer models. The relocated events allowed us to identify the active faults in this area of the north Aegean Sea by locating two, NE–SW linear patterns of aftershocks. The aftershock distribution of the mainshock event clearly reveals a NE–SW striking fault about 40 km offshore Lemnos Island that extends from 2 km up to a depth of 14 km. After the mainshock most of the seismic activity migrated to the east and to the north of the hypocenter due to (a) rupture directivity towards the NE and (b) Coulomb stress transfer. A stress inversion analysis based on 14 focal mechanisms of aftershocks showed that the maximum horizontal stress is compressional at N84°E. The static stress transfer analysis for all post-1943 major events in the North Aegean shows no evidence for triggering of the 2013 event. We suggest that the 2013 event occurred due to tectonic loading of the North Aegean crust.     

RELOCATIONS AND FOCAL MECHANISM SOLOTIONS OF MS5.5 QIANGUO EARTHQUAKE SWARM IN JILIN PROVINCE IN 2013

We relocated a seismic swarm, which started in a mass from 31 October, 2013 in Qianguo County of Jilin Province, by using double difference location method, based on the phase data of regional digital seismic network and the crustal velocity model of Sunliao Basin. The characteristics of seismogenic fault have been investigated based on the spatial distribution image of the seismic swarm and the geophysical data near the epicenter area. The relocated epicenters of the swarm earthquakes have a precision of 0.9km in E-W, 0.7km in N-S and 1.2km in U-D direction, and show an apparent concentrated seismic belt trending N-W, with a length and width of 12km and 6km, respectively. The source depths of all events are shallow, with 80%in a range of 6~8km, and the events are apparently crowded together on the depth cross section. According to the relocated spatial distribution characteristics of the seismic swarm, the features that the medium size events happened successively, and the focal mechanism of the large size events in the swarm, we infer that the seismogenic tectonics of Qianguo seismic swarm is the thrust nappe structure inside the Keshan-Da'an fault zone. The fault plane inclines to the East direction, and is steep when close to the ground surface, which shows the typical characteristics of a listric thrust fault. The longitudinal length of the rupture plane is greater than the transverse length. According to the features of seismogenic tectonics, we infer that the three M_S ≥ 5.0 earthquakes occurred at the lower layer of the thrust rupture surface of the fault, while the aftershocks were triggered by the three events and occurred mainly at the upper layer of the rupture surface.     

Fault complexity associated with the 14 August 2003 M<Subscript>w</Subscript>6.2 Lefkada,Greece, aftershock sequence

The M _w 6.2 Lefkada earthquake occurred on 14 August 2003 beneath the western coastline of Lefkada Island. The main shock was followed by an intense aftershock activity, which formed a narrow band extending over the western coast of the Island and the submarine area between Lefkada and Kefalonia Islands, whereas additional off fault aftershocks formed spatial clusters on the central and northwestern part of the Island. The aftershock spatial distribution revealed the activation of along-strike adjacent fault segment as well as of secondary faults close to the main rupture. The properties of the activated segments were illuminated by the precisely located aftershocks, fault plane solutions determination and the cross sections performed parallel and normal to their strike. The aftershock focal mechanisms exhibited mainly strike slip faulting throughout the activated area, although deviation of the dominant stress pattern is also observed. The results help to emphasize the importance of the identification of activated nearby fault segments possibly triggered by the main rupture. Because such segments are capable to produce moderate events causing appreciable damage, they should be viewed with caution in seismic hazard assessment in addition to the major regional faults.     

Relocation of the M 8.0 Wenchuan earthquake and its aftershock sequence

We relocated M8.0 Wenchuan earthquake and 2706 aftershocks with M⩾2.0 using double-difference algorithm and obtained relocations of 2553 events. To reduce the influence of lateral variation in crustal and upper mantle velocity structure, we used different velocity models for the east and west side of Longmenshan fault zone. In the relocation process, we added seismic data from portable seismic stations close to the shocks to constrain focal depths. The precisions in E-W, N-S, and U-D directions after relocation are 0.6, 0.7, and 2.5 km respectively. The relocation results show that the aftershock epi-centers of Wenchuan earthquake were distributed in NE-SW direction, with a total length of about 330 km. The aftershocks were concentrated on the west side of the central fault of Longmenshan fault zone, excluding those on the north of Qingchuan, which obviously deviated from the surface fault and passed through Pingwu-Qingchuan fault in the north. The dominant focal depths of the aftershocks are between 5 and 20 km, the average depth is 13.3 km, and the depth of the relocated main shock is 16.0 km. The depth profile reveals that focal depth distribution in some of the areas is characterized by high-angle westward dipping. The rupture mode of the main shock features reverse faulting in the south, with a large strike-slip component in the north. Supported by the Basic Research Project of Institute of Geophysics, China Earthquake Administration (Grant No. DQJB08Z03)     

芦山MS7.0级地震余震序列重新定位及构造意义

利用四川省地震台网的震相数据和双差定位方法对芦山M_S7.0级地震及其余震序列进行了精确定位,根据余震分布确定了发震断层的位置和断层面的几何特征,并对余震活动进行了分析.结果显示,芦山M_S7.0级地震的震中位于30.28°N、102.99°E,震源深度为16.33 km.余震沿发震断层向主震两侧延伸,主要分布在长约32 km、宽约15~20 km、深度为5~24 km的范围内.地震破裂带朝西南方向扩展范围较大,东北方向略小,余震震级随时间迅速衰减.震源深度剖面清晰地显示出发震断层的逆冲破裂特征,推测发震断层为大川—双石断裂东侧约10 km的隐伏断层.该断层走向217°、倾向北西,倾角约45°,产状与大川—双石断裂相比略缓,它们同属龙门山前山断裂带的叠瓦状逆冲断层系.受发震断裂影响,部分余震沿大川—双石断裂分布,西北方向的余震延伸至宝兴杂岩体的东南缘,与汶川地震的破裂带之间存在50 km左右的地震空区,有可能成为未来发生强震的潜在危险区.     

2015年7月3日皮山6.5级地震发震构造初步研究

基于新疆区域数字地震台网记录,采用CAP（Cut and Paste）方法反演了2015年7月3日皮山6.5级主震和部分M_S3.6以上余震的震源机制解和震源深度;采用HypoDD方法重新定位了序列中M_L2.5以上地震序列的震源位置,并利用小震分布和区域应力场拟合了可能存在的发震断层面参数.基于上述研究,综合分析了皮山6.5级地震序列的震源深度、震源机制和震源破裂面特征,探讨可能的发震构造.结果显示,利用CAP方法得到的最佳双力偶机制解节面I：走向280°/倾角60°/滑动角90°;节面Ⅱ：走向100°/倾角30°/滑动角90°,矩心深度19 km,表明该地震为一次逆冲型地震事件.大部分M_S3.6以上余震震源机制与主震具有一定的相似性.双差定位结果显示,M_L2.5以上的余震序列主要分布在主震的西南方向,深度主要分布在0～15 km范围内,余震分布显示出与发震构造泽普隐伏断裂一致的倾向南西的特征.利用小震分布和区域应力场拟合得到发震断层参数为走向104°/倾角34°/滑动角94°,该结果与主震震源机制解中节面Ⅱ的滑动角较为接近,绝大多数余震发生在断层面附近10 km左右的区域.根据本研究得到的震源机制、精定位结果以及利用小震分布和区域应力场拟合得到的断层面的参数,结合震源区地质构造情况,初步给出了此次皮山6.5级地震的发震模式.     

2008年8月30日攀枝花-会理6.1级地震序列 ML≥4.0事件的震源机制解

针对2008年8月30日在四川攀枝花-会理发生的M_s6.1地震序列,本研究基于四川和云南两省数字地震台网的宽频带波形记录,采用CAP方法反演了该序列主震及M_L≥4.0余震的震源机制解.结果显示:主震震源机制解的两个主应力轴仰角小于10°,其中,主压力轴方位为140°;节面之一走向185°、西倾83°、滑动角5°,显示左旋走滑略兼逆冲分量的断层作用性质.结合余震、烈度分布以及震区的活动构造,判定该节面代表了主震的发震断层面,相应的发震断层应是穿越震区的近南北向红格断裂(南段).本研究还获得主震震源机制解的最佳拟合误差深度为10 km,与该事件的定位结果相一致.该序列中6次M_L≥4.0余震也具有与主震类似的震源机制解.分析初步表明:空间上,2008年攀枝花-会理M_s6.1地震序列的震源机制解与研究区内更早地震的震源机制解具有良好的协调性,反映了该序列是在川滇地块SE-SSE向水平运动的背景下、沿近S-N向红格断裂发生左旋走滑略兼逆冲运动的结果.     

华北地区震源断层与深浅构造关系的初步研究

使用地震资料（等震线、地表破裂带、地形变、余震平面和剖面分布及震源机制）描述强震震源断层的立体特征。发现震源断层的上界即为由地震测深资料得到的Ｇ界面，并与该界面之下的高倾角深断裂吻合，证明地震时深断裂发生了剪切错动。震源断层与Ｇ界面之上的铲形正断层构成不连接的两套断裂系统，它们之间存在有对应和无对应的复杂关系。     

Rupture process and aftershock focal mechanisms of the 2022 M6.8 Luding earthquake in Sichuan

According to the China Earthquake Networks Center, a strong earthquake of M6.8 occurred in Luding County, Ganzi Tibetan Autonomous Prefecture, Sichuan Province, China (102.08°E, 29.59°N), on September 5, 2022, with a focal depth of 16 km. Rapid determination of the source parameters of the earthquake sequence is vital for post-earthquake rescue, disaster assessment, and scientific research. Near-field seismic observations play a key role in the fast and reliable determination of earthquake source parameters. The numerous broadband seismic stations and strong-motion stations recently deployed by the National Earthquake Intensity Rapid Report and Early Warning project have provided valuable real-time near-field observation data. Using these near-field observations and conventional mid- and far-field seismic waveform records, we obtained the focal mechanism solutions of the mainshock and M ≥ 3.0 aftershocks through the waveform fitting method. We were further able to rapidly invert the rupture process of the mainshock. Based on the evaluation of the focal mechanism solution of the mainshock and the regional tectonic setting, we speculate that the Xianshuihe fault formed the seismogenic structure of the M6.8 strong earthquake. The aftershocks formed three spatially separated clusters with distinctly different focal mechanisms, reflecting the segmented nature of the Xianshuihe fault. As more high-frequency information has been applied in this study, the absolute location of the fault rupture is better constrained by the near-field strong-motion data. The rupture process of the mainshock correlates well with the spatial distribution of aftershocks, i.e., aftershock activities were relatively weak in the maximum slip area, and strong aftershock activities were distributed in the peripheral regions.     

2010年以来四川地区中强地震震源机制反演及深度确定

2008年5月12日四川龙门山断裂带发生了汶川8.0级地震,之后四川境内发生了两次7.0级地震(其中一个是芦山地震),为了研究汶川地震之后龙门山断裂带及周边区域的地震活动性,本研究收集了国家地震台网和四川区域地震台网2010年1月1日—2017年12月31日四川地区发生的17次M≥5.0地震以及120多次5.0>M≥4.0地震的波形资料,利用波形拟合法反演了震源机制解及区域应力场.反演结果显示,位于龙门山断裂带上的地震,震源机制以逆冲型为主,鲜水河断裂带地震震源机制以走滑型为主,而川滇块体西南部的理塘断裂、金沙江断裂附近,震源机制解以正断层为主.根据震源机制解反演得到的龙门山地区、鲜水河地区的主压应力场方向为WNW、近EW向.川滇块体的巴塘、理塘等地区,其主压应力轴方向为12°左右,接近SN向,且仰角接近40°左右.本研究利用面波振幅谱特征对震源深度进行了精确定位,定位结果与中国地震台网中心(CENC),美国地震调查局(USGS),国际地震中心(ISC)等机构地震目录进行了对比.结果显示,四川地区强震震源深度主要分布在20km以上的中上地壳.龙门山地区震源优势分布在10～20km,鲜水河断裂地震震源深度在10km左右,川滇块体西南部的理塘断裂,巴塘断裂,金沙江断裂等地区,震源深度一般在5～10km范围.