2021年青海玛多MW7.3地震同震和早期震后断层滑动模型及发震构造摩擦属性

地震同震和早期震后断层滑动是研究发震断层深浅部孕震形态及摩擦性质的重要手段.本文基于雷达干涉测量(InSAR)数据获得了玛多地震同震和早期震后形变场,并反演得到了发震断层的滑动分布模型.研究结果表明同震滑动主要分布在上地壳浅部,并且存在多个滑动亏损区域.发震断层在东端分叉处的倾向与主断裂和西端倾向相反.基于N-SBAS方法获得的震后4.5个月形变场显示,断层近场震中区域的累积形变量达5 cm,远场区域累积形变达2 cm.震后早期余滑分布在断层浅部和深部以及两侧的横向延展区域;部分余滑区域与同震区域重叠.其中,上地壳浅部的余滑达20 cm且填充了同震滑动亏损.时序地震矩显示上地壳浅部区域在震后快速余滑而深部稳定滑动,表明了发震断层区域的复杂摩擦属性.     

基于GPS观测研究2010年青海玉树MS7.1地震震后地壳形变特征及其机制

2010年4月14日青海玉树M_S7.1地震发生在青藏高原东南部甘孜-玉树地震带,在震后7~10天内,我们快速建立了由15个GPS测站组成的跨地震破裂带观测剖面,包括1个连续站,3个半连续站和11个流动站,对所有站进行了240多天的观测,获取了该次地震的震后形变时空特征.采用欧拉矢量和位错模型解算了背景速度场,并从GPS观测的形变场中扣除该分量.采用分层黏弹性位错模型计算余震引起的地表形变,结果表明余震对部分测站的位移造成不可忽视的影响.采用对数模型拟合位移时间序列,表明特征衰减时间为6.7±1.2天.利用最速下降法反演震后余滑时空分布,反演结果表明震后断层活动以左旋滑动为主,断层南盘具有少量的抬升.在空间分布上,余滑主要位于同震破裂区的两侧,西北侧的余滑几乎达到地表,而东南区的余滑基本在同震破裂区的下方,余滑最大的区域位于结古镇东南下方10~20 km的深度范围.随着震后离逝时间的增加,2个余滑区在空间上保持不变,余滑区的面积逐渐扩大.余滑的矩释放为（1.5~5.1）×10¹⁸Nm,相当于1个M_W6.1~6.4地震释放的能量.分层岩石圈黏弹性模型计算的地壳孔隙弹性反弹形变与地表观测值相差较大,不能解释观测到的震后变形.采用麦克斯维尔流变体模型计算下地壳和上地幔松弛引起的地表形变,显示出其对地表形变的贡献较小.GPS观测得到的震后形变所具有的快速衰减特征,以及余滑模型能够较好地拟合GPS地表形变,表明2010年玉树M_S7.1地震后早期阶段的地壳形变主要是由余滑机制决定的.     

基于D-InSAR技术的玉树Ms7.1地震同震形变场提取与分析

2010年4月14日青海玉树县发生了Ms7.1地震.本文利用差分干涉合成孔径雷达(D-InSAR)测量技术,通过对日本先进陆地观测卫星搭载的ALOS PALSAR获得的雷达数据进行处理,提取了此次地震的同震干涉形变场.结果表明,形变场分布特征与其发震断层甘孜-玉树断裂带的左旋走滑特征一致;最大视线向形变量为-61.4c...     

2017年伊朗MW7.3地震的震后余滑分布及其对2018年MW6.0地震的触发影响

收集了覆盖监测区域的Sentinel-1卫星雷达影像,利用短基线集干涉测量技术提取了2017年伊朗萨波尔扎哈布M_W7.3地震后283天的地表时序形变,通过二步法反演得到其震后余滑分布,之后采用差分干涉技术获取了2018年发生于同一地区的贾万鲁德M_W6.0地震的同震形变场,并将反演所得的发震断层参数作为应力计算的接收断层参数,来分析2017年M_W7.3强震及其震后活动对2018年贾万鲁德M_W6.0地震的触发影响。结果表明:萨波尔扎哈布地震的震后形变主要由孕震断层面的余滑运动所致,震后283天余滑模型的累积滑移量达到0.7 m;2018年贾万鲁德地震的发震断层走向为355.6°,倾角为89.4°,同震断层破裂以右旋走滑为主,兼具部分正断层运动。本文所得的贾万鲁德地震断层平面上的库仑应力变化表明,2017年M_W7.3主震及其震后余滑对2018年M_W6.0地震的发生具有一定的触发效应,M_W6.0地震的发生可能与区域板块的活动性相关。      

