沿着横断板块边界地震滑动和无震滑动之间的相互作用：地震序列的一…

本文提出了或许能解释在同一孕震区经常观测到的地震多次发生的模型。这种地震序列或震群可能包含在空间上相距几十公里、在时间上相隔几天至数月的一些事件。本模型考虑了一条埋在弹性半空间内的长垂直走滑断层。断层的强度和一滑动类型（地震或无震）是非均匀的，并受到缓慢增加的周围剪切应力的作用。我们假设，断层面的一部分存在有稳定无震滑动的条件，并伴有非稳定的凹凸体。假设产生该序列中一次地震的凹凸体破裂，引起了远离     

1946年南海地震的同震滑动分布及该地震引起的无震滑动

由海啸波形的反演估算了1946年南海地震(Mw8．3)断层面上的同震滑动分布。比以前的研究(Satake，1993)有以下3个方面的进展：(1)使用了大量的较小次级断层；(2)次级断层更拟合板块的几何形状；(3)应用了更详尽的探测数据。反演的结果表明，实测波形和合成波形间的一致性比以前的研究有极大改善。在四国近海震源区西半部，在闭锁区下倾端附近发生了约6m的大滑动。上倾或浅部的滑动很小，说明该地区为弱的地震耦合。在纪伊半岛近海震源区东半部的约3m的大滑动延伸到整个闭锁区。由大地测量数据(Sagiyaand Thatcher，1999)估计了上板块八字形断层上的大滑动，对于解释海啸波形并不需要这样大的滑动，表明这些大的滑动是无震滑动。板块界面下倾端上的两个滑动分布，一个是由大地测量数据得到的，而另一个是由海啸波形得到的，除了在四国室户角下的滑动外，两个滑动分布都相当一致。据此来看，无震滑动也发生在室户角下部的板块交界处。     

用DDA模拟华北地区的无震断层滑动和断块形变

华北地区的地质构造具有典型的断块系统特征。用非连续形变分析（DDA）方法模拟了该地区的长期形变，包括无震断层滑动和断块形变。通过模拟我们发现：（1）断层沿不同方向的滑动通常与ENE－WSW向的区域构造挤压一致。（2）如果沿主要活动断层的平均无震滑动速度与跨断层大地测量观测的结果一样是每年零点几毫米的量级，那么断块内的典型应变速率是10^-8/a或更小量级，所以一些小区域的观测报告所得的应变率10^-6/a就不能代表该地区的形变速率。（3）对一条断层，由相邻断块内刚体旋转引起的断层滑动方向可能与区域构造挤压引起的滑动方向相反。但由于由构造挤压引起的滑动量级远大于由断块旋转引起的断层滑动，因此，一般来讲，断层滑动图像总体上与这一地区构造挤压方向一致，也就是说，由区域挤压引起的滑动控制了整个滑动情况。（4）根据（3），某些观测到的与ENE构造挤压方向相反的滑动可能是由某些更为局部化的因素引起的，没有构造意义。     

地震前无震形变速率变化的证据

在无岩浆环境下发生的10次地震中,已公布的可靠数据表明,震前形变率变化所持续的时间能够从几百秒钟到长达十多年。尽管大部分M>7.5级的地震是在没有探测到震前形变的情况下发生的,但无震形变的探测阈值仍然很高,地震矩相当于M5级地震的无震滑动在大部分情况下(尽管不是全部)会被忽视。震间形变率产生变化时并不伴随地震,而且1989年加利福尼亚洛马普列塔M6.9级地震前出现的应变率变化在尺度上同自那时起所发生的其他速率变化相近似。大部分持续几百秒钟或者更长时间的震前无震形变的例子也许起源于地震破裂面附近或者下方,而不是起源于正在发生地震破坏的区域内部。     

尼泊尔M_W7.8地震InSAR同震形变场及断层滑动分布

