THE STATIC STRESS TRIGGERING INFLUENCES OF THE 2015 MW6.4 PISHAN,XINJIANG EARTHQUAKE ON THE NEIGHBORING AREAS

Based on the rupture models of the 2015 Pishan M_W6.4 earthquake and half space homogeneous elastic model, the Coulomb stress changes generated by the earthquake are calculated on the active faults near the earthquake region. The horizontal stress changes and the displacement field are estimated on the area around the epicenter. Results show that:(1)The Coulomb stress is decreased in the west of the western Kunlun frontal thrust fault(9.5×10³Pa), and increased in the east of the western Kunlun frontal thrust fault and the middle of the Kangxiwa faults. More attention should be taken to the seismic rick of the east of the western Kunlun frontal thrust fault; (2)Based on the analysis on the location of the aftershocks, it is found that most of the aftershocks are triggered by the earthquake. In the region of increased Coulomb attraction, the aftershock distribution is more intensive, and in the area of the Coulomb stress reduction, the distribution of aftershocks is relatively sparse; (3)The horizontal area stress increases in the north and south of the earthquake(most part of the Qaidam Basin and the northwest of the Qinghai-Tibet plateau), and decreases in the east and west of the earthquake(northern part of the Qinghai-Tibet plateau and eastern part of the Pamir Mountains). In the epicenter area, the principal compressive stress presents nearly NS direction and the principal extensional stress presents nearly EW direction. The principal compressive stress shows an outward radiation pattern centered on the epicenter with the principal extensional stress along the direction of concentric circles. The principal compressive stress presents NW direction to the west of the epicenter, and NE to the east of the epicenter. With the increase of radius, the stress level gradually decays with 10⁷Pa in the epicenter and hundreds Pa in the Maidan Fault which is in the north of the Qaidam Basin.     

2017年九寨沟MS7.0地震对周围地区的静态应力影响

基于九寨沟M_S7.0地震的破裂模型及均匀弹性半空间模型,本文计算了该地震在周围主要活动断层上产生的库仑应力变化、在周围地区产生的应力场和位移场和同震库仑应力变化对余震的触发.结果表明：（1）九寨沟地震造成虎牙断裂中段库仑破裂应力有较大增加,已经超过0.01 MPa的阈值,虎牙断裂北段、塔藏断裂中段和岷江断裂北段北部的库仑破裂应力有较大降低,因此尤其要注意虎牙断裂中段的危险性.（2）水平面应力场在该地震震中东西两侧增加（拉张）,张应力起主要作用.在震中南北两侧降低（压缩）,压应力起主要作用.从水平主压和主张应力方向来看,均呈现出条形磁铁的磁场形态.从剖面上的应力场来看,在上盘的面膨胀区域内,大部分点的主张应力方向与地表是垂直的,在其他区域内,主张应力和主压应力均以震中为中心,向外呈辐射状.（3）从地表水平位移场来看,震中东西两侧物质朝震中位置汇聚,南北两侧物质向外流出,在震中处的最大水平位移量达43 mm.从地表垂直位移场来看,震中南北两侧出现明显的隆升,隆升最大值达56.8 mm.震中东西两侧出现明显的沉降,沉降最大值达74.5 mm.从剖面的位移场来看,九寨沟地震为左旋走滑地震,且有一定的正断成分.由分析可以推测该断层破裂在大致22~26 km的深度上就截止了.并推测下盘物质运动的动力来自震源北东东方向（四川块体）深度在6~30 km的下盘下层物质,上盘物质运动的动力来自震源北西西方向（巴颜喀拉块体）深度在0~6 km的上盘上层物质.（4）通过计算不同深度上主震对余震的触发作用可知,主震后的最大余震受到了主震的触发作用,多数其他余震也受到主震的触发作用.主震的发生促使了库仑应力增加地区余震的发生,抑制了一部分库仑应力减少地区余震的发生.     

2015年尼泊尔强震序列对中国大陆的应力影响

