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.     

A discussion on the application ofb-value to the prediction of seismic tendency

In the estimation of seismic tendency, using Gutenberg-Richter’sb-value and using Hurst exponent are two commonly used methods. Based on the fractal geometry of earthquake time series, we point out that these two methods correlate to each other. In the perspective of fractional Brownian motion (FBM), an earthquake sequence withb>3/4 and that withb<3/4 have different dynamic properties. Foundation item: The MOST Project G1998040705. Contribution No. 01FE2010, Institute of Geophysics, China Seismological Bureau.     

Frequency-Magnitude Relationships for seismic areas around the capital cities of Solomon, Vanuatu, Tonga and Fiji Islands

We studied the relationships between the frequency ofoccurrence and the magnitudes in seismogenic areasthat will affect four capital cities in the SouthPacific, namely: Honiara in Solomon Islands, Port Vilain Vanuatu, Suva in Fiji and Nuku'alofa in Tonga. Weused the NEIS catalogues for the period 1973–1997, formagnitudes greater than about 5, in the proximity tothe city under investigation. The definitions of thegeographic boundaries of the seismogenic zones arebased on the classical concepts of the distribution ofthe tectonic plates, the overwhelming number ofstudies that describe the seismotectonics in thoseactive regions and on the observed seismicity andstudies of the local people. Completeness of thecatalogue is assumed because of the high seismicity inthese areas, despite the relatively short time span.We have used a hybrid process combining least squaresfitting and Newtonian search process to find the bestfit of the statistical parameters. The characteristicb value is 1.27, and seems to be irrespective of thedepth and/or region.     

2016年新疆呼图壁6.2级地震前b值异常特征研究

利用2009年以来中国地震台网中心地震编目资料,以2016年新疆呼图壁6.2级地震作为研究对象,采用最大似然法进行b值空间扫描计算,获取震前震中及其邻区地震b值的空间图像。研究结果显示,2016年呼图壁6.2级地震发生在震前显著低b值区域,该异常特征可能反映了地震孕育的应力积累过程,印证了b值的物理意义。     

A fracture-mechanical cellular automaton model of seismicity

We present a cellular automaton model which simulates the process of seismogenesis using rules for evolution which are derived from the field of fracture mechanics, and include an interplay of positive and negative feedbacks. We describe the implementation of this model, and its analysis, in a massively parallel environment using the Connection Machine. Starting from a lattice with a fractal distribution of fracture toughnesses, theb value evolves in a way which closely mimics both the evolutions ofb value observed in the laboratory and derived from earthquake catalogues, reaching a broad and irregular maximum in the period preceding a major event, and declining rapidly during catastrophic failure. We conclude that the processes modelled are a reasonable representation of those occurring in Nature, and that the cellular automaton paradigm is a valuable way of simulating these processes on a large scale in an economical manner.     

Anomalously high b-values in the South Flank of Kilauea volcano, Hawaii: evidence for the distribution of magma below Kilauea's East rift zone

The pattern of b-value of the frequency–magnitude relation, or mean magnitude, varies little in the Kaoiki-Hilea area of Hawaii, and the b-values are normal, with b=0.8 in the top 10 km and somewhat lower values below that depth. We interpret the Kaoiki-Hilea area as relatively stable, normal Hawaiian crust. In contrast, the b-values beneath Kilauea's South Flank are anomalously high (b=1.3–1.7) at depths between 4 and 8 km, with the highest values near the East Rift zone, but extending 5–8 km away from the rift. Also, the anomalously high b-values vary along strike, parallel to the rift zone. The highest b-values are observed near Hiiaka and Pauahi craters at the bend in the rift, the next highest are near Makaopuhi and also near Puu Kaliu. The mildest anomalies occur adjacent to the central section of the rift. The locations of the three major and two minor b-value anomalies correspond to places where shallow magma reservoirs have been proposed based on analyses of seismicity, geodetic data and differentiated lava chemistry. The existence of the magma reservoirs is also supported by magnetic anomalies, which may be areas of dike concentration, and self-potential anomalies, which are areas of thermal upwelling above a hot source. The simplest explanation of these anomalously high b-values is that they are due to the presence of active magma bodies beneath the East Rift zone at depths down to 8 km. In other volcanoes, anomalously high b-values correlate with volumes adjacent to active magma chambers. This supports a model of a magma body beneath the East Rift zone, which may widen and thin along strike, and which may reach 8 km depth and extend from Kilauea's summit to a distance of at least 40 km down rift. The anomalously high b-values at the center of the South Flank, several kilometers away from the rift, may be explained by unusually high pore pressure throughout the South Flank, or by anomalously strong heterogeneity due to extensive cracking, or by both phenomena. The major b-value anomalies are located SSE of their parent reservoirs, in the direction of motion of the flank, suggesting that magma reservoirs leave an imprint in the mobile flank. We hypothesize that the extensive cracking may have been acquired when the anomalous parts of the South Flank, now several kilometers distant from the rift zone, were generated at the rift zone near persistent reservoirs. Since their generation, these volumes may have moved seaward, away from the rift, but earthquakes occurring in them still use the preexisting complex crack distribution. Along the decollement plane at 10 km depth, the b-values are exceptionally low (b=0.5), suggesting faulting in a more homogeneous medium.     