采用DInSAR技术和欧空局2014年新发射的Sentinel-1A/IW数据,获取了2015年4月25日尼泊尔M_W7.8地震的InSAR同震形变场.所用InSAR数据扫描范围东西长约500 km,南北宽约250 km,覆盖了整个变形区域,揭示了形变场的全貌及其空间连续变化形态.此次地震造成的地表形变场总体呈现为中部宽两端窄的纺锤形,从震中向东偏南约20°方向延伸,主要形变区东西长约160 km,南北宽约110 km,由规模较大的南部隆升区和规模较小的北部沉降区组成,南部最大LOS向隆升量达1.1 m,北部最大LOS向沉降量约在0.55 m.在隆升和沉降区之间干涉纹图连续变化,没有出现由于形变梯度过大或地表破裂而导致的失相干现象,表明地震断层未破裂到地表.基于InSAR形变场和部分GPS观测数据,利用弹性半空间低倾角单一断层面模型进行了滑动分布单独反演和联合反演,三种反演结果均显示出一个明显的位于主震震中以东的滑动分布集中区,向外围衰减很快,主要滑动发生于地下7~23 km的深度范围内.InSAR单独反演的破裂范围,特别是东西向破裂长度大于GPS单独反演的破裂长度,而InSAR单独反演的最大滑动量则低于GPS单独反演的滑动量.因此认为联合反演结果更为可靠.联合反演的破裂面长约150 km,沿断层倾向宽约70 km,最大滑移量达到4.39 m,矩震级为M_W7.84,与之前用地震波数据和GPS数据反演的结果一致.     

2021年云南漾濞MS6.4地震同震地表形变与断层滑动分布

2021年5月21日云南漾濞县城西侧发生MS6.4地震.文中以欧洲航空局升、降轨Senti-nel-1 SAR为数据源,基于D-InSAR技术获取了此次地震的InSAR同震形变场:升轨LOS(Line of Sight)向最大形变量约为0.07m,降轨最大LOS向形变量约为0.08m.以升、降轨同震形变场数据和GNSS...     

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年以上.     

2022年1月8日青海门源MS6.9地震GNSS同震形变场及其断层滑动分布

2022年1月8日在青海省海北州门源县发生M_S6.9地震,本次地震是继2016年门源M_S6.4地震后冷龙岭断裂周边发生的又一次强震。确定本次地震的破裂分布对分析该地区震害风险具有重要意义。通过收集震中及周边12个GNSS连续站点和震后加密观测的17个流动站点观测资料,获取了震中100 km范围内29个测站的GNSS静态同震形变场,并以此为约束反演了本次地震同震滑动分布。结果显示,近场GNSS观测到的最大形变量可达1.3 m。反演的最优破裂模型显示该地震主破裂区深度位于0～10 km,滑动破裂出露地表,最大滑动量为4.07 m,地震矩释放能量约1.1×10¹⁹ N·m,对应矩震级M_W6.7。门源地震破裂至地表是造成该地区基础设施破坏的直接原因。     

利用GPS和水准测量资料反演2001年昆仑山口西8.1级地震的同震滑动分布

通过分析 2 0 0 1年昆仑山口西MS8.1地震前后的GPS资料和 1 979— 2 0 0 2年的水准测量资料 ,获得了地表同震位移场。利用这些同震位移数据 ,以震后详细野外地质调查破裂数据为约束 ,反演了破裂断层面上的同震滑动分布。结果表明 ,破裂下界深度为 1 4 .2～ 2 1km(70 %置信区间 ) ,最佳破裂深度 1 7km。虽然在太阳湖段和主破裂带西端的中间区域未观测到地表破裂 ,但反演的结果表明此区域存在 2～ 3m左右的左旋水平走滑 ,与InSAR资料分析得到的结果一致。地表以下的破裂西端止于太阳湖段 ,而东端似乎在地表破裂迹线以东 30km范围内仍存在 1 .5～ 2 .0m的左旋滑动。反演的垂直位移表明断层在东经 93°以西部分大体表现为南盘上升 ,而东经 93°以东部分表现为北盘上升。由大地测量和地表破裂调查得到的地震矩释放为 6 .1× 1 0 2 0 N·m ,与地震波资料的反演结果基本一致     

