Temporal variations of seismicity parameters in the Central Alborz,Iran

Temporal changes of b-value, fractal (correlation) dimensions of epicenters (D _e ) and occurrence time of earthquakes (D _t ) and relations between these parameters were calculated to investigate precursory changes before 28 May 2004, Baladeh-Kojour earthquake (M _w = 6.3) of Central Alborz, Iran. 2086 events with M _N ≥ 1.7 were selected for our analyses. A wide range of variation was seen in these parameters: b-value ～ 0.6–1.11, D _e ～ 0.97–1.64, and D _t ～ 0.13–0.93. The results showed decreases in all fractal parameters several months before the main shock. This decrease, which might have arisen due to clusters of events occurred between 2002–2003, was followed by a systematic increase, corresponding to the increased level of low-magnitude seismicity. It seems that changes in fractal parameters may be precursors of Baladeh-Kojour earthquake which was caused by seismic activation and quiescence. Furthermore, a positive correlation between b-value and D _e was detected before the main shock (D _e = 0.87 + 0.7b) and during aftershock sequences (D _e = 2b ± 0.09), which was further on changed to a negative one (D _e = 2.56–1.32b).     

Joint analysis of b-value and apparent stress before the 2011 MW9.0 Tohoku-Oki,Japan earthquake

Detecting tempo-spatial changes of crust stress associated with major earthquakes has implications for understanding earthquake seismogenic processes. We conducted a joint analysis of b-value and apparent stress in the source region before the March 11, 2011 M_W9.0 Tohoku-Oki, Japan earthquake. Earthquakes that occurred between January 1, 2000 and March 8, 2011 were used to estimate b-values, while source parameters of events with magnitudes of Ms5.0–6.9 between January 1, 1997 and March 8, 2011 were used to calculate the apparent stresses. Our results show that the average b-value decreased steadily from 1.26 in 2003 to 0.99 before the Tohoku-Oki mainshock. This b-value decrease coincided with an increase in the apparent stress from 0.65 MPa to 1.64 MPa. Our results reveal a clear negative correlation between the decrease in b-value and increase in apparent stress, which lasted for approximately eight years prior to the 2011 mainshock. Additionally, spatial pattern results of the relative change in b-value show that the area associated with drastic b-value decreases (25% or greater) was concentrated near the 2011 mainshock epicenter. The joint analysis of b-value and apparent stress provides a promising method for detecting anomalies that could serve as potential indicators of large earthquakes.     

The b-value and fractal dimension of local seismicity around Koyna Dam (India)

Earthquakes began to occur in Koyna region (India) soon after the filling of Koyna Dam in 1962. In the present study, three datasets 1964–1993, 1993–1995, and 1996–1997 are analyzed to study the b-value and fractal dimension. The b-value is calculated using the Gutenberg–Richter relationship and fractal dimension D _corr. using correlation integral method. The estimated b-value and D _corr. of this region before 1993 are found to be in good agreement with previously reported studies. In the subsequent years after 1995, the b-value shows an increase. The estimated b-values of this region are found within the limits of global average. Also, the pattern of spatial clustering of earthquakes show increase in clustering and migration along the three zones called North-East Zone, South-East Zone (SEZ), and Warna Seismic Zone. The earthquake events having depth ≤5 km are largely confined to SEZ. After 1993, the D _corr. shows decrease, implying that earthquake activity gets clustered. This seismic clustering could be helpful for earthquake forecasting.     

Properties of the Aftershock Sequences of the 2003 Bingöl, M D = 6.4, (Turkey) Earthquake

