2016年射阳MS 4.4地震前后重力场变化

利用江苏重力测网2014-2017年重力场观测资料,采用绝对重力控制与相对重力联测相结合的平差方法,获取2016年射阳M_S 4.4地震前后重力场变化图。根据射阳M_S 4.4地震前后射阳地区各测线重力段差变化特点,绘制重力场等值线并进行对比分析,结合相关机理,探讨重力场变化与该地震的内在联系。同震观测数据显示：射阳M_S 4.4地震发生在重力异常值高梯度带附近,发震时震中地区位于NS挤压正异常、EW张拉负异常状态,震后区域重力梯度变化量开始减小,是一种典型的重力异常调整现象。     

全球及中国大陆周边强震活动状态研究

本文应用中国地震台网给出的全球大震目录, 结合前人的研究, 分析了1900年以来全球M_S≥8.0巨大地震活动的时间变化与M_S≥7.8大地震的空间迁移特征, 同时还探讨了包含中国大陆西部及其周边的“大三角”地区M_S≥7.8地震的发展态势。 20世纪全球百年的大震活动表现为20世纪前半叶活动强而后半叶活动弱的前高后低的活动特点, 而进入21世纪以来有增强活动显示, 呈现百年左右起伏的长期变化。 在这百年尺度的起伏之中还含有20年左右的韵律变化, 主要显示为地中海—喜马拉雅地震带巨大地震的平静—活跃交替活动特征。 根据大震活动的起伏特征和地中海—喜马拉雅地震带幕式活动特性, 进入21世纪以来, 地中海—喜马拉雅带巨震频繁发生, 显示其地震活动显著增强的高活动态势。 通过计算中国大陆西部及周边地区的不同震级地震的变异系数C_v发现, 7.0级以上强震与7.8级以上大地震活动类型不同; 采用对数正态分布对M_S≥7.8地震活动进程进行拟合分析发现, 到2015年的累积发震概率达0.80, 到2020年的累积发震概率达0.89, 表明与中国大陆关系紧密的“大三角区”已进入大地震活跃时段。     

威远MS 5.4地震序列活动特征与成因分析

通过对威远M_S 5.4地震背景与活动特征研究,表明自贡地震监测中心加强流动监测提高了威远M_S 5.4地震序列地震参数的精度和地震信息报送的准确性。威远M_S 5.4地震前经历了区域地震平静至活跃、再发震的过程,威远M_S 5.4地震为一般主—余型地震,其ρ、b值异常可能反映了其后将发生资中M_S 5.2地震。威远M_S 5.4地震序列及其附近的华蓥山断裂地震带和马边地震带震后时间扫描地震学参数均正常,威远M_S 5.4地震区近期再次发生更大地震的可能性较小;并认为威远M_S 5.4发震断裂为荣县—威远基底断裂,华蓥山断裂带及其附近荣县—威远—资中地区2018年开始的4—6级地震活跃与川西地区东昆仑断裂带2017年8月8日九寨沟M_S 7.0地震的发生存在呼应关系。     

GRACE卫星观测到的与汶川Ms8.0地震有关的重力变化

利用GRACE卫星重力资料,计算了中国大陆及周边的卫星重力时变场和地表密度变化分布,获取了具有代表性的点位区域的每月重力变化时间序列.同时获得了WUSH、LHAS、KUNM、LUZH站相对于区域参考框架的GPS位移时间序列.卫星重力观测结果显示喜马拉雅弧形带的重力在2004年苏门答腊M_w9.3地震后快速下降, 2006~2008年尤为明显,西域地块西北边界带上震后重力下降也较为显著;而沿青藏高原北至东边界2007年出现明显的重力上升沿构造边界的弧形分布,且2008年南北地震带中南段重力上升变化显著.这些苏门答腊地震后的重力变化趋势到汶川地震发生后才开始改变.GPS位移结果显示四个台站均记录到苏门答腊大地震的同震信号,震后WUSH、LHAS、KUNM站水平位移向量出现明显的运动趋势改变,且一直持续到2008年汶川M_s8.0地震的发生.GRACE卫星揭示的青藏高原及周边地表质量的变化为解释汶川地震的动力机制提供了新的观测途径和资料.本文结合区域构造运动的特点和GPS位移,对GRACE观测的时变重力场特征及汶川地震的动力机制进行了初步解释和讨论.     

