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.     

Physical background on estimating b value

The G-R relation lgN=a-bM (1954) is an empirical formula used widely in the seismicity research. But the linearity of b curves has great difference in different time and space domains. An interested question in this paper is that in how large a space-time-strength domain the b value has certain physical connotation. This study told us that we can get optimal statistical results of b value in those space-time domains which can develop correspondent strong shocks with magnitude interval (M _s≥8.5, 8.0≤M _s<8.5, 7.0≤M _s<8.0). Thus, the possible seismogenic areas in which strong shocks with different magnitude intervals develop can be inferred in different regions of the mainland of China. Finally, some new problems are proposed, such as the delimitation of seismic province, the seismicity parameter determination in seismic hazard analysis and in earthquake predictions by using b value. Contribution No. 96A-0074, Institute of Geophysics, SSB, China.     

Recent seismicity in Northeast India and its adjoining region

Recent seismicity in the northeast India and its adjoining region exhibits different earthquake mechanisms – predominantly thrust faulting on the eastern boundary, normal faulting in the upper Himalaya, and strike slip in the remaining areas. A homogenized catalogue in moment magnitude, M _W, covering a period from 1906 to 2006 is derived from International Seismological Center (ISC) catalogue, and Global Centroid Moment Tensor (GCMT) database. Owing to significant and stable earthquake recordings as seen from 1964 onwards, the seismicity in the region is analyzed for the period with spatial distribution of magnitude of completeness m _t, b value, a value, and correlation fractal dimension D _C. The estimated value of m _t is found to vary between 4.0 and 4.8. The a value is seen to vary from 4.47 to 8.59 while b value ranges from 0.61 to 1.36. Thrust zones are seen to exhibit predominantly lower b value distribution while strike-slip and normal faulting regimes are associated with moderate to higher b value distribution. D _C is found to vary from 0.70 to 1.66. Although the correlation between spatial distribution of b value and D _C is seen predominantly negative, positive correlations can also be observed in some parts of this territory. A major observation is the strikingly negative correlation with low b value in the eastern boundary thrust region implying a possible case of extending asperity. Incidentally, application of box counting method on fault segments of the study region indicates comparatively higher fractal dimension, D, suggesting an inclination towards a planar geometrical coverage in the 2D spatial extent. Finally, four broad seismic source zones are demarcated based on the estimated spatial seismicity patterns in collaboration with the underlying active fault networks. The present work appraises the seismicity scenario in fulfillment of a basic groundwork for seismic hazard assessment in this earthquake province of the country.     

古登堡-里希特定律中的b值统计样本量研究

b值是研究地震活动的重要指标,其广泛应用于地震危险性分析和地震预测研究之中,与实际资料的完整性、样本量的大小、计算方法等因素有着重要的关系。常见的b值计算方法有最小二乘法和最大似然法,样本量的大小对这2种方法影响很大。本文利用蒙特卡罗模拟地震目录和汾渭地震带实际目录作为样本,从中抽取不同大小的样本量进行计算,研究不同样本量下这2种方法计算得到的b值与设定值或真实值之间的差别。结果表明,最小二乘法需要的最低样本量为1000,最大似然法为200;当样本量达不到要求时,计算出的b值是不可靠的;由于对样本量的要求不同,前者适用于计算区域的整体b值,而后者在研究某区域b值在时间轴上的变化方面更有优势。本研究为确定2种b值计算方法对样本量的最低要求提供了参考依据。     

Analysis of the Microseismicity Induced by Fluid Injections at the EGS Site of Soultz-sous-Forêts (Alsace, France): Implications for the Characterization of the Geothermal Reservoir Properties

In June/July 2000, a hydraulic stimulation experiment took place at the geothermal EGS site of Soultz-sous-Forêts (Alsace, France) in order to enhance the permeability of the fractured granitic massif at 5 km depth. As it is well known that fluid injections tend to induce microseismic events, a downhole and a surface seismological network have been installed to monitor the seismic activity during the stimulation test. 23400 m³ of fluid have been injected in the rock volume through the open-hole section (4400 m–5000 m) of the well GPK2 at increasing rates of 30 l.s⁻¹, 40 l.s⁻¹ and then 50 l.s⁻¹. More than 7200 microseismic events in the magnitude range –0.9 to 2.6 have been precisely located through a simultaneous inversion of the seismic velocity structure and location parameters. The analysis of the behavior of the seismicity relative to the hydraulic parameters gives important information about the geothermal reservoir. It appears that the evolution of the seismicity strongly depends on the variations of the injection rate: An increase or a decrease leads to changes of the velocity structure, the number and magnitude of microseismic events. This involves different hydro-mechanical processes between the fluid flow and the fracture planes, which will control the final shape of the microseismic cloud. Moreover, the study of the variations of the b-value with time suggests that the stimulation experiment produces a large proportion of small earthquakes, but records of events of magnitude higher than 2 indicate that fluid injection could reactivate structures whose dimensions allow the generation of such earthquakes.     

