小滇西地区M≥5地震前b值时空变化特征

同时使用完整地震目录和去余震目录,采用最大似然法对小滇西地区进行b值时空扫描,研究该地区M≥5.0地震前b值的时空变化特征。结果显示:2种地震目录b值时间曲线的整体变化趋势无明显的差异,去余震目录的空间扫描结果低值异常更加突出;小滇西地区M≥5.0地震前,b值随时间变化趋势呈现出多样化的特征,但是75%的地震发生之前,该地区处于低b值状态;该地区M≥5.0地震均发生在低b值区域或者高低b值过渡地区。另外,在不同空间尺度的b值时间曲线中,均发现了较长时间的低b值现象,尽管起始时间不同,但是结束时间相同,分析认为该现象可能是2014年盈江M6.1地震的前兆异常。因此,小滇西地区低b值是本区M≥5.0地震的一个异常指标,可以根据空间扫描的低b值地区确定未来可能的震源位置,但是b值无法单独作为一个临震指标,需要与其他地震活动性参数联合使用;b值时间曲线上的长期低值,可能是较大地震的前兆异常。     

历史有感地震目录编制原则方法的研究——以《中国历史有感地震目录》编制为例

在吸收历史地震研究、历史政区地理研究等多方面成果基础上,基于历史有感地震记载特点,探索了历史有感地震目录的编制原则,给出了地震日期、震中位置、震级的确定方法。以《中国地震历史资料汇编》等为依据,编制了公元前618年-公元1949年期间中国境内3≤M<443/4级历史有感地震目录9 121条。目录数量为中国同期文字记载的43/4级以上破坏性地震的7倍多,比全球其他国家或地区最早的同类目录早1 695年,约占全球同期M≤4.5级有感地震目录的2/3。     

地震活动性的统计分析:由过去推测将来的可能性研究

为了由过去的地震活动性推测将来的地震活动性,引入了地震（震级≥m）的期望年发生率v（≥m）来描述一个地区的地震活动性．根据全球地震目录（1964-1994年）以及南加州（1932-1995年）和华北（1970-1994年）两个区域地震目录资料,以统计样本量作为目录记录时间长短的相对量度,对由不同的统计样本量计算得出的地震实际年发生率v（≥m,T,t）进行了统计分析,得到三点结沦：①在统计样本量n足够大的情况下,地震实际年发生率表现出准平稳时间过程的特征,可近似地看作地震期望年发生率,本文给出了这种近似的误差（离差系数）与统计样本量之间的定量关系;②离差系数与统计样本量之间的关系与震级无关,表现出不同震级层次的相似性;③统计样本量相同时,不同震级的地震期望年发生率之间满足logv（≥）=a－bm的关系,形式上与G－R关系相似,但它给出了由小地震的统计特征估计大地震的期望年发生率及其统计误差的方法．基于上述结论,进一步讨论了地震活动性的统计特征在地震危险性分析中的潜在应用。     

Induced seismicity: Realia and myths

The problems of initiation of earthquakes by other strong earthquakes and by underground nuclear explosions are considered. For this purpose, special software was designed that made it possible to process all of the world catalog data on earthquakes of 1950–2005. A more effective and reliable method was applied to the analysis of seismic activity, taking into account its background state. As a result, it is unambiguously established that the effect of strong earthquakes or explosions is negligibly small. In any case, their effect cannot be recognized from catalog data that served as a basis for many studies of induced seismicity problems. This conclusion is of great significance for assessment of the possibility of real-time prediction of earthquakes. If this effect were to exist, such prediction would become impossible due to unpredictability of the triggering action.     

Preliminary analysis of earthquake probability based on the synthetic seismic catalog

