Comparison between earthquake magnitudes determined by China seismograph network and US seismograph networks (I): Body wave magnitude

Introduction Gutenberg (1945a, b) introduced body wave magnitude based on P, PP and S waves (with a period of 0.5~12.0 s) of teleseismic events. Body wave magnitude includes mb determined with short-period seismograph and mB determined with middle- and long-period seismographs. Some-times it is written as m, which is referred to as unified earthquake magnitude. mb represents earth-quake magnitude measured with body wave amplitude around 1 s, while mB represents earthquake magnitude measured …     

Comparison between earthquake magnitudes determined by China seismograph network and US seismograph network （Ⅱ）： Surface wave magnitude

By using orthogonal regression method, a systematic comparison is made between surface wave magnitudes determined by Institute of Geophysics of China Earthquake Administration （IGCEA） and National Earthquake Information Center of US Geological Survey （USGS/NEIC） on the basis of observation data collected by the two institutions between 1983 and 2004. A formula is obtained which reveals the relationship between surface wave magnitudes determined by China seismograph network and US seismograph network. The result shows that, as different calculation formulae and observational instruments are used, surface wave magnitude determined by IGCEA is generally greater by 0.2 than that determined by NEIC： for M=3.5-4.5 earthquakes, it is greater by 0.3; for M=5.0-6.5 earthquakes, it is greater by 0.2; and for M≥7.0 earthquakes, it is greater by no more than 0.1.     

Empirical Global Relations Converting M S and m b to Moment Magnitude

The existence of several magnitude scales used by seismological centers all over the world and the compilation of earthquake catalogs by many authors have rendered globally valid relations connecting magnitude scales a necessity. This would allow the creation of a homogeneous global earthquake catalog, a useful tool for earthquake research. Of special interest is the definition of global relations converting different magnitude scales to the most reliable and useful scale of magnitude, the moment magnitude, M _W. In order to accomplish this, a very large sample of data from international seismological sources (ISC, NEIC, HRVD, etc.) has been collected and processed. The magnitude scales tested against M _W are the surface wave magnitude, M _S, the body wave magnitude, m _b, and the local magnitude, M _L. The moment magnitudes adopted have been taken from the CMT solutions of HRVD and USGS. The data set used in this study contains 20,407 earthquakes, which occurred all over the world during the time period 1.1.1976–31.5.2003, for which moment magnitudes are available. It is shown that well-defined relations hold between M _W and m _b and M _S and that these relations can be reliably used for compiling homogeneous, with respect to magnitude, earthquake catalogs.     

Source process of the 14 November 2001 western Kunlun Mountain M S=8.1 earthquake

Based on digital teleseismic P-wave seismograms recorded by 28 long-period seismograph stations of the global seismic network, source process of the November 14, 2001 western Kunlun Mountain M _S=8.1 (M _W=7.8) earthquake is estimated by a new inversion method. The result shows that the earthquake is a very complex rupture event. The source rupture initiated at the hypocenter (35.95°N, 90.54°E, focal depth 10 km, by USGS NEIC), and propagated to the west at first. Then, in several minutes to a hundred minutes and over a large spatial range, several rupture growth points emerged in succession at the eastern end and in the central part of the finite fault. And then the source rupture propagated from these rupture growth points successively and, finally, stopped in the area within 50 km to the east of the centroid position (35.80°N, 92.91°E, focal depth 15 km, by Harvard CMT). The entire rupture lasted for 142 s, and the source process could be roughly separated into three stages: The first stage started at the 0 s and ended at the 52 s, lasting for 52 s and releasing approximately 24.4% of the total moment; The second stage started at the 55 s and ended at the 113 s, lasting for 58 s and releasing approximately 56.5% of the total moment; The third stage started at the 122 s and ended at the 142 s, lasting for 20 s and releasing approximately 19.1% of the total moment. The length of the ruptured fault plane is about 490 km. The maximum width of the ruptured fault plane is about 45 km. The rupture mainly occurred within 30 km in depth under the surface of the Earth. The average static slip in the underground rocky crust is about 1.2 m with the maximum static slip 3.6 m. The average static stress drop is about 5 MPa with the maximum static stress drop 18 MPa. The maximum static slip and the maximum stress drop occurred in an area within 50 km to the east of the centroid position. Foundation item: Joint Seismological Science Foundation of China (103066) and Foundation of the Seismic Pattern and Digital Seismic Data Application Research Office of Institute of Earthquake Science of the China Earthquake Administration.     

