全球定位系统测定的尼泊尔M_w7.8级地震同震位移

据尼泊尔境内的连续GPS观测资料显示,2015年4月25日的尼泊尔M_w7.8级地震造成尼泊尔向南移动,最大位移量达到了1.88 m.利用国家重大科技基础设施项目"中国大陆构造环境监测网络"、中国气象局中日合作JICA项目及加州理工学院在尼泊尔境内所建的连续GPS观测资料,分析了此次地震对中国大陆构造形变场的震时动态响应及同震影响.结果显示,地震造成中国西藏地区毫米至厘米级的同震水平位移,最大值在震中北东方向,震中距220 km的珠峰站,约为30 mm,波及范围远至1500 km之外的高原东北部的祁连山地区.1 Hz的高频GPS数据分析表明,震中距1200 km内的站点记录了弹性震波和永久变形的叠加和破裂方向效应等丰富信息.通过分析中国境内地震活动频度并对比2011年日本宫城地震发现,本次尼泊尔地震并未造成中国境内微震活动显著增加,同震过程对中国西藏部分地区、川滇地区及南北地震带断裂带的应力扰动和构造加载有限,与其震级相对较小有关.     

GPS测定的尼泊尔M_w7.9和M_w7.3级地震同震形变场

2015年尼泊尔M w7.9和M w7.3级地震致灾范围包括尼泊尔、印度北部、巴基斯坦、孟加拉和中国藏南地区,地震应变应力调整将对震区和邻区的地震活动产生不同程度的影响.尼泊尔境内的GPS连续观测数据、"中国大陆构造环境监测网络"GPS基准站和地震应急流动站观测数据计算结果揭示了尼泊尔地震的静态和动态同震形变场.对M w7.9级地震,尼泊尔境内近场静态同震GPS水平位移最大为1.89 m,距震中100~400 km的中国藏南地区观测到几毫米到几十厘米的静态水平位移,毗邻尼泊尔的聂拉木县最大形变54.0 cm.M w7.3级地震静态同震形变最大为2 cm,局限在震中附近200 km范围内.距震中约2000 km内的GPS基准站均记录到M w7.9大震明显的动态形变信号,高频GPS动态形变幅度与地震破裂方向有关,位于破裂方向上的测站动态形变幅度明显大于其他方向的测站.用弹性半空间位错模型正演模拟了震区和青藏高原南部格网点上的同震形变,并分析了地震应变影响,认为尼泊尔地震对中国藏南地区产生一定程度的拉张型应力变化,需要持续关注.     

全球定位系统(GPS)测定的印尼苏门达腊巨震的远场同震地表位移

以中国的连续GPS观测资料为基础, 结合其他全球GPS观测资料, 计算了2004年12月和2005年3月发生在苏门达腊岛9.3和8.7级巨大地震的同震位移场. 2004年地震的同震永久位移的影响范围达6000~7000 km, 并呈现波状起伏的特征. 2005年地震的同震位移主要局限在震中区附近, 可能只是2004年地震诱发的强余震.     

GPS观测的2015年尼泊尔M_S 8.1级地震震前应变积累及同震变形特征

2015年4月25日,在印度板块与欧亚板块交界区的喜马拉雅地震带上发生了尼泊尔MS 8.1级大地震.震前GPS速度场和应变率场显示,喜马拉雅地震带整体表现为15.94±1.82 mm/a的压缩特征,同时还具有分段活动特征.此次地震发生在速度场顺时针旋转和逆时针旋转的分界带,该处最大主压应变率的量值在喜马拉雅地震带并非最大.GPS观测的同震位移场揭示了尼泊尔MS 8.1级地震引起的地壳变形特征,分别有9和6个测站观测到明显的水平向和垂向同震位移,其水平分量的运动方向整体表现为南向运动,位于震中东南侧的3个测站垂直分量表现为上升,其余测站为下降.中国境内距离震中最近的5个测站的垂向同震位移显示,此次地震造成珠穆朗玛峰的沉降量约为4mm.依据GPS观测到的同震位移场,利用非负最小二乘方法反演震源断层面上的滑动分布.反演结果表明最大滑动量为6.84 m,滑动量较大的区域分布在加德满都附近及其以北区域的下方,这可能是造成加德满都地区具有较大破坏的原因之一,该滑动分布模型能够很好地解释GPS观测到的同震位移.利用此滑动分布模型计算的地震矩为8.21×1020 N m,对应的矩震级为MW7.9.     