Seismicity and Seismic Hazard in Alexandria (Egypt) and its Surroundings

— Alexandria City has suffered great damage due to earthquakes from near and distant sources, both in historical and recent times. Sometimes the source of such damages is not well known. Seismogenic zones such as the Red Sea, Gulf of Aqaba-Dead Sea Hellenic Arc, Suez-Cairo-Alexandria, Eastern-Mediterranean-Cairo-Faiyoum and the Egyptian costal area are located in the vicinity of this city. The Egyptian coastal zone has the lowest seismicity, and therefore, its tectonic setting is not well known. The 1998 Egyptian costal zone earthquake is a moderate complex source. It is composed of two subevents separated by 4 sec. The first subevent initiated at a depth of 28 km and caused a rupture of strike (347°), dip (29°) and slip (125°). The second subevent occurred at a shallower depth (24 km) and has a relatively different focal parameter (strike 334°, dip 60° and slip 60°). The available focal mechanisms strongly support the manifestation of a complex stress regime from the Hellenic Arc into the Alexandria offshore area. In the present study a numerical modeling technique is applied to estimate quantitative seismic hazard in Alexandria. In terms of seismic hazard, both local and remote earthquakes have a tremendous affect on this city. A local earthquake with magnitude M_s = 6.7 at the offshore area gives peak ground acceleration up to 300 cm/sec². The total duration of shaking expected from such an earthquake is about three seconds. The Fourier amplitude spectra of the ground acceleration reveals that the maximum energy is carried by the low frequency (1–3 Hz), part of the seismic waves. The largest response spectra at Alexandria city is within this frequency band. The computed ground accelerations due to strong earthquakes in the Hellenic Arc, Red Sea and Gulf of Aqaba are very small (less than 10 cm/sec²) although with long duration (up to 3 minutes).     

Seismicity parameters for important urban agglomerations in India

India’s urban population has increased in the recent times. An earthquake near an urban agglomeration has the potential to cause severe damage. In this article, seismicity parameters for region surrounding important urban agglomerations in India are estimated. A comprehensive earthquake catalogue for the region (6°E–42°E latitude and 60°N–100°N longitude) including historic and pre-historic events has been compiled from various sources. To estimate the parameters, past earthquake data in a control region of radius 300 km has been assembled to quantify the seismicity around each urban agglomeration. The collected earthquake data is first evaluated for its completeness. From combined (historical and instrumental) data, the seismicity parameters b-value, seismic activity rate, λ and maximum expected magnitude (m _max ) have been obtained from the methodology proposed by Kijko and Graham (1998). The obtained activity rates indicate that region surrounding Guwahati urban agglomeration is the most seismically active region followed by Srinagar, Patna, Amritsar and Chandigarh.     

Predictive ground motion scaling in Switzerland: Best estimates and uncertainties

Spectral ground motion (1 to 15 Hz) as a function of distance is modeled for events spanning 3.0 <M_w ≤ 7.0 in Switzerland. The parameters required to simulate ground motion with a stochastic approach are inverted from 2958 horizontal and vertical component waveforms of small to moderate size events (2.0 ≤ M_{L} ≤ 5.2) in the distance range 10 to 300 km recorded on hard rock sites. Using a Monte Carlo simulation, we establish a significantly different amplification of about a factor of 1.9 between the Alpine Foreland and the Alps. To assess the trade-off between the free parameters of our stochastic model and their influence on the predictive ground motion relationship, we perform a grid search over the five-dimensional solution space. The uncertainties are separated into epistemic and aleatory parts; the main epistemic uncertainty is attributed to the lack of data forM > 5. To constrain the viable models at large magnitudes, results from worldwide scaling studies are evaluated in light of the Swiss data. The model that explains best the low observed stress drops at small magnitudes (Δσ ≅ 3 bar) yet matches observed intensities of historical earthquakes assumes a stress drop increasing with moment asM₀^0.25. For three sites in Switzerland we evaluate the sensitivity of the epistemic uncertainty by computing probabilistic hazard curves. Our model offers the most comprehensive and detailed study of spectral ground motion for Switzerland to date.     