基于2015年尼泊尔地震序列的破裂模型及均匀弹性半空间模型,计算了该地震序列传递到中国西藏境内发生在定日县地震和聂拉木县地震的应力.2015年尼泊尔地震序列导致定日县地震和聂拉木地震节面和滑动方向的库仑应力增加(2~3)×103 Pa和(2.4~3.1)×105 Pa,表明这两个地震受到尼泊尔地震序列的触发.其次,我们计算了2015年尼泊尔地震序列在中国大陆及其附近主要活动断层上产生的库仑应力变化.喜马拉雅主山前逆冲断裂和青藏高原内部的拉张正断层上的库仑应力有较大的增加,而青藏高原的走滑断裂,如阿尔金断裂、东昆仑断裂、玉树玛曲断裂、班公错断裂西部、嘉黎断裂的库仑应力有较大的降低.天山南北两侧的断裂库仑应力降低.而华北及东北、华南地区的库仑应力变化几乎可以忽略不计.最后,计算了该地震序列造成的水平应力变化.水平面应力在2015年尼泊尔地震序列北向(青藏高原大部和新疆区域)增加(拉张),而在地震序列东侧的西藏南部和川滇地区南部降低(压缩),在华北和东北仅有少许增加,在华南地区有少许降低.在中国西部,主压应力表现为以2015年地震序列为圆心的向外辐射状,而主张应力方向与同心圆切线方向大体一致.水平主压应力方向在东北地区为北东向,在华北地区为北东东向,在华南地区为南东东向.这种模式与现今构造应力场方向相似,表现了2015尼泊尔地震序列所代表的印度板块和欧亚板块的碰撞是中国大陆构造变形的主要动力来源.     

基于GNSS与InSAR约束的九寨沟MS7.0地震滑动模型及其库仑应力研究

2017年8月8日的九寨沟M_S7.0地震发生在岷江断裂、塔藏断裂及虎牙断裂交汇地区,地处青藏高原东北部的川甘交界地区,位于巴颜喀拉地块的东缘,地质构造复杂,对于九寨沟地震震中位置和发震断层的确定,存在不同意见.本文利用GNSS及升降轨InSAR观测,在获取九寨沟地震同震形变场的基础上,基于均匀弹性半无限位错模型,联合反演了发震断层的滑动分布模型,并计算了同震库仑应力变化.InSAR同震形变场显示,视线向最大沉降量和抬升量分别为0.21 m和0.16 m,形变场长轴为NW向,形变主要集中在断层西侧.距震中40 km和65 km的九寨和松潘两县,水平向的GNSS同震位移分别达14.31 mm和8.22 mm.联合GNSS和InSAR同震形变场反演得到的滑动分布主要集中在沿走向5~33 km,倾向2~20 km的范围内,平均滑动量为0.18 m,最大滑动量为0.91 m.发震断层长40 km,宽30 km,走向155°,倾角81°,滑动角-9.56°.同震位移场及滑移分布模型表明此次地震为一次左旋走滑为主的地震事件,地震破裂并未完全到达地表,与虎牙断裂北段的几何产状和运动学性质更为接近,结合精定位余震的分布,我们确定虎牙断裂北段为此次地震的发震断层,震中位于北纬33.25°,东经103.82°,震源深度10.86 km,矩震量为7.754×10¹⁸ Nm,相应的矩震级为M_W6.5,与美国地调局和哈佛大学给出的震源机制解基本一致.同震库仑应力导致了虎牙断裂北段延长线的东北和西南两端应力增强,其中塔藏断裂的罗叉段和马磨段未来强震的危险性值得关注.     

云南景谷M_S6.6地震对南汀河断裂带地震危险性的影响