2011 年3月11日日本宫城Mw 9.0 级地震对其周边地区火山活动的影响

2011年3月11日日本宫城Mw 9.0 级大震后的两天九州岛上新燃岳火山发生大规模喷发. 此次火山喷发是否是受地震触发所致, 以及地震是否同样也会造成富士山火山和长白山火山的喷发成为人们关注的热点. 通过对地震造成的区域应力应变场的计算, 发现富士山火山、新燃岳火山和长白山火山均处于体膨胀区, 在地下10 km 处体膨胀分别达到~220, ~8 和~14 nano-strain. 由同震GPS 位移场计算得到的3 个火山地区地表面应变也表现为膨胀, 符合模型计算结果. 考虑到拉张应力可能对岩浆通道产生的扩张作用、对包容在岩浆内二氧化碳气体的分异作用以及由此引发的一系列正反馈过程, 此次地震可能造成这些火山活动的增加, 应当警惕触发影响可能造成的火山喷发.     

利用GPS 位移反演日本Mw 9.0 仙台地震及Mw 7.9 强余震静态位错模型

基于GPS 观测的M_w 9.0 日本仙台地震同震形变及震后8 h 的地表形变, 考虑分层介质模型的影响, 反演了主震及震后8 h 的断层的位错分布. 结果显示, 主震位错对应的矩震级为M_w 8.98, 最大同震位错量可达23.3 m, 在震中区显示出纯逆冲特性, 震中两侧断层错动显示出一定的走滑分量. 另外, 约90%的地震能量发生在40 km 深度以上. 而由震后8 h 的地表形变反演的结果显示, 这一期间释放能量约相当于一个M_w 8.13 的地震, 位错主要分布在主震破裂的西南区域, 最大位错量约1.5 m, 破裂峰值区与M_w 7.9 强余震有明显的对应关系, 暗示震后8 h的断层错动主要是由M_w 7.9 强余震引起的. 另外, 主震破裂的深部延伸面在震后0.2~0.4 m 的滑动可能主要是无震余滑.     

张北-尚义地震同震形变场雷达差分干涉测量

用合成孔径雷达差分干涉测量技术获取1998年1月10日张北-尚义地震(Ms = 6.2级)的同震形变图, 利用地震前后欧洲遥感卫星1和2号合成孔径雷达(ERS-1/2 SAR)的三景数据做重复轨道差分处理得到震前震后的干涉纹图. 由视向形变图可知, 此次地震造成隆起形变, 形变中心位于114°20′E, 40°57′N, 最大视向位移量达25 cm, 形变集中在300 km2范围内, 根据地震同震形变场的空间分布特征分析了震源机理和孕震构造.     

2015年智利8.3级地震构造背景及对俯冲带地震的启示

2015年发生在智利Coquimbo的M_w8.3级地震是一个典型的俯冲带地震,其大小、破裂滑移分布与利用震间期GPS观测得到的智利南部俯冲带闭锁模型一致.利用反投影方法得到的破裂过程初步解显示,这次地震的破裂由震中(约25 km深)向浅部传播,并且这次地震破裂的能量辐射与频率明显相关,与2010年发生在Coquimbo地震破裂南部约50 km的M_w8.8级地震的观测结果一致,表明这部分俯冲带的断层性质随深度存在着系统变化.这次2015年8.3级地震的震前及同震观测资料将对地震破裂动力学的数值模拟提供一个难得的模型约束信息.     

中国大陆地震震源深度及其构造含义