Aftershock sequences of the magnitude M _{W =}6.4 Bingöl earthquake of 1 May, 2003 (Turkey) are studied to analyze the spatial and temporal variability of seismicity parameters of the b value of the frequency-magnitude distribution and the p value describing the temporal decay rate of aftershocks. The catalog taken from the KOERI contains 516 events and one month’s time interval. The b value is found as 1.49 ± 0.07 with Mc =3.2. Considering the error limits, b value is very close to the maximum b value stated in the literature. This larger value may be caused by the paucity of the larger aftershocks with magnitude M _D ≥ 5.0. Also, the aftershock area is divided into four parts in order to detect the differences in b value and the changes illustrate the heterogeneity of the aftershock region. The p value is calculated as 0.86 ± 0.11, relatively small. This small p value may be a result of the slow decay rate of the aftershock activity and the small number of aftershocks. For the fitting of a suitable model and estimation of correct values of decay parameters, the sequence is also modeled as a background seismicty rate model. Constant background activity does not appear to be important during the first month of the Bingöl aftershock sequences and this result is coherent with an average estimation of pre-existing seismicity. The results show that usage of simple modified Omori law is reasonable for the analysis. The spatial variability in b value is between 1.2 and 1.8 and p value varies from 0.6 to 1.2. Although the physical interpretation of the spatial variability of these seismicity parameters is not straightforward, the variation of b and p values can be related to the stress and slip distribution after the mainshock, respectively. The lower b values are observed in the high stress regions and to a certain extent, the largest b values are related to Holocene alluvium. The larger p values are found in some part of the aftershock area although no slip occurred after the main shock and it is interpreted that this situation may be caused by the alluvium structure of the region. These results indicate that the spatial distribution in b and p values are generally related to the rupture mechanism and material properties of an aftershock area.     

b-value and fractal dimension variations in Iran

Iran sits on a region with a high intrinsic level of seismic activity due to its tectonic setting. Through statistical examination of the earthquakes listed in the catalogue from International Institute of Earthquake Engineering and Seismology (IIEES), this research attempted to calculate some seismicity factors and find correlation between them. A preliminary analysis indicated changes in the b-value of the Gutenberg-Richter relationship over the study region. Thus, the study area was divided into five zones (Alborz, Zagros, Azerbaijan, Central and East) and b-value was computed for each zone. Considering faulting mechanism styles and the b-values in the region, it was found that the lowest b-values belong to the thrust events and strike-slip faulting earthquakes have intermediate values. These findings support previous studies. Furthermore, results of b-value calculation were used for the estimation of accumulated differential stresses (σ₁–σ₂) over each zone. Overall, the b-value for Iran is averagely low which signifies the high stress tectonic regime in this region. Also, by having calculated fractal dimension (D) in each zone, a correlation obtained showing that in Iran region, the b-value correlates to fractal dimension by D = 4.2b–2 relation which does not support Aki's (1981) speculation of D = 3b/c.     

Preliminary report of the September 5, 2022 MS 6.8 Luding earthquake,Sichuan, China

The 2022 M_S 6.8 Luding earthquake is the strongest earthquake in Sichuan Province, Western China, since the 2017 M_S 7.0 Jiuzhaigou earthquake. It occurred on the Moxi fault in the southeastern segment of the Xianshuihe fault, a tectonically active and mountainous region with severe secondary earthquake disasters. To better understand the seismogenic mechanism and provide scientific support for future hazard mitigation, we summarize the preliminary results of the Luding earthquake, including seismotectonic background, seismicity and mainshock source characteristics and aftershock properties, and direct and secondary damage associated with the mainshock. The peak ground displacements in the NS and EW directions observed by the nearest GNSS station SCCM are ～35 mm and ～55 mm, respectively, resulting in the maximum coseismic dislocation of 20 mm along the NWW direction, which is consistent with the sinistral slip on the Xianshuihe fault. Back-projection of teleseismic P waves suggest that the mainshock rupture propagated toward south-southeast. The seismic intensity of the mainshock estimated from the back-projection results indicates a Mercalli scale of VIII or above near the ruptured area, consistent with the results from instrumental measurements and field surveys. Numerous aftershocks were reported, with the largest being M_S 4.5. Aftershock locations (up to September 18, 2022) exhibit 3 clusters spanning an area of 100 km long and 30 km wide. The magnitude and rate of aftershocks decreased as expected, and the depths became shallower with time. The mainshock and two aftershocks show left-lateral strike-slip focal mechanisms. For the aftershock sequence, the b-value from the Gutenberg-Richter frequency-magnitude relationship, h-value, and p-value for Omori’s law for aftershock decay are 0.81, 1.4, and 1.21, respectively, indicating that this is a typical mainshock-aftershock sequence. The low b-value implies high background stress in the hypocenter region. Analysis from remote sensing satellite images and UAV data shows that the distribution of earthquake-triggered landslides was consistent with the aftershock area. Numerous small-size landslides with limited volumes were revealed, which damaged or buried the roads and severely hindered the rescue process.     

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值存在缓慢升高的趋势,未来较长时期内余震序列仍将处于持续衰减的状态。     

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

