Discrimination of earthquakes and explosions using multi-fractal singularity spectrums properties

A new method of discrimination of seismic records from earthquakes and explosions is proposed which is based on using properties of their multi-fractal singularity spectrums. The efficiency of the method is illustrated by analysing seismic records from the region of Aswan Dam in Egypt. Current pattern of seismicity in the Upper Egypt is composed by three types of records: tectonic earthquakes, reservoir-induced earthquakes and seismic events generated by quarry blasts. To discriminate quarry blasts from earthquakes of both natures, multi-fractal analysis of records were used. Singularity spectrum support width and multi-fractal generalised Hurst exponents were calculated for all seismic events in the selected data set from a given area. The linear Bayesian discriminator using these characteristics of seismic records provides correct classification for 93 % of earthquakes records and for 99 % of quarry blasts records.     

Site characterizations for the Tehran network (TDMMO) in Tehran Region Using Micro-Earthquake,Microtremor and Quarry Blast Data

We did a Study of Horizontal-to-Vertical Component Spectral Ratio in the Tehran seismic zone. Micro-earthquakes, microtremors and quarry blasts data were used as an estimation of the site response in the Tehran zone. Site effects were studied based on horizontal to vertical ratios by the Nakamura׳s technique. Also, we used the spectra of signals for three components with the lowest noise levels for spectral slope studies. The analysis used seismic events from a network of 13 seismic stations by the permanent local seismological network of the Tehran Disaster Mitigation and Management Organization (TDMMO) from 2004 to 2007. The number of events used were different for each station. Quarry blast events were with 1.2≤M_L≤2.2 and micro-earthquakes were with 1.1≤M_L≤4.1.By comparing results for earthquake, microtremor and quarry blast, we could see that there is a significant difference between them. The data showed clear observations, especially in high-frequencies. The H/V spectral ratios indicate dominant frequency for rock/soft site with a higher ratio level for quarry blast ratios, which are comparable to the earthquake results due to their difference sources. The results derived by spectral H/V ratios and spectral analysis may be used to distinguish between local earthquakes and quarry blasts.     

A GIS-based time-dependent seismic source modeling of Northern Iran

The first step in any seismic hazard study is the definition of seismogenic sources and the estimation of magnitude-frequency relationships for each source. There is as yet no standard methodology for source modeling and many researchers have worked on this topic. This study is an effort to define linear and area seismic sources for Northern Iran. The linear or fault sources are developed based on tectonic features and characteristic earthquakes while the area sources are developed based on spatial distribution of small to moderate earthquakes. Time-dependent recurrence relationships are developed for fault sources using renewal approach while time-independent frequency-magnitude relationships are proposed for area sources based on Poisson process. GIS functionalities are used in this study to introduce and incorporate spatial-temporal and geostatistical indices in delineating area seismic sources. The proposed methodology is used to model seismic sources for an area of about 500 by 400 square kilometers around Tehran. Previous researches and reports are studied to compile an earthquake/fault catalog that is as complete as possible. All events are transformed to uniform magnitude scale; duplicate events and dependent shocks are removed. Completeness and time distribution of the compiled catalog is taken into account. The proposed area and linear seismic sources in conjunction with defined recurrence relationships can be used to develop time-dependent probabilistic seismic hazard analysis of Northern Iran.     

结合天然地震推导水库地震的建筑物易损性

通过类比的方法,选择有现场调查资料并有建筑物破坏比结果的天然地震,基于地震中建筑物震害表现的一致性,以震级、震源深度和极震区烈度为标准,聚类得到与水库地震震害相似的天然地震震例。考虑到宏观经济水平与地震经济损失已有的统计模型,采用了人口密度、人均GDP和三产比例等3个宏观经济指标来反映不同地区的建筑物总体抗震水平,以加权海明(Hamming)距离来定义已知矩阵与目标矩阵的近似度,最终的加权结果即为待求地区的易损性矩阵。通过实际检验,所得的建筑物易损性关系能够反映水库地震的破坏特点,较中强天然地震的易损性要高。这种工作思路也可以用于水库地震其它方面的研究中     

