2022年青海门源MS6.9地震余震的空间迁移特征

精确识别和定位中强地震序列中的微震事件对于准确判定发震构造具有重要的意义。文章对2022年1月8日发生在青藏高原东北缘的青海门源M_S6.9地震序列进行精定位及微震检测研究。首先运用双差定位法对2022年1月8—16日由中国地震台网中心记录的1 010个门源地震序列原始目录进行重定位,得到404个精定位地震目录。分别采用原始地震目录(CENC)和双差重定位目录(HypoDD),对震源区150 km范围内9个台站的连续波形数据进行微震检测。结果表明,基于CENC目录识别的余震个数是原始目录地震数量的3.0倍,基于HypoDD目录识别的地震个数是原始目录的2.1倍,是HypoDD目录的5.8倍;两种地震目录的微震检测均使得M_L震级完备性从1.7级降低至1.1级。新的地震目录空间位置显示,主震发生后余震主要沿托莱山断裂向西侧扩展,8分钟以后,在托莱山断裂和冷龙岭断裂均发生破裂。根据本研究获取的更高空间分辨率的地震序列,同时结合震源机制解,认为2022年门源M_S6.9地震初始破裂位于近E-W向的托莱山断裂,并触发了NW-SE向的冷...     

2020年新疆伽师MS6.4地震震源区 微震检测与目录完备性分析

中强地震发生后,地震检测因受到尾波的干扰可能会遗漏部分微震事件,影响地震目录的完备性。文章利用波形模板匹配方法对2020年新疆伽师M_S6.4地震序列开展微震检测,相比原始的中国地震台网中心统一地震目录,新检测出1 756个微震事件,地震数量增加了1.3倍。基于检测后的余震目录计算最小完备震级为M_L1.2,地震活动性b值为0.76,较原始目录的M_L1.6和0.77均有所降低。通过伽师震源区地震序列活动特征分析,结果表明前震序列在主震前短时间内(前36小时)出现地震活动的密集增强,相应的b值显示为低值;主震发生后地震序列完备震级较高,随着时间的推移,完备震级缓慢降低并趋于稳定,并且呈周期性的波动。本研究提高了伽师震源区地震目录的完备性,为精细化描述该地区地震序列时空演化特征提供了关键数据基础。     

2014年金寨ML3.9震群序列遗漏地震检测及发震构造分析

针对2014年8月-2015年1月安徽金寨发生的ML3.9震群,利用匹配滤波技术补充台网目录遗漏的地震事件,再利用波形互相关震相检测技术标定P波和S波到时,进一步采用双差定位方法对震群进行重定位,结合震源机制解等分析此次震群活动可能的发震构造.计算结果显示,通过互相关扫描检测到1376个地震台网常规分析遗漏的地震,数量...     

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，地震目录最小完整性震级的减小有利于地震工作者对区域地震活动性作出更准确全面的结论，并使地震危险性分析更可靠。     

Efficient testing of earthquake forecasting models

Computationally efficient alternatives are proposed to the likelihood-based tests employed by the Collaboratory for the Study of Earthquake Predictability for assessing the performance of earthquake likelihood models in the earthquake forecast testing centers. For the conditional L-test, which tests the consistency of the earthquake catalogue with a model, an exact test using convolutions of distributions is available when the number of earthquakes in the test period is small, and the central limit theorem provides an approximate test when the number of earthquakes is large. Similar methods are available for the R-test, which compares the likelihoods of two competing models. However, the R-test, like the N-test and L-test, is fundamentally a test of consistency of data with a model. We propose an alternative test, based on the classical paired t-test, to more directly compare the likelihoods of two models. Although approximate and predicated on a normality assumption, this new T-test is not computer-intensive, is easier to interpret than the R-test, and becomes increasingly dependable as the number of earthquakes increases.     

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.     

Ergodicity examined by the Thirumalai-Mountain metric for Taiwanese seismicity

Ergodicity is a behavior generally limited to equilibrium states and is here defined as the equivalence of ensemble and temporal averages. In recent years, effective ergodicity is identified in simulated earthquakes generated by numerical fault models and in real seismicity of natural fault networks by using the Thirumalai-Mountain metric. Although the effective ergodicity is already reported for Taiwanese seismicity, an immediate doubt is the unrealistic gridded sizes for discretizing the seismic data. In this study, we re-examined the effective ergodicity in Taiwanese seismicity by using reasonable gridded sizes which corresponded with the location errors in the real earthquake catalogue. Initial time and magnitude cut-off were examined for the validity of ergodic behavior. We found that several subsets extracted from Taiwanese seismicity possessed effectively ergodic intervals and all terminations of these ergodic intervals temporally coincided with the occurrences of large earthquakes (M _L < 6.5). We thus confirm the ergodicity in the crustal seismicity by using the Thirumalai-Mountain metric.     