垂直摆倾斜仪地震效应与倾斜场模拟

在整理垂直摆地震响应的基础上,采用USGS提供的有限断层模型和速度结构模型,运用德国地学研究中心（GFZ）王荣江研究员开发的（PSGRN/PSCMP）计算程序,分析日本M_W 9.0地震和中俄交界M_S 6.6地震倾斜场变化,并与垂直摆记录数据进行对比。通过整理现有地质资料和断层数据,模拟依舒断裂带北段鹤岗周边发生M_S 6.0地震时,该地区同震及震后倾斜场的变化,以了解不同震源性质的地震产生的地倾斜变化及对鹤岗地区产生的影响。     

利用同震和震后位移数据联合反演2011年日本MW9.0地震同震断层滑动

基于黏弹性球体位错理论, 联合陆地和海底同震GPS数据以及日本本岛330个陆地GPS站点5~10年的震后数据, 反演了日本M_W9.0地震的断层滑动模型, 提升了断层滑动分布在细节上的合理性。 首先, 基于日本本岛330个陆地GPS站点震前2年和震后10年的连续观测数据, 获取了日本M_W9.0地震震后5~10年的年平均位移, 该时段的位移几乎完全由地幔黏弹性松弛效应引起; 接着, 利用黏弹性球体位错理论对震后5~10年的位移进行反复拟合, 确定了日本M_W9.0地震震源及周边地区的地幔黏滞性系数最优解(9.0×10¹⁸ Pa·s)。 然后, 联合同震和震后位移数据, 引入黏弹性位错格林函数, 反演了2011年日本M_W9.0地震的断层滑动分布。 结果表明, 该地震同震破裂的最大值达到了62.72 m, 同震滑动的总地震矩为4.48×10²² Nm, 相应的矩震级为M_W9.03。 由于黏弹性松弛效应引起的震后位移中包含了同震破裂的信息, 基于黏弹性球体地震位错理论, 联合同震和震后位移数据反演断层同震破裂, 有效提高了日本M_W9.0地震断层滑动分布的可靠性。 最后, 本文提出的反演方法为同震观测结果缺乏的大地震震后科考提供了理论支撑: 在大地震发生之后, 即使在同震期间没有足够的观测数据, 也可以在震后通过对震源区的加密观测积累的震后数据, 使用本文提出的反演方法优化同震断层滑动模型。     

2016年3月2日苏门答腊MS7.8地震同震位移和应力场数值模拟研究

苏门答腊地区发育多条左旋走滑性质断层,地震活动活跃.2006年3月2日该区西南海域发生了M_S7.8大地震.大地震的发生常常会引发区域位移场和应力场及周围断层应力状态发生变化.本文建立全球PREM有限元地球模型,据已有的断层滑动模型计算了此次苏门答腊地震引发的同震位移和应力及库仑应力变化,并进一步讨论了此次地震对周围断层的影响以及区域构造应力场对库仑应力变化计算的影响.初步结果表明此次苏门答腊M_S7.8地震造成较大的南北向水平位移且集中在探测者破裂区（Investigator Fracture Zone）,最大水平位移量约6.74 m,断层倾角接近垂直,下盘向北运动而上盘向南,进一步表明M_S7.8地震为典型的左旋走滑为主的地震,发生海啸的可能性较低;库仑应力变化达MPa量级的区域集中在震中,但近场大部分余震分布在库仑应力减小区域,有效摩擦系数变化和区域构造应力场的耦合作用可能是其原因;利用改进的库仑应力变化计算方法和最优破裂方向计算得出的结果显示库仑应力触发理论可较好地解释余震分布.     

2015年尼泊尔MS8.1大地震孕育的深层过程与发生的动力学响应