历史强震震中精度统计特征及其对地震危险性研究的影响

以《中国历史强震目录(公元前23世纪——公元1911年)》为基础,将目录中给出的震中精度参数在分段时间和分段震级上进行统计分析,并给出各类精度地震的地理分布．研究结果表明, 58.4%的历史强震震中精度都在3类及以上,其震中不确定范围大于25 km;一般来说,年代越久远的历史地震,其震中精度越差,而且公元1500年以后的历史强震低精度地震所占比例仍然不小;MS6.0以上的历史强震, 3类、 4类和5类精度地震占了60.1 %;不同精度历史强震的地理分布各有特点,云南省1类、 2类精度地震最多,而新疆和西藏的历史强震绝大部分都是低精度地震事件．以震中精度表述的历史强震震中位置不确定性会对潜在震源区划分、近场区地震活动性研究以及地震构造环境评价等工作造成影响,因此在地震危险性研究与地震地质工作中,都需对低精度历史地震重新考证或适当取舍．      

The 2008 West Bohemia earthquake swarm in the light of the WEBNET network

A swarm of earthquakes of magnitudes up to M _L = 3.8 stroke the region of West Bohemia/Vogtland (border area between Czechia and Germany) in October 2008. It occurred in the Novy Kostel focal zone, where also all recent earthquake swarms (1985/1986, 1997, and 2000) took place, and was striking by a fast sequence of macroseismically observed earthquakes. We present the basic characteristics of this swarm based on the observations of a local network WEBNET (West Bohemia seismic network), which has been operated in the epicentral area, on the Czech territory. The swarm was recorded by 13 to 23 permanent and mobile WEBNET stations surrounding the swarm epicenters. In addition, a part of the swarm was also recorded by strong-motion accelerometers, which represent the first true accelerograms of the swarm earthquakes in the region. The peak ground acceleration reached 0.65 m/s². A comparison with previous earthquake swarms indicates that the total seismic moments released during the 1985/1986 and 2008 swarms are similar, of about 4E16 Nm, and that they represent the two largest swarms that occurred in the West Bohemia/ Vogtland region since the M _L = 5.0 swarm of 1908. Characteristic features of the 2008 swarm are its short duration (4 weeks) and rapidity and, consequently, the fastest seismic moment release compared to previous swarms. Up to 25,000 events in the magnitude range of 0.5 < M _L < 3.8 were detected using an automatic picker. A total of nine swarm phases can be distinguished in the swarm, five of them exceeding the magnitude level of 2.5. The magnitude–frequency distribution of the complete 2008 swarm activity shows a b value close to 1. The swarm hypocenters fall precisely on the same fault portion of the Novy Kostel focal zone that was activated by the 2000 swarm (M _L ≤ 3.2) in a depth interval from 6 to 11 km and also by the 1985/1986 swarm (M _L ≤ 4.6). The steeply dipping fault planes of the 2000 and 2008 swarms seem to be identical considering the location error of about 100 m. Furthermore, focal mechanisms of the 2008 swarm are identical with those of the 2000 swarm, both matching an average strike of 170° and dip of 80° of the activated fault segment. An overall upward migration of activity is observed with first events at the bottom and last events at the top of the of the activated fault patch. Similarities in the activated fault area and in the seismic moments released during the three largest recent swarms enable to estimate the seismic potential of the focal zone. If the whole segment of the fault plane was activated simultaneously, it would represent an earthquake of M _L ~5. This is in good agreement with the estimates of the maximum magnitudes of earthquakes that occurred in the West Bohemia/Vogtland region in the past.     

Evaluation of seismicity and seismic hazard parameters in Turkey and surrounding area using a new approach to the Gutenberg–Richter relation

