利用InSAR资料反演2008年西藏改则Mw6．4和Mw5．9地震的断层参数

2008年1月9日在我国西藏改则发生了一次MW6.4地震,随后有40次3.5级以上余震发生,其中最大的一次为1月16日的MW5.9余震.本文处理了ENVISAT ASAR两轨(升轨和降轨)同震资料,精确确定了同震地表位移的空间分布;随后利用弹性半空间的位错模拟确定了上述事件的断层面参数;最后,基于非均匀滑动模型反演确定了两次地震断面上的滑动分布.结果表明,MW6.4主震断层为走向218°、倾角52°的西倾断层,最大滑动量约1.9 m,出现在地表以下约7.6 km处;而MW5.9余震发生在主震断层西3.2 km的地方,发震断层为走向200°、倾角59°的西倾断层,最大滑动量约1.0 m,出现在地表以下约3.9 km处.     

利用InSAR资料反演2008年西藏改则MW6.4和MW5.9地震的断层参数

2008年1月9日在我国西藏改则发生了一次M_W6.4地震,随后有40次3.5级以上余震发生,其中最大的一次为1月16日的M_W5.9余震.本文处理了ENVISAT ASAR两轨（升轨和降轨）同震资料,精确确定了同震地表位移的空间分布;随后利用弹性半空间的位错模拟确定了上述事件的断层面参数;最后,基于非均匀滑动模型反演确定了两次地震断面上的滑动分布.结果表明,M_W6.4主震断层为走向218°、倾角52°的西倾断层,最大滑动量约1.9 m,出现在地表以下约7.6 km处;而M_W5.9余震发生在主震断层西3.2 km的地方,发震断层为走向200°、倾角59°的西倾断层,最大滑动量约1.0 m,出现在地表以下约3.9 km处.     

玉树MS7.1级地震部分余震重新定位及发震构造分析

综合利用玉树震区应急流动台站观测数据和青海地震台网固定台站观测数据,依据最新的人工地震宽角反射/折射剖面的速度模型,采用Hypo2000地震定位法,对2010年4月18日至4月29日期间玉树震区发生的部分余震进行了重新定位.重新定位后,震源位置的水平和垂直方向平均误差分别为1.35 km和4.68 km,走时残差为0.49 s.震源深度分布范围为1.48～19.85 km,平均震源深度为10.28 km.定位研究结果表明:玉树地震余震沿北西-南东向的甘孜-玉树断裂带的北支,即玉树-隆宝断裂分布,长约97 km.余震分布特征在主震(微观震中)两侧存在差异,可能反映了两侧构造特征存在差异.截止到4月29日,主震东南仍是应力的主要释放区域,余震强度大且活动密集的区域位于主震东南距主震约5 km、横向范围约20 km.主震破裂区的大部分应力在主震过程中得以释放,主震时应力未释放的区域成为主要的余震分布区.余震的连续发生可能已造成主震破裂区相互连通,且破裂范围向西北方向扩展.玉树主震及余震的发震构造为甘孜-玉树断裂的北支,即玉树-隆宝断裂段,断层性质为北东倾向的高角度左旋走滑断层.发震断层的倾角和宽度在帮洞两侧有所不同,帮洞以东发震断层宽度约为12 km,倾角约为83°;而帮洞以西发震断层宽度约为6.5 km,断层倾角约减缓为63°.     

由初至P震相重新定位2017年九寨沟地震序列的主震与ML≥3.0余震并分析发震构造