南汀河断裂带是滇西南块体内部的一条左旋走滑断裂带,几乎横切整个块体,总长度达380km,是块体内NE向断裂中最长的一条,也是次级块体的边界断裂.南汀河断裂带晚第四纪活动性非常强烈,但仅南段在1941年记录到一次约7级的地震,其余段落还没有5级以上地震的记录,目前可能正处于应力积累和孕震阶段.2014年10月7日发生的景谷MS6.6地震位于断裂带南东约94km,其地震烈度等震线长轴与余震皆呈北西展布,指向南汀河断裂带.为了解景谷地震对周边构造特别是南汀河断裂带的影响,本文通过数值模拟方法计算了地震触发的同震静态库仑应力变化.利用两种同震滑动分布模型计算获得的结果显示,景谷地震对震中附近的断裂,如澜沧江断裂和景谷断裂影响较大,局部应力增加可达90kPa;对较远的断裂,如南汀河断裂带、龙陵—澜沧断裂带和无量山断裂带的影响较小,应力变化值均小于10kPa.通过设置不同断层参数进一步计算,南汀河断裂带北段两支断裂断层面上的静态库仑应力扰动呈半圆形分布,应力增加的最大值位于24.15°N附近的地表,沿断层的走向和深度都逐渐减小.其中西支断裂上应力变化最大值为0.89kPa,东支断裂上为1.18kPa.此外,在南汀河断裂带北段的古地震研究结果显示,该断裂段全新世以来发生过产生地表破裂的大地震,震级应当不低于7级.放射性碳测年将该次古地震事件的发震时间限定在900—1480AD,离逝时间为535—1115年.结合古地震事件的离逝时间和断裂带的滑动速率,本文计算得到南汀河断裂带北段已经积累的水平滑动量为2.8+1.5/-1.0m,进一步利用滑动量与震级的经验公式可估算出该断裂段目前积累的滑动量如果完全释放将会产生一个7.5+0.1/-0.2级的地震.虽然景谷地震在南汀河断裂带上触发的静态库仑应力变化值表明,该地震可能不会引起南汀河断裂带地震危险性的突变,但仍起到一定的加速作用.再考虑到断裂带北段目前已经积累了约7.5级地震所需的能量,该断裂段在未来具有较高的地震危险性.     

基于InSAR同震形变观测反演2010年新西兰南岛Mw7.1 Darfield地震同震破裂分布

2010年9月4日新西兰南岛Canterbury平原发生了M_w7.1地震,震源深度约为10 km.本次地震发生在一条震前不为人所知的断层上.我们利用覆盖整个震区的合成孔径雷达(SAR)观测资料,通过干涉处理分析获得雷达视线向(LOS)同震形变场;以此资料为约束反演了断层的几何参数以及同震破裂分布.结果显示,该地震造成四条相对独立断层的破裂.大部分的地震矩释放发生在Greendale断层(编号1-4),其错动以右旋走滑为主,最大破裂约为8.5 m.其它三条断层中,经过震源的逆冲断层最大破裂为5.1 m (编号6),位于Greendale断层以西的逆冲断层最大破裂为3.5 m (编号5),位于Greendale断层北面的走滑断层最大破裂为1.9 m(编号7).反演的Greendale断层地表滑动与地质调查得到的地表破裂在形态和数值上均吻合较好.本次地震释放的地震矩为5.0×10¹⁹N·m,矩震级为7.1.板块边界带形变场分析表明,Darfield地震的发生受边界带应变分配在该地区残留构造应力场控制,其复杂性体现了区域构造应力场的特点.地震对其周围地区的应力场影响较大,库仑应力增加区与余震分布有一定对应关系,并在2011年Christchurch 6.3级地震发震断层区域造成约0.1bar的库仑应力增加,对此地震有一定的触发作用.     

新疆精河M_W6.3地震产生的静态应力变化研究

利用精河M_W6.3地震有限断层破裂模型,计算了精河地震产生的位移场、应力场、周围主要断层上的静态库仑应力变化以及主震对余震的触发作用。结果表明:(1)精河地震产生的地表隆升最大值约为6.6cm,沉降最大值约为1.8cm;水平位移方向呈现震中南北侧向震中汇聚、震中东西侧向外"流出"的特点。(2)精河地震产生的水平面应力场展布南北侧物质主要受到指向震中的拉张力作用,东西两侧物质主要受到因震中过剩物质东西向排出而导致的东西向挤压力作用。(3)震中西侧距震中约20km的库松木契克山前断裂中段和震中东北部距震中约50km的四棵树-古尔图南断裂西段的库仑应力加载均大于0.01MPa,即2处为地震危险区。(4)在震源深度为8~12km的余震事件中,约有85.5%处于库仑应力加载区,即受到主震的的触发作用;在深度为4~8km的余震事件中,约有87%受到主震的应力触发作用。     

2017年8月8日四川九寨沟7.0级地震InSAR同震形变场及断层滑动分布反演