利用全国统一目录和流动台站目录,研究了四川盆地东南部长宁地区的地震活动特征和b值的空间分布特征。研究结果显示,长宁地区的地震活动在时间上呈现明显的分段特征,地震活动在2015年后明显增强;在空间上,长宁地区的地震活动主要集中在以28.3°N为界限的南、北2个地区,对于这2个区域的b值演化,计算结果显示出不同的分段特征。此外,研究了不同震中距下,长宁地区2018年兴文5.7级、2019年珙县5.3级和2019年长宁6.0级3次主要地震事件震前与震后的b值演化,结果表明统计范围会显著影响b值的变化特征,当计算范围为20km时,3次主震震前短期内均出现了低于背景b值90%的低b值异常,这表明在地震目录相对完备的情况下,利用低b值异常分析区域地震危险性,应采用较小范围的地震数据。     

Analysis of aftershock sequences in South and Southeastern Spain

A probabilistic modeling is used to analyze the spatio-temporal behavior of eleven aftershock sequences occurred in South and Southeastern Spain. This study focuses on the analysis of two seismicity parameters: the b-value of the frequency-magnitude distribution, and the p-value, explaining the temporal decay rate of aftershocks. The estimated b values range between 0.77 ± 0.05 and 1.18 ± 0.10 close to the typical b-values of the aftershock frequency-magnitude relationship b ≈ 1.0. The estimated p-values range between 0.75 ± 0.03 and 1.43 ± 0.10 showing broad regimes of the temporal decay of aftershocks. The modified Bath’s law used to analyze the energy partitioning, suggests that a large fraction of the accumulated energy is released in the mainshock and relatively small fraction of energy is released during aftershock sequence, for example 80% of the total energy is released during the Mula 1999 mainshock, 88% during Bullas 2002 mainshock and 87% during La Paca 2005 mainshock. The fractal dimension D₂ is estimated using the correlation integral, and then used to derive the slip ratio, as the ratio of the slip occurred on primary fault segment to the total slip. For example, we obtained a slip ratio equal to 71% for the Mula 1999 aftershock sequence, 61% for the Bullas 2002 event, 58% for the La Paca 2005 aftershock, 50% for the Lorca 2011 sequence and 63% for the sequence triggered by the Gador 2002 mainshock.Finally, the correlations between the fractal dimension, the b-value and the p-value is analyzed, and the Aki’s relation D = 3b/c is discussed as well.     

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…     

Dynamics of Induced Seismicity during the Filling of the Nurek Reservoir

The seismological data in the area of induced seismicity in the region of the Nurek reservoir are analyzed. The analysis is based on the developed database for the earthquakes that occurred from 1955 to 1989 and is aimed at finding the regularities in the variations of the parameters of the transitional seismic regime caused by filling a reservoir. These parameters include the b-value—the slope of the graph of the Gutenberg–Richter frequency–magnitude relationship, the fractal dimension d of the set of the epicenters, and fracture cycle parameter q = αb ? d, where coefficient α determines the ratio between the magnitude and source size M = α log l + β. It is shown that during the filling of a reservoir, these parameters undergo statistically reliable variations: at the initial stages, the b-value increases, the fractal dimension of the set of epicenters decreases, and the fracture cycle parameter q grows and becomes positive in the middle of the time interval of reservoir filling. After a reservoir is filled, these parameters recover their background values. The aftershock sequences of the three strongest earthquakes—before, in the beginning, and in the middle of the reservoir filling period—are studied. It is confirmed that the Omori parameter p for the aftershock sequences during filling is smaller than for the earthquake before filling. Based on the dynamics of the studied parameters, it is conjectured that the relaxation time of the transitional seismic regime after the emergence of induced seismicity is about 10 years.     

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级地震的发生。     

Analysis of the 2012–2013 Torreperogil-Sabiote seismic swarm