Space-time variations of the frequency-energy relation for mining-induced seismicity in the Ostrava-Karviná mining district

The values of parameters of the relation logN=a–b logE characterize in general the level of seismicity of individual regions. In the present paper, a detailed analysis ofb values of the frequency-energy distributions was carried out with the aid of a database of seismic events recorded by the seismological networks in the Ostrava-Karviná Coal Mine District (Czech Republic), using the least squares regression and maximum likelihood method, as well. The determination ofb value was performed for a relatively large number of regions investigated and for different time series of observations. Special attention has been paid to induced seismic events statistics for the vicinity of one of the coalfaces in the Lazy Mine. It could be generally stated that lowerb values correspond to a higher level of induced seismic activity, while the higher ones correspond to a low and a moderate seismic activity.     

2017年九寨沟MS7.0地震序列活动特征分析

中强地震序列的主震发生后,短时间内受台站距震中较远、尾波干扰和波形重叠等因素的影响,往往会遗漏大量的地震,而地震目录的完整性会直接影响到震后趋势判定和余震序列特征分析的科学性和可靠性。本文利用基于GPU加速的模板匹配方法对2017年8月1～12日的连续波形进行扫描计算,检测九寨沟M_S7.0地震前后遗漏的地震事件,选取台网目录中信噪比较高的1033个地震事件作为模板,在主震前7天至震后5天期间识别出4854个检测地震事件,为台网可定位目录的3.3倍,除去对台网单台地震事件的修正外,还检测到1797个遗漏地震事件,将完备震级从1.6级降低到1.4级。基于补充了遗漏地震的完整地震目录,对2017年8月8日九寨沟M_S7.0地震序列活动特征进行分析。结果表明,前震序列在主震前短时间内出现了地震活动的密集增强,b值也显示为低值状态,可能是深部断层发生破裂之前的加速蠕动的结果。随着时间的推移,余震序列的完备震级逐渐下降并趋于稳定,b值存在缓慢升高的趋势,未来较长时期内余震序列仍将处于持续衰减的状态。     

The Study on Seismicity and Relocation of the July 3,2015 Pishan MS6.5 Earthquake

The July 3, 2015 Pishan M_S6.5 earthquake occurred in the intersection area of the Tarim block and West Kunlun block where the moderate-strong earthquakes have become active in recent years. This paper has studied the seismicity parameters of the earthquake sequences such as the b-value in the Pishan region and its vicinity. In addition, we also relocated the aftershocks of the Pishan M_S6.5 earthquake using the seismic phase report by the double-difference method. The temporal and spatial variation characteristics of the Pishan earthquake sequence in the rupture zone are analyzed. The study is of great significance in the seismic hazard assessment in this region.     

甘东南重力、b值异常与岷县-漳县M6.6地震关系探讨