基于InSAR技术,利用欧空局升降轨Sentinel-1A/IW宽幅数据,获取了2017年8月8日四川九寨沟7.0级地震InSAR同震形变场,并以升降轨InSAR观测结果为约束,反演了断层滑动分布,基于三种不同接收断层计算了同震库仑应力变化.结果表明,同震形变场发生在塔藏断裂、岷江断裂和虎牙断裂交汇的三角地带,升降轨干涉位移均显示本次地震的形变场影响范围约为50 km×50 km,形变场长轴方向为NW向,升降轨观测的形变量相反,反映断层运动性质以走滑运动为主,升降轨数据观测得到的最大LOS （Line of Sight,视线向）形变量分别为~22 cm和~14 cm.非对称形变场反映出断层两侧的运动差异.反演结果显示,最大滑动量约为1 m,平均滑动角为-9°,矩震级为M_W6.5,地震破裂主要集中在地下1~15 km深度范围内,但整体而言本次地震破裂较为充分,基本将该区域1973年及1976年4次 > M_W6.0地震的破裂空区完全破裂.考虑到塔藏断裂和虎牙断裂的运动性质,可初步判定发震断层为虎牙断裂北侧延伸分支.基于三种不同接收断层模型的同震库仑应力变化计算结果反映出该区域以应力释放为主,进一步触发较大走滑型余震的可能性不大.     

帕米尔—兴都库什地区构造应力场反演及拆离板片应力形因子特征研究

从Global CMT目录搜集了1976年1月至2016年6月之间的震源深度大于70km的255个震源机制解,用阻尼应力反演方法,分70～160km和170～310km两个深度,计算了帕米尔—兴都库什地区的构造应力场;同时以10km为间隔计算了兴都库什地区深度介于70～310km之间的应力形因子.得到以下初步结论:兴都库什板片向下俯冲和帕米尔地区断裂带的横向拉张,可能是导致应力场不同的原因.兴都库什俯冲带与帕米尔俯冲带碰撞,导致交汇地区(37°N—37.5°N)的应力场参数突变.兴都库什俯冲板片受到深部温度、压力等因素,出现薄弱面进而形成拆离板片.其脱离了主俯冲板片的束缚后,重力的上下拉张作用导致空区附近张轴倾伏角接近90°,拆离板片俯冲至上地幔不连续面,导致板片部分熔融进而应力形因子随着深度变小.而拆离板片受到地幔挤压其内部发生破碎,其压应力轴由西部的NS到东部NW-SE方向偏转及纵向张应力轴倾伏角变小.     

九寨沟M_s7.0地震的InSAR观测及发震构造分析

2017年8月8日四川省九寨沟县发生M_s7.0地震.本文基于Sentinel-1 SAR影像,利用InSAR技术获取了此次地震的同震形变场,反演获得同震滑动分布,计算了同震位错对余震分布和周边断层的静态库仑应力变化,并对发震构造进行了分析讨论.结果表明:①InSAR同震形变场显示,九寨沟地震造成地表形变最大量级约为20 cm(雷达视线方向),同震形变存在非对称性分布特征.②同震位错以左旋走滑为主,主要发生在4~16 km深度,最大滑动量约为77 cm,位于9 km深处.反演得到的矩震级为Mw6.46.同震错动未破裂到地表.③大部分余震发生在库仑应力增加区.此次地震增加了震中周边地区一些断裂的库仑应力,如东昆仑断裂带东段、龙日坝断裂、虎牙断裂等.④东昆仑断裂东段的未来地震危险性值得关注.⑤九寨沟地震的发震断层为树正断裂,可能是虎牙断裂的北西延伸隐伏部分,此次地震是巴颜喀拉块体南东向运动受到华南块体的强烈阻挡过程中发生的一次典型构造事件.     

STUDY ON DISPLACEMENT OF THE PEAKS OF THE HIMALAYA GENERATED BY THE 2015 NEPAL EARTHQUAKE SEQUENCE

