2018年11月25日新疆博乐MS4.9地震余震检测与目录完整性研究

本文基于匹配滤波技术，通过SEPD（Seismic Events and Phase Detection）对2018年11月25日新疆博乐M_S4.9地震序列进行检测，检测出遗漏地震32条，84.4%地震为M_L0.0—1.0，9.4%地震小于M_L0.0，较地震目录中原有15条地震多213%，检测出的遗漏地震事件使地震目录更加完整。检测后的最小完整性震级由检测前的M_L1.6减至M_L0.8，地震目录最小完整性震级的减小有利于地震工作者对区域地震活动性作出更准确全面的结论，并使地震危险性分析更可靠。     

Am improved source mechanism for the 1935 Timiskaming,Quebec earthquake from regional waveforms

The Timiskaming earthquake, which occurred near the Quebec-Ontario border at the northwest end of the Western Quebec seismic zone in 1935, is one of the five largest instrumentally recorded southeastern Canadian earthquakes. Previous studies of this earthquake concentrated on modeling teismograms recorded at regional distances, a better constrained focal mechanism is obtained. The waveforms indicate thrust faulting on a moderately dipping northwest striking plane at a depth of 10 km. TheM _w of 6.1 determined in this study is in good agreement with previous magnitude estimates (m _b 6.1,M _s 6.0, andm _bLg 6.2–6.3). The focal mechanism is similar to those of many recent small to moderate earthquakes in the region, and the inferred (from theP axis) acting stress of northeast compression is consistent with the overall eastern North American stress field. The Lake Timiskaming Rift Valley in which the earthquake occurred, comprises several northwest striking faults consistent with the strike of the 1935 event. Thus, the 1935 earthquake appears to be a result of faulting on the reactivated Timiskaming graben.     

Updated earthquake catalogue for seismic hazard analysis in Pakistan

A reliable and homogenized earthquake catalogue is essential for seismic hazard assessment in any area. This article describes the compilation and processing of an updated earthquake catalogue for Pakistan. The earthquake catalogue compiled in this study for the region (quadrangle bounded by the geographical limits 40–83° N and 20–40° E) includes 36,563 earthquake events, which are reported as 4.0–8.3 moment magnitude (M_W) and span from 25 AD to 2016. Relationships are developed between the moment magnitude and body, and surface wave magnitude scales to unify the catalogue in terms of magnitude M_W. The catalogue includes earthquakes from Pakistan and neighbouring countries to minimize the effects of geopolitical boundaries in seismic hazard assessment studies. Earthquakes reported by local and international agencies as well as individual catalogues are included. The proposed catalogue is further used to obtain magnitude of completeness after removal of dependent events by using four different algorithms. Finally, seismicity parameters of the seismic sources are reported, and recommendations are made for seismic hazard assessment studies in Pakistan.     

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.     

芦山地震和九寨沟地震台站实测震级差异性分析

本文选取2013年芦山地震和2017年九寨沟地震波形,重新量取垂直向振幅,计算宽频带面波震级M_S（BB）,分析各台站实测震级出现方向性差异的原因。其中,通过572个宽频带台站实测芦山地震震级M_S（BB）7.1,通过603个宽频带台站实测九寨沟地震震级M_S（BB）6.9。芦山地震实测震级大于M_S（BB）7.3的台站呈现WN-ES向分布,与断层倾向一致;实测震级小于M_S（BB）7.0的台站呈现NE-WS向分布,与其所在断层走向一致。九寨沟地震实测震级大于M_S（BB）7.0的台站分布呈现NE向分布,与断层倾向一致;实测震级小于M_S（BB）6.8的台站总体分布较为离散,大体呈现NW-SE向分布,与树正断裂走向一致。实测震级偏大的台站方向性分布与多普勒效应和P波辐射花样联系不明显。对比分析芦山地震和九寨沟地震,去除场地响应和仪器自身影响,台站实测震级差异性仍然存在,因此,台站实测震级差异性是由于受到了多普勒效应、辐射花样、仪器和场地响应之外的因素影响。综合考虑地震震级涉及的影响因素,芦山地震和九寨沟地震的台站实测震级差异性可能与地震波的传播路径有关。     

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

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

