四川芦山Ms7.0级强烈地震震源运动学特征

使用中国数字地震台网记录的区域宽频带波形,通过频率域和时间域多步反演,研究了2013年四川芦山“4·20”7.0级强烈地震的震源运动学特征.基于点源的震源机制解揭示：地震发震断层面参数分别为走向214°/倾角47°/滑动角96°,表现为一次高倾角的逆冲型事件.矩心在水平方向上位于震中（30.303°N/102.988°E）西南向约4.5 km,矩心深度约17 km.平均总标量地震矩M₀为1.16×10¹⁹ N·m,矩震级M_w约6.6.进一步模拟高达0.5 Hz高频波形,获得了芦山地震破裂过程图像,结果显示：此地震为一次不对称双侧破裂事件.破裂半径约15 km,整个破裂面积为706.7 km²,平均滑动量约0.231 m.破裂在8 s内释放了大多数能量.震后0~3 s内,破裂以孕震点为中心向四周同时扩展,3 s后,破裂表现出明显的方向性,主要向北北东扩展,导致位于震中北东向多数台站视破裂持续时间总体偏小,最小值为4 s.破裂约8 s后基本停止.     

Rupture pattern of the Oct 23, 2011 Van-Merkez,Eastern Turkey earthquake

High-frequency rupture process of the Oct 23, 2011 Van-Merkez earthquake is imaged by back-projection method using high-quality teleseismic P wave data from the US Array, and prestack Kirchhoff migration using P wave data from a subarray of global seismic networks. The rupture model with two asperities is confirmed by previous two methods. In low-frequency imaging, a large asperity derived from the migration method corresponds to the second one from the high-frequency P waves. The consistency of the locations of asperities from datasets with different frequency bands indicates that there is possible insignificance of the frequency-dependent feature for the earthquake. The resultant images illustrate the spatial and temporal evolution of the rupture, which mainly propagated WSW over a length of 33 km during the first 18 s, accompanying with bursts of two asperities at 3 and 11–13 s. The rupture direction is confirmed by the S wave corner frequency variations of strong ground accelerations. The rupture fronts are mainly located at the updip of the causative fault. Based on polarities of the P waveforms and focal mechanisms of the mainshock and aftershocks, the failure of these two asperities is determined to have occurred on a reverse fault with a dip angle of 47°. Hence, the rupture pattern of the 2011 Van-Merkez earthquake was dominated by a unilateral rupture toward the west-southwest direction.     

1999年台湾集集地震震源破裂过程

使用GPS同震位移资料和远场P波记录,研究了1999年台湾集集地震震源破裂过程．根据地质构造和余震分布引入了一个由弯曲断层面构造的三段“铲状”断层模型．在使用静态GPS位移资料反演集集地震的断层破裂滑动分布时,由于集集地震断层北部近地表破裂的复杂性,在位错模型中考虑拉张分量对地表同震位移的贡献,可更好地同时拟合GPS观测资料的水平和垂向分量．而纯剪切位错弹性半空间模型和分层地壳模型都无法同时拟合水平和垂向GPS观测资料．在此基础上,同时使用静力学同震位移资料和远场地震波形记录,反演集集地震的震源破裂过程．结果表明,一种垂直于断层面的“挤压性”(负)拉张分量几乎集中分布于地震断层的浅部和北部转折处,而这一带地表破裂远较没有(负)拉张分量出现的南部断层复杂．“冒起构造”的数字模拟表明,这种在集集地震破裂转折处及北部断裂带广为出现的典型破裂造成的地表位移可以用具有负拉张分量(挤压)的逆冲断层更好地模拟．而这种负拉张分量(挤压)的分布正是地震破裂性质和几何复杂性的综合反映,震源破裂过程也显示北部转折处破裂在空间和时间上的复杂性.高滑区域与余震分布表现为负相关.     

1976年MS7.8唐山地震断层动态破裂及近断层强地面运动特征