对中国大陆地震的震源深度资料作了系统研究.震源深度的研究,对于探索地震孕育和发生的深部环境,地震能量集结、释放的活动构造背景,以及地壳内部构造变形及其力学属性等都有非常重要的意义.以1970年1月～2000年5月期间中国大陆ML≥2.0级并给出深度的31282次1和2类精度的浅源地震为基础,研究了深度分布特征,并采用网格滑动平均方法,统计了网格内地震的平均深度.结果表明,中国大陆平均深度为（16±7）km,东部地区为（13±6）km,西部为（18±8）km,东部比西部平均偏浅5km.我国大陆地震深度最大的地区是新疆西南部地区,即塔里木地块的西端和西南缘.震源深度与构造分区的关系密切,青藏活动地块震源深度平均为（33±12）km.新疆活动地块为（21±10）km,华北为（14±7）km;东北为（11±5）km;华南为（10±4）km.完整地块边界上的地震较深,其中最明显的是塔里木盆地的西南缘、北缘;准噶尔盆地的南边缘,阿拉善地块的南边缘;鄂尔多斯地块的东、西两侧以及四川盆地的西边缘.在新生性的破裂带上地震偏浅,如滇西南地震带及张渤地震带.另外,还根据我国震源深度分布特征,讨论了壳幔（主要是地壳）力学行为、变形属性和破坏方式.     

2008年汶川8.0级特大地震孕育和发生的多单元组合模式

2008年5月12日四川省汶川县境内发生8.0级特大地震. 这次逆冲型地震发生在大陆内部的高角度逆冲断裂之上, 与有历史记载以来所发生的逆冲型特大地震是不同的. 通过对汶川地震的地表破裂、震源机制、余震定位、地震破裂过程、同震地壳形变、强地面运动等的综合研究, 认为汶川特大地震的孕育和发生是3个地质单元共同作用的结果. 川西高原作为变形单元震前发生长期持续的变形, 并且将变形转换为积累在龙门山断裂带的应力; 龙门山断裂带作为闭锁单元震前变形缓慢但积累很大的应力, 当其超过断裂的摩擦强度或岩体的破裂强度时就突发破裂, 形成地震, 释放出巨大的能量; 四川盆地作为支撑单元对川西高原和龙门山的向东运动产生阻挡, 是汶川地震孕育不可缺少的元素. 汶川地震的孕育和发生可以用多单元组合模型来理解.     

基于水准数据的芦山7.0级地震震间期和同震位移场特征

利用1985~2010年跨龙门山断裂带南段的精密水准观测,并消除2008年汶川地震同震变形的影响,获取了龙门山断裂带南段震间期的地壳垂直形变场.结果表明,双石-大川断裂(龙门山前山断裂)和新开店断裂(大邑断裂南段)之间的芦山地块处于正、负垂直形变梯度带相交的区域,说明该区域在芦山地震前25年内已闭锁.根据2010~2013年穿过震区的水准观测数据,同时消除震前近3年的地壳垂直变形,得到了芦山地震同震垂向位移.同震位移场揭示距离震中最近的达多35水准点最大同震位移达到198.4 mm(相对马雅165基),西北盘(上盘)相对东南盘(下盘)向上运动,显示发震断层以逆冲错动为主.对比芦山地震垂向同震变形和震间期变形,可发现该地震符合弹性回跳学说,芦山地块震间期积累的弹性应变能(位移亏损)在震时得到弹性回跳释放.     

2022年泸定MS6.8地震同震形变特征及周边强震危险性

2022年四川泸定6.8级地震发生于鲜水河断裂磨西段,为揭示此次地震的同震形变特征,收集了震中200 km范围内的连续全球定位系统(global positioning system, GPS)测站和震中50 km范围内的强震动数据,进行了高精度处理以提取测站的同震水平位移.此外,还收集了覆盖震中区域的Sentinel和ALOS-2降轨数据,通过干涉差分技术处理并获得了卫星视线向同震形变场.结果显示:水平同震形变呈四象限分布,表明泸定地震的震源机制为左旋走滑;距震中16 km的强震动台站记录到了明显的位移波形,其峰值位移达14 cm、永久变形约为12 cm;合成孔径雷达干涉测量(interferometric synthetic aperture radar, InSAR)显示断层近场30 km×30 km区域具有明显的同震变形,断层西侧卫星视线向最大位移达15 cm,挖角乡附近的局部形变超过15 cm;野外调查推测在“二台子-爱国村”间的断层段发育地表破裂,但从InSAR形变场上难以判定,同震是否破裂到地表仍需详细的野外考察予以确定.反演了川滇块体北东边界主干断裂在泸定地震前的闭锁分...     