In this study, the spatial distributions of seismicity and seismic hazard were assessed for Turkey and its surrounding area. For this purpose, earthquakes that occurred between 1964 and 2004 with magnitudes of M ≥ 4 were used in the region (30–42°N and 20–45°E). For the estimation of seismicity parameters and its mapping, Turkey and surrounding area are divided into 1,275 circular subregions. The b-value from the Gutenberg–Richter frequency–magnitude distributions is calculated by the classic way and the new alternative method both using the least-squares approach. The a-value in the Gutenberg–Richter frequency–magnitude distributions is taken as a constant value in the new alternative method. The b-values calculated by the new method were mapped. These results obtained from both methods are compared. The b-value shows different distributions along Turkey for both techniques. The b-values map prepared with new technique presents a better consistency with regional tectonics, earthquake activities, and epicenter distributions. Finally, the return period and occurrence hazard probability of M ≥ 6.5 earthquakes in 75 years were calculated by using the Poisson model for both techniques. The return period and occurrence hazard probability maps determined from both techniques showed a better consistency with each other. Moreover, maps of the occurrence hazard probability and return period showed better consistency with the b-parameter seismicity maps calculated from the new method. The occurrence hazard probability and return period of M ≥ 6.5 earthquakes were calculated as 90–99% and 5–10 years, respectively, from the Poisson model in the western part of the studying region.     

基于深度学习的地震震级分类

为了探索地震加速度时程记录的震级信息,训练卷积神经网络基于地震震级大小对地震记录进行分类,将K-NET和KiK-net中将近12万个地震记录作为样本,对其进行信息筛选和归一化,之后将地震加速度时程记录用作输入,训练卷积神经网络模型以M5.5为分类界限来区分大震和小震。结果显示,在训练集中基于该模型的分类准确率为93.6％,在测试集中的准确率为92.3％,具有良好的分类效果,这表明大震记录与小震记录之间存在一些根本的区别,即可通过地震动加速度时程记录获取一定的震级信息。      

The 20-second regional magnitude <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">S</Emphasis></Subscript>(20R) for the Russian Far East

A new modified magnitude scale M _S (20R) is elaborated. It permits us to extend the teleseismic magnitude scale M _S (20) to the regional epicenter distances. The data set used in this study contains digital records at 12 seismic stations of 392 earthquakes that occured in the northwest Pacific Ocean in the period of 1993–2008. The new scale is based on amplitudes of surface waves of a narrow range of the periods (16–25 s) close to the period of 20 s, for distances of 80–3000 km. The digital Butterworth filter is used for processing. On the basis of the found regional features concerning distance dependence for seismic wave attenuation, all the stations of the region have been subdivided into two groups, namely, “continental” and “island-arc.” For each group of stations, its own calibration function is proposed. Individual station corrections are used to compensate for the local features.     

Discrimination of earthquakes and quarry blasts in the eastern Black Sea region of Turkey

In recent years, a large number of quarry blasts have been detonated in the eastern Black Sea region. When these blasts are recorded by seismic stations, they contaminate the regional earthquake catalog. It is necessary to discriminate quarry blast records from the earthquake catalogs in order to determine the real seismicity of the region. Earthquakes and quarry blasts can be separated through different methods. These methods should be applied concurrently in order to safely distinguish these events. In this study, we discriminated quarry blasts from earthquakes in the eastern Black Sea region of Turkey. We used 186 seismic events recorded by the Karadeniz Technical University and Bogaziçi University Kandilli Observatory Earthquake Research Institute stations which are Trabzon, Espiye, Pazar, Borçka, Ayd?ntepe, and Gümü?hane between years of 2002 and 2010. For the discrimination of quarry blasts from earthquakes, we used both, statistical methods (calculation of the maximum ratio of S to P waves (S/P), complexity (C)) and spectral methods (spectrogram calculation). These methods included measuring the maximum amplitude S/P, C, spectral ratio, and time-frequency analysis. We especially relied on two-dimensional time-frequency analysis methods to discriminate quarry blasts from earthquakes in Turkey. As a result of this study, 68 % of the examined seismic events were determined to be quarry blasts and 32 % to be earthquakes. The earthquakes occurring on land are related to small faults and the blasts are concentrated in large quarries. Nearly 40 % of the earthquakes occurred in the Black Sea, most of them are related to the Black Sea thrust belt, where the largest earthquake was observed in the time period studied. The areas with the largest earthquake potential in the eastern Black Sea region are in the sea.     