采用美国南加州地震委员会（SCEC）Steven Day博士提供的三维有限差分断层瞬态破裂动力学模型（3D-FDM）,以1976年唐山M_S7.8地震为例,从简化的断层双侧破裂模式出发,对该地震发震断层的动态破裂过程及近断层地表运动特征进行了仿真模拟和计算.研究区域为围绕发震断层200 km×140 km×40 km（深度）的长方形块体组成,模拟计算的空间分辨率和时间分辨率分别为200 m和0.012 s,形成的空间网格节点数为1051×701×201.在DELL小型工作站上,我们实现了对源程序的移植和并行计算.同时,通过引进计算机可视化技术,对模拟数据进行了3D/4D解释分析.另外,在对源程序修改过程中,实现了对京津唐地区三维地壳速度结构的嵌入,在一定程度上增强了对地震波传播以及地面运动模拟的真实性,并讨论了地震破裂的方向性对近断层地表运动的影响.最后根据初步研究结果结合京津唐地区活动断层构造特征,对唐山M_S7.8级主震后随之而来的1976滦县M_S7.1级余震及宁河M_S6.9级余震的动态触发机制提出了新的解释.由于受主震破裂方向性作用的影响,使得主震对后续两个较大余震产生的动态应力变化的峰值在断层的走滑方向上较大,为2~3 MPa,在逆冲方向上较小,为0.1~0.2 MPa.即唐山主震的发生使得其周边的应力场有一个瞬态的应力调整,唐山主震对后续余震的发生有促发作用.     

玉树地震地表变形InSAR观测及初步分析

采用玉树M_S7.1级地震前后两期PALSAR雷达数据(震前2010年1月15日, 震后4月17日)进行了“两轨+DEM”的InSAR处理, 获得了高质量的差分干涉雷达条纹图像和同震变形场。参考该区的基本构造格局, 根据干涉图像的变形范围、变形量和变形梯度可以初步判断:(1)玉树地震诱发了总体上NWW走向, 全长约70km地表陡变带, 陡变带南段位错及陡变梯度较大, 会在地表产生地表破裂;而西北部4段位错及陡变梯度较小, 不易在地表诱发破裂, 但可能在地下一定层位产生了隐伏破裂带;(2)陡变带两侧的雷达视线向运动方向预示发震断裂以左旋走滑运动为主;(3)宏观震中位于玉树县城西北约16km的地表陡变带上。D-InSAR解译结果与中国地震台网中心震源机制解、野外发震断裂调查结果及地貌特征吻合较好, 证明了干涉雷达解译成果的可靠性, 可以为准确定位玉树地震发震断裂地表行迹和快速评定震害损失提供有力的技术支持。      

2010年玉树地震(Ms7.1)地表破裂特征、破裂机制与破裂过程

2010年4月14日,青海省玉树地区发生Ms7.1级地震,造成大量人员伤亡和财产损失.地震发生后,我们对地震地表破裂带进行了详细的考察,并对同震位移量进行了精确的测量.根据野外考察和测量的结果,对玉树地震的地表破裂特征、同震位移量及其分布特征进行了分析,并对地震的破裂机制和破裂过程进行了探讨,取得如下认识:(1)玉树地震形成了沿鲜水河断裂带西北段(甘孜-玉树断裂)分布的东、西两条地表破裂带,西段破裂带分布在微观震中附近的隆宝湖拉分盆地中,长约19km;东段破裂带沿扎曲河南岸及巴塘河西岸山坡展布,长度约31km;上述两条破裂带之间存在约15km的地表破裂空区;(2)野外测量获得玉树地震的最大同震位移量为2.3m,位于东段地表破裂带中部郭央烟宋多附近;(3)地表破裂和野外构造地貌特征均反映了发震断层处于走滑伸展环境,断层左旋走滑过程中伴随正断作用;(4)地震波反演结果和地表破裂分布特征表明,玉树地震的破裂过程包括两次子事件,分别在地表形成了隆宝湖破裂带和扎曲河、巴塘河破裂带,隆宝湖及玉树县城西侧的山间谷地是在甘孜-玉树断裂长期活动的破裂带阶区转换拉张过程中形成的两个拉分盆地.     

利用近场高频GPS、 强地面运动和远场地震波形数据联合反演2008年汶川MS8.0地震的震源时空破裂过程

