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

利用2009年以来中国地震台网中心地震编目资料,以2016年新疆呼图壁6.2级地震作为研究对象,采用最大似然法进行b值空间扫描计算,获取震前震中及其邻区地震b值的空间图像。研究结果显示,2016年呼图壁6.2级地震发生在震前显著低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值异常分析区域地震危险性,应采用较小范围的地震数据。     

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.     

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

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.     

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.

     

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.     

Study of the Q -value in and around Minle Basin

The digital seimograph network set up by China and France in Zhangye of China had been operated in 1988. The Zhangye network is situated in the middle segment of Hexi corridor and Qilina mountain, which was regarded as a monitoring earthquake area. Using the records of Zhangye digital network theQ-values in and around Minle basin have been measured. The results of this study showed that theQp-values range from 500 to 780, andQs-values range from 230 to 460. TheQ-values of inside of Minle basin are higher than that around the basin. The greater parts of moderate and strong earthquakes occurred along tectonic belts around the Minle basin. Moreover, TheQ-values increased with the depth of penctration of wave ray. The attenuation of S wave is stronger than P wave in shallow layer of crust. Some problems ofQ-value change versus time before and after Sunan eathquake (M _s=5.7) have been also disscussed. These results can be applied to study and to monitor seismic danger of the Minle monitoring area. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 296–302, 1993.     

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.     

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.     

Homogenization of Earthquake Catalog for Northeast India and Adjoining Region

A catalog for northeast India and the adjoining region for the period 1897–2009 with 4,497 earthquakes events is compiled for homogenization to moment magnitude M _w,GCMT in the magnitude range 3–8.7. Relations for conversion of m _b and M _s magnitudes to M _w,GCMT are derived using three different methods, namely, linear standard regression, inverted standard regression (ISR) and orthogonal standard regression (OSR), for different magnitude ranges based on events data for the catalog period 1976–2006. The OSR relations for M _s to M _w,GCMT conversion derived in this paper have significantly lower errors in regression parameters compared to the relations reported in other studies. Since the number of events with magnitude ≥7 for this region is scanty, we, therefore, considered whole India region to obtain the regression relationships between M _w,GCMT and M _s,ISC. A relationship between M _w,GCMT and M _w,NEIC is also obtained based on 17 events for the range 5.2 ≤ magnitude ≤ 6.6. A unified homogeneous catalog prepared using the conversion relations derived in this paper can serve as a reference catalog for seismic hazard assessment studies in northeast India and the adjoining region.     

河北平原地震带b值时空变化特征

利用华北地区近44年地震资料,在区域地震序列完整性分析的基础上,用最小二乘法进行b值时间扫描计算,用最大似然法进行b值空间扫描计算。时间扫描中的b值为研究区内每个扫描窗口的平均b值,因此其变化幅度不大,基本保持在0.62~1.05之间。研究区b值空间分布范围基本维系在0.5~1.4,低b值区域为昌平—宝坻断裂段和唐山—迁安断裂段,变化范围为0.5~0.7,表明该区域地壳介质正处于相对高应力或闭锁状态,存在未来可能发生中强以上地震的潜在危险。     

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.     

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 regional characteristics of seismic activity in China

ＴｈｅｒｅｇｉｏｎａｌｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｓｅｉｓｍｉｃａｃｔｉｖｉｔｙｉｎＣｈｉｎａＺｈｅｎ－ＬｉａｎｇＳＨＩ，ＪｉａｎＷＡＮＧａｎｄＸｉａｏ－ＤｏｎｇＺＨＡＮＧ（时振梁，王健，张晓东）（ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｐｈｙｓｉｃｓ，Ｓ...     

Study of theQ-value in and around Minle Basin

The digital seimograph network set up by China and France in Zhangye of China had been operated in 1988. The Zhangye network is situated in the middle segment of Hexi corridor and Qilina mountain, which was regarded as a monitoring earthquake area. Using the records of Zhangye digital network theQ-values in and around Minle basin have been measured. The results of this study showed that theQp-values range from 500 to 780, andQs-values range from 230 to 460. TheQ-values of inside of Minle basin are higher than that around the basin. The greater parts of moderate and strong earthquakes occurred along tectonic belts around the Minle basin. Moreover, TheQ-values increased with the depth of penctration of wave ray. The attenuation of S wave is stronger than P wave in shallow layer of crust. Some problems ofQ-value change versus time before and after Sunan eathquake (M _s=5.7) have been also disscussed. These results can be applied to study and to monitor seismic danger of the Minle monitoring area.     

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.     

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

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

Monitoring the West Bohemian earthquake swarm in 2008/2009 by a temporary small-aperture seismic array

The most recent intense earthquake swarm in West Bohemia lasted from 6 October 2008 to January 2009. Starting 12 days after the onset, the University of Potsdam monitored the swarm by a temporary small-aperture seismic array at 10 km epicentral distance. The purpose of the installation was a complete monitoring of the swarm including micro-earthquakes (M _L < 0). We identify earthquakes using a conventional short-term average/long-term average trigger combined with sliding-window frequency-wavenumber and polarisation analyses. The resulting earthquake catalogue consists of 14,530 earthquakes between 19 October 2008 and 18 March 2009 with magnitudes in the range of − 1.2 ≤ M _L ≤ 2.7. The small-aperture seismic array substantially lowers the detection threshold to about M _c = − 0.4, when compared to the regional networks operating in West Bohemia (M _c > 0.0). In the course of this work, the main temporal features (frequency–magnitude distribution, propagation of back azimuth and horizontal slowness, occurrence rate of aftershock sequences and interevent-time distribution) of the recent 2008/2009 earthquake swarm are presented and discussed. Temporal changes of the coefficient of variation (based on interevent times) suggest that the swarm earthquake activity of the 2008/2009 swarm terminates by 12 January 2009. During the main phase in our studied swarm period after 19 October, the b value of the Gutenberg–Richter relation decreases from 1.2 to 0.8. This trend is also reflected in the power-law behavior of the seismic moment release. The corresponding total seismic moment release of 1.02×10¹⁷ Nm is equivalent to M _L,max = 5.4.     

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.