Based on the rupture models of the 2015 Nepal earthquake sequence and half space homogeneous elastic model, the displacement field near the epicenters is estimated. The horizontal components converge to the epicenters from north and south with maximum value of 871~962mm. The farther the epicenter distance is, the smaller of the horizontal displacement occurred. The displacement on the south side of the epicenters decreases more rapidly than that on the north side as the distance from the epicenter increased. Significant settlement occurred on the north side of the epicenters with maximum of 376~474mm, while large uplift occurred on the epicenters and its south side with maximum value of 626~677mm. Then, the displacement of the peaks of the Himalaya near the epicenters is estimated. The largest displacement occurred at the peak of Shishapangma with 393mm horizontal component and 36mm settlement. Mt. Everest, the world's highest peak, moves 36mm in nearly southward direction with 9mm settlement. The displacements of other peaks of the Himalaya are different with the epicentral distance and azimuth of the 2015 Nepal earthquake sequence.     

COULOMB STRESS EVOLUTION AND SEISMIC HAZARD ALONG THE EASTERN BOUNDARY OF THE SICHUAN-YUNNAN BLOCK

Using a more realistic model of multi-layered viscoelastic media, and considering the effects of the coseismic dislocation and the postseismic viscoelastic relaxation caused by the 34 great earthquakes occurring along the eastern boundary of the Sichuan-Yunnan block since 1480 and the interseismic stress accumulation caused by the tectonic loading generated by plate motions which were modeled by introducing "virtual negative displacements" along the major fault segment in the region under study, we calculated the evolution of the Coulomb stress change in each fault plane of 18 major fault segments along the eastern boundary caused by the coseismic, postseismic and interseismic effects. We studied the interactions of the Xianshuihe, Anninghe, Zemuhe and Xiaojiang fault zones on the eastern boundary of the Sichuan-Yunnan block. By evaluating if the previous earthquake could bring another earthquake closer to or farther from failure, we analyzed the interactions of the earthquakes which occurred in the different segments in the same fault zone, or in the different fault zones respectively. And further based on the calculation results of the Coulomb stress change on the fault planes, we analyzed the seismic hazard of each fault segment.The results show that the previous earthquake may trigger another earthquake which can occur in the same fault zone or in the different fault zone. And the calculation results on the evolution of the cumulative Coulomb stress change in the each fault segment show that, the Coulomb stress increases significantly in the middle section and the Moxi segment of the Xianshuihe fault zone, the Mianning-Xichang segment of the Anninghe fault zone, the Qiaojia-dongchuan segment and the Jianshui segment of the Xiaojiang fault zone, and the seismic hazard in these fault segments is worthy paying attention to.     

Coseismic deformation of the 2020 Yutian MW 6.4 earthquake from Sentinel-1A and the slip inversion

A remarkable earthquake struck Yutian, China on June 26th, 2020. Here, we use Sentinel-1 images to investigate the deformation induced by this event. We invert the InSAR observations using a two-step approach: a nonlinear inversion to constrain fault geometries with uniform slip based on the rectangular plane dislocation in an elastic half-space, followed by a linear inversion to retrieve the slip distribution on the fault plane. The results show that the maximum LOS displacement is 22.6 ​cm, and the fault accessed to the ruptured characteristics of normal faults with the minor left-lateral strike-slip component. The fault model indicates a 210° strike. The main rupture zone concentrates in the depth of 5–15 ​km, and the fault slip peaks at 0.89 ​m at the depth of 9 ​km. Then, we calculate the variation of the static Coulomb stress based on the optimal fault model, the results suggest that the Coulomb stress of the Altyn Tagh fault and other neighboring faults has increased and more attention should be paid to possible seismic risks.     

2008年新疆于田7.3级地震对周围断层的影响及其正断层机制的区域构造解释