2022年泸定6.8级地震GNSS同震形变场及其约束反演的破裂滑动分布

2022年9月5日四川省甘孜州泸定县发生6.8级地震,地震发生在川滇菱形块体东边界鲜水河断裂磨西段西侧附近.本文基于震中350 km范围内中国大陆构造环境网络和中国地震科学实验场118个GNSS连续站数据,观测获得了精细的同震形变场,结果显示：同震形变跨发震断层呈空间对称分布,表明本次地震具有显著的左旋走滑特征;记录到的最大同震形变发生在震中距40 km的SYD5(石棉安顺场)站,东西向和南北向形变量分别达到-22.0±1.2 mm和11.6±0.9 mm,震中距100 km以外测站的水平同震形变均小于5 mm;垂向同震形变不显著.结合走向163°和倾角77°的发震断层模型,本文对发震断层面的同震破裂滑动分布进行了反演,结果显示：同震滑动主要集中在主震东南侧余震空区0~15 km深度范围内,且破裂达到了地表,同震释放的地震矩与矩震级为M_W6.57地震相当.结合理论同震形变场、主应变场和邻近区域主要活动断裂库仑应力变化的空间分布特征,本次可能导致处于强闭锁和地震空区的安宁河断裂石棉—冕宁段未来发震风险性增强.

     

基于Sentinel-1A的新疆阿克陶MS6.7地震同震形变与滑动分布特征

文中利用Sentinel-1A升、降轨数据获取了2016年新疆阿克陶地震的同震形变场。形变场以水平运动为主,形变主要发生于断层南盘;升、降轨最大形变量分别约为12cm和-21cm;基于雷达影像观测右视成像的特点可知阿克陶地震具有右旋走滑的破裂特征,结合形变场形态特征与余震剖面推断认为,阿克陶地震的发震断层为S倾的木吉断裂。基于均匀弹性半空间模型反演双断层面的静态滑动分布,结果显示:2个断层面上各存在1个椭圆状的滑动集中区,破裂基本到达地表,滑动主要发生在沿断层面走向长约50km、沿断层面倾向方向0～20km的范围内,最大滑动量位于约9km深处,量级约为0.7m;西段主要以走滑为主,东段为走滑兼具少量正断性质,反演得到的地震矩约为8.81×10¹⁸N·m,相当于M_W6.57。综合分析认为,本次地震的震源特征为右旋走滑兼有少量正断分量,是发生在公格尔拉张系内拉张环境下的构造地震事件。     

基于GNSS和InSAR约束的2021年玛多MS7.4地震同震滑动分布及应用

2021年5月22日青海省果洛藏族自治州玛多县发生M_S7.4地震,震中位于青藏高原中部的巴颜喀拉块体,这是近20多年来在巴颜喀拉块体周边发生8次M≥7级强震后,块体内部的一次强震,也是汶川地震以来中国大陆发生的最大一次地震,因此该地震的成因及周边地区未来的地震危险性值得重点关注.本文利用震后及时获取的39个近场流动GNSS观测,联合61个GNSS连续观测、Sentinel-1和ALOS-2 InSAR观测获取了本次地震精细的同震形变场,以此为约束,基于均匀弹性半无限位错模型,反演了发震断层的滑动分布,并计算了同震库仑应力变化.GNSS水平同震形变十分显著,断层南北两侧的GNSS点位,最大水平形变分别达0.7 m和-1.2 m,距震中200 km的测点仍有1 cm左右的同震形变.Sentinel-1和ALOS-2的升降轨InSAR同震形变场显示此次地震造成了约160 km长的地表破裂,最大视线向形变分别达0.9 m和1.2 m.同震滑动分布模型显示,发震断层由主段和次段组成,长度分别为170 km和20 km,主段倾向北,倾角85°,平均滑动角为-4.36°,表明玛多地震是一次典型的走滑型地震.次段倾向南,倾角68°,平均滑动角为-11.84°.地震破裂主要集中在0~15 km深度范围,最大滑动量为4.4 m,对应深度6.97 km.反演给出的矩震量为1.61×10²⁰N·m,对应矩震级M_W7.4.主发震断层上存在4个凹凸体,玛多地震是一次不对称双侧破裂事件.结合余震精定位、野外调查及地质资料,我们认为主发震断裂为昆仑山口—江错断裂,东部的次级破裂与主破裂机制不同.同震库仑应力结果显示,东昆仑断裂玛沁段应力有所增加(>0.01 MPa),处于应力加载状态,未来发生强震的危险性较高.     

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地震的破裂空区完全破裂.考虑到塔藏断裂和虎牙断裂的运动性质,可初步判定发震断层为虎牙断裂北侧延伸分支.基于三种不同接收断层模型的同震库仑应力变化计算结果反映出该区域以应力释放为主,进一步触发较大走滑型余震的可能性不大.

     