四川地区强震前弱震空区的变化特征

本文使用四川台网地震目录全部资料,适当参考邻省地震目录,对于四川省从1970年下半年以来发生的15组5级以上地震(另有5组缺乏完整资料无法讨论),作 R-t 图、D-t 图寻找时间空区,并作相应的平面空区图.通过对比分析,发现6组6级以上和4组5级地震前震中附近出现空区现象.根据形成空区的变化过程可分为先扩(大)后缩(小)型和逐渐扩大的开放型两种类型.平面空区的位置和形状往往与构造展布有关.有5组5级震例,或因震前震中区小震不断,或因震中及其附近地区小震太少,看不出明显的空区.本文利用出现空区的10个震例,总结出空区参数与未来地震的震级的经验关系:M=5.50 lg Rmax-2.331.32 M=2.98 lgT+2.941.06 M=2.52lg(RmaxT)-0.510.96文中将上述结果与大竹政和的结果作了对比.空区不同变化形态可用金森博雄概括的突起点模式加以解释.      

Classification of seismic events using linear discriminant function (LDF) in the Sakarya region,Turkey

The Sakarya prefecture is an interesting area with various seismicity types. This activity comes from earthquakes occurring at the North Anatolian Fault Zone and from a few quarry blast areas in the region. These quarry blast recordings produce errors in the determination of active faults and mapping of the microearthquake activity. Therefore, to recognize the tectonic activity in the region, we need to be able to discriminate between earthquakes and quarry blasts in the catalogues. In this study, a statistical analysis method (linear discriminant function) has been applied to classify seismic events occurring in the Sakarya region. We used 110 seismic events that were recorded by Sakarya University Seismic Station between 2012 and 2014. Time and frequency variant parameters, maximum S wave and maximum P wave amplitude ratio (S/P), the spectral ratio (Sr), maximum frequency (f_max), and total signal duration of the waveform were used for discrimination analyses. The maximum frequency (f_max) versus time duration of the seismic signal gives a higher classification percentage (94%) than the other discriminants. At the end of this study, 41 out of 110 events (44%) are determined as quarry blasts, and 62 (56%) are considered as earthquakes.     

利用历史有感地震目录（M≥3（1／2））研究华北地区地震活动性特征

对新编的华北地区历史有感地震目录（３（１／２）≤Ｍ＜４（１／２））的精度和可信度进行了讨论。认为该地震目录从１４８４年后的记载率较高，它的发震时间、震中位置和震级的精度基本可以满足研究历史地震活动性的需要。它与中国地震目录（Ｍ≥４（３／４））比较，补充了相当数量的３（１／２）≤Ｍ≤４（１／２）的地震。利用历史有感地震目录分析了华北北部部分６级以上的历史地震，结果表明，这些地震在强震发生前具有明显的空区、条带、地震活动性增强、平静等前兆特征，同时，也显示了某些强＊发生前的前震活动和主震之后起伏的强余震活动     

1976年松潘大震前中强震的活动特征