综合利用区域台站和流动台站(近台)的记录,基于初至P震相重新测定了2017年九寨沟序列M_S7.0主震和M_L≥3.0余震的震源位置,并利用较高精度的定位结果分析余震分布与地震构造的关系,解释发震断裂带的结构.获得的新认识有:(1)九寨沟主震震源深度为16km,位于余震带中段的南缘;余震主要分布深度为4～17km.(2)沿余震带的走向,余震分布与主震同震位错大小的分布明显相关.余震带中段8～16km深度存在的余震稀疏区与同震位错的高值区相吻合,应是发震断裂带主凹凸体的部位,也是主震时应变释放较充分的部位;余震带南东段10～18km深度的余震密集区对应了同震位错的亏损区之一,三次M_L≥5.0余震都发生于此;余震带西北段在5～10km之下既缺少余震,又属同震位错的亏损区,可能与那里多条断裂的交汇或合并造成的构造复杂性有关;余震带中-北西段3～5km深度的也缺少余震,也对应了浅部的同震位错亏损区.(3)证实了九寨沟地震的发震构造为虎牙断裂带北段,同时新揭示出发震断裂带表现为由主断裂和分支断裂构成的、向上分叉的花状结构,尺度约为4.5km宽(最大)、35km长,主断裂朝SW陡倾.这些反映主震破裂可能不只受控于单一的断裂,而有可能是沿主断裂发生主破裂,而沿分支断裂发生次要破裂.另外,本文对发震断裂带结构的分段解释,是遵循构造地质学原理去综合震源排列、震源机制解、地表断层已知位置、相邻剖面断层解释结果等信息的分析结果,而不仅仅依据余震的密集分布进行推断.     

1973年炉霍大地震(Ms=7.6)最大余震(Ms=6.3)的库仑破裂应力触发

1973年2月6日在四川省炉霍发生Ms7.6大地震，震后在震中周围出现丰富的余震，其最大余震（Ms6.3）发生在鲜水河断裂和玉树－甘孜断裂之间的正断层活动区内。本文根据震源机制解，地表破裂，同震位错分布，地震地质构造等资料，建立炉霍走大地震破裂的弹性位错模型，计算该大地震在周围正断层滑动方向上引起的库仑破裂应力变化（△CFS），结果表明，最大余震发生在库仑破裂应力增加（△CFS＞0）的地区，△CFS＝4.5MPa,因此认为，最大余震可能是由主震引起的库仑破裂应力变化触发的。     

基于余震重定位和震源机制解研究青海玛多MS7.4地震序列的发震构造和断裂形态

2021年5月22日02时04分(北京时间),青海果洛州玛多县发生M_S7.4地震,震后余震不断.地质调查和卫星观测对地表断裂痕迹有较好的约束.然而,对于理解区域应力场、地震的产生、传播和终止具有重要意义的地下断层几何结构的约束精度略显不足.利用国家地震台网的连续波形记录,本研究首先基于双差定位法对玛多地震震后25天的余震序列进行重定位,结果显示余震序列大致沿NWW向的江错断裂呈线性分布,位于主震震中两侧,延伸总长~170 km.主震东南侧存在一余震稀疏区,在断裂带东西两端余震分布转向且出现分叉现象,反映出发震断层的复杂几何形态,这与前人研究结果基本一致.进一步采用波形反演方法和P波初动极性反演方法,获得了玛多震源区132个中小余震的震源机制解与震源矩心深度,并基于此对该主余型地震的发震构造与断裂形态进行了初步分析.震源机制解结果表明,玛多M_S7.4主震的发震断裂主要为左旋走滑性质,余震与主震性质整体相同,在断裂带东段存在部分逆冲型余震.震源机制解约束的区域主应力方向约N60°E,与区内整体走滑断裂作用相一致.余震震源深度略微起伏,主要集中在10~12 km,且浅部余震较少,表明浅部应力可能主要通过主震释放,余震深度分布可能限定了主震同震破裂的下边界.玛多主震破裂起始于断裂带走向和倾向发生明显变化的位置,表明断裂带的复杂几何结构可能是此次玛多M_S7.4地震初始破裂空间分布的决定因素.主震破裂结束的两端都有"马尾状"构造(或次级断层),表明这种分叉断层复杂的几何形态可能控制着主震破裂的最终位置.     

玉树地震震源区速度结构与余震分布的关系