两个自由度的滑块——弹簧模型中的幕式无震滑动

设想了一种两个自由度的滑块——弹簧模型，其中的2个滑块(滑块1和滑块2)由一根弹簧连接，并由一个慢慢移动的传动装置驱动。假设采用与速率和状态相关的摩擦定律，设定摩擦参数，使滑块1发生动态失稳。当摩擦参数接近稳定性转换点时，滑块2就会发生幕式无震滑动。在地震间隔期，当应力累积到大约稳态水平时，滑块2就开始缓慢滑动。随着振幅波动的衰减而成为稳态值，应力和滑动速度会发生准静态振荡。振荡衰减可能是自2001年以来日本东海地区观测到的幕式无震滑动的合理的产生机制。我们还在不同摩擦性质相互作用的双滑块模型的基础上，设想了板块边界上的各种滑动模式。     

2020年伽师MS6.4地震InSAR同震形变与滑动分布特征分析

利用Sentinel-1A升轨和降轨数据,基于D-InSAR技术,获取2020年1月19日伽师M_S6.4地震同震形变场,并结合其他研究机构给出的震源机制解参数和已有研究成果,反演得到伽师地震的发震断层几何特征和滑动分布。研究结果表明,伽师地震同震形变在地表有明显差异;升轨同震形变在卫星视线方向北侧抬升55 mm,南侧下降42 mm;降轨同震形变在卫星视线方面北侧抬升63 mm,南侧下降23 mm。通过反演得到发震断层走向为275°,倾角为20°,地震滑动主要分布在地下5 km处,最大滑动量约为0.32 m,平均滑动角为89.3°,累积地震矩为1.46×1018 N·m,合矩震级M_W6.1,发震构造为具有少量走滑性质的逆冲断裂。从发震构造特征、同震滑动分布推测,伽师地震发震构造是柯坪塔格褶皱带滑脱面以上沉积盖层内的逆冲断裂,支持了柯坪推覆体的薄皮构造模型观点。     

基于InSAR技术的伽师MS 6.4地震同震形变与滑动反演研究

<正>1研究背景地震震级较大或震源较浅时,地表会产生一定程度的变形,通常称之为同震形变。及时、准确和全面地获取同震形变场,能够快速认识地震地表形变整体宏观特征,甚至有助于研判发震断层位置及其走向变化等,这对政府应急响应和灾后救援部署等至关重要,具有重大的社会效益和经济效益（温少妍等,2020）。合成孔径雷达干涉测量技术（InSAR）是在干涉雷达基础上发展起来的一种微波遥感技术,具有形变灵敏度高、高空间分辨率高、覆盖范围大、几乎不受云雨天气制约等突出的技术优势,在对同震形变等大范围缓慢形变的测量研究中迅速得到广泛应用。     

利用GPS和InSAR数据反演2011年日本东北MW9.0地震断层的同震滑动分布