2015年4月25日发生在尼泊尔博克拉M_S8.1大地震的深层动力过程与“地中海-喜马拉雅-南亚地震带”的中段喜马拉雅地震活动带密切相关.这次大地震是该地震强烈活动带上长期以来深部物质与能量强烈交换、运动,并导致构造活动和应力积累的产物.综合分析与研究提出：（1）博克拉M_S8.1大地震的孕育、发生和发展具有长期活动和近年来相邻地带地震活动频繁的背景;这一地带自1505年-2015年,即500多年来相继发生多次M_S≥8.0的大地震.（2）这一地带具有特异的地球物理边界场响应和深层动力过程,显示深部物质的重新分异、调整与能量交换.（3）大地震发生与周边地带应力场分布特异,壳、幔结构与介质属性变异及破裂响应与断层面解的属性相关.（4）喜马拉雅地带的三条北倾断裂带以不同角度向深部延伸、震源位置及浅表层的变形特征尚应深化理解.（5）M_S8.1大地震的发生对相邻地带的波场影响强烈,故应强化高精度地球物理场的观测和探测,以“捕捉”未来可能大地震的孕育与发生.     

2017年伊拉克MW7.3地震的类型界定及其震后趋势分析

2017年伊拉克地震发生在我们划分的巴格达地震区,鉴于不同机构提供的该震震级参数相差较大,本文利用孕震断层多锁固段脆性破裂理论,分情况讨论了该震所属地震类型,并分析了巴格达地震区地震趋势.结果表明：若2017年伊拉克地震为M_W7.3,则该震为第3锁固段向峰值强度点演化过程中发生的1次显著前震,该区未来将发生M_W7.7~8.2（双震型为M_W7.5~8.0）标志性地震,目前已接近临界状态;若2017年伊拉克地震为M_W7.5,除可能为显著前震外,还可能为标志性地震（双震）之一,若如此两年内该区将发生另一次M_W7.5地震;若2017年伊拉克地震为M_S7.8,则该震为第3锁固段峰值强度点对应的标志性地震,与我们对该震的前瞻性中长期预测结果相符.我们判断该震不为主震,预测该区未来还将发生M_W7.8~8.3（双震型为M_W7.6~8.1）标志性地震,目前该区远离临界状态.     

2015年尼泊尔地震三维发震构造及地震危险性研究

2015年4月25日尼泊尔M_W7.8地震发生在喜马拉雅造山带上,人们普遍认为该地震不足以释放该造山带上累积的能量,但对该地区后续地震危险性评价多基于二维或是假三维的形变反演计算结果.本研究从2015年尼泊尔地震主震与其最大余震M_W7.3地震之间的关系出发,着重分析讨论了尼泊尔地震的时间和空间的非均匀性,结合震源机制解、地壳速度结构、精定位后的余震分布及InSAR反演结果,建立了三维发震构造模型,利用非线性摩擦有限元方法,对一个地震周期内断层摩擦行为和块体变形进行了模拟,将计算结果和地表同震形变、形变反演的同震破裂、历史地震时空演化进行对比,在确认该三维模型可靠性的基础上,讨论了该地区后续地震的可能位置,认为在1934年比哈—尼泊尔M_W~8.1地震以东区域,还存在发生大地震的可能,在最大余震M_W7.3地震东南部位,还存在发生中大地震的可能.     

2010智利莫尔M_W8.8地震对其破裂区及附近地区地震活动性的影响

强震的发生与地壳中应力增加有关。库仑破裂应力被用于研究由于某次强震发生所产生的应力变化。2010年智利莫尔M_W8.8地震发生之后,其破裂区以南地区没有发生M≥6.0地震,而在其以北地区于2015年在智利发生了伊亚佩尔M_W8.3级地震。在这两个地区地震活动性悬殊。本文利用库仑破裂应力研究了2010年智利莫尔M_W8.8地震在这两个地区产生的应力变化。结果表明,2010年智利莫尔M_W8.8地震在其破裂区以北地区产生的应力变化为正值,有利于断层错动发震,而在其破裂区以南地区产生的应力变化为负值,具有抑制断层错动发震的作用。     

九寨沟地震发震断层属性及青藏高原东南缘现今应变状态讨论