利用玉树震区21个应急流动地震台站和青海省地震台网固定地震台站的观测数据,采用双差层析成像方法,对2010年4月14日至6月15期间发生的地震进行了重定位,并反演得到了玉树地震震源区的三维速度结构.重定位结果揭示余震主要沿NW向成窄带状分布在断层的两侧,表明脆性破裂应力释放主要集中于一个狭窄的区域内.在西北端,余震偏离玉树—甘孜断裂分布,在SW向也有分布,推测可能与南西向次级断裂有关.双差层析成像得到的速度结构在浅部与地表地质构造相一致,中上地壳的速度结构显示巴颜喀拉地块为高速异常,羌塘地块为低速异常.玉树地震余震分布与特定的速度结构存在相关性：主震发生在高低速过渡带偏高速体的一侧,余震主要分布在高速体外围,高速体内部几乎没有余震分布.一般说来,中上地壳的高速体通常具有较高的强度,可以积累较强的孕震能量.主震发生后,高速体内积累的弹性能量向周边释放,可能是导致高速体周边余震发生的主要原因.     

2010年MS7.1级玉树地震同震库仑应力变化以及对2011年MS5.2级囊谦地震的影响

2010年4月14日青海省玉树藏族自治州发生M_S7.1级地震.和传统的板内地震相比,玉树M_S7.1级地震的余震具有数量少、震级大的特点.研究玉树地震主震与余震之间的关系,对于我们了解余震的发震机理具有十分重要的参考价值.本文利用弹性位错理论和分层岩石圈模型,计算玉树地震引起的同震及震后黏弹松弛应力场变化,讨论M_S7.1级玉树地震对余震分布的影响以及与2011年囊谦M_S5.2级地震之间的触发关系.结果显示,玉树地震导致了四处明显的库仑应力增强的扇区,2010年4月13日至6月17日的870次M_L>1.0级余震主要分布于主震破裂面附近区域以及破裂面东北端的应力增强扇区.分析玉树地震对余震分布的影响时,有效摩擦系数以及计算深度的选取对计算结果的影响较小,是否考虑区域构造应力场的影响较大.考虑区域构造应力场时,占总数86.7%的余震位于库仑应力增强区,地震应力触发理论较好地解释了余震的分布.选取囊谦地震震源机制解的两个节面作为库仑应力计算中的接收断层参数,并且考虑不同黏滞系数下的玉树地震同震及震后黏弹松弛效应,模型计算结果均表明囊谦地震位于玉树地震所导致应力影区,仅依靠地震的静态、震后黏弹松弛应力触发理论,无法解释囊谦地震的发生,说明该次地震可能是一次独立的事件.     

利用D-InSAR技术研究西藏改则地震同震形变场

针对2008年1月9日MW6.4西藏改则地震和2008年1月16日的MW5.9余震,通过两通(2-pass)加外部DEM差分干涉处理技术(D-InSAR),提取了地震区域2次地震累积的视线向(LOS)同震形变场。结果表明:发震断层均为正断层,位于依布茶卡-日干配错断裂端点附近。主震发震断层走向为N30°E,余震发震断层走向为N21°E,两断层距离约7km;在影像上主震发震断层有造成地表破裂的痕迹,余震未见地表破裂的痕迹;这次地震造成的同震形变场长约30km,宽约20km,主震断层上盘和下盘视线向最大形变量分别为39.2cm和11.2cm,两盘相对位错达50.4cm,余震造成的视线向形变量为9.4cm     

2015尼泊尔M_S8.1地震中等余震震源机制研究

本研究利用西藏台网记录的波形数据,采用gCAP方法反演了2015年4月25日尼泊尔MS8.1大震5次中等余震(5.0≤MS≤6.5)及西藏定日MS5.9地震震源机制解.结果显示,6次地震包含2个正断、2个走滑及2个逆冲型地震.其中2个正断型地震位于主震的东北方向,即发震断层的上盘,表明该区域受到主震同震位移的影响,表现出应力拉张的变化特征;2个走滑型地震在主震破裂的东南方向上,说明随着破裂往东南方向延伸,余震的走滑分量增强;另外2个逆冲型地震位于5月12日MS7.5强余震区域,与MS7.5地震的滑移状态一致,可能与主震同震位移引起该区域处于应力挤压状态密切相关.这些结果表明,尼泊尔MS8.1主震发生后,由于同震位移的影响,不同区域处于不同的应力状态,从而使中等余震表现出不同的震源类型.     

