Earth's rotation-triggered earthquakes before the 2018 MS5.7 Xingwen earthquake

For earthquakes (M_L≥2.0) that occurred from January 2006 to October 2018 around the M_S5.7 Xingwen earthquake occurred on 16 December 2018 in Xingwen, Sichuan province, China, we statistically investigated the correlation between the phase of Earth's rotation and the occurrence of earthquakes via Schuster's test to determine the signals that triggered earthquakes before the M_S5.7 Xingwen event. The results were evaluated based on the P-value where a smaller P-value corresponded to a higher correlation between the occurrence of an earthquake and Earth's rotation. We investigated the spatial distribution of P-values in the region around the epicenter of the M_S5.7 Xingwen event, and obtained a result exhibiting a extremely low-P-value region. The M_S5.7 event occurred inside near the northern boundary of this region. Furthermore, we analyzed the temporal evolution of P-values for earthquakes that occurred within the extremely low-P-value region and found that some extremely low P-values (less that 0.1%), i.e., significant correlation, were calculated for earthquakes that occurred before the M_S5.7 Xingwen earthquake. Among sixty-one earthquakes with the lowest P-value, occurred from May 2014 to April 2018, a vast majority of them occurred during the acceleration of Earth's rotation. The lower P-value obtained in this study reveals that the Xingwen source body probably was extremely unstable prior to the occurrence of the M_S5.7 Xingwen earthquake.     

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 depth variations in the b-value of frequency-magnitude distribution of the earthquakes in the Garm region of Tajikistan

The depth changes in the b-value and density of the number of earthquakes in different magnitude bins (M ≥ 1.8, M ≥ 3.0, M ≥ 3.5) are analyzed using highly accurate seismological observations carried out in 1955–1991 at the Garm prognostic area in Tadjikistan. It is found that the observed b-values are controlled by the variations in the proportion between weak and strong earthquakes. Two horizons with different patterns of the b-value are identified in the Earth’s crust above and below a depth of 15–16 km. The b-value in the upper and lower horizons is close to 0.8 and 1.2, respectively. The lower horizon is marked by almost complete absence of relatively strong earthquakes with M ≥ 3.0. The observed changes in the b-value with increasing depth could probably be due to the increase in the strength of crustal material caused by the growth in temperature and confining pressure in the depth interval from 0 to 15 km. The transitional interval between the upper and lower crustal horizons (～13–18 km), which is characterized by a sharp drop in seismic activity, can probably be associated with the zone of the phase transition of crustal material from an elastic brittle state to a plastic state, as suggested by some authors. Typically, the top of this zone hosts the hypocenters of the strongest earthquakes in a given territory. The correlation is established between the crustal areas with low b-values and the locations of the strongest earthquakes in the region. It is suggested that the three-dimensional mapping of the b-value can be helpful for estimating the location, depth, and maximal magnitude of the probable strong earthquakes in seismically active regions and can be used to assess seismic risks.     

An estimation of the maximum b-value in the Gutenberg-Richter relation

An estimation of the maximum b-value in the Gutenberg-Richter relation is achieved from empirical laws of earthquake occurrence and from dimensional considerations of power law distributions of earthquake size. The deduced maximum b-value, 1.5 (or with error limits 1.30⩽b_max⩽1.64) agrees very well with maximum b-values computed from earthquake catalogues or published elsewhere in the literature. © 1999 Elsevier Science Ltd. All rights reserved.     

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.     

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

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

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 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.     

Probabilistic Assessment of Earthquake Hazard in Gujarat and Adjoining Region of India