Estimation of Seismic Hazard Parameters for the Himalayas and its Vicinity from Complete Data Files

—The maximum likelihood estimation of earthquake hazard parameters has been made in the Himalayas and its surrounding areas on the basis of a procedure which utilizes data containing complete files of the most recent earthquakes. The entire earthquake catalogue used covers the period from 1900–1990. The maximum regional magnitude M _max?, the activity rate of the seismic event λ, the mean return period R of earthquakes with a certain lower magnitude M _max≥ m along with their probability of occurrence, as well as the parameter b of of Gutenberg Richter magnitude-frequency relationship, have been determined for six different seismic zones of the Himalayas and its vicinity. It is shown that in general the hazard is higher in the zone NEI and BAN than the other four zones. The high difference of the b parameter and the hazard level from zone to zone reflect the high seismotectonic complexity and crustal heterogeneity.     

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.     

青海祁连MS5.2地震微震检测与目录完备性研究

利用模板匹配方法对2015年11月23日青海省祁连县M_S5.2地震进行遗漏地震检测研究,由于主震后短时间内目录中遗漏事件较多,故对主震后1天的连续波形进行检测。主震后1天内青海测震台网记录到的余震个数(包括单台)共62个,选取主震后M_L1.0以上余震30个作为模板事件,通过匹配滤波的方式扫描出遗漏地震31个,约为台网目录给出的0.5倍。基于包络差峰值振幅与震级的线性关系估测检测事件的震级参数,最后将检测后的余震目录与台网余震目录在主震后1天内的最小完备震级进行对比分析,结果发现检测后最小完备震级从M_L1.2降到了M_L0.7,得到青海测震台网在祁连地区最小完整性震级为M_L0.7。     

Estimating the size of an earthquake using short-period seismograms of the first three seconds: A simulated experiment using the 1999 Chi-Chi earthquake sequence

Most of the present earthquake early warning systems are based on broadband or strong motion recordings. How-ever, the short-period instruments are still deployed. It is well-known that short-period recordings have saturation problems for large earthquakes when estimating the size of an earthquake. Thus, it is necessary to make clear the magnitude at which saturation starts to occur for the commonly used τc and Pd measurements, respectively. To investigate the possibility of using short-period seismic recordings for earthquake early warning, we conducted a simulated experiment using the strong motion data of the 1999 Chi-Chi earthquake sequence including its main shock and 31 aftershocks, with magnitude span from 4 to 7.6. The strong motion acceleration recordings were convolved with the instrument response of short-period seismographs in northern China to simulate short-period seismograms. Parameters τc and Pd from the first-three-second seismograms were calculated for the simulated short-period recordings and compared with that obtained by the original strong ground motion recordings. The result showed that to some extent, short-period recordings can be used for threshold earthquake early warning, while the magnitude saturation of Pd estimation can be up to 6.5, better than τc estimation.     

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

An evaluation of seismic hazard parameters in southern Turkey

To investigate the characteristics of earthquake hazard parameters as a means of identifying different zones of seismicity, we have compiled a catalogue of about 1850 moderate to large-sized earthquakes with magnitudes m4.0 or greater in southern Turkey for the time period from 1900 to 1990. Several methods have been applied to the earthquake catalogue to assess seismic hazard. The study area is divided into 77 overlapping cells of 2° size. Theoretical calculations were made for the prediction of maximum magnitude, intensity, b-values, strain energy release and corresponding m3 and peak ground acceleration levels for a given period of time. The resultant seismic hazard for each parameter is depicted as a contour map to indicate lateral variations in areas of seismic source. A combination and evaluation of various hazard parameters resulted in more reasonable estimates of hazard. It is found that the most hazardous seismic zones are the Rhodes and Burdur zones where the level of peak ground acceleration reaches up to 280 cm s-2 for an average return period of 100 years.     