The analysis of seismic hazards relies on the statistical analysis of historical seismic data and the instrumental seismic catalog to obtain the regional earthquake recurrence interval and earthquake probability. The accuracy of analysis thus depends strongly on the completeness of the seismic data used. However, available seismic catalogs are too short or incomplete for the reliable analysis of the statistical characteristics of earthquakes. If a long-term synthetic seismic catalog can be generated using a physics-based numerical simulation, and the simulation results match the crustal deformation, seismicity, and other observations,then such a synthetic catalog helps us to further understand the characteristics of seismic activity and analyze the regional seismic hazard. In this paper, taking the northeastern Tibetan Plateau as a case study, we establish a three-dimensional visco-elastoplastic finite-element model to simulate earthquake cycles and the spatiotemporal evolution of earthquakes on the model fault system and obtain a seismic catalog on a time scale of tens of thousands of years. On the basis that the model satisfies the regional geodynamics of the northeastern Tibetan Plateau, we analyze seismicity on the northeastern Tibetan Plateau using the simulated synthetic earthquake catalog. The characteristics of earthquake recurrence at different locations and different magnitudes, and the long-term average probability of earthquake occurrence within the fault system on the northeastern Tibetan plateau are studied. The results are a reference for regional seismic hazard assessment and provide a basis for the physics-based numerical prediction of earthquakes.     

基于深度学习到时拾取自动构建长宁地震前震目录

将深度学习到时拾取、震相关联技术与传统定位方法联系起来,构建一套连续波形自动化处理与地震目录自动构建流程,对于高效充分利用地震资料,提升微震检测能力具有十分重要的意义.我们应用最新发展的迁移学习震相识别技术、震相自动关联技术,对长宁M_S6.0地震震中附近21个台站震前半个月(6月1日—6月17日)的连续记录波形进行P、S震相识别、震相自动关联和初步定位,并应用传统绝对定位和相对定位技术得到了长宁地震震前微震活动的绝对和相对定位目录.其中绝对定位目录能在较小的误差范围匹配85%的人工处理目录,其发震时刻平均误差为0.36±0.07 s,震级平均误差为0.15±0.024级,水平定位平均误差为1.45±0.028 km,其识别的1.0级以下微震数目是人工的8倍以上,将长宁地震震前微震目录的检测下限提升至M_L-1左右,证明了基于深度学习到时识取和REAL (Rapid Earthquake Association and Location,快速震相关联和定位技术)震相自动关联来构建微震目录具有较好的实用性.我们的自动地震目录揭示了长宁M_S...     

Assessment of Creep Events as Potential Earthquake Precursors: Application to the Creeping Section of the San Andreas Fault, California

—We report the analysis of over 16 years of fault creep and seismicity data from part of the creeping section of the San Andreas fault to examine and assess the temporal association between creep events and subsequent earthquakes. The goal is to make a long-term evaluation of creep events as a potential earthquake precursor. We constructed a catalog of creep events from available digital creepmeter data and compared it to a declustered seismicity catalog for the area between San Juan Bautista and San Benito, California, for 1980 to 1996. For magnitude thresholds of 3.8 and above and time windows of 5 to 10 days, we find relatively high success rates (40% to 55% 'hits') but also very high false alarm rates (generally above 90%). These success rates are statistically significant (0.0007 < P < 0.04). We also tested the actual creep event catalog against two different types of synthetic seismicity catalogs, and found that creep events are followed closely in time by earthquakes from the real catalog far more frequently than the average for the synthetic catalogs, generally by more than two standard deviations. We find no identifiable spatial pattern between the creep events and earthquakes that are hit or missed. We conclude that there is a significant temporal correlation between creep events and subsequent small to moderate earthquakes, however that additional information (such as from other potential precursory phenomena) is required to reduce the false alarm rate to an acceptable level.     

2013年3月涿鹿微震群遗漏地震事件检测和发震构造分析