The Gujarat and adjoining region falls under all four seismic zones V, IV, III and II of the seismic zoning map of India, and is one of the most seismically prone intracontinental regions of the world. It has experienced two large earthquakes of magnitude M _w 7.8 and 7.7 in 1819 and 2001, respectively and several moderate earthquakes during the past two centuries. In the present study, the probability of occurrence of earthquakes of M ≥ 5.0 has been estimated during a specified time interval for different elapsed times on the basis of observed time intervals between earthquakes using three stochastic models namely, Weibull, Gamma and Lognormal. A complete earthquake catalogue has been used covering the time interval of 1819 to 2006. The whole region has been divided into three major seismic regions (Saurashtra, Mainland Gujarat and Kachchh) on the basis of seismotectonics and geomorphology of the region. The earthquake hazard parameters have been estimated using the method of maximum likelihood. The logarithmic of likelihood function (ln L) is estimated and used to test the suitability of models in three different regions. It was found that the Weibull model fits well with the actual data in Saurashtra and Kachchh regions, whereas Lognormal model fits well in Mainland Gujarat. The mean intervals of occurrence of earthquakes are estimated as 40.455, 20.249 and 13.338 years in the Saurashtra, Mainland Gujarat and Kachchh region, respectively. The estimated cumulative probability (probability that the next earthquake will occur at a time later than some specific time from the last earthquake) for the earthquakes of M ≥ 5.0 reaches 0.9 after about 64 years from the last earthquake (1993) in Saurashtra, about 49 years from the last earthquake (1969) in Mainland Gujarat and about 29 years from the last earthquake (2006) in the Kachchh region. The conditional probability (probability that the next earthquake will occur during some specific time interval after a certain elapsed time from last earthquake) is also estimated and it reaches about 0.8 to 0.9 during the time interval of about 57 to 66 years from the last earthquake (1993) in Saurashtra region, 31 to 51 years from the last earthquake (1969) in Mainland Gujarat and about 21 to 28 years from the last earthquake (2006) in Kachchh region.     

Seismic Hazard Evaluation in Western Turkey as Revealed by Stress Transfer and Time-dependent Probability Calculations

Western Turkey has a long history of destructive earthquakes that are responsible for the death of thousands of people and which caused devastating damage to the existing infrastructures, and cultural and historical monuments. The recent earthquakes of Izmit (Kocaeli) on 17 August, 1999 (M _w  = 7.4) and Düzce (M _w  = 7.2) on 12 November, 1999, which occurred in the neighboring fault segments along the North Anatolian Fault (NAF), were catastrophic ones for the Marmara region and surroundings in NW Turkey. Stress transfer between the two adjacent fault segments successfully explained the temporal proximity of these events. Similar evidence is also provided from recent studies dealing with successive strong events occurrence along the NAF and parts of the Aegean Sea; in that changes in the stress field due to the coseismic displacement of the stronger events influence the occurrence of the next events of comparable size by advancing their occurrence time and delimiting their occurrence place. In the present study the evolution of the stress field since the beginning of the twentieth century in the territory of the eastern Aegean Sea and western Turkey is examined, in an attempt to test whether the history of cumulative changes in stress can explain the spatial and temporal occurrence patterns of large earthquakes in this area. Coulomb stress changes are calculated assuming that earthquakes can be modeled as static dislocations in elastic half space, taking into account both the coseismic slip in large (M ≥ 6.5) earthquakes and the slow tectonic stress buildup along the major fault segments. The stress change calculations were performed for strike-slip and normal faults. In each stage of the evolutionary model the stress field is calculated according to the strike, dip, and rake angles of the next large event, whose triggering is inspected, and the possible sites for future strong earthquakes can be assessed. A new insight on the evaluation of future seismic hazards is given by translating the calculated stress changes into earthquake probability using an earthquake nucleation constitutive relation, which includes permanent and transient effects of the sudden stress changes.     

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

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

Depth and region dependence of b-value for micro-aftershocks of the May 12th, 2008 Wenchuan earthquake and its tectonic implications