2008年10月西藏当雄MW6.3地震震后形变提取与余滑反演

2008年10月6日西藏当雄发生M_W6.3地震.本文利用震后2008年10月26日至2010年8月22日的16期ENVISAT ASAR数据,通过小基线集干涉测量、误差校正与MInTS（Multiscale InSAR Time Series）技术提取高精度的震后形变场,利用SDM（Steepest Descent Method）方法反演断层震后余滑演化过程,并分析震后余滑与同震滑动的关系.结果表明：当雄M_W6.3地震的近场震后形变场主要位于断层西侧,在时间演化上具有明显的对数函数衰减规律;震后余滑主要集中于断层中南段深0~15 km区间,最大的余滑量约0.07 m,位于断层深约9.28 km处,滑动角约-103°;震后余滑引起的地震矩能量M₀与矩震级M_W在时间演化上具有指数函数递增规律;当剪切模量μ=32 GPa,震后665天余滑释放的地震矩能量约为1.92×10¹⁷N·m,约占同震滑动释放地震矩的4.8%,相当于矩震级M_W5.46;虽然震后余滑已经延伸到断层浅部0~5 km区间,但由于余滑量相对较小,没有改变同震滑动在断层浅部区域的滑动亏损现象,这可能是2010年11月30日该区域又发生M_W5.3级余震的主要原因之一.     

基于InSAR分析2015年尼泊尔MW7.8地震形变场特征及断层滑动分布反演

2015年4月25日尼泊尔发生了M_W7.8地震, 本文基于震前、 震后两景Sentinel-1A雷达影像, 采用D-InSAR两轨差分干涉法提取了此次地震的同震形变场。 结果显示, 同震形变场位于喜马拉雅造山带—主边界逆冲断裂(MBT)和主前锋逆冲断裂(MFT)附近, 形变场整体表现为自西北向往东南方向延伸近150 km的纺锤形包络状, 以大面积隆起抬升形变为主, 视线向最大隆升形变达1.18 m, 抬升区北侧存在一小沉陷区, 以InSAR观测值定位同震最大形变中心。 基于均匀介质弹性半空间模型(Okada模型)与InSAR观测数据反演了断层滑动分布。 反演结果表明该地震属于典型逆冲型地震, 发震断层为主喜马拉雅逆冲断裂(MHT), 同震破裂从主喜马拉雅逆冲断裂(MHT)向上沿着主前锋逆冲断裂(MFT)传递。 基于InSAR同震形变场局部形变细节, 结合震区地质背景、 断裂分布及断层运动特征, 获得了同震破裂拟出露地表迹线。     

利用GPS观测资料分析2008年于田Ms7.3地震的同震位移及震后形变

2008年3月21日新疆于田发生Ms7.3级地震.本文通过处理、分析GPS数据,得到破裂断层北侧100 km附近的同震位移及震后形变信息.在观测区域GPS点监测到10 mm左右的同震位移,其中最大为南向14 mm,东向5 mm.同震位移呈现一致性的东南向运动特征,证实于田地震存在显著的左旋走滑分量.震后台站向西南方向运...     

2001年昆仑山口西8.1级地震同震形变场特征的初步分析a

利用差分干涉雷达测量技术获取的宏观震中区的同震形变场，结合对地震活动性、震源机制、野外考察等资料分析，对昆仑山口西8.1级地震同震形变场特征进行了研究. 结果表明:宏观震中位于库赛湖东北侧，宏观震中区发震断层可分为两个形变中心区域，其中西段长约42 km，东段长约48 km，整个发震断层主破裂段长90 km；由干涉形变条纹分布格局可清楚地判断出发震断层的左旋走滑特征；断层两盘变形特征不同，南盘变形程度明显大于北盘；宏观震中附近最大斜距向位移量为288.4 cm，最小斜距向位移量为224.0 cm，宏观震中发震断层最大左旋水平位错为738.1 cm，最小地面左旋水平位错为551.8 cm.      