本文首先对Envisat/ASAR数据进行干涉处理, 获取2011年日本东北M_W9.0地震的地表InSAR同震形变场; 然后通过对InSAR同震形变数据重采样方法的深入分析, 选择条纹率法结合干涉图的空间相干性对InSAR同震形变数据进行重采样; 最后基于弹性半空间位错模型, 联合InSAR与GPS形变数据, 采用最小二乘法反演发震断层的滑动分布. 研究结果表明: 考虑相干性的条纹率重采样方法, 更适用于形变场中存在除断层外的有限边界、 且形变场范围较大的InSAR数据重采样处理; 断层滑动主要发生在地表以下50 km范围内, 最大滑动量为49.9 m, 矩张量为4.89×1022 N·m, 所对应的矩震级为M_W9.1, 与地震学反演的结果比较吻合.      

2022年青海门源MW6.9地震同震形变及断层滑动分布反演

2022年1月8日,青海海北州门源县发生M_W 6.9地震,震中位于青藏高原东北缘祁连地震带上冷龙岭断裂和托莱山断裂的交汇处.门源地震活动强烈,造成地表破裂明显,因此研究门源地震的发震机理,对评估周边区域及主要断裂的地震危险性具有重要意义.本文基于D-InSAR技术,利用升降轨Sentinel-1A SAR数据,获取门源地震的同震形变场.结果显示,同震形变主要集中在冷龙岭断裂和托莱山断裂的交汇处,形变长轴整体呈NWW-SEE,同震引发的升轨隆升形变量0.40 m,沉降量0.65 m,降轨隆升形变量0.80 m,沉降量0.70 m,升降轨视线向(LOS)形变呈现符号相反大小相近的特征,断裂运动以左旋走滑为主.以升降轨同震形变信息为约束条件,基于弹性半空间位错理论,采用两步法进行反演,获取门源地震断层的几何参数.结合该区域的地质构造得到断层面上的精细滑动分布,结果表明,断层破裂延伸至地表,破裂迹线长度达22 km.断层滑动主要集中在地下2～12 km,最大滑动量为4.2 m,位于地下7 km处.断层走向约为109°,倾角约为82°,释放的地震矩为2.67×10^...     

2020年西藏尼玛MS6.6地震InSAR同震形变场及断层滑动分布反演

2020年7月23日,西藏尼玛县发生MS 6.6地震,研究本地区的发震断层对认识周边活动断裂有着重要意义.由于发震地区位于青藏高原腹地,地震监测台站分布稀疏,监测能力较弱.故本文基于升、降轨Sentinel-1A SAR数据,通过差分干涉技术获取了西藏尼玛MS 6.6地震InSAR同震形变场,显示本次地震所引起的地表形...     

基拉韦厄火山南翼的突发无震断层滑动

与海底火山相关的最大灾害之一不在于火山本身，而在于潜在的火山侧翼的塌陷（McGuire,1996;Moore,etal,1989)。这种侧翼塌陷可以酿成毁灭性的海啸。不仅会波及邻近地区，而且对整个大洋盆地沿岸城市都会造成巨大威胁（Ward,2001)。位于美国夏威夷岛的基拉韦厄火山是一座具有潜在测翼塌陷危险的火山（Ward,2001;Lipman,etal,1988)，因此该火山一直被大地测量台网密切监测，该台网包括全球定位系统（GPS）的接收仪，倾斜仪和应变仪，我们发现，2000年11月上旬，该台网记录到短瞬的东南向位移，我们将其解释为发生了一次无震断层滑动，这一滑动持续了约36h，等同于一次矩震级为5.7的地震，并且最大滑动速度为6cm/d，对GPS数据的反演表明，4.5km深处存在一个小倾角逆冲断层，我们将它解释为希利纳巴利-霍莱伊巴利正断层系的下倾部分，这表明连续记录的大地测量台网可以监测到加速滑动，也许可以对火山侧翼塌陷作出预警。