Rapid evaluation of radiated seismic energy for great shallow earthquakes from 2014 to 2019
-
摘要: 根据地震波衰减特性,采用一维速度模型开展了快速测定辐射能量ES和能量震级Me的方法研究。利用全球地震台网和国家数字测震台网提供的宽频带资料,测定了2014—2019年间MW≥6.0的115次浅源地震的辐射能量和能量震级,将计算结果与其它机构的结果进行对比。结果表明:利用本文方法可在得到地震数据半小时内计算出稳定的Me,且本文的测定结果与美国地震学研究联合会的结果基本一致。地震造成的灾害与能量震级的大小密切相关,当Me>MW时,地震灾害较为严重;在所有的地震类型中,发震断层类型为走滑型时,其地震辐射能量的效率高,Me明显大于MW。通过分析2018年2月4日和2019年4月18日台湾花莲两次MW6.1地震能量释放的差异得出,发生在相似位置且具有相同的震源机制的两次地震,尽管它们的MW相同,但Me相差很大,接近0.5。Me与MW的差异表明,MW只能获得有关震源的静态特征,它与地震引起的断裂面积、断裂平均位错等静态构造效应密切相关,而Me可以提供震源的动态信息,从而客观评价地震的破坏强度。因此,本文使用的方法既能准确测定能量震级Me又能极大提高其测定速度,非常适用于快速反应系统。本研究可以为未来地震台网将能量震级Me作为日常产出震级提供参考,为快速评估大地震造成的灾害提供更多信息。Abstract: Based on seismic waves’ attenuation characteristics and the one-dimensional velocity model, we conducted a study for the rapid determination of radiated seismic energy ES and energy magnitude Me. The seismic recordings obtained from the Global Seismographic Network and China Seismological Digital Network were used to calculate the ES and Me of 115 shallow earthquakes with MW≥6.0 from 2014 to 2019, and the results were compared with MW and Me produced by other institutions. The results show the stable Me can be obtained within half an hour after obtaining the seismic data, and the Me values of IRIS were consistent with the earthquakes analyzed in this study. The earthquake damage is closely related to the size of energy magnitude Me, and the devastation that an earthquake can cause is more serious when the energy magnitude Me is far greater than the moment magnitude MW. In all types of earthquakes, the efficiency of radiated energy for earthquakes with strike-slip faults is high, with Me being significantly greater than MW. Next, by analyzing two MW6.1 earthquakes occurred in Hualien, Taiwan region on February 4, 2018 and April 18, 2019, it is concluded that two earthquakes with the same source mechanism in similar locations, although they have the same moment magnitude MW, the energy magnitude Me vary greatly, close to 0.5. MW can only obtain the source’s static characteristic, which is closely related to the static tectonic effects such as fault area and average dislocation of rupture caused by earthquakes. In contrast, Me can provide the source’s dynamic information and objectively evaluate the damage intensity of an earthquake. Therefore, the approach applied in this study can ensure the accuracy of the results and greatly improve the measurement speed of energy magnitude Me, which is very suitable for a rapid response system. Our research results can provide a reference for the future seismic network to take energy magnitude Me as daily output magnitude and provide more information for rapid assessment of disasters caused by large earthquakes.
-
Keywords:
- radiated seismic energy /
- energy magnitude /
- broadband recording /
- moment magnitude
-
引言
非火山型微震颤(non-volcanic tremor,简写为NVT)是位于板块隐没、远离火山地区的一种类似火山地区微震颤的地震活动。该类活动首次在日本西南部被发现,不同于一般的地震活动,NVT信号没有明显的体波到时;振幅小;能量频段集中在1—10 Hz,较普通地震低;活动持续时间较长,约数分钟至数天甚或数周。大地测量学证据表明非火山型微震颤常伴随着缓慢的板块滑移事件(Obara,2002;Rogers,Dragert,2003)。全球许多地方均观测到了NVT信号,例如卡斯凯迪亚断层区(Ide,2012)、圣安得烈斯断层区(Nadeau,Dolenc,2005;Shelly,2010)、中国台湾(Peng,Chao,2008;Tang et al,2010)、阿拉斯加(Peterson,Christensen,2009)、智利南部及新西兰(Fry et al,2011;Kim et al,2011)等地区。常规地震通常发生在上地壳脆性带,而NVT通常发生在下地壳至莫霍面之间的区域,因此观测NVT信号成为监测下地壳破裂及研究变形机理的一种方法(Daub et al,2011;Chao et al,2012)。
Peng和Chao (2008)首次在台湾观察到NVT事件,此次事件由2001年11月14日我国青海省昆仑山口西MS8.1地震所触发,NVT震源位于中央山脉深处的滑脱断层上,深度约20 km;出现信号的时间点与远场地震的勒夫波到时几乎一致,因此此次事件被认为是由昆仑山口西地震勒夫波产生的动态应力触发。Tang等(2010)分析得出中央山脉南部的部分NVT信号由一些P波和S波信号较微弱的低频地震(low-frequency earthquakes,缩写为LFEs)组成,这与观测到的位于其它地区板块边界的NVT非常相似。Chao等(2012)分析了1998—2009年45次MW≥7.5的远场地震事件,其中有9次触发了NVT事件,其定位结果显示,这9次NVT事件的震源位于台湾中央山脉的南段,震源深度介于15—25 km之间,分析认为远场地震面波的振幅大小是触发的关键条件,并得出触发NVT的动态应力阈值为7—8 kPa。由于Chao等(2012)主要针对切向分量的分析,确定的触发阀值不够全面,因此,本文拟针对波形的三分量进行综合考量,以确定更全面的触发阈值。
虽然NVT已在全球很多地区被观测到,但对于不同构造环境中NVT的活动特性尚未明确。因此,本文选择对NVT事件有一定先验条件、区内构造复杂、台站分布密集、地震活动活跃的我国台湾地区为研究区,讨论NVT触发所需的条件和潜在机理,并对其空间分布特征予以分析,以期了解NVT震源附近的孕震构造,为地震学的进一步研究提供更多的参考。
1. 研究区域概况与研究方法
1.1 研究区域概况
中国台湾地处地球上板块运动最活跃的区域—环太平洋构造带上,位于欧亚大陆板块与菲律宾海板块的碰撞边界。古生代以来,台湾地区就位于华夏古陆的边缘,也因如此,长期受板块作用的影响远大于板块内部环境的影响。中新世至今,菲律宾海板块每年以7—8 cm的速度向西北方向移动(Yu et al,1997),在两个板块相互作用之下,菲律宾海板块向北隐没于欧亚大陆板块之下(Teng,1996),而南中国海板块则向东俯冲于菲律宾海板块之下(Chai,1972)(图1b)。台湾岛位在两个俯冲带之间,两个俯冲带与左移兼具逆冲的花东纵谷断层相连接。复杂的构造环境导致了台湾地区的地震活动频繁(Wu et al,2007)。本文在台湾南部地区(21°30'N—24°00'N,120°00'E—121°30'E)展开研究,研究区内台站和由4次远震触发的NVT事件震中和震源深度分布如图1所示。
图 1 研究区域的构造背景及台站和非火山型微震颤(NVT)事件分布(a) 研究区台站和NVT震中分布;(b) 区域构造背景;(c) 触发NVT的4次地震位置示意图;(d) AA′ 剖面上的NVT事件的震源深度分布Figure 1. Tectonic settings and distribution of stations and non-volcanic tremor (NVT) locations(a) Distribution of seismic stations and locations of NVTs;(b) Regional tectonic settings;(c) Schematic diagram for location of four earthquakes that triggered NVT;(d) Distribution of NVT depths along the cross-section AA′1.2 数据与方法
本研究首先从全球地震矩张量解(global centroid-moment-tensor,缩写为GCMT)目录中选取2011年1月1日至2016年12月31日的38次距离台湾岛大于1 000 km且MW≥7.5的地震事件。采用台湾宽频地震观测网(Broadband Array in Taiwan for Seismology,缩写为BATS)、台湾气象局(Weather Bureau Seismic Network,缩写为WBSN)及美国地震学研究联合会(Incorporated Research Institutions for Seismology,缩写为IRIS)发布的地震发生后3小时内的连续波形资料进行分析,波形数据的采样率为20 Hz或100 Hz。
参考Peng和Chao (2008)一文中的方法,去除仪器响应后,对波形进行频谱分析,如图2所示。在远震信号中,频率最高的能量通常来自P波,若在频谱分析过程中面波经过的时段内存在比P波频率高的信号,则可以初步认定此信号为区域地震信号。由于受波形数据采样率限制,参考Chao等(2012)的信号处理流程,本文统一对波形进行2—8 Hz带通滤波,若5个以上台站均接收到高频、无脉冲且持续时段大于100 s的连续波包信号时,即可认定为NVT事件发生。
检测到NVT信号后,计算每个信号的包络函数,并使用包络函数互相关法进行定位(Obara,2002;Peng,Chao,2008;Chao et al,2013)。此方法的原理是将台站对之间的观测走时差与理论走时差进行比较,以网格搜寻的方式将所有台站对的走时残差均方根(root mean square,缩写为RMS)最小的位置确定为震源位置,本文中4次NVT事件的定位结果如图3所示。假设使用N个台站的包络函数进行定位,则有ni, j=N(N-1)/2个台站对,网格点(经度x,纬度y,深度z)的均方差如下:
图 3 4次NVT事件的定位示意图(a) 东日本MW9.1地震触发事件;(b) 苏门答腊MW8.6地震触发事件;(c) 所罗门群岛MW7.6地震触发事件;(d) 尼泊尔MW7.9地震触发事件Figure 3. The schematic diagrams of locating four NVT events(a) NVT triggered by MW9.1 Tohoku earthquake;(b) NVT triggered by MW8.6 Sumatra earthauake; (c) NVT triggered by MW7.6 Solomon earthquake;(d) NVT triggered by MW7.9 Nepal earthquake${\rm{RM}}{{\rm{S}}_{\left( {x{\text{,}}\!\!\!\! y{\text{,}}\!\!\!\!z} \right)}} {\text{=}} \sqrt {\frac{{\displaystyle\mathop \sum \limits_{i {\text{=}} 1,j {\text{=}} 2}^n {{(\Delta {T_{i,j}})}^2}}}{n}} {\text{=}} \sqrt {\frac{{\displaystyle\mathop \sum \limits_{i {\text{=}} 1,j {\text{=}} 2}^n {{\left( {\Delta t_{i{\text{,}}\!\!\!\!j}^{{\rm{theo}}} - \Delta t_{i{\text{,}}\!\!\!\!j}^{{\rm{obs}}}} \right)}^2}}}{n}}{\text{,}} $
(1) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
表 1 台湾地区一维速度模型(Wu et al,2007)Table 1. 1-D velocity model in Taiwan region (Wu et al,2007)vP/(km·s−1) vP/vS vS/km·s−1 深度范围/km vP/(km·s−1) vP/vS vS/km·s−1 深度范围/km 3.82 1.74 2.20 0—2 6.13 1.75 3.50 21—25 4.93 1.68 2.93 2—4 6.25 1.71 3.65 25—30 5.51 1.70 3.24 4—6 6.54 1.73 3.78 30—35 5.63 1.70 3.31 6—9 7.03 1.75 4.02 35—50 5.83 1.76 3.31 9—13 8.00 1.74 4.60 50—70 5.99 1.75 3.42 13—17 8.34 1.73 4.82 70—80 6.06 1.75 3.46 17—21 This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
2. 远场大规模地震触发的NVT事件
在38次远场大规模地震中,触发研究区内NVT事件的4次地震分别为:2011年3月11日东日本MW9.1地震,2012年4月11日苏门答腊岛MW8.6地震,2014年4月12日所罗门群岛MW7.6地震和2015年4月25日尼泊尔MW7.9地震,这4次地震的相对发震位置如图1c所示,具体信息列于表2。
表 2 触发NVT的四次远震事件信息概要Table 2. General information of the four NVT-triggered teleseismic events发震时间 地点 震中位置 震源深度/km MW 震中距/° TPUB台站
反方位角/°年−月−日 时:分:秒 北纬/° 东经/° 2011−03−11 5:47:32.8 日本本州岛东岸近海 37.52 143.05 20.0 9.1 23.893 48 2012−04−11 8:39:31.4 苏门答腊岛北部西岸远海 2.35 92.82 45.6 8.6 34.049 236 2014−04−12 20:14:49.9 所罗门群岛 11.35 162.24 27.3 7.6 53.327 126 2015−04−25 6:11:58.6 尼泊尔 27.91 85.33 12.0 7.9 32.092 286 注:表中震中距是以TPUB台站为标准计算得到的。 2.1 2011年3月11日东日本MW9.1地震触发的NVT事件
2011年3月11日东日本地震距研究区2 652 km,勒夫波理论到时为647 s,瑞雷波理论到时为758 s。此次地震在4次触发NVT的地震事件中震中距最小、震级最大,而且在全球很多地区均触发了浅层微地震和NVT事件(Miyazawa,2011;Gonzalez-Huizar et al,2012;Hill et al,2013;Peng et al,2013),在日本四国岛、美国阿留申岛、美国阿拉斯加、美国卡斯凯迪亚、美国加州中部的帕克菲尔德、美国加州南部的圣哈辛托、新西兰北岛等地区均观测到NVT信号(Chao et al,2013)。而台湾地区观测到该地震触发的较为明显的NVT信号从震后650 s左右开始,持续了约300 s (图4a)。通过对明显NVT的出现时间及面波理论到时(图4b)的对比,推断本次NVT事件由东日本MW9.1地震的勒夫波触发,触发的NVT震中位置为(23.05°N,120.90°E)。
图 4 2011年3月11日东日本大地震触发的NVT事件波形图(a) 各台站南北向分量2—8 Hz带通滤波后的波形图,波形左侧为台站名及其距NVT震源距离;(b) 仪器校正后距离NVT震源最近台站的径向、垂向、切向的波形及滤波后波形,红色、黑色箭头分别指勒夫波(4.1 km/s)和瑞雷波(3.5 km/s)理论到时,下同Figure 4. Tremor triggered by the Tohoku earthquake on 11 March 2011(a) The 2−8 Hz bandpass-filtered seismograms of N-S component, The hypocentral distance between each station and tremor source is shown next to the station name on left of each trace;(b) The radial,vertical and transverse component seismograms,in which instrument response has been removed,recorded at the nearest station and bandpass-filtered seismograms The red and black vertical arrows indicate the predicted arrivals of the Love and Rayleigh waves with the apparent velocity of 4.1 and 3.5 km/s,respectively,the same below2.2 2012年4月11日苏门答腊岛MW8.6地震触发的NVT事件
2012年4月11日印度尼西亚苏门答腊岛北部西岸远海发生MW8.6地震,距本文研究区约3 779 km,勒夫波与瑞雷波理论到时分别为922 s和1 080 s。该地震在日本多个位置触发了NVT事件,这些震源位于不同的断层系统中,如:日本九州岛西部的走滑断层、日本关东的逆冲断层、琉球群岛西南部的俯冲带区域等(Chao,Obara,2016)。由图5可见:南北分量原始波形经过2—8 Hz带通滤波后可以看到NVT信号,从1 150 s左右开始,持续了大约500 s,被YULB,FULB,TWGB,ALSB和TPUB等5个台站记录到;而一些微弱的NVT信号在920 s左右出现,被YULB,FULB,TWGB和ALSB等4个台站记录到。使用网格搜寻的方法对此次事件进行定位,结果显示其震中位置为(23.05°N,120.90°E)。
图 5 2012年4月11日苏门答腊岛地震触发NVT事件波形图(a) 各台站南北向分量2—8 Hz带通滤波后的波形图;(b) 仪器校正后距离NVT震源最近台站的径向、垂向、切向的波形及滤波后波形Figure 5. Tremor triggered by the Sumatra earthquake occurred on 11 April 2012(a) The 2−8 Hz bandpass-filtered seismograms of N-S component;(b) The radial,vertical and transverse component seismograms,in which instrument response has been removed,recorded at the nearest station and bandpass-filtered seismograms2.3 2014年4月12日所罗门群岛MW7.6地震触发的NVT事件
2014年4月12日所罗门群岛MW7.6地震,距本文研究区5 913 km,勒夫波理论到时为1 444 s,瑞雷波理论到时为1 691 s。与其它3次触发NVT的地震事件相比,此次地震的震中距最大,震级最小。由图6波形可以看出,明显NVT信号从1 450 s左右开始,持续了大约450 s。此NVT事件的震中位置为(23.05°N,121.00°E)。
图 6 2014年4月12日所罗门群岛地震触发的NVT事件波形图(a) 各台站南北向分量2—8 Hz带通滤波后的波形图;(b) 仪器校正后距离NVT震源最近台站的径向、垂向、切向的波形及滤波后波形Figure 6. Tremor triggered by the Solomon earthquake on 12 April 2014(a) The 2−8 Hz bandpass-filtered seismograms of N-S component;(b) The radial,vertical and transverse component seismograms,in which instrument response has been removed,recorded at the nearest station and bandpass-filtered seismograms2.4 2015年4月25日 尼泊尔MW7.9地震触发的NVT事件
2015年4月25日尼泊尔MW7.9大地震,距本文研究区3 562 km,勒夫波理论到时为869 s,瑞雷波理论到时为1 018 s。此次地震与其它3次地震事件相比,其震源深度最浅,为12 km,且其在重庆綦江触发了ML3.8地震(Han et al,2017)。明显的NVT信号从1 150 s左右开始,持续了大约420 s,一些微弱的NVT信号在900 s左右出现,被SYNB,TPUB,TWGB和FULB等4个台站记录到(图7),推断本次事件由尼泊尔地震的勒夫波触发,其震中位置为(23.15°N,120.90°E)。
图 7 2015年4月25日尼泊尔地震触发的NVT事件波形图(a) 各台站南北向分量2—8 Hz带通滤波后的波形图;(b) 仪器校正后距离NVT震源最近台站的径向、垂向、切向的波形及滤波后波形Figure 7. Tremor triggered by the Nepal earthquake on 25 April 2015(a) The 2−8 Hz bandpass-filtered seismograms of N-S component;(b) The radial,vertical and transverse component seismograms,in which instrument response has been removed,recorded at the nearest station and bandpass-filtered seismograms3. 讨论与结论
3.1 NVT位置
本研究中,4次远场大地震在台湾南部触发了明显的NVT事件,震源位于中央山脉南段,潮州—梨山断层与花东纵谷之间,即欧亚大陆板块与菲律宾海板块交界处,深度分布在21—53 km。Peng和Chao (2008)以及Tang等(2010)定位到NVT震中在TPUB台站与TWGB台站之间,震源深度范围为15—25 km,本研究结果与其相比,震中位置大致相同,震源深度较深。Peng和Gomberg (2010)研究表明,NVT的出现与断层的蠕滑行为有密切关联,而蠕滑产生的重要因素则是液体的存在。中央山脉南段深部区域(15—50 km)的地下介质品质因子Q值较周边地区低(陈光荣,1993),电导率较周边地区高(Chen,Chen,1998),且液体孔隙压力较高(Lee et al,2010),表明该地区流体含量较周边地区大,导致断层面的耦合度降低,在受到来自面波的瞬时动态应力扰动时产生蠕滑(Chao et al,2012;Tang et al,2010),进而产生NVT信号。此外,NVT震源区域几乎无地震发生(图1a),而在此区域之内发生的地震震源机制多以正断层为主(Kao,Jian,2001;Hsu et al,2010),因此推测NVT发生在拉张型的地质环境中。
3.2 触发NVT的地震特征
为了检测动态应力是否为触发NVT的主要因素,本文利用PGV计算了38次远震的切向、垂向与径向动态应力,结果如图8所示,可见:切向、垂向和径向的动态应力范围分别为0.207—33.6 kPa,0.279—50.5 kPa,0.257—38.8 kPa;除了所罗门群岛MW7.6地震,其余3次触发NVT的地震动态应力均高于未触发地震;如图8c所示,在尼泊尔MW7.9地震附近的灰色正方形代表东日本MW9.1地震的一次MW7.9余震,其后十分钟左右又有一次MW7.6余震,余震MW7.9可能触发NVT,但由于其面波受到后面发生地震的干扰,无法辨别出NVT事件。而所罗门群岛MW7.6地震的震中距较其它3次地震大,且震级最小,这使得其地震波传播到研究区时产生的地表运动速度较小,动态应力值未超过14 kPa,但也触发了NVT事件。Chao等(2012)提出切向分量的触发阈值为7—8 kPa,本研究中仅所罗门群岛地震的阈值不满足此条件。除此之外,我们还计算了波形资料的信噪比,结果显示信噪比结果并无规律可循,不能利用其对是否触发NVT事件进行有效区分。
图 8 TPUB台站记录到的大地震径向 (a)、垂向 (b) 和切向 (c) 分量的PGV、动态应力及信噪比分布图背景噪声分析时窗长度为远震事件发震前600 sFigure 8. The scattergram of PGV,dynamic stress and signal-to-noise ratio for the radial (a),vertical (b) and tangential (c) components of the earthquakes recorded by the station TPUBThe ambient noise level of each event is calculated from a 600 s time window before the occurrence of each main shock通过NVT出现时间和面波理论到时的对比,分析得出本文研究的4次远震触发NVT事件均由勒夫波触发,但在瑞雷波到达后振幅增大,这与Chao等(2012)的研究结果有差异。该研究认为当远震面波垂直于中央山脉走向(N16°E)时,勒夫波触发NVT,但除了垂直中央山脉走向入射的尼泊尔地震的勒夫波触发NVT以外,其余3次地震的勒夫波到达时均有NVT信号出现,故推断勒夫波到达时NVT事件被触发,随后瑞雷波使NVT震源振动幅度增大,起到了调节的作用。
3.3 结论
本研究对2011年至2016年间发生的38次MW≥7.5地震(震中距离台湾大于1 000 km)在台湾南部触发的NVT行为进行了调查,结果显示共有4次远震在台湾南部触发了NVT事件,主要频率范围为2—8 Hz,持续时间为300 s至500 s不等,震中位置位于台湾中央山脉南段,深度范围为21—53 km,震源深度的定位结果仍有较大的误差,希望日后能够改进定位方法,得到更加精确的深部位置。本文还讨论了触发NVT事件的远震特征,包括远震产生的径向、垂向和切向动态应力,其中:东日本MW9.1地震、苏门答腊岛MW8.6地震、尼泊尔MW7.9地震产生的动态应力均大于14 kPa,其它未触发NVT事件的地震产生的应力均小于该值;所罗门群岛MW7.6地震产生的动态应力较小,但也触发了NVT事件,究其原因还需今后对大地震的触发行为进行统计,同时需要对远震进行更多空间域及频率域的讨论。
台湾宽频地震观测网、台湾气象局和美国地震学研究联合会为本文提供了地震目录和地震波形资料,审稿专家对本文提出了宝贵建议,作者在此一并表示感谢。
-
表 1 根据震源机制修正后的
$ M_{\rm{e}}^{{\rm{rev}}}$ 与未修正的Me的比较Table 1 Comparison of Me corrected for focal mechanisms with uncorrected Me
发震日期 MW $M_{\rm{e}}^{{\rm{IRIS}}} $ Me $M_{\rm{e}}^{{\rm{rev}}} $ 震源机制解 地点 2018-12-20 7.4 7.7 7.5 7.6 正断型 科曼多尔群岛 2019-09-29 6.7 6.9 6.9 6.9 逆冲型 智利中部沿海 表 2 台湾地区花莲县两次地震震源参数对比
Table 2 Source parameters of two earthquakes occurred in Hualien County,Taiwan region
发震日期 东经/° 北纬/° 震源深度/km 震源机制 MW Me MeIRIS 2018-02-04 121.72 24.20 7.8 逆断型 6.1 6.2±0.23 5.7 2019-04-18 121.69 23.99 20 逆断型 6.1 6.6±0.29 6.4 -
刘瑞丰,陈运泰,薛峰. 2018. 测定的震级之间不应相互换算[J]. 地震地磁观测与研究,39(3):1–9. doi: 10.3969/j.issn.1003-3246.2018.03.001 Liu R F,Chen Y T,Xue F. 2018. The measured magnitude should not be converted to each other[J]. Seismological and Geomagnetic Observation and Research,39(3):1–9 (in Chinese).
Boatwright J,Choy G L. 1986. Teleseismic estimates of the energy radiated by shallow earthquakes[J]. J Geophys Res:Solid Earth,91(B2):2095–2112. doi: 10.1029/JB091iB02p02095
Bormann P,Saul J. 2008. The new IASPEI standard broadband magnitude Mb[J]. Seismol Res Lett,79(5):698–705. doi: 10.1785/gssrl.79.5.698
Bormann P. 2009. New Manual of Seismological Observatory Practice (NMSOP-2)[M]. Potsdam: Deutsches GeoForschungszentrum GFZ: 124–133.
Bormann P,Di Giacomo D. 2010. The moment magnitude MW and the energy magnitude Me:Common roots and differences[J]. J Seismol,15(2):411–427. doi: 10.1007/s10950-010-9219-2
Choy G L,Boatwright J L. 1995. Global patterns of radiated seismic energy and apparent stress[J]. J Geophys Res:Solid Earth,100(B9):18205–18228. doi: 10.1029/95jb01969
Convers J A,Newman A V. 2011. Global evaluation of large earthquake energy from 1997 through mid‐2010[J]. J Geophys Res:Solid Earth,116(B8):B08304. doi: 10.1029/2010JB00728
Di Giacomo D,Grosser H,Parolai S,Bormann P,Wang R J. 2008. Rapid determination of Me for strong to great shallow earthquakes[J]. Geophys Res Lett,35(10):L10308. doi: 10.1029/2008GL033505
Di Giacomo D,Parolai S,Bormann P,Grosser H,Saul J,Wang R J,Zschau J. 2010. Suitability of rapid energy magnitude determinations for emergency response purposes[J]. Geophys J Int,180(1):361–374. doi: 10.1111/j.1365-246X.2009.04416.x
Gutenberg B,Richter C F. 1956. Earthquake magnitude,intensity,energy,and acceleration (Second paper)[J]. Bull Seismol Soc Am,46(2):105–145.
Kanamori H. 1977. The energy release in great earthquakes[J]. J Geophys Res,82(20):2981–2987. doi: 10.1029/JB082i020p02981
Lomax A. 2005. Rapid estimation of rupture extent for large earthquakes:Application to the 2004,M9 Sumatra-Andaman mega-thrust[J]. Geophys Res Lett,32(10):L10314. doi: 10.1029/2005gl022437
Picozzi M,Bindi D,Brondi P,Di Giacomo D,Parolai S,Zollo A. 2017. Rapid determination of P wave-based energy magnitude:Insights on source parameter scaling of the 2016 Central Italy earthquake sequence[J]. Geophys Res Lett,44(9):4036–4045. doi: 10.1002/2017gl073228
Stockwell R G,Mansinha L,Lowe R P. 1996. Localization of the complex spectrum:The S transform[J]. IEEE Trans Signal Process,44(4):998–1001. doi: 10.1109/78.492555
Venkataraman A,Kanamori H. 2004. Effect of directivity on estimates of radiated seismic energy[J]. J Geophys Res:Solid Earth,109(B4):B04301. doi: 10.1029/2003jb002548
Wang R J. 1999. A simple orthonormalization method for stable and efficient computation of Green’s functions[J]. Bull Seismol Soc Am,89(3):733–741.