近断层应力变化快速计算方法及对强余震空间分布的指示意义:以M_w7.9汶川地震为例

构造地震一般由断层摩擦失稳所致.断层内部及周边所累积的剪切形变则通过同震滑动位移得到局部释放.因此,震后断层内部及近断层周边的静态剪切应力变化量的空间分布可通过断层面上的滑动位移分布计算得到.本文采用傅氏变换方法(FTM)计算单一有限断层同震滑移场所形成的静态剪切应力变化分布,近断层两侧的应力变化计算可由波数域内应力近似算法获得.结果表明,FTM快速有效、易于实现,有效地避免了常规应力计算中奇异值的出现.以2008年Mw7.9中国汶川大地震为例,采用前人所得有限断层滑动位移模型,得到了断层面和近断层周边准3D剪切应力分布解,并同主震后中强余震的空间分布特征作了比较.结果表明,大部分的中强余震震源位置处于剪切应力变化值为正的区域,由同震位移所产生的局部应力降峰值和均值大小同应力变化的正值大致相当,从而表明了快速且有效地计算断层内部及近断层附近的应力变化分布可以为主震后强余震发生的潜在区域提供指示意义.需要强调的是,应力变化空间分布特征的获取强烈地依赖于有限断层滑移模型解.有关滑动位移反演解的多解性对应力变化计算结果的影响,本文作了必要的讨论.     

Fault complexity associated with the 14 August 2003 M<Subscript>w</Subscript>6.2 Lefkada,Greece, aftershock sequence

The M _w 6.2 Lefkada earthquake occurred on 14 August 2003 beneath the western coastline of Lefkada Island. The main shock was followed by an intense aftershock activity, which formed a narrow band extending over the western coast of the Island and the submarine area between Lefkada and Kefalonia Islands, whereas additional off fault aftershocks formed spatial clusters on the central and northwestern part of the Island. The aftershock spatial distribution revealed the activation of along-strike adjacent fault segment as well as of secondary faults close to the main rupture. The properties of the activated segments were illuminated by the precisely located aftershocks, fault plane solutions determination and the cross sections performed parallel and normal to their strike. The aftershock focal mechanisms exhibited mainly strike slip faulting throughout the activated area, although deviation of the dominant stress pattern is also observed. The results help to emphasize the importance of the identification of activated nearby fault segments possibly triggered by the main rupture. Because such segments are capable to produce moderate events causing appreciable damage, they should be viewed with caution in seismic hazard assessment in addition to the major regional faults.     

Numerical simulation on the influences of Wenchuan earthquake on the surrounding faults

On 12 May 2008, the devastating Wenchuan earthquake struck the Longmenshan fault zone, which comprised the eastern margin of the Tibetan Plateau, and this fault zone was predominantly a convergent boundary with a right-lateral strike-slip component. After such a large-magnitude earthquake, it was crucial to analyze the influences of the earthquake on the surrounding faults and the potential seismic activity. In this paper, a complex viscoelastic model of western Sichuan and eastern Tibet regions was constructed including the topography. Based on the findings of co-seismic static slip distribution, we calculated the stress change caused by the Wenchuan earthquake with the post-seismic relaxation into consideration. Our preliminary results indicated that: (1) The tectonic stressing rate was relatively high in Kunlun mountain pass-Jiangcuo, Ganzi-Yushu, Xianshuihe and Zemuhe faults; while in the east Kunlun and Longriba was medium; also the value was less in the Minjiang, Longmenshan, Anninghe and Huya faults. As to the Longmenshan fault, the value was 0.28×10^-3 MPa/a to 0.35×10^-3 MPa/a, which is coincident with the previous long recurrence interval of Wenchuan earthquake; (2) The Wenchuan earthquake not only caused the Coulomb stress decrease in the source region, but also the stress increase in the two terminals, especially the northeastern segment, which is comparatively consistent with the aftershock distribution. Meanwhile, the high concentration areas of the static slip distribution were corresponding to the Coulomb stress reductions; (3) The Coulomb stress change caused by Wenchuan earthquake showed significant increase on five major faults, which were northwestern segment of Xianshuihe fault, eastern Kunlun fault, Longriba fault, Minjiang fault and Huya fault respectively; also the Coulomb stress on the fault plane of the Yushu earthquake was faintly increased; (4) We defined the recurrence interval as the time needed to accumulate the magnitude of the stress drop, and the recurrence interval of Wenchuan earthquake was estimated about 1 714 a to 2 143 a correspondingly.     