1976年8月四川松潘,平武发生了强烈地震,震级分别为7.2,6.7和7.2,这是我国南北地震带中段自1933年以来的最大地震。1933—1976年间该段上发生了一系列5.0—6.5级中强地震,本文讨论了这些中强震在震中跳迁,能量释放,大震源区小震频度,震中分布以及中强震余震的强度衰减等方面的地震活动特征。指出1973—1974年间的中强震活动,应当是1976年松潘平武地震的前兆地震活动,並且依据上述特征有可能将这类中强震从一般的中强震中区分出来。主要的活动特征可以归结为:自1972年开始,地震能量释放一反常态由原来源区和外围区的同步成比例变化改为截然反向的变化,这些中强震余震的扩散不明显,相反,主震却很清楚;这几个中强震和他们自已的最大强余震之间的时间间隔都大于4天,这明显大于附近地区其他中强震的这一间隔,也大于国内近年来大多数5.0—6.0级中强震的这一时间隔;1976年初,即大震前半年,震源区的地震活动异常平静,一直持续到大震发生。此外结合该区历史强震的迁移特点有可能预示1976年松潘平武地震的时间和地点,而其孕育过程似乎可由1949,1972及1976年初分为几个阶段。     

The 1992 Tafilalt seismic crisis (Anti-Atlas, Morocco)

The Tafilalt region, located at the eastern end of the Anti-Atlas chain in Morocco, was shaken on 23 and 30 October 1992 by two moderate earthquakes of magnitude mb ∼ 5 and intensity ∼ VI MSK64, which caused two deaths and great damage in the area between Erfoud and Rissani. The review of data available on the seismic crisis allowed us to improve the knowledge on the macroseismic, instrumental and source parameters of the earthquakes. The main results of the present study are: (1) location of the epicentres with the help of data from a close portable network allowed us to propose new epicentral coordinates at 31.361° N, 4.182° W (23 October) and 31.286° N, 4.347° W (30 October); both events have focal depths of 2 km; (2) the shock of 30 October was followed by a series of 305 aftershocks, most of which were located west of Rissani; the 61 best-constrained events had focal depths of 5 to 19 km and magnitudes 0.7 to 3; (3) the largest damage was located in an area between the two epicentres within the Tafilalt valley and was probably amplified by site effects due to the proximity of the water table within the Quaternary sediments; (4) focal mechanisms of the main events correspond to strike-slip faulting with fault planes oriented N–S (left lateral) and E–W (right lateral); the only mechanism available for the aftershocks also corresponds to strike-slip faulting; (5) spectral analysis shows that the scalar seismic moment (Mo) of the first event is slightly larger than the second; the corresponding values of Mw are 5.1 and 5.0, respectively; (6) the dimensions of the faults for a circular fault model are 7.7 ± 1.4 and 7.4 ± 1.2 km, respectively; the average displacement is 4 cm for the first event and 3.7 cm for the second; the stress drop is 0.4 and 0.3 MPa, respectively, in agreement with standard values; (7) the Coulomb Stress test performed for both earthquakes suggests a relationship between both events only when the used location is at the limit of the horizontal uncertainty; (8) finally, the occurrence of these shocks suggests that the Anti-Atlas is undergoing tectonic deformation in addition to thermal uplift as suggested by recent publications.     

2017年9月临城震群遗漏地震检测及发震构造初探

以河北测震台网记录的连续波形数据为基础,采用模板匹配滤波方法,检测2017年9月4日—6日河北临城震群遗漏地震。利用波形互相关扫描,检测到19次被测震台网常规分析遗漏的地震,震级范围为M_(L )0.0—1.2,并标定新检测地震事件的P波和S波到时,估计其震中位置及震级。通过地震精定位和最大地震震源机制解,认为此次临城地震序列的主要发震构造应为近EW向活动断层。     

Spectral classification methods in monitoring small local events by the Israel seismic network