Micro-aftershocks with magnitude range of 1.5?4 around the Wenchuan earthquake epicenter, the southern part of the Longmenshan fault zone, exhibit good frequency-magnitude linear relationships, thus enabling b-value analysis. The average b-value for micro-aftershocks of M1.5?4 from July to December of 2008 in our local study region is about 0.88, similar to the b-value for all aftershocks of M3.0?5.5 from May, 2008 to May, 2009 along the whole Longmenshan fault zone. The similarity between the local and regional b-values possibly indicates that the southern part of the Longmenshan fault zone has similar seismogenic environment to the whole Longmenshan fault zone. Alternatively, it may also imply that b-values derived from all events without consideration of structural variation can not discriminate local-scale tectonic information. The present study shows that the b-value for the Wenchuan earthquake micro-aftershocks varies with different regions. The b-value in southwest of the Yingxiu town is higher than that in the northeast of the Yingxiu town. The high b-value in the southwest part where the Wenchuan earthquake main shock hypocenter located indicates that the current stress around the hypocenter region is much lower than its surrounding area. The b-values are also dependent on depth. At shallow depths of < 5 km, the b-values are very small (~0.4), possibly being related to strong wave attenuation or strong heterogeneity in shallow layers with high content of porosity and fractures. At depths of ~5?11 km, where most aftershocks concentrated, the b-values become as high as ~0.9?1.0. At the depth below ~11 km, the b-values decrease with the depth increasing, being consistent with increasing tectonic homogeneity and increasing stress with depth.

     

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.     

Gorkha earthquake (MW7.8) and aftershock sequence: A fractal approach

On April 25, 2015, Nepal was struck by the M_W7.8 Gorkha earthquake followed by an intense aftershock sequence. It was one of the most destructive earthquakes in the Himalayan arc, causing more than 8900 fatalities. In this study, we analyzed the dataset (429 events, magnitude of completeness (M_c) ≥ 4.2 local magnitude) of the first 45 days after the Gorkha earthquake to estimate the seismicity parameters b-value, D-value, and p-value. We used the maximum likelihood method to estimate the b-value and Omori-Utsu parameters, whereas the correlation integral method was applied to estimate the fractal dimension (D-value). The analysis was carried out using running and sliding window techniques. The lowest b-value (0.57 ± 0.04) and the highest D-value (1.65 ± 0.02) were computed at the time of the Gorkha earthquake, after which the b-value significantly increased to a maximum of 1.57. It again dropped to 0.93 at the time of the major aftershock on May 12, 2015. The D-value showed an initial quick drop and then decreased in a wavy pattern until the end of the study period, indicating the clustering and scattering of earthquakes in a fault region. The b-value contour map identified the eastern part of the study area as a high stress region (b = ～0.8), implying that the stress shifted to that region. The D-value contour map reveals that the seismogenic structure shifted from linear to planar in the region. The rate of aftershock decay (p = 0.86 ± 0.04) for a short period reflects that the level of stress decreased rapidly. This study helps to understand the level of stress and seismicity pattern of a region, which could be useful for aftershock studies.     

Large Earthquake Hazard of the San Jacinto Fault Zone,CA, from Long Record of Simulated Seismicity Assimilating the Available Instrumental and Paleoseismic Data

We investigate spatio-temporal properties of earthquake patterns in the San Jacinto fault zone (SJFZ), California, between Cajon Pass and the Superstition Hill Fault, using a long record of simulated seismicity constrained by available seismological and geological data. The model provides an effective realization of a large segmented strike-slip fault zone in a 3D elastic half-space, with heterogeneous distribution of static friction chosen to represent several clear step-overs at the surface. The simulated synthetic catalog reproduces well the basic statistical features of the instrumental seismicity recorded at the SJFZ area since 1981. The model also produces events larger than those included in the short instrumental record, consistent with paleo-earthquakes documented at sites along the SJFZ for the last 1,400 years. The general agreement between the synthetic and observed data allows us to address with the long-simulated seismicity questions related to large earthquakes and expected seismic hazard. The interaction between m ≥ 7 events on different sections of the SJFZ is found to be close to random. The hazard associated with m ≥ 7 events on the SJFZ increases significantly if the long record of simulated seismicity is taken into account. The model simulations indicate that the recent increased number of observed intermediate SJFZ earthquakes is a robust statistical feature heralding the occurrence of m ≥ 7 earthquakes. The hypocenters of the m ≥ 5 events in the simulation results move progressively towards the hypocenter of the upcoming m ≥ 7 earthquake.     