2008年于田MW7.2地震滑动分布反演及三维同震形变场解算

基于ENVISAT ASAR升降轨数据, 利用InSAR获取2008年于田M_W7.2地震同震形变场; 采用SDM反演本次地震断层滑动分布; 使用PSCMP正演获取于田M_W7.2地震南北向模拟形变量, 并结合升降轨同震形变场, 解算三维同震形变场。 同震形变场分析表明, 2008年于田M_W7.2地震以正断为主, 且带有走滑运动特征, 破裂带走向为NNE向。 同震滑动分布反演结果显示, 断层沿走向被分为4段F1、 F2、 F3、 F4, 其滑动分布集中在0~14 km区间, 以F2、 F3段为主, 最大滑动量约5.31 m, 位于F2段深部2.76 km处; 沿破裂带走向, 左旋走滑位移与垂直位移比值有增大的趋势; 反演获得的地震矩M₀=5.58×10¹⁹N·m, 相当于矩震级M_W7.1。 三维同震形变场解算结果显示, 断层上盘整体表现为沉降, 断层下盘整体表现为隆升, 且沉降量明显大于隆升量, 表明地震以正断破裂为主; 除靠近断裂带中上部表现为向东南运动外, 上盘整体上表现为向西南运动; 断层下盘则整体表现为向东北运动, 证明破裂兼有左旋走滑运动。 滑动分布反演、 正演与三维同震形变场解算结果皆表明, 于田M_W7.2地震破裂以正断为主, 且带有一定的左旋走滑。     

The Rupture Process of the Manjil,Iran Earthquake of 20 June 1990 and Implications for Intraplate Strike-Slip Earthquakes

In terms of seismically radiated energy or moment release, the earthquake of 20 January 1990 in the Manjil Basin-Alborz Mountain region of Iran is the second largest strike-slip earthquake to have occurred in an intracontinental setting in the past decade. It caused enormous loss of life and the virtual destruction of several cities. Despite a very large meizoseismal area, the identification of the causative faults has been hampered by the lack of reliable earthquake locations and conflicting field reports of surface displacement. Using broadband data from global networks of digitally recording seismographs, we analyse broadband seismic waveforms to derive characteristics of the rupture process. Complexities in waveforms generated by the earthquake indicate that the main shock consisted of a tiny precursory subevent followed in the next 20 seconds by a series of four major subevents with depths ranging from 10 to 15 km. The focal mechanisms of the major subevents, which are predominantly strike-slip, have a common nodal plane striking about 285°–295°. Based on the coincidence of this strike with the dominant tectonic fabric of the region we presume that the EW striking planes are the fault planes. The first major subevent nucleated slightly south of the initial precursor. The second subevent occurred northwest of the initial precursor. The last two subevents moved progressively southeastward of the first subevent in a direction collinear with the predominant strike of the fault planes. The offsets in the relative locations and the temporal delays of the rupture subevents indicate heterogeneous distribution of fracture strength and the involvement of multiple faults. The spatial distribution of teleseismic aftershocks, which at first appears uncorrelated with meizoseismal contours, can be decomposed into stages. The initial activity, being within and on the periphery of the rupture zone, correlates in shape and length with meizoseismal lines. In the second stage of activity the aftershock zone expands and appears to cluster about the geomorphic and geologic features several tens of kilometres from the rupture zone. The activity is interpreted as a regional response to quasistatic stress migration along zones of tectonic weakness. The radiated energy of the main shock and the estimate of seismic moment yields an apparent stress of 20 bars. High apparent stress may be typical of strike slip earthquakes occurring in intracontinental environments undergoing continental collision.     