高频GPS测定的汶川Ms8.0级地震震时近场地表变形过程

以龙门山断裂带下盘的四川GPS连续观测网高频(1 Hz)GPS观测资料为基础, 得出了2008年5月12日发生在四川汶川的Ms8.0级大地震的震时近场地表形变过程. 结果显示: 地震在近场产生的最大形变量明显大于震后位移, 破裂带北段各站水平分量震时先向震中方向运动, 后转折垂直于破裂带方向运动, 南段各站水平分量形变相对较小且基本为可恢复性变形. 各站垂向均先下沉, 然后呈周期性上下起伏波动. 将高频GPS所得位移与强震仪所得结果进行对比分析, 发现震动初期两者具有较好的一致性, 但中后期虽然相位基本同步, 但其振幅存在10 cm左右的差异, 其具体原因还需进一步研究论证. 此次记录到的汶川大地震近场地表变形过程, 可为进一步研究地表破裂过程和地震波的传播方式提供非常有价值的基础资料.     

汶川MS8地震最大地表同震垂直位移量及其地表变形样式

汶川地震的同震变形量及其分布是揭示汶川地震孕育机制、破裂扩展特征的重要科学依据. 已有研究显示汶川地震致使龙门山中央断裂中-北段和前山断裂中段分别产生长度为240和90 km的地表破裂带. 在对北川县沙坝村一带地质地貌调查和广泛走访老乡, 获得震前地貌、建筑物信息和照片实证的基础上, 对主要的地物标志、被断错建筑物位差等进行全站仪和差分GPS实测, 得到沙坝村邹家院一带最大垂直位移为(9±0.5) m, 右旋位移量约为(2±0.5) m. 尽管沙坝村一带地表表现为正断层, 但与其他地段的同震变形一致, 表现为断层西北盘上升, 不存在一般正断层上盘重力下滑迹象; 这些说明, (9±0.5) m的位移值应为汶川5.12地震地表最大同震垂直位移量.     

汶川M_S8地震最大地表同震垂直位移量及其地表变形样式

汶川地震的同震变形量及其分布是揭示汶川地震孕育机制、破裂扩展特征的重要科学依据.已有研究显示汶川地震致使龙门山中央断裂中-北段和前山断裂中段分别产生长度为240和90km的地表破裂带.在对北川县沙坝村一带地质地貌调查和广泛走访老乡,获得震前地貌、建筑物信息和照片实证的基础上,对主要的地物标志、被断错建筑物位差等进行全站仪和差分GPS实测,得到沙坝村邹家院一带最大垂直位移为(9±0.5)m,右旋位移量约为(2±0.5)m.尽管沙坝村一带地表表现为正断层,但与其他地段的同震变形一致,表现为断层西北盘上升,不存在一般正断层上盘重力下滑迹象;这些说明,(9±0.5)m的位移值应为汶川5.12地震地表最大同震垂直位移量.     

汶川地震破裂近场震后变形观测及其意义

汶川地震破裂带上的近场震后变形, 是一种震后的蠕滑行为, 而且大部分表现为与同震滑移的方向相反. 地震断层的活动通常划分成震前、同震、震后和震间4个阶段, 反映了一次地震从孕育、发生到结束的整个演化过程, 而不同活动阶段的变形特征反映出不同的应力状态和力学性质. 汶川地震沿龙门山断裂带形成了两条长分别为250和72 km的地表破裂带. 为了了解汶川地震地表破裂带的震后变形特征, 对中央主破裂带上人为破坏较轻的断层崖或断层挠曲崖进行反复测量, 结果显示在19个观测点中, 13个观测点(68%)的断层崖或断层挠曲崖的高度(垂直同震位移)震后降低回落, 平均降低了9.7%; 5个观测点(26%)没有发生变化; 1个观测点(6%)的在震后继续抬升, 抬升了12.8%, 而且该观测点位于中央主破裂带的南西端部. 尽管这种变化中存在着上冲断层盘虚假抬升后压实回落的影响, 但主要是沿汶川地震破裂带发生的震后滑移造成的, 而且大部分震后滑移(68%)与同震滑移的方向相反. 除破裂端部存在同震位移亏损, 弹性能释放不完全外, 其他部位同震位移要么与震间累积达到平衡, 要么过冲产生能量亏损, 揭示了汶川地震的能量可能基本释放完全, 发生7级以上强余震的可能性不大. 此外, 这种震后变形特征还告诉我们, 在进行活动断层构造地貌研究时, 特别是通过断层崖高度(或其他水平位错量)判断断层运动速率和估计古地震事件大小时, 除侵蚀作用产生的误差外, 还需要考虑10%左右来自震后滑移的系统误差.     