中国海域及邻区基于BP神经网络的矩震级转换模型

海域地震对我国海洋资源开发和沿海地区的经济建设形成严重威胁,开展相关地震活动性研究的重要前提是编译我国海域及邻区的地震目录。我国常用的震级标度为面波震级(MS),而国际上最新的地震活动模型多采用矩震级(MW),因此在应用这些模型时需要拟合面波震级与矩震级之间的转换关系。本文以中国海域及邻区为研究区,收集了1988—2020年中国地震台网的面波震级和全球矩心矩张量(GCMT)项目的矩震级数据,从中提取年份、深度、经度、纬度、面波震级作为影响因子,以实际记录的矩震级值作为标记,训练BP神经网络建立以GCMT的矩震级为目标的震级转换模型。同时,使用最小二乘回归和正交回归建立线性模型作为对比。结果显示,最小二乘回归和正交回归的平均绝对误差和均方根误差比BP神经网络高40%左右。此外,BP神经网络的残差绝对值更小、分布更集中。     

Comparison between earthquake magnitudes determined by China seismograph network and US seismograph network (Ⅱ):Surface wave magnitude

By using orthogonal regression method, a systematic comparison is made between surface wave magnitudes determined by Institute of Geophysics of China Earthquake Administration (IGCEA) and National Earthquake Information Center of US Geological Survey (USGS/NEIC) on the basis of observation data collected by the two institutions between 1983 and 2004. A formula is obtained which reveals the relationship between surface wave magnitudes determined by China seismograph network and US seismograph network. The result shows that, as different calculation formulae and observational instruments are used, surface wave magnitude determined by IGCEA is generally greater by 0.2 than that determined by NEIC: for M=3.5～4.5 earthquakes, it is greater by 0.3;for M=5.0～6.5 earthquakes, it is greater by 0.2;and for M≥7.0 earthquakes, it is greater by no more than 0.1.     

The December 26, 2004, off the west coast of northern Sumatra, Indonesia, M W=9.0, earthquake and the critical-point-like model of earthquake preparation

Long-term seismic activity prior to the December 26, 2004, off the west coast of northern Sumatra, Indonesia, M _W=9.0 earthquake was investigated using the Harvard CMT catalogue. It is observed that before this great earthquake, there exists an accelerating moment release (AMR) process with the temporal scale of a quarter century and the spatial scale of 1 500 km. Within this spatial range, the M _W=9.0 event falls into the piece-wise power-law-like frequency-magnitude distribution. Therefore, in the perspective of the critical-point-like model of earthquake preparation, the failure to forecast/predict the approaching and/or the size of this earthquake is not due to the physically intrinsic unpredictability of earthquakes. Foundation item: Ministry of Science and Technology Project (2004CB418406). Contribution No. 05FE3010, Institute of Geophysics, China Earthquake Administration.     

ML of the 1989 Loma Prieta earthquake

The use of regional attenuation in computing the local magnitude, M_L, from strong motion data gathered at distances less than 100 km may lead to systematic underestimates approaching 0·5 magnitude units (Trifunac & Herak, Soil Dynamics and Earthquake Engineering, 1992, 18, 229-41). The use of the attenuation law Att(Δ), for example, with synthetic estimates of Wood-Anderson seismometer response, during the Loma Preita earthquake, leads to estimates of M_L which agree with the surface wave and moment magnitudes, and which are essentially distance-independent.     

2001年11月14日昆仑山库赛湖地震(M_S8.1)地表破裂带的基本特征

震后野外考察表明 ,2 0 0 1年 11月 14日昆仑山库赛湖地震 (MS8 1)发生在青藏高原北部东昆仑断裂带库赛湖段上 ,发震断层具有高速率左旋滑动的基本特征 ,晚更新世晚期以来的平均滑动速率达 (14 8± 2 4 )mm/a ;地震地表破裂带沿库赛湖段西起布喀达板峰东缘 (91°0 8′E) ,向东经库赛湖北缘、青藏公路 2 894里程碑、玉珠峰南麓 ,东止于青藏公路东 70km附近 (94°4 8′E) ,地震地表破裂带沿N70°～ 90°W走向线状展布 ,全长约 35 0km ,由一系列走向N4 5°～ 5 0°E拉开状张裂缝、走向N6 0°～ 75°E张剪切裂缝、走向N80°W剪切裂缝以及隆起鼓包或开裂陷坑等斜列状组合而成 ,显示出纯剪切走滑的破裂特征 ,最大左旋水平位移 6m ;宏观震中位于昆仑山口西 80～ 90km附近的库赛湖东北角山麓地带 ,地震地表破裂带宽度 30 0m ,在库赛湖北岸至山麓地带的地震地表破裂带和由地震动或重力效应引起的次生破裂带总宽度可达 2km。库赛湖地震地表破裂的左旋走滑特征表明 ,青藏高原物质确实存在着向东的滑移或流动 ,东昆仑断裂带东部与库赛湖段斜列的东大     