计算了2008年3月21日新疆于田M_s7.3地震导致周围断层的库仑破裂应力变化. 结果表明,库仑应力在贡嘎错断裂的西南段和东北段、康西瓦断裂东段和西段、普鲁断裂西段、龙木错-邦达错西段、阿尔金断裂和西昆仑山前逆冲断裂增加分别达6.9×10⁴ Pa、4.4×10³ Pa、8.0×10⁵ Pa、2.8×10³ Pa、2.3×10⁴ Pa、5.7×10³ Pa、1.6×10³ Pa和1.1×10³ Pa,而使得贡嘎错断裂中段、康西瓦断裂中段、龙木错-邦达错东段、普鲁断裂中段和喀喇昆仑北段断层上的库伦破裂应力减小分别达3.9×10⁵ Pa、5.8×10⁵ Pa、1.5×10⁴ Pa、1.5×10⁴ Pa和1.1×10³ Pa. 根据康西瓦的断层滑动速率估计该地震使得中段特征地震延迟1~127年,使得东段特征地震提前2~138年. 计算还得到2005年10月8日巴基斯坦M_w7.6地震在于田地震破裂面和滑动方向上产生的库仑破裂应力增加为10~20 Pa,虽然有利于于田地震的发生,但数值太小,难以看到明显的触发作用. 分析了该地震及其周围已发生地震出现正断层机制的原因. 认为在印度板块向欧亚板块北北东向推挤作用下,相对塔里木块体,青藏高原北缘的柴达木块体沿阿尔金断裂东向运动,而兴都库什块体沿喀喇昆仑断裂整体北西向运动. 地震发生区位于此两块体之间,在双边动力相互作用下导致了东西向拉张的形变模式.     

苍山地震库仑破裂应力变化与地震活动相关性探讨

本文收集了苍山地震周围断层的几何参数和运动学参数, 选择不同的接受断层, 研究了苍山地震库仑应力变化对余震及周边断层活动的影响。 结果表明, 苍山地震正的库仑应力变化对余震活动有促进作用, 大部分苍山余震分布于库仑应力变化增加区域, 应力增加量值为0.003~0.005 MPa。 通过对比该地区苍山地震前后16年的地震资料, 发现苍山地震正的库仑应力变化对苍山断裂附近地震活动增强亦有一定的促进作用, 而对于周边其他断层, 由于距离较远, 库仑应力变化量较小, 苍山地震库仑应力变化的影响可以忽略。     

利用差分干涉雷达测量技术(D-InSAR)提取同震形变场h

简要介绍了合成孔径雷达干涉测量技术、差分干涉雷达测量技术，并对干涉测量精度进行了简单讨论．以西藏玛尼地区为例，通过三通差分干涉处理，获取了玛尼地震同震形变场.结果表明:形变场长200 km、宽115 km．干涉条纹以北东东向发震断层——玛尔盖茶卡断层为中心分布，且基本与发震断层平行；通过对干涉形变图进行分析，发震断层可分为3段，其中西段长约23 km，中段长约60 km，东段长约26 km，整个发震断层共长110 km；震中附近最大隆起斜距向位移量为162.4 cm，断层西侧最大沉降斜距向位移量为103.6 cm，震中最大地面水平位错为7.96 m．      

2021年云南漾濞MS6.4地震序列静态应力触发研究

为深度剖析2021年云南漾濞M_S6.4地震对周围断层的影响以及地震序列间的静态库仑应力影响关系,文章基于主震的破裂模型和Okada给出位移对空间偏导数的解析式,首先计算了主震在周围断层上产生的静态库仑破裂应力,结果表明维西—乔后断裂带中段、澜沧江断裂带北端、红河断裂带北端以及怒江断裂带中段库仑应力均有千帕量级的增加。其次,计算了此次地震对周围地区产生的水平应力场及位移场,发现震中东西两侧物质向外流出,南北两侧向震中汇聚;震中南北两侧沉降,东西两侧隆升;产生的应力场呈EW向挤压,NS向拉张,在一定程度上抵消了该区域背景构造应力场。最后计算了前震-主震-余震序列间的静态库仑应力影响,结果表明前震产生的静态库仑应力促进了主震的发生;在2 km、13 km和18.5 km深度附近,触发的余震(前震和主震产生的库仑应力变化为正)比例很高,但在7 km深度处(同震破裂模型中滑移量最大)大部分余震分布在库仑应力负值区(应力影区),考虑到该深度余震与主震震源机制相差较大,因此通过模拟最易错动的断层面作为余震接受断层面,从而计算出最大静态库仑破裂应力,发现应力影区的余震仍有被触发的...     

采用有限元等效体力方法数值分析火山区岩浆压力变形源——以长白山火山为例