基于分区解缠方式研究2010年玉树M_W6.9地震同震形变

2010年玉树MW6.9地震发生后,已有学者利用DInSAR技术得到了该次地震的同震形变场,并在此基础上反演了其震源破裂滑动分布。本文以提高玉树地震同震形变场准确度作为出发点,设计了一种新的解缠方式获取相应的同震形变场:首先沿地表同震破裂迹线将InSAR干涉图像分割成上下两部分,然后利用网络流解缠算法对其分别进行相位解缠,最后通过统计重叠部分的相位一致性信息将两部分结果进行拼接。利用震中区GPS同震位移数据对分区解缠结果进行验证。结果表明:相对于整体解缠结果,分区解缠方式得到的最大视线向沉降值由42.6cm增加到48.1cm,所得结果与GPS观测结果更加接近,采用分区解缠方式提高了断层附近形变场的准确性。     

鲁甸地震同震与震后形变、重力时空分布的理论模拟

基于位错理论,考虑重力和黏弹性的影响,在分层介质模型下计算鲁甸地震引起的同震、震后形变和重力变化.结果表明形变和重力的显著变化主要发生于断层在地表投影附近区域.同震形变场显示发震断层有明显的走滑性质.考虑黏弹松弛效应,随着时间的推移,震后形变和重力有了明显改变,同震效应为正的区域得到加强,为负区域进一步减弱.震后松弛效应的影响范围相比同震明显增加.在靠近断层的GPS观测台站处,计算了由黏弹松弛效应引起的震后形变和重力时间序列.震后松弛效应引起的重力变化在50年之后均达到同震水平,除了NJ13的纬向、垂直位移,NJ16的垂直位移,NJ15的径向位移,其余台站的所有震后形变都超过1mm.观测台站的震后重力和垂直位移时间序列在震后100年趋于稳定,纬向位移和经向位移在震后50年趋于稳定.     

Numerical Simulation of Displacement and Stress Fields Associated with the 1993 Kushiro-oki, Japan, Earthquake

—We constructed a three-dimensional finite element model to simulate coseismic and postseismic displacement and stress fields associated with the 1993 Kushiro-oki earthquake, which was a very large intermediate-depth earthquake that occurred within the subducted Pacific plate at a depth of 107 km beneath the southeastern part of Hokkaido, Japan. Taking the configuration of the subducted Pacific plate into account, we constructed a realistic model with lateral heterogeneity of viscoelastic structure. We assigned a variable slip distribution to the fault plane, which was obtained from inversion analysis of near-field seismic waveforms. The result shows that elastic deformation associated with the faulting reflects the assigned inhomogeneous slip distribution on the fault plane near the fault region, while it does not reflect the distribution on the free surface of the model. The calculated postseismic deformation does not reflect the slip distribution, but shows symmetric spatial patterns concerning the dipping direction of the fault both near the fault region and on the model surface. For the next 20 years following the earthquake, the amount of the calculated deformation is a fraction of the coseismic deformation. The calculated coeseismic deformation is large just above and below the fault plane, reaching 1 m, while the postseismic deformation is dominant near the upper and lower material boundaries between the subducted plate and the surrounding asthenosphere. The spatial distribution of maximum shear stress near the fault plane corresponds to the assigned slip distribution, amounting to 32 MPa. The directions of principal stress-change axes represent reverse fault type in the SSE region of the fault, whereas normal fault type is dominant in the NNW region with the exception of some asymmetrical spatial patterns of the principal stress-change axes on the fault due to the inhomogeneous slip distribution. Time variations both in the amount and the directions of stresses are minor, suggesting that the coseismic state of the stress would remain unchanged for two decades after the event.     

基于Sentinel-1A数据的四川长宁MS6.0地震同震形变场分析及断层滑动分布反演

利用大地测量数据研究2019年6月17日四川长宁M_S6.0地震同震形变场特征和发震断层参数, 基于DInSAR技术处理升降轨Sentinel-1A数据获取干涉相位图, 并考虑大气折射效应和余震形变误差实现同震形变场改正。四叉树采样后的形变数据作为反演数据源, 采用弹性半空间分层模型反演发震断层几何面滑动分布。结果表明本次地震发震机制为兼具逆冲和左旋走滑, 矩震级为M_W5.9, 断层破裂尺度达28 km×20 km, 震源深度约9.4 km。升降轨视线向同震形变场在断层两侧呈现形变特征差异, 最大沉降量分别是8.34 cm(升轨)和4.23 cm(降轨), 最大抬升量分别是5.5 cm(升轨)和7.5 cm(降轨); 发震断层走向为302°, 倾角为43°, 平均滑动角为50°, 断层面最大滑动量达到0.28 m。     