利用甘东南测网2010—2014年的流动重力观测数据,分析2013年7月22日甘肃岷县—漳县M6.6地震前后1a尺度的重力场变化,并结合地震b值图像和GNSS结果进一步分析甘东南重力、b值异常与岷县—漳县M6.6地震的关系,结果表明:(1)测区重力异常经历了由"区域重力负值变化→持续负值变化→转折上升正变化—发震→震后恢复变化"的时空演化过程;(2)岷县—漳县M6.6地震前后b值曲线出现"低值—高值—回落发震"的变化特征;(3)岷县—漳县M6.6地震发生在重力变化高梯度带上、"0"等值线附近和低b值区域。     

2017年九寨沟7.0级地震前应力状态及b值异常特征研究

2017年8月8日青藏高原东缘四川九寨沟地区发生7.0级强震,依据前人研究结果分析九寨沟7.0级地震发震构造,并计算震前应力状态。结果显示:本次地震受到构造和历史强震的影响,是发生在历史强震引起的应力加载区域。另外,采用中国地震台网1990年以来的地震目录,在评估目录完整性的基础上,利用最大似然法计算得到2017年8月8日九寨沟7.0级地震前震源区及邻区地震b值空间图像。结果显示,九寨沟7.0级地震发生在四川北部地区显著低b值高应力异常区域内部(0.82b0.75)。所以,研究区域内外历史强震可能促进了九寨沟7.0级地震的发生。     

Seismicity anomalies before the great earth—quake of Ms=8.1 in the Kunlun Pass and its significance to earthquake prediction

A great earthquake of M _S=8.1 took place in the west of Kunlun Pass on November 14, 2001. The epicenter is located at 36.2°N and 90.9°E. The analysis shows that some main precursory seismic patterns appear before the great earthquake, e.g., seismic gap, seismic band, increased activity, seismicity quiet and swarm activity. The evolution of the seismic patterns before the earthquake of M _S=8.1 exhibits a course very similar to that found for earthquake cases with M _S≥7. The difference is that anomalous seismicity before the earthquake of M _S=8.1 involves in the larger area coverage and higher seismic magnitude. This provides an evidence for recognizing precursor and forecasting of very large earthquake. Finally, we review the rough prediction of the great earthquake and discuss some problems related to the prediction of great earthquakes.     

The experimental research on b value of AE for the rock specimens with pre-existing crack or notch under uniaxial compression

A lot of researches onb value have been made in seismology. Since the 1960’s Mogi, Scholz and others have studied AE of rock specimens in laboratory and discovered that it is related to natural earthquakes. All former researchers used integral specimens to studyb value in the laboratory. However a major earthquake is usually related to a existing seismic-fault in that area. For this reason, a series of fracture experiments with rock and glass specimens having pre-existing crack or notch is performed in order to examine the effect of preexisting crack tob value. The experimental results show that theb value begins to decrease as soon as the initiation of the crack and finally drop to a very low value when the specimen breaks unstably. Based on these, a brief discussion on the possible mechanism ofb value change for natural earthquakes is given. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,9, 393–400, 1987. Projects sponsored by the Chinese Joint Seismological Science Foundation.     

Monitoring of fault healing after the 1999 Kocaeli, Turkey, earthquake

The North Anatolian fault zone that ruptured during the mainshock of theM 7.4 Kocaeli (Izmit) earthquake of 17 August 1999 has beenmonitored using S wave splitting, in order to test a hypothesisproposed by Tadokoro et al. (1999). This idea is based on the observationof the M 7.2 1995 Hyogo-ken Nanbu (Kobe) earthquake, Japan.After the Hyogo-ken Nanbu earthquake, a temporal change was detectedin the direction of faster shear wave polarization in 2–3 years after the mainshock (Tadokoro, 1999). Four seismic stations were installed within andnear the fault zone at Kizanlik where the fault offset was 1.5 m, about80 km to the east of the epicenter of the Kocaeli earthquake. Theobservation period was from August 30 to October 27, 1999. Preliminaryresult shows that the average directions of faster shear wave polarization attwo stations were roughly parallel to the fault strike. We expect that thedirection of faster shear wave polarization will change to the same directionas the regional tectonic stress reflecting fault healing process. We havealready carried out a repeated aftershock observation at the same site in2000 for monitoring the fault healing process.     

中、强震前川滇块体遥感红外辐射异常分析——以宁蒗M5.7地震和鲁甸M6.5地震为例