联合近场GPS测站1-Hz运动学位移、 强震仪加速度波形和全球台站P震相波形作为约束, 以时空滑动分布约束条件和ABIC模型参数选择方法, 结合先验的滑动方向变化范围, 反演2008年汶川M_S8.0地震的震源时空破裂过程, 给出了能够综合反映震源破裂过程的统一模型。 结果表明, 汶川地震总体上存在4个主要的破裂区, 最主要的一个破裂区位于震源东北40~120 km, 断层面上的最大位错量约为10 m, 主体滑动分布在2~20 km深度范围, 破裂达到地表; 第二个主体破裂区位于断层破裂带南段, 最大滑动量达到6 m; 另外2个主体滑动区位于断层破裂带北段, 但滑动破裂量小于断层南段破裂区的滑动量, 滑动破裂值最大值为4 m, 超过1 m的区域在走向上超过70 km。 反演得到的断层滑动模型的地震矩为9.5×10²¹ Nm, 相应的矩震级为M_W7.95。 汶川地震破裂表现为单侧破裂, 起始破裂在汶川下方16 km深度, 向东北方向一致性地传播, 过程持续~120 s。 在地震发生后0~10 s内, 破裂集中在震源起始破裂区, 滑动破裂值为~1.0 m, 之后破裂向东北方向扩展, 震后20~40 s是主要的破裂时段。 在40~60 s, 破裂跨越断层南段和北段。 在80~90 s破裂最大值开始下降, 在100~110 s时, 下降为~0.5 m, 在110~120 s时, 下降为~0.1 m。 加入近场GPS测站1-Hz 波形数据与近场强震仪波形和远场长周期体波联合反演, 提高了震源破裂模型的空间分辨率, 特别是浅部滑动破裂区的分辨率, 反演的最大滑动破裂值比不用1-Hz 波形数据反演的结果增大, 表明近场1-Hz GPS波形数据对于揭示汶川地震的时空破裂过程具有重要的作用。     

2018年斐济MW8.2深震:高频辐射过程和发生原因

2018年8月19日,在斐济东部海域563 km深处发生了一次M_W8.2地震.我们首先挑选位于美国阿拉斯加地区的131个宽频带台站构成台阵,选用垂直分量0.5~2 Hz的高频信号,利用广义台阵反投影技术对这次地震的破裂过程进行了成像,然后基于破裂速度对地震的辐射效率进行了估计.结果表明,这次地震总体上呈单侧破裂,破裂方位在3.0°左右,破裂总长度约51 km,持续时间22 s,平均破裂速度为2.5 km·s^-1.但能量释放有2次高峰,形成两次子事件.第一次为前10 s,峰值在7 s左右,破裂速度为2.9 km·s^-1,辐射效率为45%.第二次为10~22 s,峰值在15 s左右,破裂速度为1.6 km·s^-1,辐射效率为26%.结合震源位置、震源机制、破裂速度以及辐射效率,我们认为这次地震是由于俯冲板块前缘受到下部地幔物质上浮阻力引起的剪切失稳所致,起初板块内部的脆性破裂表现突出,致使辐射效率较高,后来震源处高温高压下的熔融耗散特征逐渐凸现,致使辐射效率下降.     

Rupture line and damage effects in Aqaba area, Egypt

Strong-motion parameters, peak ground acceleration, peak ground velocity, and peak ground displacement depend on several factors, such as the source of earthquake, distance between the source and site, and the characteristics of that site. Five seismographs and two accelographs were installed by the Egyptian Geological Survey team along the western side of the Gulf of Aqaba from Taba to Sharm El-Sheikh to record seismic events during the period from September 1995 to June 1996. During this period, two events were recorded by the accelographs. The two accelographs were located on the surface of the basement rocks, the first near the epicenter and the second at a farther distance. However, the farthest accelograph recorded higher values compared to the nearer one. Fault mechanics are an important factor in determining the values of strong motion parameters, where the direction of the rupture line plays an important role in detecting the values of strong motion parameters, the strong motion parameters and damage effects seem to be attenuated very fast in the direction perpendicular to the rupture line. This can be interpreted by the fact that the farthest accelograph lies at the extension of the fault rupture (azimuth = 30°), while the nearer one was perpendicular to the strike of the fault rupture.