We use the dense Israel Seismic Network (ISN) to discriminate between low magnitude earthquakes and explosions in the Middle East region. This issue is important for CTBT monitoring, especially when considering small nuclear tests which may be conducted under evasive conditions. We explore the performance of efficient discriminants based on spectral features of seismograms using waveforms of 50 earthquakes and 114 quarry and underwater blasts with magnitudes 1.0–2.8, recorded by ISN short-period stations at distances up to 200 km. The single-station spectral ratio of the low and high-frequency seismic energy shows an overlap between explosions and earthquakes. After averaging over a subnet of stations, the resolving power is enhanced and the two classes of events are separated. Different frequency bands were tested; the (1–3 Hz)/(6–8 Hz) ratio provided the best discriminant performance. We also estimated normalized r.m.s. spectral amplitudes in several sequential equal frequency windows within the 1–12 Hz band and applied multiparametric automatic classification procedures (Linear Discrimination Function and Artificial Neural Network) to the amplitudes averaged over a subnetwork. A leave-one-out test showed a low rate of error for the multiparametric procedures. An innovative multi-station discriminant is proposed, based on spectral modulation associated with ripple-firing in quarry blasts and with the bubbling effect in underwater explosions. It utilizes a distinct azimuth-invariant coherency of spectral shapes for different stations in the frequency range (1–12 Hz). The coherency is measured by semblance statistics commonly used in seismic prospecting for phase correlation in the time domain. After modification, the statistics applied to the network spectra provided event separation. A new feature of all the above mentioned procedures is that they are based on smoothed (0.5 Hz window), instrument-corrected FFT spectra of the whole signal; they are robust to the accuracy of onset time estimation and, thus well suited to automatic event identification.     

地震定位方法适用性分析

选取辽宁区域三种不同类型的六次地震事件,运用单纯型、Hyposat、HYP2000三种定位方法对其进行重新定位,产出的结果进行对比分析。结果表明,Hyposat方法在定位这三种类型的地震时,残差范围最小最为稳定,定位网内地震时震中差也相对比较小,较适用于内陆地震和矿震的地震定位。     

Analysis and interpretation of the October 21, 1766 earthquake in the Southeastern Caribbean

The October 21, 1766 earthquake is the most widely felt event in the seismic history of Trinidad and Venezuela. Previous works diverged on the interpretation of the historical data available for this event. They associated the earthquake either with the Lesser Antilles subduction zone, with strike-slip motion along El Pilar fault, or with intraplate deformation at the edge of Guyana shield. Isoseismal areas are proposed after a new search and analysis of primary and secondary sources of historical information. Two of the largest earthquakes of the twentieth century which occurred in the region, the 1968 (M _S 6.4, h = 103 km), and the 1997 (M _W 6.9, h = 25 km) events, for which both intensity data and instrumentally determined source parameters are available, are used to calibrate the isoseismal areas and to interpret them in terms of source depth and magnitude. It is concluded that the large extent of intensity values higher than V is diagnostic of the depth (85 ± 20 km) of the 1766 source, and of local amplifications of ground motion due to soft soil conditions and to strong contrasts of impedance at the edge of Guyana shield. It is proposed that the event occurred either in slab, or close to the bottom lithospheric interface between the Caribbean and South American plates (∼11°N; ∼62.5°W). The value of the magnitude is estimated at 6.5 < M _S < 7.5 depending on the source depth and on the decay of ground motion as a function of distance. Deep and intermediate depth earthquakes can induce important casualties in Trinidad, Venezuela, and Guyana, possibly more damaging than those induced by shallower earthquakes along the strike of El Pilar Fault.     

Evolution of the seismic activity from the 1982 eruption of El Chichon Volcano, Chiapas, Mexico

 We analyzed more than 1700 earthquakes related to the 1982 eruption of El Chichon volcano in southern Mexico. The data were recorded at specific periods throughout the whole eruptive interval of March to April 1982, by three different networks. The seismic activity began several months before the first eruption on 28 March. During this period the seismicity consisted of hybrid and long-period shallow earthquakes most likely related to processes of faulting, fracturing, and fluid movement underneath the volcano. The foci of events occurring before the eruption circumscribe an aseismic zone from approximately 7 to 13 km below the volcano. After the eruption, the seismic activity consisted of tectonic-type earthquakes that peaked at 1200 events/h. This later activity occurred over a wide range of depths, mostly between 5 and 20 km, that includes the former aseismic zone and is roughly limited by the major tectonic faults in the area. Received: 19 May 1998 / Accepted: 13 June 1999