川东南长宁地区地震活动及b值演化特征

利用全国统一目录和流动台站目录,研究了四川盆地东南部长宁地区的地震活动特征和b值的空间分布特征.研究结果显示,长宁地区的地震活动在时间上呈现明显的分段特征,地震活动在2015年后明显增强;在空间上,长宁地区的地震活动主要集中在以28.3°N为界限的南、北2个地区,对于这2个区域的b值演化,计算结果显示出不同的分段特征....     

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

2022年9月5日四川泸定发生M6.8地震,为研究泸定地震孕震区的应力变化,选取b值、小震调制比和丛集率这3个参数,对泸定地震前的区域地震活动状态进行计算研究。结果显示:泸定及周边区域几次强震发生前,区域地震活动均存在持续时间较长的低b值时段,且在低b值状态下震前短期内出现小震高丛集、高调制比的现象;鲜水河断裂带的地震活动状态分析显示,此次泸定地震前该断裂带存在持续时间近10个月的低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.     

Probabilistic Seismic Hazard Assessment for Iraq Using Complete Earthquake Catalogue Files

Probabilistic seismic hazard analysis (PSHA) has been carried out for Iraq. The earthquake catalogue used in the present study covers an area between latitude 29°–38.5° N and longitude 39°–50° E containing more than a thousand events for the period 1905–2000. The entire Iraq region has been divided into thirteen seismogenic sources based on their seismic characteristics, geological setting and tectonic framework. The completeness of the seismicity catalogue has been checked using the method proposed by Stepp (1972). The analysis of completeness shows that the earthquake catalogue is not complete below M_s=4.8 for all of Iraq and seismic source zones S1, S4, S5, and S8, while it varies for the other seismic zones. A statistical treatment of completeness of the data file was carried out in each of the magnitude classes. The Frequency Magnitude Distributions (FMD) for the study area including all seismic source zones were established and the minimum magnitude of complete reporting (M_c) were then estimated. For the entire Iraq the M_c was estimated to be about M_s=4.0 while S11 shows the lowest M_c to be about M_s=3.5 and the highest M_c of about M_s=4.2 was observed for S4. The earthquake activity parameters (activity rate , b value, maximum regional magnitude m_max) as well as the mean return period (R) with a certain lower magnitude m_min m along with their probability of occurrence have been determined for all thirteen seismic source zones of Iraq. The maximum regional magnitude m_max was estimated as 7.87 ± 0.86 for entire Iraq. The return period for magnitude 6.0 is largest for source zone S3 which is estimated to be 705 years while the smallest value is estimated as 9.9 years for all of Iraq.The large variation of the b parameter and the hazard level from zone to zone reflects crustal heterogeneity and the high seismotectonic complexity. The seismic hazard near the source boundaries is directly and strongly affected by the change in the delineation of these boundaries. The forces, through which the geological structure along the plate boundary in Eastern and Northeastern Iraq are evolved, are still active causing stress-strain accumulation, deformation and in turn producing higher probabilities of earthquake activity. Thus, relatively large destructive earthquakes are expected in this region. The study is intended to serve as a reference for more advanced approaches and to pave the path for the probabilistic assessment of seismic hazard in this region.     

Lateral variations of the modal (a/b) values for the different regions of the world

Several catalogues of global earthquakes reported for the time period from 1900 to 2000 have been compiled to examine lateral variations of the modal (a/b) values as derived from the Gutenberg–Richter empirical law. For this purpose, the world was divided into 27 different seismic regions in terms of tectonic environments. The parameters a and b were calculated using the least-squares method. The modal values computed for each region were used to produce a global map of the modal values using a grid space of 3°. The results show that a and b-values do not always supply much information about tectonic environments of the different regions. It is observed that the modal values estimated for different tectonic regions are consistent with seismicity of the world and represent global seismic sources better than a or b values. The highest modal values have been found in the oceanic subduction zones, and the lowest values in the oceanic ridges. Lowest b values are observed in trenches. These observations suggest that there is a correlation between apparent stresses and b values. Mapping of the modal values provides detailed images of the zones presenting low and high seismic activity and it may be used as a measure of seismic potential sources and relative hazard levels.     