九寨沟地震(M_s7.0或M_w6.5)震中位于青藏高原巴颜喀拉块体东缘东昆仑断裂带东端塔藏断裂、岷江断裂和虎牙断裂交汇部位,中国地震局相关科研机构的研究人员曾将该震中区判定为玛沁—玛曲高震级地震危险区.地震应急科学考察期间没有发现地震地表破裂带,但地震烈度等震线长轴方位、极震区基岩崩塌和滑坡集中带、重新定位余震空间展布和震源机制解等显示出发震断层为NNW向虎牙断裂北段,左旋走滑性质,属东昆仑断裂带东端分支断层之一.此外,汶川地震后,在青藏高原东缘和东南缘次级活动断层上发生了包括2017年九寨沟地震(Mw6.5)、2014年鲁甸(M_w6.2)、景谷(M_w6.2)、康定(M_w6.0)等多次中强地震,显示出青藏高原东缘至东南缘各块体主干边界活动断层现今处于中等偏高的应变积累状态,即在巴颜喀拉、川滇等块体主干边界活动断层上具备了发生高震级(M_w≥7.0)地震的构造应力-应变条件,未来发生高震级地震的危险性不容忽视.     

川滇地区重力场动态变化及其强震危险含义

基于川滇地区2011—2014年的重力复测资料,系统分析了区域重力场时-空动态变化及其与2012年云南彝良MS5.7、2013年四川芦山MS7.0、2014年云南鲁甸MS6.5和四川康定MS6.3地震发生的关系.结合GPS、水准观测成果和区域地质构造动力环境,进一步研究了区域重力场变化的时空分布特征及其机理,讨论了近期区域重力场动态变化的强震危险含义.结果表明:1重力变化与川滇地区断裂构造活动存在密切空间联系,重力变化较好地反映了伴随活动断层的物质迁移和构造变形引起的地表重力变化效应.2重力资料对测区内2012年以来发生的4次MS5.7以上强震均有较好反映,地震前震中区及其附近观测到明显的区域性重力异常及重力变化高梯度带,可能是地震孕育过程中观测到的重力前兆信息.3区域重力场动态演化大体反映了青藏高原物质东流的动态效应,龙门山断裂带地壳受挤压隆起、面压缩率和重力上升变化的特征最为显著.4重力场的空间分布及其随时间变化与地壳垂直与水平运动及地质构造活动等观测结果有一定的对应关系,强震易发生在重力变化四象限分布中心地带或正、负异常区过渡的高梯度带上,研究区的一些重力异常部位仍存在中-长期大震危险背景.     

2017年四川九寨沟MS7.0地震前区域重力场变化

利用南北地震带2014-2017年期间的流动重力观测资料,系统分析了区域重力场变化及其与2017年8月8日四川九寨沟7.0级地震发生的关系.结果表明：①区域重力场异常变化与北西西向塔藏断裂和南北向岷江断裂带在空间上关系密切,反映了沿控震断裂在2013-2017年期间发生了引起地表重力变化效应的地壳变形和构造活动.②九寨沟地震前,测区内出现了大空间范围的区域性重力异常,而震源区附近产生了局部重力异常,沿塔藏断裂带形成了重力变化高梯度带,其中,甘肃玛曲、迭部、青海河南蒙古族自治县、四川若尔盖、九寨沟一带重力差异变化达100×10^-8m·s^-2以上;这些可能反映九寨沟地震前,区域及震源区附近均产生与该地震孕育、发生有关的构造运动或应力增强作用.③九寨沟地震震中位于重力差异运动剧烈的鞍部等值线附近,与断裂走向基本一致的重力变化高梯度带零值线上,这一观测事实进一步佐证了重力场动态变化图像对强震地点预测具有重要的指示意义.     

2016年青海门源MS6.4地震前重力变化

利用青藏高原东北缘2011-2015年期间的流动重力观测资料,系统分析了区域重力场变化及其与2016年1月21日青海门源M_S6.4地震发生的关系,结合GNSS、水准观测成果和区域地质构造动力环境,进一步研究了区域重力场变化的时空分布特征及其机理.结果表明：（1）测区内重力场异常变化与祁连山断裂带在空间上关系密切,反映沿祁连山断裂带（段）在2011-2015年期间发生了引起地表重力变化效应的构造活动或变形.（2）门源M_S6.4地震前,测区内先出现了较大空间范围的区域性重力异常,到临近发震前显示出相对闭锁的现象,且围绕震中区周围出现四象限分布特征的局部重力变化,地震发生在重力反向变化过程中,并出现显著的四象限分布特征的重力异常变化,其中,青海门源与甘肃天祝一带重力差异变化达100×10^-8m·s^-2以上.（3）区域重力场动态演化大体反映了青藏高原东北缘物质北东流的动态效应,门源震中附近区域地壳受挤压变形显著、面压缩率和重力剧烈变化的特征最为显著.（4）重力场的空间分布及其随时间变化与地壳垂直与水平运动及地质构造活动等观测结果有一定的对应关系,强震易发生在重力变化四象限分布中心地带或正、负异常区过渡的高梯度带上.     