Source time functions of the 1999, Jiji (Chi-Chi) earthquake from GDSN long period waveform data using aftershocks as empirical Green’s functions

A large earthquake (M _W=7.6) occurred in Jiji (Chi-Chi), Taiwan, China on September 20, 1999, and was followed by many moderate-size shocks in the following days. Two of the largest aftershocks with the magnitudes of M _W=6.1 and M _W=6.2, respectively, were used as empirical Green’s functions (EGFs) to obtain the source time functions (STFs) of the main shock from long-period waveform data of the Global Digital Seismograph Network (GDSN) including IRIS, GEOSCOPE and CDSN. For the M _W=6.1 aftershock of September 22, there were 97 pairs of phases clear enough from 78 recordings of 26 stations; for the M _W=6.2 aftershock of September 25, there were 81 pairs of phases clear enough from 72 recordings of 24 stations. For each station, 2 types of STFs were retrieved, which are called P-STF and S-STF due to being from P and S phases, respectively. Totally, 178 STF individuals were obtained for source-process analysis of the main shock. It was noticed that, in general, STFs from most of the stations had similarities except that those in special azimuths looked different or odd due to the mechanism difference between the main shock and the aftershocks; and in detail, the shapes of the STFs varied with azimuth. Both of them reflected the stability and reliability of the retrieved STFs. The comprehensive analysis of those STFs suggested that this event consisted of two sub-events, the total duration time was about 26 s, and on the average, the second event was about 7 s later than the first one, and the moment-rate amplitude of the first event was about 15% larger than that of the second one. Foundation item: State Natural Science Foundation of China (49904004) and IPGP of France. Contribution No. 02FE2007, Institute of Geophysics, China Seismological Bureau.     

Near-realtime source analysis of the 20 March 2012 Ometepec-Pinotepa Nacional,Mexico earthquake

We apply a single-step, finite-fault analysis procedure to derive a coseismic slip model for the large M_W 7.4 Ometepec-Pinotepa Nacional, Mexico earthquake of 20 March 2012, using teleseismic P waveforms recorded by the Global Seismographic Network. The inversion is conducted in near-realtime using source parameters available from the USGS/NEIC and the Global Centroid Moment Tensor (gCMT) project. The fault orientation and slip angle are obtained from the gCMT mechanism assuming that the fault coincides with the shallow-dipping nodal plane. The fault dimensions and maximum rise time are based on the magnitude reported for the event. Teleseismic data from the USGS/NEIC Continuous Waveform Buffer database are used in the inversion with record start times set to the P-wave arrivals used to compute the earthquake hypocenter. The inversion is stabilized by requiring a smooth transition of slip across the fault while minimizing the seismic moment. These constraints are applied using a smoothing weight that is estimated from the inverse problem, allowing the recovery of the least-complicated rupture history in a single step. Inversion of the deconvolved, ground-displacement waveforms reveals a simple, circular rupture similar in extent to the source identified by the USGS/NEIC using body-and surface-wave data, indicating that the teleseismic P waves can provide a first-order source model for the event in near-realtime. Additional inversions conducted using velocity records identify a more-detailed rupture model characterized by an elliptical 2500 km² source region extending updip and downdip from the hypocenter. This elliptical source preserves the orientation and overall dimensions of a dual-source slip model obtained recently by other investigators using local strong motions and global seismic waveforms. The results indicate that velocity waveforms could provide additional details of the earthquake rupture in near-realtime, finite-fault inversions using teleseismic P waves.     

Spatial damage distribution of August 16, 2003, Inner Mongolia, China, M S=5.9 earthquake and analysis