许多研究表明中强地震前出现了红外辐射异常。近年来,川滇块体中、强地震频发,利用遥感数据对其红外辐射异常变化开展深入分析十分必要。结合13年MODIS/Terra遥感红外数据,采用空间距平分析法、辐射增强阈值面积指数法及辐射亮温均值低频时间序列法,对2004年4月—2017年6月红外亮温数据进行深入分析,研究川滇块体及周边构造断裂M≥5.5地震的震前异常信息。结果表明:(1)M≥6.0地震在震前3～6个月川滇块体会出现较为明显的红外辐射增强异常现象;(2)震前红外异常并非随机分布,而可能是受到相关块体的控制;(3)以宁蒗M5.7地震和鲁甸M6.5地震为例研究发现,震前辐射增强较高的小范围区域与震中位置存在一定关联。     

2022年泸定M6.8地震前地震活动参数变化研究

2022年9月5日四川泸定发生M6.8地震,为研究泸定地震孕震区的应力变化,选取b值、小震调制比和丛集率这3个参数,对泸定地震前的区域地震活动状态进行计算研究。结果显示:泸定及周边区域几次强震发生前,区域地震活动均存在持续时间较长的低b值时段,且在低b值状态下震前短期内出现小震高丛集、高调制比的现象;鲜水河断裂带的地震活动状态分析显示,此次泸定地震前该断裂带存在持续时间近10个月的低b值状态,且短期内出现丛集率升高、调制比高值现象。通过对比分析,认为泸定地震是鲜水河断裂带构造运动的结果。综合分析认为,结合应力场背景和构造条件研究地震活动b值、固体潮调制比和丛集率的时空变化有助于理解大地震的孕育演化过程。     

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.     

九寨沟7.0级地震松潘至九寨沟震害特征及分析

对九寨沟7.0级地震松潘至九寨沟各烈度区的房屋结构特点和震害特征进行了分析。该地区新建建筑(尤其是框架结构)抗震性能良好,砖混、砖木等结构整体性较差,地震中损坏普遍。通过计算分析震害指数后发现,框架结构震害指数远小于土木、砖木结构;在Ⅵ度区框架结构震害指数普遍为0,随烈度增加,框架结构的震害指数快速增加,而增加幅度逐渐减小。此外,探讨了缺乏建筑物的地区以及建筑物结构为单一框架结构地区的烈度评定问题并给出了相关的建议,以期对今后的地震现场工作起到参考借鉴作用。     

Simulating strong ground motion with the `k -2 ' kinematic source model: An application to the seismic hazard in the Erzincan basin, Turkey

We present a seismic hazard application of a kinematic broad-bandrupture model. This model is based on the k-square dislocation distribution concept (Herrero and Bernard, 1994).Synthetic seismograms are calculated in the far-field approximation with alayered velocity medium for the 13 March 1992 Erzincan earthquake.With a parametrization of the source constrained by other studies,the far-field contribution correctly fits the recorded strong ground motion, which presents a 0.5 to 2 Hz dominant frequency range.As the k-square model is particularly well adapted to synthetize realistic strong-motion at short distances from the fault,it is a reliable tool for calculating seismic hazard maps around active faults. We thus present a synthetic peak ground acceleration map associatedwith the 13 March 1992 activated fault, for a 60 km × 60 km regionaround the epicenter taking into account a smoothed velocity structure ofthe basin in agreement with the absence of significant site effects related to1D resonance deduced from the aftershock records study.This map is compared with several post seismic reports: macroseismicintensities, detailed distribution of damage, and soil cracking andliquefaction. Our model shows that the values of the peak acceleration andvelocity can explain the dominant spatial distributionof these effects, which concentrates in a narrow band along theactivated segment fault, and in particular at its southernextremity. These results enable us to present such maps forhypothetical future earthquake ruptures, located on the major visible orinferred active fault segments in and around the basin. The effects ofthese potential sources are analyzed in relation to the 1992 eventeffects in order to eliminate unknown site responses. We show that thesouthern part of the basin is particularly exposed because of thepresence of strike-slip faults,and that the western part of the basin would suffer a significantlyhigher strong motion levelthan during the 1992 event with the activation of moderate sizednormal faults evidenced on thesouthwestern edge of the basin.