New Zealand Earthquake Forecast Testing Centre

The New Zealand Earthquake Forecast Testing Centre is being established as one of several similar regional testing centres under the umbrella of the Collaboratory for the Study of Earthquake Predictability (CSEP). The Centre aims to encourage the development of testable models of time-varying earthquake occurrence in the New Zealand region, and to conduct verifiable prospective tests of their performance over a period of five or more years. The test region, data-collection region and requirements for testing are described herein. Models must specify in advance the expected number of earthquakes with epicentral depths h ≤ 40 km in bins of time, magnitude and location within the test region. Short-term models will be tested using 24-h time bins at magnitude M ≥ 4. Intermediate-term models and long-term models will be tested at M ≥ 5 using 3-month, 6-month and 5-year bins, respectively. The tests applied will be the same as at other CSEP testing centres: the so-called N test of the total number of earthquakes expected over the test period; the L test of the likelihood of the earthquake catalogue under the model; and the R test of the ratio of the likelihoods under alternative models. Four long-term, three intermediate-term and two short-term models have been installed to date in the testing centre, with tests of these models commencing on the New Zealand earthquake catalogue from the beginning of 2008. Submission of models is open to researchers worldwide. New models can be submitted at any time. The New Zealand testing centre makes extensive use of software produced by the CSEP testing centre in California. It is envisaged that, in time, the scope of the testing centre will be expanded to include new testing methods and differently-specified models, nonetheless that the New Zealand testing centre will develop in parallel with other regional testing centres through the CSEP international collaborative process.     

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

Temporal and spatial distributions of precursory seismicity rate changes in the Thailand-Laos-Myanmar border region: implication for upcoming hazardous earthquakes

To study the prospective areas of upcoming strong-to-major earthquakes, i.e., M _w  ≥ 6.0, a catalog of seismicity in the vicinity of the Thailand-Laos-Myanmar border region was generated and then investigated statistically. Based on the successful investigations of previous works, the seismicity rate change (Z value) technique was applied in this study. According to the completeness earthquake dataset, eight available case studies of strong-to-major earthquakes were investigated retrospectively. After iterative tests of the characteristic parameters concerning the number of earthquakes (N) and time window (T _w ), the values of 50 and 1.2 years, respectively, were found to reveal an anomalous high Z-value peak (seismic quiescence) prior to the occurrence of six out of the eight major earthquake events studied. In addition, the location of the Z-value anomalies conformed fairly well to the epicenters of those earthquakes. Based on the investigation of correlation coefficient and the stochastic test of the Z values, the parameters used here (N = 50 events and T _w  = 1.2 years) were suitable to determine the precursory Z value and not random phenomena. The Z values of this study and the frequency-magnitude distribution b values of a previous work both highlighted the same prospective areas that might generate an upcoming major earthquake: (i) some areas in the northern part of Laos and (ii) the eastern part of Myanmar.     

Temporal variation of b value associated with M ��4 earthquakes in the reservoir-triggered seismic environment of the Koyna�CWarna region, Western India

It is generally found that the b values associated with reservoir-triggered seismicity (RTS) are higher than the regional b values in the frequency magnitude relation of earthquakes. In the present study, temporal and spatial variation of b value is investigated using a catalog of 3,000 earthquakes from August 2005 through December 2010 for the Koyna?CWarna region in Western India, which is a classical site of RTS globally. It is an isolated (30?×?20?km²) zone of seismicity where earthquakes of up to M ??5 are found to occur during phases of loading and unloading of the Koyna and Warna reservoirs situated 25?km apart. For the Warna region, it is found that low b values of 0.6?C0.9 are associated with earthquakes of M ??4 during the loading phase. The percentage correlation of the occurrence of an M????4 earthquake with a low b value outside the 1?? or 2?? level is as high as 78?%. A drastic drop in the b value of about 50?% being reported for an RTS site may be an important precursory parameter for short-term earthquake forecast in the future.     

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.