Design Earthquakes in the U.K.

Three sites in the UK are taken, representative of low, medium and high hazard levels (by UK standards). For each site, the hazard value at 10⁻⁴ annual probability is computed using a generic seismic source model, and a variety of ground motion parameters: peak ground acceleration (PGA), spectral acceleration at 10 Hz and 1 Hz, and intensity. Disaggregation is used to determine the nature of the earthquakes most likely to generate these hazard values. It is found (as might be expected) that the populations are quite different according to which ground motion parameter is used. When PGA is used, the result is a rather flat magnitude distribution with a tendency to low magnitude events (\le 4.5 ML) which are probably not really hazardous. Hazard-consistent scenario earthquakes computed using intensity are found to be in the range 5.8–5.9 ML, which is more in accord with the type of earthquake that one expects to be a worst-case event in the UK. This revised version was published online in July 2006 with corrections to the Cover Date.     

First Order Seismic Microzonation of Haldia, Bengal Basin (India) Using a GIS Platform

The seismic microzonation of the Bengal Basin, Haldia region, India is carried out using the Analytical Hierarchy Process (AHP) on the Geographic Information System (GIS). Three themes are used for the seismic microzonation, namely Peak Ground Acceleration (PGA), predominant frequency and elevation map. An analysis of the maximum magnitude (m _max) and the b value is carried out after preparing the earthquake catalogue from various sources. On the basis of the tectonic set up and seismicity of the region, five seismic zones are delineated which can be a threat to Haldia. They are broadly classified as Zone 1: Arakan-Yoma Zone (AYZ), Zone 2: Himalayan Zone (HZ), Zone 3: Shillong Plateau Zone (SPZ), Zone 4: Bay of Bengal Zone (BBZ) and Zone 5: Shield Zone (SZ). The m _max for Zones 1, 2, 3, 4 and 5 are 8.30 ± 0.51, 9.09 ± 0.58, 9.20 ± 0.51, 6.62 ± 0.43 and 6.61 ± 0.43, respectively. The PGA value is computed for Haldia following the attenuation relationship taking the m _max of each source zone. The expected PGA at Haldia varies from 0.09–0.19 g. The predominant frequency of Haldia is also calculated using the H/V ratio with a frequency ranging from 0.1–3.0 Hz. The elevation map of Haldia is also generated using the Shuttle Radar Topography Mission (STRM) data. A first-order seismic microzonation map of Haldia is prepared in which four zones of hazard have been broadly classified for Haldia as very high seismic hazard zone, high seismic hazard zone, moderate seismic hazard zone and less seismic hazard zone. The very high seismic hazard zone is observed along the southern part of Haldia where there are major industrial and port facilities. The PGA for the four hazard zones are: 0.09–0.13 g for low hazard zone, > 0.13–0.15 g for moderate hazard zone, > 0.15–0.16 g for high hazard zone and > 0.16–0.19 g for very high hazard zone.     

昆仑山口西8.1级地震前的地震活动异常及预测意义的讨论

昆仑山口西8．1地震前明显地出现大区域地震活动平静、高6值、高调制、低锅值的异常配套现象．回顾性分析了巨大地震前不同时问和空间尺度的地震活动性异常,探讨了巨大地震的预测可能性。     

The reservoir-associated earthquakes of April 1983 in western Thailand: Source modeling and implications for induced seismicity

On 22 April 1983, a very large area of Thailand and part of Burma were strongly shaken by a rare earthquake (m _b=5.8,M _s=5.9). The epicenter was located at the Srinagarind reservoir about 190 km northwest of Bangkok, a relatively stable continental region that experienced little previous seismicity. The main shock was preceded by some foreshocks and followed by numerous aftershocks. The largest foreshock ofm _b=5.2 occurred 1 week before the main shock, and the largest aftershock ofm _b=5.3 took place about 3 hours after the main shock. Focal mechanisms of the three largest events in this earthquake sequence have been studied by other seismologists using firts-motion data. However, the solutions for the main shock and the largest aftershock showed significant inconsistency with known surface geology and regional tectonics. We reexamined the mechanisms of these three events by using teleseismicP-andS-waveforms and through careful readings ofP-wave first motions. The directions of theP axes in our study range from NNW-SSE to NNE-SSW, and nodal planes strike in the NW-SE to about E-W in agreement with regional tectonics and surface geology. The main shock mechanism strikes 255°, dips 48°, and slips 63.5°. The fault motions during the main shock and the foreshock are mainly thrust, while the largest aftershock has a large strike-slip component. The seismic moment and the stress drop of the mainshock are determined to be 3.86×10²⁴ dyne-cm and 180 bars, respectively. The occurrence of these thrust events appears to correlate with the unloading of the Srinagarind reservoir. The focal depths of the largest foreshock, the main shock, and the largest aftershock are determined to be 5.4 km, 8 km, and 22.7 km, respectively, from waveform modeling and relative location showing a downward migration of hypocenters of the three largest events during the earthquake sequence. Other characteristics of this reservoir-induced earthquake sequence are also discussed.     