震群活动时,短时间发生大量地震,不同地震事件的记录波形相互交叠影响,易造成地震目录的遗漏,对震群发震构造分析等研究带来不利的影响.本文针对2013年3月3日至5日在河北涿鹿发生的微震震群,利用匹配滤波技术,以地震台网观测目录所记录地震事件的波形为模板,在连续波形记录中搜索与模板相似的信号,从而检测台网目录遗漏的地震.利用波形互相关标定新检测到地震事件的P波和S波到时,进而对其震中位置和震级做出估计.计算结果显示,通过互相关扫描检测到52个地震台网常规分析遗漏的地震,约为地震目录给出的45个事件的1.16倍.检测到的遗漏地震震级估算为ML0.1~0.9,通过震级-频次统计分析,加入遗漏地震后地震目录的完整性在ML0.3~0.8范围内有较明显的改善.根据地震事件精定位结果,推测此次震群的发震构造为北西走向倾角较大的断层,施庄断裂为发震构造的可能性较大.     

2012年9月7日彝良地震及余震序列双差定位研究

本文提出了时域多通道相关检测函数并用其计算波形互相关走时差数据,采用双差定位法对2012年9月7日云南彝良地震和余震序列共944个地震进行重定位,得到652个重定位事件,并与目录数据的结果进行了对比.本文采用了多个准则对走时差数据进行筛选,确保定位结果稳定可靠.得到MS5.7主震的震中为27.516°N,103.951°E,震源深度6.9 km;MS5.6主震的震中为27.543°N,104.023°E,震源深度7.27 km;重定位结果显示,地震序列紧缩为条带状并沿附近断裂走向分布,深度总体分布较重定位前变浅,集中分布在5~8 km,地震群出现轻微倾斜.东西向、南北向、深度和发震时刻的平均相对误差分别为55.2 m,43.0 m,186.7 m和0.01 s,走时残差16ms.研究表明:互相关数据的结果要优于目录数据;震源深度与速度模型存在较大的相关性;确定彝良—会泽断裂为本次彝良地震序列的发震构造.     

2013年甘肃岷县漳县M_S6.6地震余震序列时域衰减特征分析

本文基于补充遗漏地震事件后的地震目录和修正的大森公式,对2013年7月22日甘肃岷县漳县6.6级地震余震序列时间域衰减特征进行了分析.鉴于主震后短时间内目录遗漏的余震较多,首先利用岷县台连续地震记录波形的高频包络差,检测主震后3h内目录遗漏的地震.经检测在主震后3h内共发现目录遗漏的ML1.0以上地震139个,最大震级为ML3.6.主震后1000s内检测到遗漏地震69个,约为目录给出余震数量的6倍.而后使用补充遗漏地震的目录,基于修正的大森公式分别拟合余震频度和余震地震矩随时间的变化.结果显示拟合p值约为1.07,表明岷县漳县地震余震序列衰减速率与全球平均水平接近,而未补充遗漏地震的频度拟合会造成余震序列衰减速率的低估.利用高频包络差直接计算地震频度曲线,通过三种衰减模式对地震频度曲线拟合参数比较,未观察到岷县漳县地震主震后存在早期余震缺失现象.分析认为,加入遗漏地震可以提高余震频度拟合估计衰减速率结果的准确度和精度,若缺少遗漏地震检测结果,则使用地震矩拟合所得衰减速率结果准确度较优,但需充分考虑其精度上的误差.在分析余震序列衰减特征的实际研究工作中,需根据地震目录完整性选择适当的拟合方法.     

Earthquake risk prediction as a stochastic process

The extrapolation in time of an earthquake sequence considered as a multidimensional stochastic point process is discussed. Estimates of seismic risk for both long- and short-term predictions are considered and an algorithm for the calculations is proposed. Several examples of short-term extrapolations are carried out by means of Monte Carlo simulations of the process. An assessment of the predictability of the seismic process shows that the catalog of strong earthquakes (M ? 7.0) contains about 0.4 bits of information per earthquake for the particular model of the process applied here.     

A catalog of aftershock sequences in Greece (1971–1997): Their spatial and temporal characteristics

