科里奥利力对断层作用的统计研究

首先给出了分解到已知断层面法向和切向上的断层错动科里奥利应力(简称为科里奥利法向应力和切向应力)表达式,然后从哈佛大学矩心矩张量目录中选取主余震资料进行分析,按断层分类研究了科里奥利法向应力和主震震级与其最大余震震级差及主震地震矩与余震地震矩总和之差的折合矩震级间的关系.研究结果显示:地震断层错动过程中虽然产生了使断层两盘相互拉离(或挤压)的科里奥利法向应力,它降低(或增加)了断层错动时断层面上的摩擦阻力,但是应力量值太小(科里奥利法向应力估计最大不到0.1MPa,即不到一个大气压),不足以对断层错动及主震能量释放产生影响,从而影响余震的最大震级和总体水平.     

2018年9月4日新疆伽师MS5.5地震序列及发震构造讨论

2018年9月4日新疆伽师发生M_S5.5地震,震中处于塔里木地块西北缘,位于1997~1998年伽师强震群震区内。此次伽师地震前发生了M_S4.7前震,截至9月30日最大余震震级为M_S4.6（M_L5.0）,初步判定为前-主-余型地震序列。序列精定位结果显示,余震沿近NE向展布,主震震源深度与1997~1998年伽师强震主震基本一致,发震断层陡立。本文从区域的构造环境、地震震源机制解和余震分布特征等方面分析认为,地震发生在伽师隐伏断裂东南端部,为1997~1998年伽师强震群震区的一次新的构造活动。序列参数、视应力等计算结果显示,伽师M_S5.5地震的预测最大余震震级与最大余震震级M_S4.6接近,表明序列最大余震已经发生。     

强震等震线、余震区形状与地震构造关系的研究

本文通过对川滇和华北地区不同地震构造环境中, 1966年以来强震的等震线长宽比、 余震区长度的对比研究, 认为强震等震线的长宽比、 余震区长度都与地震构造环境和震级相关。 对于强震等震线的长宽比和余震区长度, 地震构造环境和震级的影响相同的是, 都与发震断裂是否位于活动地块边界, 位于哪一级活动地块边界, 是与边界方向一致的大断裂, 还是配套小断裂等有关。 最重要的是与活动断层的规模, 尤其是与发震断裂的长度有关, 也与震级有关。 不同的是, 强震等震线的长宽比明显受地表覆盖层的影响, 而余震区长度与地表覆盖层关系不明显, 似乎与震源破裂面长度的关系更密切。     

台湾集集和华北大同地震序列的定量模型分析

介绍了传染型余震序列(ETAS)模型的基本原理和发展概况，并利用该模型分析了1999年开始的台湾集集地震序列(主震Ms7．6)和1989年开始的大同地震序列(主震为Ms5．7)。通过分析计算，得到了二个序列的ETAS参数向量，并结合各自的孕震环境和机制分析了二个地震序列的异同和序列衰减规律。     

2015年尼泊尔MW7.9地震应力状态与余震空间分布规律

本文利用2015年尼泊尔M_W7.9地震断层面滑动位移分布的运动学反演结果,通过傅里叶变换法得到了主断层面上的两分量应力状态,并研究了余震的空间分布和断层面上应力状态之间的关系.发现滑动位移分布与应力状态分布都相对较为集中,大约70%的余震分布在应力变化为正的区域,而其余发生在应力降区域的余震,又大多发生在应力变化梯度较大的地区.为了得到一个更符合实际的滑动模型来解释余震的触发机制,我们计算了波数域中滑动位移和应力状态的傅里叶谱,发现此次地震的滑动位移和应力状态近似满足k^-3和k^-2衰减.我们利用简化的圆盘模型说明了非均匀应力变化下的衰减过程,计算了圆盘模型的有效半径r_e约等于0.7倍的圆盘半径.这就说明圆盘模型中应力增加的部分应该占整个圆盘破裂面积的51%.在本次尼泊尔M_W7.9地震实例中,断层面上应力状态为负的区域比滑动位移为正的区域有了明显地缩小.事实表明,余震可以发生在有滑动位移的区域,非均匀应力降模型比均匀应力降模型更加接近真实的震源破裂过程.     

云南禄劝地震对余震有触发作用吗--静态应力触发地震的讨论

采用文献[1]计算地震产生的静态应变场的解析表达式,对云南禄劝地震的15次余震,计算了主震在余震节面上产生的正应力和静态库仑破裂应力变化。分析了正应力变化和静态库仑破裂应力变化是否对余震有触发作用。结果表明,主震产生的正应力变化和库仑破裂应力变化并没有明显的证据“有利于”余震的发生。另外,还比较了运用与主震相同的余震机制和反演得到的余震机制在库仑破裂模型中计算结果的差异。使用反演的余震机制得出的结果稍优于假定余震与主震有相同震源机制得到的结果。     