The spatial damage distribution of August 16, 2003, Inner Mongolia, China, M _S=5.9 earthquake is summarized through field investigation. The moment tensor solution and focal mechanism are inverted using the digital long-period waveform records of China Digital Seismograph Network (CDSN). The relation between the spatial damage distribution and focal mechanism is analyzed according to the focal mechanism, the aftershock distribution and the spatial damage distribution. The possible relation between the characteristics of ground motion and the tectonic background of the source region is discussed in terms of the global ground motion records, historical earthquake documents and the damage distribution. Investigation reveals that the meizoseismal region is in east-west direction, which is consistent with the nodal plane of focal mechanism inversion. The meizoseismal area is relatively large and the damage of single-story adobe houses or masonry houses is more severe. This may have relations with local seismotectonic environment. Foundation item: The Special Funds for Major State Basic Research Project (2002CB412706). Contribution No. 05FE3009, Institute of Geophysics, China Earthquake Administration.     

内蒙古中西部地区中小地震矩震级研究

基于S震相"S窗"内的波形信号识别、品质因子Q(f)和22个台站场地响应,利用2009~2016年3月内蒙古中西部地区地震的波形资料,反演了182次中小地震的震源波谱参数,得到这些小震的零频幅值及其拐角频率,据此计算了这些地震的地震矩M_0、矩震级M_W和应力降Δσ。利用回归分析方法得到了近震震级与矩震级、矩震级与应力降的关系式。分析表明,近震震级与矩震级、矩震级与应力降呈线性关系。可见,将矩震级纳入地震的快报与正式目录中,可以丰富地震观测报告内容,更好地为地震应急和地震科研服务。     

Comparison between earthquake magnitudes determined by China seismograph network and US seismograph networks (Ⅰ): Body wave magnitude

By using orthogonal regression method, a systematic comparison is made between body wave magnitudes determined by Institute of Geophysics of China Earthquake Administration (IGCEA) and National Earthquake Information Center of US Geological Survey (USGS/NEIC) on the basis of observation data from China and US seismograph networks between 1983 and 2004. The result of orthogonal regression shows no systematic error between body wave magnitude mb determined by IGCEA and mb (NEIC). Provided that mb (NEIC) is taken as the benchmark, body wave magnitude determined by IGCEA is greater by 0.2～0.1 than the magnitude determined by NEIC for M=3.5～4.5 earthquakes; for M=5.0～5.5 earthquakes, there is no difference; and for M≥6.0 earthquakes, it is smaller by no more than 0.2. This is consistent with the result of comparison by IDC (International Data Center).     

Inversion of long-period body-wave data for the source process of the Gonghe, Qinghai, China earthquake

Long period body waves data recorded by the China Digital Seismograph Network (CDSN) are inverted for the seismic moment tensors of the April 26, 1990, Gonghe, QinghaiM _S=6.9 earthquake and itsM _S=5.0 after-shock occurred on May 7, 1990. In the inversion, the generalized reflection-transmission coefficient matrix method is used to generate Green’s function. From the inversion it is obtained that the rupture process of theM _S=5.0 aftershock is relatively simple, and that of the main shock is rather complex. There are at least two events during main shock rupture process with an interval about 35 seconds. The focal mechanisms of two events are roughly the same as that of the aftershock, all of them were mainly reverse dip-slipping faulting with minor left-lateral strike-slip motion. These results indicate that the Gonghe earthquake was the result of the farther extension of one NWW-SEE striking buried fault on the southern margin of Gonghe basin from shallower depth to deeper depth and from NW to SE under the action of a nearly horizontal NE direction compressive stress. Contribution No. 95A0111, Institute of Geophysics, SSB, China.     

Comparison between different earthquake magnitudes determined by China Seismograph Network