基于Gutenberg-Richter定律快速估算最大余震震级:以2017年九寨沟MS7.0地震为例

针对九寨沟M_S7.0地震之后不同时间段的余震序列目录,利用推定最大余震震级,给出了实际最大余震震级的估计值。结果表明,推定最大余震震级随主震后时间尺度的延长而趋于稳定,且该值与实际发生的最大余震的震级一致。需要强调的是,就九寨沟地震序列而言,当余震数据较为完备时,采用主震后较短时间段内（1～2天）的余震目录就可以较准确地估算出主震区域内可能发生的最大余震震级。实际上,主震后12h（0.5天）的余震数据已完全可以给出最大余震震级的有效下限。此外,计算中我们采用了里氏震级M_L和面波震级M_S的余震目录,结果显示,2种震级类型目录的估算结果完全一致,表明利用推定最大余震震级估算实际最大余震震级的方法不受震级类型的影响。据此,该最大余震震级快速评估方法可进一步推广应用于我国大陆地区中强震后强余震灾害分析评估中。目前的拟合技术也显示出随着测震技术的不断进步以及余震识别能力的提高,快速评估方法可以在主震后短时间（<1天）内准确地预测可能发生的最大余震震级。     

Moment Magnitude and Its Calculation

In this paper, we give a brief introduction to the proposal and development history of the earthquake magnitude concept. Moment magnitude M_W is the best physical quantity for measuring earthquakes. Compared with other magnitude scales used traditionally, moment magnitude is not saturated for all earthquakes, regardless of big and small earthquakes, deep and shallow earthquakes, far field and near field seismic data, geodetic and geological data, moment magnitude can be measured, and can be connected with well-known magnitude scales such as surface wave magnitude M_S. Moment magnitude is a uniform magnitude scale, which is suitable for statistics with wide magnitude range. Moment magnitude is the preferred magnitude selected by the International Seismological community, and it is preferred by the departments responsible for publishing seismic information to the public.Moment magnitude is a uniform magnitude scale, which is suitable for statistics with wide magnitude range. Moment magnitude is a preferred magnitude for international seismology, it is preferred by the agency responsible for providing information about earthquakes to the public. We provide all formulas used in the calculation of moment magnitude, and the calculation steps in detail. We also analyzed some problems and rules to solve these problems by using different formulas and numerical value calculation steps.     

Determining the time and space domains for estimating b value in seismic zoning

To take the seismic zone that includes the great shock with M _s≥8.2 as the statistical unit of estimating b value can often lead to more large variance, because the seismogenic zone of the great shock with M _s≥8.2 are larger than that delineated in general seismic zone. Two-level statistical units are considered in this paper. The seismic province is the first level unit that is suitable for group of earthquakes including the great shock of M _s≥8.5. A seismic province can be divided into several seismic zones. They can be taken as the second level unit for group of quakes in which the super magnitude of the greatest shock do not exceed 8. Because of the nonstationarity in time of seismic activity, the unbalancedness of data and differential of seismic temporal series feature in different areas need to be considered when we select the time period for estimating b value. According to local conditions, the time period is selected at one’s discretion in order to reflect seismicity level of this statistical unit in future 100 years. Contribution No. 98A02039, Institute of Geophysics, China Seismological Bureau, China.     

Statistical testing of the relation between burst of aftershocks and strong earthquakes