2020年5月6日新疆乌恰5.0级和5月9日柯坪5.2级地震总结

2020年5月6日、5月9日,新疆地区南天山西段先后发生乌恰5.0级和柯坪5.2级地震,系统总结2次地震发生前出现的地震活动和地球物理观测异常,其中：①地震活动：震前存在调制地震集中、地震窗、5级以上地震成组等中短期异常;②地球物理观测：2次地震震中附近震前出现形变、电磁和流体观测异常,其中形变异常3项、电磁异常4项、流体异常1项,主要分布在柯坪5.2级地震震中附近。通过对2次地震序列进行跟踪,发现：乌恰5.0级地震余震较少,震后60天内共记录M_L 3.0以上余震4次,最大震级为M_L 4.5;柯坪5.2级地震后余震较丰富,震后60天内共记录M_L 3.0以上余震10次,最大震级为M_L 4.7,计算得到序列h值为1.6,b值为0.73。综合分析认为,2020年5月新疆地区2次5级以上地震前存在的地震活动异常较少,但区域地震活动水平较强,主要存在具有中短期指示意义的地球物理观测异常。     

地震大小的度量

本文首先对地震大小度量的研究历程作了简要回顾。在此基础上着重指出以下几点: ① 近几十年广泛使用的M_L、m_b(m_B)、M_S震级标度,作为对地震大小的度量,不仅都存在着“以偏概全”和“震级饱和”的问题,而且由于未充分顾及地震波衰减的区域差异,尤其台站场地对地震动放大作用的差异,因此都只是对地震大小的不精确的度量。 ② 地震矩M₀不仅物理含义明确,而且克服了M_L、m_b(m_B)、M_S震级标度所存在的各种问题,是最适合科学度量地震大小的物理量。为了继续延用“震级”这一术语,Kanamori定义了矩震级标度M_W,虽然其假定的前提条件仍有待研究,但也是对地震相对大小的较合理的度量。 ③ 为更加科学地度量地震的大小,必须充分利用现代区域数字地震台网产出的大量波形数据,加强地震波衰减特征、场地效应和震源参数测定及有关定标关系的研究,在完善M_L、m_b(m_B)、M_S震级测定方法的基础上,重点改进、完善M_W标度,逐渐推进以M_W作为统一度量地震相对大小的物理量,为地震科学研究和地震预测研究奠定更加坚实的基础。     

利用GRACE重力和重力梯度变化估计2012年苏门答腊地震断层参数

2012年4月11日苏门答腊北部附近海域发生M_W=8.6地震,国际上采用不同数据获得了该地震的断层滑动模型,但断层走滑性质存在较大差异.卫星重力GRACE（Gravity Recovery and Climate Experiment）观测覆盖震中区域,可以提供很好的断层参数估计约束.本文采用GRACE月重力场模型数据提取了此次地震同震重力和重力梯度变化,表明了北向分量的优越性.并与断层模型理论计算结果进行了比较,分析了T_xx（北-北向重力梯度）分量对断层参数的敏感性,以及基于粒子群算法以及Okubo位错模型反演了该地震断层参数.结果表明GRACE观测到的同震重力梯度变化在空间形态分布上与断层模型模拟结果符合较好,但振幅差异较大.重力和重力梯度的北向分量可以很好地压制条带误差,其中g_N（北向重力变化）和T_xx的误差最小,其次是T_xz（北-上向重力梯度）,误差最大和对条带误差抑制效果最差的是T_xy（北-东向重力梯度）.T_xx对断层的深度不敏感,对其余断层参数较敏感.基于GRACE反演得到的断层的走向角为113.63°,倾角为89.99°,滑移角为175.26°,平均滑移量为28.18 m,相应的矩震为8.71级,且此次地震的走滑性质为右旋走滑.