This study analyses the temporal clustering, spatial clustering, and statistics of the 2012–2013 Torreperogil-Sabiote (southern Spain) seismic swarm. During the swarm, more than 2200 events were located, mostly at depths of 2–5 km, with magnitude event up to m_bLg 3.9 (M_w 3.7). On the basis of daily activity rate, three main temporal phases are identified and analysed. The analysis combines different seismological relationships to improve our understanding of the physical processes related to the swarm's occurrence. Each temporal phase is characterized by its cumulative seismic moment. Using several different approaches, we estimate a catalog completeness magnitude of m_c≅ 1.5. The maximum likelihood b-value estimates for each swarm phase are 1.11 ± 0.09, 1.04 ± 0.04, and 0.90 ± 0.04, respectively. To test the hypothesis that a b-value decrease is a precursor to a large event, we study temporal variations in b-value using overlapping moving windows. A relationship can be inferred between change in b-value and the regime style of the rupture. b-values are indicators of the stress regime, and influence the size of ruptures. The fractal dimension D₂ is used to perform spatial analysis. Cumulative gamma and beta functions are used to analyse the behaviour of inter-event distances during the earthquake sequence.     

On the seismic regime of the Japan region before the Tohoku mega-earthquake (M w = 9)

The seismic regime taking place before the Tohoku mega-earthquake was studied using the catalog of the Japan Meteorological Agency (JMA). We show that the Tohoku earthquake was preceded by a 6–7-year period of regional reduction in the b-value and in the rate of main shocks. The space-time regions that involved precursory activation were nearly identical with the predictive phenomena that were previously detected by A.A. Lyubushin from an analysis of seismic noise based on data from the Japanese F-net. We discovered a previously unknown effect of correlation between the number of main shocks and the b-value. Both the ordinary foreshock activation and the longer weaker tendency, which consist in a precursory increase in the seismicity rate, were identified in the vicinities of M ≈ 7 Japanese earthquakes (similarly to the seismicity in the Generalized Vicinities of large earthquakes based on worldwide data).     

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.     

Continuous-cyclic variations in the b-value of the earthquake frequency-magnitude distribution

Seismicity of the Earth (M ≥ 4.5) was compiled from NEIC, IRIS and ISC catalogues and used to compute b-value based on various time windows. It is found that continuous cyclic b-variations occur on both long and short time scales, the latter being of much higher value and sometimes in excess of 0.7 of the absolute b-value. These variations occur not only yearly or monthly, but also daily. Before the occurrence of large earthquakes, b-values start increasing with variable gradients that are affected by foreshocks. In some cases, the gradient is reduced to zero or to a negative value a few days before the earthquake occurrence. In general, calculated b-values attain maxima 1 day before large earthquakes and minima soon after their occurrence. Both linear regression and maximum likelihood methods give correlatable, but variable results. It is found that an expanding time window technique from a fixed starting point is more effective in the study of b-variations. The calculated b-variations for the whole Earth, its hemispheres, quadrants and the epicentral regions of some large earthquakes are of both local and regional character, which may indicate that in such cases, the geodynamic processes acting within a certain region have a much regional effect within the Earth. The b-variations have long been known to vary with a number of local and regional factors including tectonic stresses. The results reported here indicate that geotectonic stress remains the most significant factor that controls b-variations. It is found that for earthquakes with M _w ≥ 7, an increase of about 0.20 in the b-value implies a stress increase that will result in an earthquake with a magnitude one unit higher.     

Probabilistic aftershock hazard analysis,two case studies in West and Northwest Iran

Aftershock hazard maps contain the essential information for search and rescue process, and re-occupation after a main-shock. Accordingly, the main purposes of this article are to study the aftershock decay parameters and to estimate the expected high-frequency ground motions (i.e., Peak Ground Acceleration (PGA)) for recent large earthquakes in the Iranian plateau. For this aim, the Ahar-Varzaghan doublet earthquake (August 11, 2012; M _N =6.5, M _N =6.3), and the Ilam (Murmuri) earthquake (August 18, 2014 ; M _N =6.2) have been selected. The earthquake catalogue has been collected based on the Gardner and Knopoff (Bull Seismol Soc Am 64(5), 1363-1367, 1974) temporal and spatial windowing technique. The magnitude of completeness and the seismicity parameters (a,??b) and the modified Omori law parameters (P,??K,??C) have been determined for these two earthquakes in the 14, 30, and 60 days after the mainshocks. Also, the temporal changes of parameters (a,??b,??P,??K,??C) have been studied. The aftershock hazard maps for the probability of exceedance (33%) have been computed in the time periods of 14, 30, and 60 days after the Ahar-Varzaghan and Ilam (Murmuri) earthquakes. For calculating the expected PGA of aftershocks, the regional and global ground motion prediction equations have been utilized. Amplification factor based on the site classes has also been implied in the calculation of PGA. These aftershock hazard maps show an agreement between the PGAs of large aftershocks and the forecasted PGAs. Also, the significant role of b parameter in the Ilam (Murmuri) probabilistic aftershock hazard maps has been investigated.     