Research on the seismotectonics of the January 17, 1995 Hanshin M 7.2 earthquake

Based on the geological tectonics, aftershock activity, earthquake surface rupture and peak ground motion, the geometric and dynamic characteristics of seismogenic tectonics about the 1995 Hanshin earthquake are analysed. Nojima fault and Rokko fault have the same trending direction, but opposite dips. Their rising and falling plates are in symmetrically diagonal distribution. The two faults can be defined as thrust-strike slip faults and constitute a pivotal strike-slip fault. The earthquake just occurred at the pivot, which is the seismotectonics for the earthquake to develop and occur. The pivotal movement along a strike-slip fault often leads to the occurrence of large earthquakes, whose dynamic process can be demonstrated by the stress analysis on the torsion of a beam with rectangle section. The displacement of earthquake surface rupture, aftershock density and peak acceleration change in a certain range of epicentral distance just similar as the shear stress changes from the center to the sides in the rectangle section. The distribution characteristics of the heaviest damage areas are also discussed in the article from the aspects of special geological tectonics and seismotectonic condition. The result obtained from the article can be applied not only to realizing the potencial earthquake sources in middle-long time, but also to build reasonably the prediction model about earthquake hazard.     

2010年4月玉树M_S7.3地震序列的断层结构

利用双差定位方法对玉树地震序列2010年4月14日至10月31日间发生的ML≥1.0地震进行双差定位,得到1545个地震的重定位结果.综合分析地震双差定位结果和玉树地震序列中强地震震源机制解,发现玉树MS7.3地震发震构造由北西向和北东东向两条相交断层组成,主震发生在北西走向的甘孜—玉树断裂带上,5月29日的MS5.9余震序列发生在北东东走向的一条隐伏断裂上,两条断裂均接近直立.甘孜—玉树断裂是羌塘地块和巴彦喀拉地块的构造边界,由于羌塘地块和巴颜喀拉地块的差异运动使甘孜—玉树断裂强耦合段应力高度积累,在应变能超过岩石强度时破裂失稳发生了MS7.3地震.主震断层的左旋滑动导致北东东向断层的正应力减小,库伦应力增加,45天后触发了MS5.9余震序列的活动.     

GPS和InSAR同震形变约束的尼泊尔M_W7.9和M_W7.3地震破裂滑动分布

2015年4月25日尼泊尔爆发MW7.9地震,继而引发5月12日MW7.3级余震,GPS、InSAR监测到震源区及周边大范围同震形变.本文以国内外的GPS和InSAR同震形变为约束,考虑喜马拉雅断裂带岩石圈垂向分层和横向差异的影响,反演主喜马拉雅逆冲断裂在这次主震和余震中破裂面形状和滑动分布.结果显示,主震从USGS确定的震中位置向东偏南延伸100km以上,破裂地面迹线与主前缘逆冲断裂迹线基本一致.破裂面倾角约7°~11°,大部分破裂集中在深度8~20km,同余震分布深度一致.主震最大滑动量约6.0~6.6m,位于14km深处.余震破裂集中在震中附近30km范围内,填补了主震东部破裂空区,最大滑动约3.6~4.6 m,位于13km深.深度20km以下基本没有破裂.地壳介质不均匀性对破裂滑动分布的影响较大,介质不均匀模型的观测值不符值比各向同性弹性半空间模型降低10%以上.本文地震破裂模型特征与地震反射剖面、以及根据震间期大地测量数据反演的喜马拉雅深部蠕滑剖面极其相似.跨喜马拉雅断裂剖面的震间形变量与地震破裂滑移量直接相关.以此推算,尼泊尔中部大震原地复发周期在300年以上.