从壳幔地震各向异性初探华北地区壳幔耦合关系

根据首都圈地震台网、国家地震台网、IRIS台站资料, 并比较了流动的华北地震台阵部分台站资料的结果, 利用地震各向异性分析, 得到减少局部构造影响后的华北地区(至少是华北地区北部)的地壳背景水平主压应力方向为北东95.1°±15.4°, 穿透壳幔的SKS分裂的快剪切波偏振平均值为北东110.2°±15.8°. 而台站分析也显示, 地壳内近场资料的快剪切波偏振与SKS快剪切波偏振总是相差十几度以上. 不同周期下的Rayleigh面波相速度方位各向异性快波优势方向也同样揭示出随深度增加快波优势方向的变化. 这个特征暗示华北地区的壳幔耦合关系可能既不是简单的壳幔解耦型, 也不是壳幔强耦合型, 可能是两种模式共存的不均匀分布, 或者是物理性质介于两者之间的渐变模式. 研究表明, 在分析穿透整个岩石圈的SKS分裂和地表的GPS观测的基础上, 增加揭示地壳各向异性的近场资料剪切波分裂, 结合面波的方位各向异性, 可以更准确地多方位刻画壳幔各向异性性质, 进而探讨壳幔的耦合关系.     

“是否存在有助于预报的地震前兆”的讨论

经过几十年地震预报探索实践,已发现了很多地震前兆异常,其中部分前兆还通过了优秀地震前兆标准的检验,并分析总结了前兆异常的群体特征.但用这些前兆异常,通过经验性统计方法进行地震预报成功率很低,其根本原因是地震孕育、发生的空间尺度较大,前兆异常的挑选较为困难;地震复发周期较长,已积累的震例资料不足于研究地震孕育发生的规律;地震的孕震环境复杂,与周边的动力、构造环境密切相关,不同地震之间前兆异常差异较大;前兆观测分布在地球表面,只是地下应力、应变的间接反映,所获取的异常与地震前兆关系有待深入研究.在目前情况下,下一次地震完全按统计得到的共性特征演化明显是不可能的.地震前兆观测应从单点观测向直接观测地下应力、应变场过渡,地震预报从经验性统计预报向基于地球物理场动态观测的物理预报过渡.     

汶川8级地震的形成机制与中国大陆近期地震概况

此文概略介绍了四川汶川8级地震的主要特征和形成机制,介绍了余震的情况,讨论了关于地震预报问题。笔者认为,以历史地震资料为依据的中长期预报是比较可信的,临震预报的难度较大,目前尚未获得满意的成果。根据近百年来中国地震活动周期性的变化规律,探讨了未来几十年内我国大陆强震活动的可能分布特点。根据上千年的历史地震资料,讨论了华北地区未来相当于唐山地震的强震活动可能将在240年以后、300年以前的期间发生。     

中国大陆现今地壳运动速度场的最新观测结果

“中国地壳运动观测网络”是国家“九五”期间投资建设的一项重大科学工程,其中GPS观测由27个连续观测的基准站、55个年复测的基本站和近1000个不定期复测的区域站组成[1,2].自1998年该网络开始观测以来,基准站一直连续观测,基本站已观测7期,区域站于1999,2001和2004年观测了3期.中国大陆地域广阔,地质构造复杂,地震活动强烈,故GPS观测点密度越高越能准确地反映中国大陆现今地壳运动的总貌和细部特征,因此,目前区域站的观测数据是刻画现今地壳运动速度场之关键所在.据1999和2001年两期GPS观测,前人已给出过中国大陆速度场的初步结果[3~5…