By linear regression and orthogonal regression methods, comparisons are made between different magnitudes (lo-cal magnitude ML, surface wave magnitudes MS and MS7, long-period body wave magnitude mB and short-period body wave magnitude mb) determined by Institute of Geophysics, China Earthquake Administration, on the basis of observation data collected by China Seismograph Network between 1983 and 2004. Empirical relations between different magnitudes have been obtained. The result shows that: 1 As different magnitude scales reflect radiated energy by seismic waves within different periods, earthquake magnitudes can be described more objectively by using different scales for earthquakes of different magnitudes. When the epicentral distance is less than 1 000 km, local magnitude ML can be a preferable scale; In case M<4.5, there is little difference between the magnitude scales; In case 4.5MS, i.e., MS underestimates magnitudes of such events, therefore, mB can be a better choice; In case M>6.0, MS>mB>mb, both mB and mb underestimate the magnitudes, so MS is a preferable scale for deter-mining magnitudes of such events (6.08.5, a saturation phenomenon appears in MS, which cannot give an accurate reflection of the magnitudes of such large events; 2 In China, when the epicentral distance is less than 1 000 km, there is almost no difference between ML and MS, and thus there is no need to convert be-tween the two magnitudes in practice; 3 Although MS and MS7 are both surface wave magnitudes, MS is in general greater than MS7 by 0.2~0.3 magnitude, because different instruments and calculation formulae are used; 4 mB is almost equal to mb for earthquakes around mB4.0, but mB is larger than mb for those of mB≥4.5, because the periods of seismic waves used for measuring mB and mb are different though the calculation formulae are the same.     

Magnetic perturbations before the strong earthquakes

The search for the magnetic precursors of earthquakes by the method of synchronous detection is described. The data of the Guam Observatory (13.6°N, 144.9°E) located in a seismically active region and the USGS/NEIC catalogue of the National Earthquake Information Center of the U.S. Geological Survey for the period from 1991 to 2009 are used. Earthquakes with magnitudes M ≥ 7 in the range of longitudes from 100° to 170°E are analyzed. The intervals of accumulation cover 40 hours before and 40 hours after an earthquake. The preliminary result reveals an enhancement of fluctuations in the Z-component within 12 hours before the earthquake.     

Precise determination of focal parameters for July 20, 1995 M_L=4.1 earthquake sequence in the Huailai basin

IntroductionSince the late 1970s, the quickly developed global digital seismograph network has been providing high quality recordings of large earthquakes in global scale, based on which digital seismology has made great progress. Compared with large earthquakes, moderate and small sized shocks have more frequent occurrence, and comprise clues to geological tectonics and tectonic stress field in a region. Preceding and following a large earthquake, usually occur numbers of small events that im…     

Updated earthquake catalogue for seismic hazard analysis in Pakistan

A reliable and homogenized earthquake catalogue is essential for seismic hazard assessment in any area. This article describes the compilation and processing of an updated earthquake catalogue for Pakistan. The earthquake catalogue compiled in this study for the region (quadrangle bounded by the geographical limits 40–83° N and 20–40° E) includes 36,563 earthquake events, which are reported as 4.0–8.3 moment magnitude (M_W) and span from 25 AD to 2016. Relationships are developed between the moment magnitude and body, and surface wave magnitude scales to unify the catalogue in terms of magnitude M_W. The catalogue includes earthquakes from Pakistan and neighbouring countries to minimize the effects of geopolitical boundaries in seismic hazard assessment studies. Earthquakes reported by local and international agencies as well as individual catalogues are included. The proposed catalogue is further used to obtain magnitude of completeness after removal of dependent events by using four different algorithms. Finally, seismicity parameters of the seismic sources are reported, and recommendations are made for seismic hazard assessment studies in Pakistan.     

芦山地震和九寨沟地震台站实测震级差异性分析

本文选取2013年芦山地震和2017年九寨沟地震波形,重新量取垂直向振幅,计算宽频带面波震级M_S（BB）,分析各台站实测震级出现方向性差异的原因。其中,通过572个宽频带台站实测芦山地震震级M_S（BB）7.1,通过603个宽频带台站实测九寨沟地震震级M_S（BB）6.9。芦山地震实测震级大于M_S（BB）7.3的台站呈现WN-ES向分布,与断层倾向一致;实测震级小于M_S（BB）7.0的台站呈现NE-WS向分布,与其所在断层走向一致。九寨沟地震实测震级大于M_S（BB）7.0的台站分布呈现NE向分布,与断层倾向一致;实测震级小于M_S（BB）6.8的台站总体分布较为离散,大体呈现NW-SE向分布,与树正断裂走向一致。实测震级偏大的台站方向性分布与多普勒效应和P波辐射花样联系不明显。对比分析芦山地震和九寨沟地震,去除场地响应和仪器自身影响,台站实测震级差异性仍然存在,因此,台站实测震级差异性是由于受到了多普勒效应、辐射花样、仪器和场地响应之外的因素影响。综合考虑地震震级涉及的影响因素,芦山地震和九寨沟地震的台站实测震级差异性可能与地震波的传播路径有关。