The large aftershocks triggered by the 2011 Mw 9.0 Tohoku-Oki earthquake,Japan

The Mw 9.0 Tohoku-Oki earthquake that occurred off the Pacific coast of Japan on March 11, 2011, was followed by thousands of aftershocks, both near the plate interface and in the crust of inland eastern Japan. In this paper, we report on two large, shallow crustal earthquakes that occurred near the Ibaraki-Fukushima prefecture border, where the background seismicity was low prior to the 2011 Tohoku-Oki earthquake. Using densely spaced geodetic observations (GPS and InSAR datasets), we found that two large aftershocks in the Iwaki and Kita-Ibarake regions (hereafter referred to as the Iwaki earthquake and the Kita-Ibarake earthquake) produced 2.1 m and 0.44 m of motion in the line-of-sight (LOS), respectively. The azimuth-offset method was used to obtain the preliminary location of the fault traces. The InSAR-based maximum offset and trace of the faults that produced the Iwaki earthquake are consistent with field observations. The fault location and geometry of these two earthquakes are constrained by a rectangular dislocation model in a multilayered elastic half-space, which indicates that the maximum slips for the two earthquakes are 3.28 m and 0.98 m, respectively. The Coulomb stress changes were calculated for the faults following the 2011 Mw 9.0 Tohoku-Oki earthquake based on the modeled slip along the fault planes. The resulting Coulomb stress changes indicate that the stresses on the faults increased by up to 1.1 MPa and 0.7 MPa in the Iwaki and Kita-Ibarake regions, respectively, suggesting that the Tohoku-Oki earthquake triggered the two aftershocks, supporting the results of seismic tomography.     

Analysis of aftershock sequences in South and Southeastern Spain

A probabilistic modeling is used to analyze the spatio-temporal behavior of eleven aftershock sequences occurred in South and Southeastern Spain. This study focuses on the analysis of two seismicity parameters: the b-value of the frequency-magnitude distribution, and the p-value, explaining the temporal decay rate of aftershocks. The estimated b values range between 0.77 ± 0.05 and 1.18 ± 0.10 close to the typical b-values of the aftershock frequency-magnitude relationship b ≈ 1.0. The estimated p-values range between 0.75 ± 0.03 and 1.43 ± 0.10 showing broad regimes of the temporal decay of aftershocks. The modified Bath’s law used to analyze the energy partitioning, suggests that a large fraction of the accumulated energy is released in the mainshock and relatively small fraction of energy is released during aftershock sequence, for example 80% of the total energy is released during the Mula 1999 mainshock, 88% during Bullas 2002 mainshock and 87% during La Paca 2005 mainshock. The fractal dimension D₂ is estimated using the correlation integral, and then used to derive the slip ratio, as the ratio of the slip occurred on primary fault segment to the total slip. For example, we obtained a slip ratio equal to 71% for the Mula 1999 aftershock sequence, 61% for the Bullas 2002 event, 58% for the La Paca 2005 aftershock, 50% for the Lorca 2011 sequence and 63% for the sequence triggered by the Gador 2002 mainshock.Finally, the correlations between the fractal dimension, the b-value and the p-value is analyzed, and the Aki’s relation D = 3b/c is discussed as well.     

Constrained edge dynamic deleting in CD-TIN based on influence domain retriangulation of virtual point

Constrained Delaunay triangulated irregular network is one kind of dynamic data structures used in geosciences. The research on point and edges insertion in CD-TIN is the basis of its application. Comparing with the algorithms of points and constrained edge insertion, there are very a few researches on constrained edge deletion in CD-TIN. Based on the analysis of the polymorphism of constrained edge, virtual points are used to describe the intersection of constrained edges. A new algorithm is presented, called as influence domain retriangulating for virtual point (IDRVP), to delete constrained edges with virtual points. The algorithm is complete in topology. Finally, the algorithm is tested by some applications cases.     

2017年九寨沟7.0级地震断层破裂尺度初步研究

本文利用2017年九寨沟7.0级地震的余震序列数据,结合震源机制解确定的地震断层面,研究了地震断层的破裂尺度。基于余震序列在断层面上的投影,分析直接余震和间接余震以及不同起始震级条件下余震区的长度和震源深度分布,估算了地震断层破裂面的破裂长度和破裂深度,并探讨了地震断层破裂对地震烈度分布的影响。研究结果显示,九寨沟7.0级地震断层的破裂长度约33—35km,破裂深度约23—26km;地震断层破裂对Ⅷ度区分布有明显的控制作用,地震断层破裂长度接近Ⅷ度区的长度。