The pattern of diurnal periodicity of weak earthquakes in Iran

The catalogue of earthquakes recorded in Iran during 2006–2010 by a dense network of digital telemetric seismic stations is analyzed. The spectrum of the time series of these earthquakes contains a sharp maximum at 24 h. The corresponding curve of diurnal periodicity constructed by the superposed epoch method has a clear double-peak maximum near noon(11 a.m.-5 p.m. local time), which exceeds the level of seismic activity observed during the rest time of the day by a factor of 3-3.5. In the same time interval, the average hourly magnitudes of seismic events sharply drop from M = 2.15 to M = 1.95. The ratio of the normalized number of earthquakes in the daytime to those at other times of the day, which was determined within a moving window half a square degree in size, has seven distinct compact spatially isolated maxima whose magnitudes attain several dozens to a hundred units. These maxima are probably caused by industrial activity, such as road building and quarry explosions. We also note the presence of the weekend effect when the daytime maximum in the weekly curve of diurnal variations in seismic activity almost completely disappears on Friday, which is the weekend in Muslim countries. At the same time, elimination of the supposed noise component from the catalogue by the approved technique changes nothing for the daytime maximum in the daily pattern of earthquakes in Iran. In order to account for this inconsistency, we suggest invoking additional information on the technogenic seismicity and considering weak earthquakes induced by quarry explosions and vibrations of industrial machines, in particular, power units of numerous hydroelectric power stations distributed over the territory of Iran.     

Test of the EEPAS Forecasting Model on the Japan Earthquake Catalogue

The EEPAS (Every Earthquake a Precursor According to Scale) model is a method of forecasting earthquakes based on the notion that the precursory scale increase () phenomenon occurs at all scales in the seismogenic process. The rate density of future earthquake occurrence is computed directly from past earthquakes in the catalogue. The EEPAS model has previously been fitted to the New Zealand earthquake catalogue and successfully tested on the California catalogue.Here we describe a further test of the EEPAS model in the Japan region spanning 1965–2001, initially on earthquakes with magnitudes exceeding the threshold value 6.75. A baseline model and the Gutenberg-Richter b-value were fitted to the JMA catalogue over the learning period 1926–1964. The performance of EEPAS, with the key parameters unchanged from the New Zealand values, was compared with that of the baseline model over the testing period, using a likelihood ratio criterion. The EEPAS model proved superior. A sensitivity analysis shows that this result is not sensitive to the choice of the learning period or b-value, but that the advantage of EEPAS over the baseline model diminishes as the magnitude threshold is lowered. When key model parameters are optimised for the Japan catalogue, the advantage of EEPAS over the baseline model is consistent for all magnitudes above 6.25, although less than in the New Zealand and California regions.These results add strength to the proposition that the EEPAS model is effective at a variety of scales and in a variety of seismically active regions.     

2016年安徽无为MS3.0震群序列发震构造与震后趋势研究

首先通过模板匹配方法检测无为震群活动期间目录遗漏的地震事件,共识别出5次遗漏地震事件,震级为M_L0.5～1.2,得到了更为完整的地震目录;然后基于波形互相关震相检测技术标定震相到时,进而采用双差定位方法进行精定位,精定位后震群分布更加集中,未见明显的优势方位分布;采用Snoke方法计算震级较大地震的震源机制,结果表明,此次震群为NEE向的水平挤压与NNE向的水平拉张应力场作用下具逆冲分量的走滑型地震活动,严家桥-枫沙湖断裂可能为其发震构造;最后,计算了震群序列的视应力,结果显示,视应力和扣除震级影响后的差视应力随着震群序列的衰减逐渐恢复,因此,分析认为,随着无为震群序列的衰减,震源区发生更大地震的可能性不大。     

中国大陆强震时空分布特征及2021年玛多7.4级地震后趋势分析