A complete catalog of aftershock sequences is provided for main earthquakes with M_L 5.0, which occurred in the area of Greece and surrounding regions the last twenty-seven years. The Monthly Bulletins of the Institute of Geodynamics (National Observatory of Athens) have been used as data source. In order to get a homogeneous catalog, several selection criteria have been applied and hence a catalog of 44 aftershock sequences is compiled. The relations between the duration of the sequence, the number of aftershocks, the magnitude of the largest aftershock and its delay time from the main shock as well as the subsurface rupture length versus the magnitude of the main shock are calculated. The results show that linearity exists between the subsurface rupture length and the magnitude of the main shock independent of the slip type, as well as between the magnitude of the main shock (M) and its largest aftershock (Ma). The mean difference M–Ma is almost one unit. In the 40% of the analyzed sequences, the largest aftershock occurred within one day after the main shock.The fact that the aftershock sequences show the same behavior for earthquakes that occur in the same region supports the theory that the spatial and temporal characteristics are strongly related to the stress distribution of the fault area.     

Spatio-temporal Scanning and Statistical Test of the Accelerating Moment Release (AMR) Model Using Australian Earthquake Data

— The Accelerating Moment Release (AMR) preceding earthquakes with magnitude above 5 in Australia that occurred during the last 20 years was analyzed to test the Critical Point Hypothesis. Twelve earthquakes in the catalog were chosen based on a criterion for the number of nearby events. Results show that seven sequences with numerous events recorded leading up to the main earthquake exhibited accelerating moment release. Two occurred near in time and space to other earthquakes preceded by AMR. The remaining three sequences had very few events in the catalog so the lack of AMR detected in the analysis may be related to catalog incompleteness. Spatio-temporal scanning of AMR parameters shows that 80% of the areas in which AMR occurred experienced large events. In areas of similar background seismicity with no large events, 10 out of 12 cases exhibit no AMR, and two others are false alarms where AMR was observed but no large event followed. The relationship between AMR and Load-Unload Response Ratio (LURR) was studied. Both methods predict similar critical region sizes, however, the critical point time using AMR is slightly earlier than the time of the critical point LURR anomaly.     

山东及邻区小震活动月频次分布浅析

采用1975-2010年4月山东区域地震台网地震目录,确定研究区地震震级下限,剔除余震,对预处理后的小震活动月频次进行Kolmogorov-Smirnov分布检验,及小震月频次(X)与累计次数之间的函数拟合.结果表明,山东地区小震活动月频次符合正态分布,月频次与其累计次数之间存在类似G-R关系.对小震月频次分布的峰度K...     

华北地区地震危险性分析和地面运动预测的一致性方法

利用华北地区地震活动性资料,建立了地震危险性计算的一致性模型.在此模型的基础上,得出了北京、天津、唐山和济南等7个城市未来2500年内地震的时空强度分布,并计算了2500年回复周期的地震动峰值加速度(PGA).结果表明,唐山和太原的PGA最大(>0.2g),石家庄和北京次之(≈0.17g).对华北地区2500年地震记录的正演计算结果表明,太原和唐山地区的潜在地震危险最有可能来源于震级在6.0～7.0、震中距离在12～15km的地震活动;而北京、天津和石家庄地区则可能来源于震级在5.5～6.0、震中距离在10km左右的地震活动.采用IBC(International Building Code)方法计算后的结果显示,太原、唐山等地区的PGA与2001年我国地震动峰值加速度值基本一致,与此地区的较高地震活动性特征相符.利用随机震源模型,还给出了影响此7个城市的最大地震记录的加速度、速度及位移时程曲线,这对本区工程建筑的抗震性设计以及对救援设施的选址等有重要作用.     

On applicability of the RTL prognostic algorithms and energy estimation to Sakhalin seismicity