余震序列分维特征的初步探讨

本文应用关联维数计算方法,分别计算了雷波5.4级地震、康定5.2级、5.8级地震和大邑6 1/4级地震的余震序列的关联维数D_2,并给出了相应的离差和相关系数。它们的关联维数分别为2.26982、1.93915、2.28738和1.5367。在这4个地震序列中,发生强余震的序列的D_2大于无强余震的序列。     

Spatial Fractal Evolutions and Hierarchies for Microearthquakes in Central Greece

—Statistical characteristics of seismicity represented by microearthquakes are examined for three regions in central Greece, in particular the fractal correlation dimension, D ₂?, and traditional b values are examined in tandem as a function of time by using the moving window technique. The Patras region contains the complicated tectonics, extending to damaging historical earthquakes of the western Corinth Gulf and the Rio Graben, yielding D ₂ values between 0.40 and 1.20 with b between 0.94 and 1.27 unusually, the temporal evolution between D ₂ and b generates a positive correlation, although the variation is mostly in D ₂?. When the whole evolution is divided into two stages then the positive correlations are even stronger than for the whole evolution. The views of Henderson and others might suggest a highly fractured, fluid-filled zone. The Pavliani region, with no known active fault, and the Volos region, containing the through-going Nea Ankhialos fault, yield D ₂ values 0.33 to 0.79 with b 0.92 to 1.30 and D ₂ 0.82 to 1.56 withb 1.02 to 1.37, respectively. Temporal evolution between D ₂ and b provides a typical negative correlation in both regions. Examination of gross seismicity (time window embracing the whole data set) in each region produces D ₂ values for the Patras and Volos regions that are both larger than that for Pavliani; there are no obvious differences amongst the b values. This accords with the knowledge that Patras and Volos are in regions with very active seismotectonic features which generate repeated strong earthquakes exceeding 6M _s ?. D ₂ for both regions is fairly close to 1, the topological dimension of a line, consistent with seismicities on leading active fault zones or through-go ing faults. These values highlight the ability for microearthquakes to illuminate the character of their parent tectonic province. Resolution and hierarchy in these data from Greece are compared with these aspects elsewhere (Japan, Turkey, South America, USA) in the sense that their banding with respect to examined magnitude, areal extent and duration of observation period for respective data sets is examined in relation to the earthquake potential of the parent seismotectonic province. Evolutions are then categorised as being macroscopic, transitional or microscopic in character.     

Universality of the Seismic Moment-frequency Relation

—We analyze the seismic moment-frequency relation in various depth ranges and for different seismic regions, using Flinn-Engdahl's regionalization of global seismicity. Three earthquake lists of centroid-moment tensor data have been used the Harvard catalog, the USGS catalog, and the Huang et al. (1997) catalog of deep earthquakes. The results confirm the universality of the β-values and the maximum moment for shallow earthquakes in continental regions, as well as at and near continental boundaries. Moreover, we show that although fluctuations in earthquake size distribution increase with depth, the β-values for earthquakes in the depth range of 0–500 km exhibit no statistically significant regional variations. The regional variations are significant only for deep events near the 660 km boundary. For declustered shallow earthquake catalogs and deeper events, we show that the worldwide β-values have the same value of 0.60 ± 0.02. This finding suggests that the β-value is a universal constant. We investigate the statistical correlations between the numbers of seismic events in different depth ranges and the correlation of the tectonic deformation rate and seismic activity (the number of earthquakes above a certain threshold level per year). The high level of these correlations suggests that seismic activity indicates tectonic deformation rate in subduction zones. Combined with the universality of the β-value, this finding implies little if any variation in maximum earthquake seismic moment among various subduction zones. If we assume that earthquakes of maximum size are similar in different depth ranges and the seismic efficiency coefficient, χ, is close to 100% for shallow seismicity, then we can estimate χ for deeper earthquakes for intermediate earthquakes χ≈ 5%, and χ≈ 1% for deep events. These results may lead to new theoretical understanding of the earthquake process and better estimates of seismic hazard.