对中国大陆强震活动的时空分布特征进行分析,有助于中国大陆强震趋势的判定。由于地震目录完整性的限制,目前对中国大陆强震时间间隔的分析多基于1900年以来地震目录。因为记录时长相对较短,难以排除当前强震时序特征基于偶然的可能性,并且可能导致过拟合而影响预测效果。针对上述问题,本文根据2021—2030年中国大陆地震重点危险区确定工作中的相关资料,对合成地震目录进行了分析,结果表明区域地震活动强弱交替是一种普遍现象,当假设中国大陆活动断层具有准周期复发特征时,合成地震目录时序特征与当前目录最为接近。基于上述认识,使用适用范围更广的单参数指数分布拟合了中国大陆强震间隔,并结合相邻强震构造关联随时间的变化对2021年5月22日玛多7.4级地震后中国大陆的强震趋势进行了分析。结果表明,继玛多7.4级地震之后,未来两年中国大陆再次发生7级以上地震的概率较大,下次强震发生在2022年年底之前的概率为61.81%。下次强震发生在2021年的概率为30.58%,最有可能的发震区域是巴颜喀拉地块;下次强震发生在2022年的概率为44.97%,重点关注南北地震带中南段。     

Analysis of Shallow and Deep Earthquake Doublets in the Fiji-Tonga-Kermadec Region

For the Fiji-Tonga-Kermadec area and for the period from January 1977 to July 2003, the Harvard CMT catalogue lists 1022 shallow, 410 intermediate and 633 deep earthquakes of moment magnitude from 4.9 to 8.0. The magnitude threshold, above which the catalogue is complete, is 5.3–5.4, and the number of earthquakes of magnitude above this value is 691 for shallow, 329 for intermediate and 476 for deep events, respectively. The proportion of earthquakes associated with doublets and multiplets against the total number of earthquakes is approximately the same in both data sets and therefore all earthquake pairs were considered regardless of their magnitude. We investigated all the pairs of earthquakes that occurred at a centroid distance of less than 40, 60 or 90 km from each other and within a time interval of 200, 300 or 450 days, depending on their magnitude. We found 208 pairs of shallow, 31 of intermediate and 92 of deep events. To ascertain whether these earthquakes in pairs are not connected by chance, the possibility of their occurrence in an uncorrelated Poissonian catalogue was considered. It was assumed that in such a catalogue the inter-event time is exponentially distributed, the earthquake magnitude follows the Gutenberg-Richter relation, and the distribution of centroid distances between the events in pairs is controlled by its non-parametric kernel estimate. The probability of the appearance of the observed proportion of doublets of shallow earthquakes in the Poissonian catalogue was found to be very low. The low probability of occurrence in a semi-random catalogue, created by randomising centroid locations in the actual data set, also indicates major importance of the distance criterion used for a doublet specification. In general, shallow earthquakes tend to form pairs at shorter distances and within shorter time intervals than deep earthquakes. Both the distance and the time intervals do not depend on the magnitude of involved events. The largest number of pairs of deep earthquakes is observed at a depth of about 600 km, and the proportion of deep events associated with doublets against the number of all events increases with depth. From comparison of the focal mechanism of earthquakes in pairs, measured by the 3-D rotation angle, it follows that deep earthquakes forming pairs have a more diverse focal mechanism than shallow events; the rotation angle for three quarters of shallow pairs and only for about one third of deep pairs is reasonably small. The azimuth between two events forming a doublet is in about 60–65% of cases close to the strike of one of nodal planes of the first or the second event.     

Tsunami Hazard Posed to New Zealand by the Kermadec and Southern New Hebrides Subduction Margins: An Assessment Based on Plate Boundary Kinematics, Interseismic Coupling, and Historical Seismicity

We assess the tsunami hazard posed to New Zealand by the Kermadec and southern New Hebrides subduction margins. Neither of these subduction zones has produced tsunami large enough to cause significant damage in New Zealand over the past 150?years of well-recorded history. However, as this time frame is short compared to the recurrence interval for major tsunamigenic earthquakes on many of the Earth’s subduction zones, it should not be assumed that what has been observed so far is representative of the long term. For each of these two subduction zones we present plate kinematic and fault-locking results from block modelling of earthquake slip vector data and GPS velocities. The results are used to estimate the current rates of strain accumulation on the plate interfaces where large tsunamigenic earthquakes typically occur. We also review data on the larger historical earthquakes that have occurred on these margins, as well as the Global CMT catalogue of events since 1976. Using this information we have developed a set of scenarios for large earthquakes which have been used as initial conditions for the COMCOT tsunami code to estimate the subsequent tsunami propagation in the southwest Pacific, and from these the potential impact on New Zealand has been evaluated. Our results demonstrate that there is a significant threat posed to the Northland and Coromandel regions of New Zealand should a large earthquake (M _w ?8.5) occur on the southern or middle regions of the Kermadec Trench, and that a similarly large earthquake on the southern New Hebrides Trench has the potential to strongly impact on the far northern parts of New Zealand close to the southern end of the submarine Three Kings Ridge. We propose logic trees for the magnitude–frequency parameters of large earthquakes originating on each trench, which are intended to form the basis for future probabilistic studies.     