Premonitory phases (seismic quiescence and foreshock activity) have been retrospectively identified before the Neftegorsk and Uglegorsk earthquakes using the RTL technique. The probabilities that these phases were accidental are less than 1 and 2%, respectively. This allows an optimistic assessment of the possibility of applying this technique to seismicity at Sakhalin. The estimates of the time and energy class for the two earthquakes, using a model of self-organized seismic criticality, proved to be unconvincing because obvious acceleration of the seismic process prior to these seismic events did not occur. The applicability of this approach to the seismicity at Sakhalin should be tested for future large earthquakes. The regional Sakhalin catalog for 1980–2000, with a lowest completely reported energy class of K = 8 (lent by the Geophysical Service, Russian Academy of Sciences) was used as the database for this study.     

Study on the Moment Magnitude of Small and Moderate Earthquakes Located in the Inner Mongolia Region

Based on digital seismic waveform data from Inner Mongolia Digital Seismic Network, the source spectrum parameters of 182 small and moderate earthquakes from January, 2009 to September, 2016 are derived, and the seismic moment M₀ and moment magnitude M_W of the earthquakes are calculated. The M_L-M_W relationship and the relationship between stress drop and magnitude are obtained using the linear regression method. It is clear that incorporating the moment magnitude into the seismic quick report catalog and the official earthquake catalog can enrich earthquake observation report content, thus providing better service for earthquake emergency and earthquake scientific research.     

Dynamic Stress Drop of Recent Earthquakes: Variations within Subduction Zones

—Stress drop is a fundamental parameter of earthquakes, but it is difficult to obtain reliable stress drop estimates for most earthquakes. Static stress drop estimates require knowledge of the seismic moment and fault area. Dynamic stress drop estimates are based entirely upon the observed source time functions. Based on analytical formulas that I derive for the crack and slip-pulse rupture models, the amplitude and time of the initial peak in source time functions can be inverted for dynamic stress drop. For multiple event earthquakes, this method only gives the dynamic stress drop of the first event. The Michigan STF catalog provides a uniform data base for all large earthquakes that have occurred in the past four years. Dynamic stress drops are calculated for the nearly 200 events in this catalog, and the resultant estimates scatter between 0.1 and 100 MPa. There is some coherent tectonic signal within this scatter. In the Sanriku (Japan) and Mexico subduction zones, underthrusting earthquakes that occur at the up-dip and down-dip edges of the seismogenic zone have correspondingly low and high values of stress drop. A speculative picture of the stress state of subduction zones emerges from these results. A previous study found that the absolute value of shear stress linearly increases down the seismogenic interface to a value of about 50 MPa at the down-dip edge. In this study, the dynamic stress drop of earthquakes at the up-dip edge is about 0.2 MPa, while large earthquakes at the down-dip edge of the seismogenic plate interface have dynamic stress drops of up to 5 MPa. These results imply that (1) large earthquakes only reduce the shear stress on the plate interface by a small fraction of the absolute level; and thus (2) most of the earthquake energy is partitioned into friction at the plate interface.     

西北地区中强震前固体潮调制比时空特征分析

基于1900—2015年中国大陆地区浅源地震阴历目录,依据地震活跃特征确定固体潮调制阴历日期和调制比阈值,并利用固体潮调制比时空扫描方法回顾性地研究近年来西北地区发生的3次中强地震,即2013年甘肃岷县漳县M_S6.6、2015年内蒙古阿拉善M_S5.8和2016年青海门源M_S6.4地震前, 区域小震的固体潮调制比时空演化过程.研究结果显示:岷县漳县M_S6.6和阿拉善M_S5.8地震前震中区域出现固体潮调制比异常,门源M_S6.4地震前震中附近的固体潮调制比异常演化时间较长; 西北地区3次中强震均发生在固体潮调制比时空扫描的高值异常集中地区及周围区域,扫描时间窗长度和固体潮调制比异常区域的大小可能与中强震的强度有一定的相关性.      

Strong earthquake activity all over the world and strong-moderate earthquake activity within and near China (June, 2001～July,2001)

~~Strong earthquake activity all over the world and strong-moderate earthquake activity within and near China (June, 2001～July,2001)@陈培善