Self-organized Fractal Seismicity of Reservoir Triggered Earthquakes in the Koyna-Warna Seismic Zone,Western India

— Analysis of the Koyna-Warna earthquake catalog (1968–1996) shows that on an average there is a positive correlation between the b value (decrease) and fractal dimensions (decrease in both D₂^s and D₂^t) of earthquake epicenters 0.5 and 2.5 years prior to 1973 (M5.2) and 1980 (M5.5) events, respectively, except a negative correlation for about five years (1988–1993) prior to the 1993/1994 sequence (M5.4). This positive correlation indicates a weaker clustering, or that the epicenters tend to fill the two-dimensional plane. While the origin of the negative correlation seems to be that during periods of large events (low b value), there is strong clustering around the main shock epicenter (high fractal dimension). Interestingly, during the last year (1995–1996) of the studied period both the b value and correlation dimensions rose significantly, suggesting that stress release occurs through increased levels of low magnitude and increasingly scattered seismicity, suggesting an increased risk of larger magnitude events. Incidentally, during 2000 three earthquakes of magnitude M 5.0, one earthquake of M 4.0, 45 earthquakes of magnitude M 3.0–3.9, and several thousand earthquakes of M < 3 have occurred in the region. Thus it can be inferred that at local scales the relationship yields both positive and negative correlation that appears to be controlled by different modes of failure within the active fault complex.Acknowledgement. The authors are grateful to Dr. B.K. Rastogi of NGRI for providing the catalog of Koyna earthquakes and for useful scientific discussions. The comments of Dr. I. G. Main have improved the quality of paper for which we extend to him our sincere thanks. One of the authors (AOM) thanks the Third World Academy of Science and the Council of Scientific and Industrial Research, India for the Postdoctoral Fellowship award under which this work was carried out.     

Seismic quiescence precursory to a past and a future Kurile island earthquake

A systematic search was made for seismicity rate changes in the segment of the Kurile island arc from 45°N to 53°N by studying the cumulative seismicity of shallow (h100 km) earthquakes within 11 overlapping volumes of radius 100 km for the time period 1960 through beginning of 1978. We found that in most parts of this island arc and most of the time the seismicity rate as obtained from the NOAA catalogue and not excluding any events is fairly constant except for increased seismicity in the mid 1960s in the southern portion due to the great 1963 mainshock there, and for seismicity quiescence during part of the time period studied within two well defined sections of the arc. The first of these is a volume of 100 km radius around a 1973 (M _s =7.3) mainshock within which the seismicity rate was demonstrated at the 99% confidence level to have been lower by 50% during 2100 days (5.75 years) before this mainshock. The second volume of seismic quiescence coincides with the 400 km long north Kuriles gap. In this gap the seismicity rate is shown (at the 99% confidence level) to be lower by 50% from 1967 to present (1978), in comparison with the rate within the gap befor 1967, as well as with the rate surrounding the gap. We propose that the anomalously low seismicity rate within the Kuriles gap is a precursor to a great earthquake, the occurrence time of which was estimated by the following preliminary relation between precursory quiescence time and source dimensionT=190L ^0.545. We predict that an earthquake with source length of 200–400 km (M>8) will occur along the north Kurile island arc between latitude 45.5°N and 49.2°N at a time between now and 1994.     

Primary component analysis method and reduction of seismicity parameters

Introduction Data mining (SHAO and YU, 2003) is a new kind of technique developed with database and artificial intelligence in recent years, which processes the data in the database to abstract the im- plied and pre-unknown, but potentially useful information and knowledge from large amounts of incomplete, noisy, blurring and stochastic data. For data mining, data purging is an important link beforehand that includes eliminating noise, making up lost domain, and deleting ineffective data, as…