山东乳山震群遗漏地震检测及震源区活动特征分析

对震群活动特征的深入研究可以为区域地震危险性分析和地震活动趋势判断提供有效的依据,但受台网布局和震群位置的影响,地震目录中往往会遗漏一些地震,而地震目录的完整性将会影响震群活动特征分析的可靠性.因此,本文利用基于GPU加速的模板匹配方法对山东乳山震群2014年5月至2015年6月期间固定台站和流动台阵记录的连续波形进行...     

中国海域及邻区地震时间分布特征研究

地震时间分布特征研究是进行地震预测和地震危险性分析的重要基础.以中国海域统一地震目录为基础资料,以指数分布模型、伽马分布模型、威布尔分布模型、对数正态分布模型以及布朗过程时间分布(BPT)模型为目标模型,采用极大似然法估算模型参数.根据赤池信息准则(AIC)、贝叶斯信息准则(BIC)以及K-S检验结果确定能够描述海域地...     

四川地区4级以上调制小震与5级以上中强地震的关系探讨

通过对四川地区1972年以来的地震目录进行挑选和分析,得到四川地区4级以上调制小震与5级以上中强地震的对应数据结果。在所发生的41组调制小震中只有14组对应发生了中强地震,这说明在四川地区调制小震与中强地震的对应关系不是太好。但是,调制小震与中强地震的发震地点间的距离对应相关却比较好,平均距离间隔为93．49千米;所对应的调制小震与中强地震的发震时间的时间间隔平均为17．5个月。     

自动任意下载网页上的地震目录资料

介绍一种能直接从网页获取地震目录的应用软件，实现对任意地震信息网站上的任意地震目录自动下载。该软件向HTTP服务器发送请求获取服务器传回包含数据的网页源文件，最后提取源文件的目录信息，并进行数据格式转换。核心功能主要是实现原来人工操作搜索地震目录，并对接收来的文件进行分析，分离其中的数据部分，以达到自动下载和整理地震目录的目的。     

The macroseismic field of the Balkan area

A catalogue of 356 macroseimic maps which are available for the Balkan area was compiled, including information on the source parameters of the corresponding earthquakes, the macroseismic parameters of their strength and their macroseismic field. The data analysis of this catalogue yields new empirical relations for attenuation, which can be applied for the calibration of historical events, modelling of isoseismals and seismic hazard assessment. An appropriate analysis allowed the separation and estimation of the average values of the geometrical spreading, n, and anelastic attenuation factor, c, for the examined area which were found equal to –3.227 ± 0.112 and –0.0033 ± 0.0010. Scaling relations for the focal macroseismic intensity, I_f, and the epicentral intensity I₀, versus the earthquake moment magnitude were also determined for each Balkan country. A gradual decrease of the order of 0.5 to 1 intensity unit is demonstrated for recent (after 1970) earthquakes in Greece. Finally the depths of the examined earthquakes as they robustly determined (error <5 km) on the basis of macroseismic data were found to have small values ( 10 km). However large magnitude earthquakes show higher focal depths ( 25 km), in accordance with an increase of the seismic fault dimensions for such events.     

The El Salvador earthquakes of January and February 2001: context, characteristics and implications for seismic risk

The small Central American republic of El Salvador has experienced, on average, one destructive earthquake per decade during the last hundred years. The latest events occurred on 13 January and 13 February 2001, with magnitudes M_w 7.7 and 6.6, respectively. The two events, which were of different tectonic origin, follow the patterns of the seismicity of the region although neither event has a known precedent in the earthquake catalogue in terms of size and location. The earthquakes caused damage to thousands of traditionally built houses and triggered hundreds of landslides, which were the main causes of fatalities. The earthquakes have clearly demonstrated trends of increasing seismic risk in El Salvador due to rapid population expansion in areas of high shaking and landslide hazard, exacerbated by deforestation and uncontrolled urbanisation. The institutional mechanisms required for the control of land use and building practice are very weak and present a major obstacle to risk mitigation.     

福建省历史强震目录的复核与评述

１９９５年出版的《中国历史强震目录》对１９８８年出版的《中国地震简目》所收录的福建历史强震目录作了全面的修订。增、删、改地震１６次,占原目录同一时段地震数的６１．５％。本文通过核查,分析地震史料,利用中国东部烈度衰减规律,参考现今闽台地震的实际烈度资料,并考虑台湾强震波及影响,对历史强震参数变动较大的１１次地震进行了复核与评述,认为其大部分地震参数的变动是需要商榷的。为此,对所增、删、改的地震分别