According to geological tectonics and seismic activites this paper devided North China (30°–45°N, 105°–130°E) into four areas. We analyzed the North China earthquake catalogue from 1970 to 1986 (from 1965 to 1986 for Huabei, the North China, plain region) and identified forty-two bursts of aftershock. Seven of them occurred in aftershock regions of strong earthquakes and seventeen of them in the seismic swarm regions. The relation between strong earthquakes with the remaining eighteen bursts of aftershocks has been studied and tested statistically in this paper. The result of statistical testing show that the random probabilityp of coincidence of bursts of aftershock with subsequent strong earthquakes is less than six percent. By Xu’sR scoring method the efficacy of predicting strong earthquake from bursts of aftershock is estimated greater than 39 percent. Following the method proposed in the paper we analyzed the earthquake catalogue of China from 1987 to June, 1988. The results show that there was only one burst of aftershock occurred on Jan. 6, 1988 withM=3.6 in Xiuyan of Northeast China. It implicates that a potential earthquake withM _S⩽5 might occur in one year afterwards in the region of Northeast China. Actually on Feb. 25, 1988 an earthquake withM _S=5.3 occurred in Zhangwu of Northeast China. Another example is Datong-Yanggao shock on October 18, 1989 which is a burst of aftershock. Three hours after an expected shock withM =6.1 took place in the same area. Two examples above have been tested in practical prediction and this shows that bursts of aftershocks are significant in predicting strong earthquakes. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 273–280, 1991. Part of earthquake catalogue is from Jinbiao Chen, Peiyan Chen and Quanlin Li.     

Probabilistic Assessment of Earthquake Recurrence in Northeast India and Adjoining Regions

Northeast India and adjoining regions (20°–32° N and 87°–100° E) are highly vulnerable to earthquake hazard in the Indian sub-continent, which fall under seismic zones V, IV and III in the seismic zoning map of India with magnitudes M exceeding 8, 7 and 6, respectively. It has experienced two devastating earthquakes, namely, the Shillong Plateau earthquake of June 12, 1897 (M _w 8.1) and the Assam earthquake of August 15, 1950 (M _w 8.5) that caused huge loss of lives and property in the Indian sub-continent. In the present study, the probabilities of the occurrences of earthquakes with magnitude M ≥ 7.0 during a specified interval of time has been estimated on the basis of three probabilistic models, namely, Weibull, Gamma and Lognormal, with the help of the earthquake catalogue spanning the period 1846 to 1995. The method of maximum likelihood has been used to estimate the earthquake hazard parameters. The logarithmic probability of likelihood function (ln L) is estimated and used to compare the suitability of models and it was found that the Gamma model fits best with the actual data. The sample mean interval of occurrence of such earthquakes is estimated as 7.82 years in the northeast India region and the expected mean values for Weibull, Gamma and Lognormal distributions are estimated as 7.837, 7.820 and 8.269 years, respectively. The estimated cumulative probability for an earthquake M ≥ 7.0 reaches 0.8 after about 15–16 (2010–2011) years and 0.9 after about 18–20 (2013–2015) years from the occurrence of the last earthquake (1995) in the region. The estimated conditional probability also reaches 0.8 to 0.9 after about 13–17 (2008–2012) years in the considered region for an earthquake M ≥ 7.0 when the elapsed time is zero years. However, the conditional probability reaches 0.8 to 0.9 after about 9–13 (2018–2022) years for earthquake M ≥ 7.0 when the elapsed time is 14 years (i.e. 2009).     

Regional characteristics of coda attenuation after the lancang-Gengma earthquake

In the light of the single scattering model of coda originating from local earthquakes, and based on the aftershock coda registered respectively at the 4 short period stations installed near the foci shortly after theM7.6 Lancang andM7.2 Gengma earthquakes, this paper has tentatively calculated the rate of amplitude attenuation and theQ _c-value of the coda in the Lancang and Gengma areas using a newly-founded synthetic determination method. Result of the study shows the rate of coda amplitude attenuation demonstrates remarkable regional differences respectively in the southern and northern areas. The southern area presents a faster attenuation (Q _c=114), whereas the northern area shows a slower attenuation (Q _c=231). The paper also discusses the reasons causing such differences. Result of the study also suggests a fairly good linear relation between the coda source factorA _o(f) and the seismic moment and the magnitude. Using the earthquake scaling law, the following formulas can be derived: lgM ₀=lgA ₀(f)+17.6,M _D=0.67lgA ₀(f)+1.21 and logM ₀=1.5M _D+15.79. In addition, the rates of amplitude attenuationβ _s andβ _m are respectively calculated using the single scattering and multiple scattering models, and the ratioβ _s/β_m=1.20−1.50 is found for the results respectively from the two models. Finally, the mean free pathL of the S-wave scattering in the southern and northern areas are determined to be 54 km and 122 km respectively by the relations which can distinguish between the inherentQ _i and scatteringQ _s, testify to this areas having lowQ-values correspond to stronger scatterings. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 71–82, 1992. This study is partly supported by the Seismological Science Foundation of the State Seismological Bureau of China, and the present English version of the paper is translated from its Chinese original by Wenyi Xia, Seismological Bureau of Yunnan Province.     