1976年唐山地震震时和震后变形的模拟

本文采用三维粘弹性有限元方法拟合唐山地区１９７６-１９８５年观测到的地震震时和震后的水平与垂直地形变，反演华北板块下方深部物质的流变学性质．模型采用多层弹性覆盖层与线性粘弹性层的有限块体，发展断层面上存在着位错运动，并用正交设计法拟合观测数据．模拟计算表明，华北板块下方软流层粘度为７．１×１０^１８Ｐａ·ｓ；上地幔粘度为２．１×１０^１９Ｐａ·ｓ．     

九寨沟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.同震错动未破裂到地表.③大部分余震发生在库仑应力增加区.此次地震增加了震中周边地区一些断裂的库仑应力,如东昆仑断裂带东段、龙日坝断裂、虎牙断裂等.④东昆仑断裂东段的未来地震危险性值得关注.⑤九寨沟地震的发震断层为树正断裂,可能是虎牙断裂的北西延伸隐伏部分,此次地震是巴颜喀拉块体南东向运动受到华南块体的强烈阻挡过程中发生的一次典型构造事件.     

Coseismic vertical deformation of the MS8.0 Wenchuan earthquake from repeated levelings and its constraint on listric fault geometry

The devastating M_S8.0 Wenchuan earthquake ruptured two large parallel thrust faults along the middle segment of the Longmenshan thrust belt. Preseismic and postseismic leveling data indicated the hanging wall of the YingxiuBeichuan-Nanba thrust fault mainly presented coseismic uplift with respect to the reference point at Pingwu county town, and the observed maximum uplift of 4.7 m is located at Beichuan county (Qushan town) which is about 100 m west of the fault scarp. The foot wall of the Yingxiu-Beichuan-Nanba thrust fault mainly showed subsidence with maximum subsidence of 0.6 m near the rupture. By employing a listric dislocation model, we found that the fault geometry model of exponential dip angle δ= 88°?×1-exp(-9/h) with depth of 18 km and uniform thrust-slip of 5.6 m could fit the observed coseismic vertical deformation very well, which verifies the listric thrust model of the Longmenshan orogenic zone.

     

2021年青海玛多Mw7.3地震地壳变形及区域孕震环境

北京时间2021年5月22日02时04分,在青海省果洛藏族自治州玛多县发生Mw7.3地震.本文利用震中及邻区2009-2021年GNSS观测资料,研究了此次地震的区域孕震环境、同震和震后初期的变形特征.首先,给出了较高空间分辨率的GNSS速度场,表明震中及其邻区的构造活动以左旋剪切为主,巴颜喀拉块体南北两侧的相对运动是...     

基于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,与美国地调局和哈佛大学给出的震源机制解基本一致.同震库仑应力导致了虎牙断裂北段延长线的东北和西南两端应力增强,其中塔藏断裂的罗叉段和马磨段未来强震的危险性值得关注.     

InSAR数据约束的2021年5月21日云南漾濞MS6.4地震发震构造研究

运用Sentinel-1A卫星数据和D-InSAR技术,获取2021-05-21云南漾濞M_S6.4地震的同震形变场。结果显示,漾濞地震同震形变场长轴近NW展布升降轨形变场符号相反,视线向最大沉降量和抬升量为0.1 m。InSAR同震形变场反演的滑动分布主要集中在沿走向2～12 km,倾向1～9 km的范围内,最大滑动量0.35 m,发震断层长9.8 km、宽4 km,滑动量主要集中在地下3～6 km范围内,滑动角-146.7°。同震位移场及滑动分布模型反映本次地震为发震断层的右旋走滑事件,地震破裂未达到地表。断层模型反演结果显示,矩震级为M_W6.1,发震断层以北西走向右旋走滑运动为主,初步认为本次M_W6.1地震发震断裂可能是一条NW向的维西—乔后断裂西侧的隐伏次生断裂。