火山区岩浆压力变形源的反演计算采用解析方法存在难以考虑地形的限制,采用传统有限元方法则存在网格依赖和计算量大的问题,反演过程中每一次正演由于岩浆房位置和大小变化都需要重新生成一次网格,耗费巨大的计算量和网格生成时间.为了克服上述问题,首次在长白山火山区使用"有限元等效体力"方法考虑地形影响反演地下岩浆压力变形源,计算岩浆应力扰动对周边断层稳定性的影响.在火山区地下压力变形源引起的地表形变计算中,地表地形影响不可忽略.埋深越浅,地表最大径向位移u_r所在的位置越靠近岩浆囊中心.当坡度达到30°时,最大垂向位移u_z所在位置不再位于岩浆囊正上方.椭球状岩浆囊压力源可以较好地模拟长白山火山地区2002—2003年间的GPS和水准测量.岩浆房扰动应力场和区域构造应力场的叠加有可能造成天池西部近EW向,天池北部以NW-NNW向为主的现今应力方向.岩浆房压力源引起的库仑应力变化有利于天池火山口NW向震群在空间上主要分布于火山口的西南和东北部.     

1997年玛尼地震对青藏川滇地区构造块体系统稳定性影响的三维DDA+FEM方法数值模拟

本文用三维非连续变形与有限元相结合（DDA+FEM）的方法,在青藏川滇地区三维构造块体相互制约的大背景中,通过用GPS资料做位移速率边界约束和震源机制约束,计算得到研究区的初始位移场和应力场与该地区GPS测量结果和震源机制分布结果基本一致.在此基础上进一步数值模拟1997年玛尼7.9级大震的发生过程,研究大震引起研究区各块体边界断层应力状态变化的特征.（1）发震断层两侧发生左旋走滑错动,最大水平位错大约7 m;（2）深部位错面上位错分布与用地震波资料震源反演的结果类似;（3）最大差应力变化等值线图与由星载D\|INSAR技术获取的地表形变场图像相似;（4）地表垂直位移表明地震断层面略向北逆冲.计算模拟得到了玛尼地震发生引起青藏川滇地区构造块体系统各边界断层上库仑破裂应力变化的分布,表明玛尼大震的发生除了使其发震断层的两端库仑破裂应力增大,应力进一步集中外,位于上地壳层上东昆仑断裂中段的2001年昆仑山8.1大震（H=11 km）发震断层段的库仑破裂应力增加约2 MPa,位于中地壳层上喀拉昆仑断裂带中的2008年改则6.9级地震（H=30 km）发震断层段的库仑破裂应力也增加约0.7 MPa,可见这两个已接近破裂强度地段的失稳对发生大震起了一定促进作用.研究结果也表明：作者发展的三维DDA+FEM方法能有效地用于大震活动与各构造块体相互作用关系的研究.     

2014年新疆于田MS7.3地震同震位移及位错反演研究

利用于田震中300 km范围内的1个GPS连续站和12个GPS流动站数据,解算得到了2014年新疆于田M_S7.3地震地表同震位移,并反演了发震断层滑动分布,探讨此次地震对周边断裂的影响.地表同震位移结果显示,GPS观测到的同震位移范围在平行发震断裂带的北东-南西向约210 km,垂直发震断裂带的北西-南东方向约为120 km,同震位移量大于10 mm的测站位于震中距约120 km以内;同震位移特征整体表现为北东-南西方向的左旋走滑和北西-南东方向的拉张特征,其中在北东-南西方向,I069测站位移最大,约为32.1 mm,在北西-南东方向,XJYT测站位移最大,约为28.1 mm;位错反演结果表明,最大滑动位于北纬36.05°,东经82.60°,位于深部约16.6 km,最大错动量为2.75 m,反演震级为M_W7.0,同震错动呈椭圆形分布,以左旋走滑为主并具有正倾滑分量,两者最大比值约为2.5：1,同震错动延伸至地表,并向北东方向延伸,总破裂长度约50 km,地表最大错动约1.0 m;同震水平位移场模拟结果显示贡嘎错断裂、康西瓦断裂和普鲁断裂等不同位置主应变特征具有差异性,这种差异特征是否影响断裂带以及周围区域的应力构造特征,值得关注.