Harmonization check of M w within the central, northern, and northwestern European earthquake catalogue (CENEC)

Large data sets covering large areas and time spans and composed of many different independent sources raise the question of the obtained degree of harmonization. The present study is an analysis of the harmonization with respect to the moment magnitude M _w within the earthquake catalogue for central, northern, and northwestern Europe (CENEC). The CENEC earthquake catalogue (Grünthal et al., J Seismol, 2009) contains parameters for over 8,000 events in the time period 1000–2004 with magnitude M _w ≥ 3.5. Only about 2% of the data used for CENEC have original M _w magnitudes derived directly from digital data. Some of the local catalogues and data files providing data give M _w, but calculated by the respective agency from other magnitude measures or intensity. About 60% of the local data give strength measures other than M _w, and these have to be transformed by us using available formulae or new regressions based on original M _w data. Although all events are thus unified to M _w magnitude, inhomogeneity in the M _w obtained from over 40 local catalogues and data files and 50 special studies is inevitable. Two different approaches have been followed to investigate the compatibility of the different M _w sets throughout CENEC. The first harmonization check is performed using M _w from moment tensor solutions from SMTS and Pondrelli et al. (Phys Earth Planet Inter 130:71–101, 2002; Phys Earth Planet Inter 164:90–112, 2007). The method to derive the SMTS is described, e.g., by Braunmiller et al. (Tectonophysics 356:5–22, 2002) and Bernardi et al. (Geophys J Int 157:703–716, 2004), and the data are available in greater extent since 1997. One check is made against the M _w given in national catalogues and another against the M _w derived by applying different empirical relations developed for CENEC. The second harmonization check concerns the vast majority of data in CENEC related to earthquakes prior to 1997 or where no moment tensor based M _w exists. In this case, an empirical relation for the M _w dependence on epicentral intensity (I ₀) and focal depth (h) was derived for 41 master events, i.e., earthquakes, located all over central Europe, with high-quality data. To include also the data lacking h, the corresponding depth-independent relation for these 41 events was also derived. These equations are compared with the different sets of data from which CENEC has been composed, and the goodness of fit is demonstrated for each set. The vast majority of the events are very well or reasonably consistent with the respective relation so that the data can be said to be harmonized with respect to M _w, but there are exceptions, which are discussed in detail.     

2017年3月27日漾濞5.1级地震前后气枪震源观测走时数据的变化

以2017年3月27日漾濞5.1级地震为例,根据区域特性和信噪比要求,选取数据较为完整的6个台站记录的2017年1月1日～6月6日期间的宾川地震信号发射台气枪震源波形资料,采用互相关检测技术提取6个台站各自稳定震相的走时数据,并对漾濞5.1级地震前后走时数据的变化情况进行分析。结果表明,漾濞5.1级地震前后6个台站各自稳定震相存在较为明显的走时变化,且短期内走时变化具有较好的同步性,相关台站异常幅度大小和异常出现时间存在细小差异。地震发生前,6个台站走时低值异常过程明显,以YUL台最为显著。地震发生前后走时变化形态特征为双“V”型,漾濞5.1级地震发生在第1个“V”型末端。地震发生后,不同方位相关台站受地震的影响程度不同,走时波动大小存在差异。