2021年青海玛多MW7.4地震GNSS同震形变场及其断层滑动分布

2021年5月22日青海玛多M_W7.4地震作为发生在巴颜喀拉块体内部的一次强震,再次引起了人们对该地区地震活动性的强烈关注.本文基于震后GNSS流动观测和区域连续GNSS站资料,解算了106个站点的同震形变及其中17个站点的高频形变波形.同震形变场显示玛多地震具有典型的左旋走滑特征,GNSS观测到的最大同震位移达到1.2 m.GNSS与InSAR数据相符度较高,GNSS提供了准确的近场形变信息.基于GNSS同震形变场,本文反演了断层滑动分布,并计算了发震断层上产生的库仑应力变化.结果表明,发震断层的滑动破裂存在多个凹凸体,破裂分段特征明显且出露地表,与野外地表破裂考察和余震分布吻合,主体破裂位于断层面0~10 km的浅部区域,最大滑动量达到4.6 m,地震矩1.63×10²⁰N·m,矩震级为M_W7.4;发震断层上静态库仑应力增加区域与余震分布具有一致性,说明余震主要是由静态库仑应力加载而触发的.     

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地震相当.结合理论同震形变场、主应变场和邻近区域主要活动断裂库仑应力变化的空间分布特征,本次可能导致处于强闭锁和地震空区的安宁河断裂石棉—冕宁段未来发震风险性增强.     

GNSS约束的2022年泸定M6.8地震滑动分布及同震应力变化

2022年9月5日四川泸定县发生M6.8级强震,震中位于鲜水河断裂东南末端,野外地质调查初步结果显示本次地震并未发现明显的地表破裂迹象.本文基于震后科考GNSS流动观测和震中周边连续站观测资料,解算并提取震中90 km范围内31个测站的静态水平向同震位移.结果显示：GNSS同震形变场空间分布呈现明显的左旋走滑特征,观测到的最大水平向同震位移达23 cm,震中距50 km范围内同震位移量普遍大于1 cm.基于GNSS观测资料反演得到的同震滑动分布显示泸定地震地表破裂主要集中在磨西至田湾之间,主破裂深度2～8 km,最大滑动量～1.96 m,地震矩9.25×10¹⁸N·m,对应矩震级M_W6.6.静态库仑应力结果显示本次地震增强了震源破裂区周边活动断层的库仑应力,并触发了大量余震,余震主要发生在库仑应力增强区域.结合震间闭锁分布、历史地震及库仑应力变化,我们认为未来需要密切关注与磨西断裂交接的安宁河、大凉山断裂以及康定—磨西段的地震危险性.     

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。门源地震破裂至地表是造成该地区基础设施破坏的直接原因。     

基于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),处于应力加载状态,未来发生强震的危险性较高.     

InSAR观测与小震分布联合分析2021年云南漾濞MS6.4地震应力卸载对周围区域的影响

为分析2021年5月21日云南省漾濞M_S6.4地震后震区应力变化对周围断层的影响,本文通过InSAR技术获得了漾濞地震的同震形变场,并联合小震分布数据建立断层破裂滑动模型,继而通过计算断层面上的同震库仑应力来评估此次地震对周边断层的影响,以便有效地分析地震破裂的时空解析度.结果显示:(1)在升降轨InSAR数据获得的精细同震形变场中,升轨最大视线向形变量约为5.00 cm,降轨最大视线向形变量约为7.80 cm;(2)余震精定位的主震震中为(99.89°E,25.67°N),震源深度为13.29 km,除主震之外震源深度主要集中在5—15 km;通过小震位置拟合出的发震断层走向为NW-SE(316.69°),断层倾角为88.56°,滑动角为177.97°;(3)基于InSAR同震形变场结果及小震拟合断层参数联合反演得到此次地震的断层滑动以右旋走滑为主,升轨断层最大滑动量为0.80 m,对应的深度为8.85 km,平均滑动量为0.22 m,矩震级为M_W6.41;降轨的断层最大滑动量为0.30m,对应的深度为6.88 km,平均滑动量为0.05 ...     

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

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,发震构造为具有少量走滑性质的逆冲断裂。从发震构造特征、同震滑动分布推测,伽师地震发震构造是柯坪塔格褶皱带滑脱面以上沉积盖层内的逆冲断裂,支持了柯坪推覆体的薄皮构造模型观点。     

GNSS观测的九寨沟7.0级地震同震位移初步结果

2017年8月8日四川省九寨沟县发生了7.0级地震,中国大陆构造环境网络与北斗地基增强系统的GNSS连续观测共同监测到了此次地震的同震位移（坐标：东向为正,北向为正）,结果显示：3个站点记录到了明显的同震位移,距离震中43 km的九寨沟台站（SCJZ）在东西向的位移为-9.8±1.5 mm,在南北向的位移为3.3±0.7 mm;距离震中65 km的松潘站（SCSP）在东西向的位移为-1.8±0.7 mm,在南北向的位移为-7.7±0.6 mm;距离震中77 km的舟曲站（GSZQ）在东西向的位移为0.4±1.2 mm,在南北向的位移为3.6±0.8 mm.通过同震位移分布特征,可以推测此次地震为一次左旋走滑型事件,引起水平向同震位移大致不超过150 km范围,地震对东南侧的龙门山断裂带影响非常小,对北侧的塔藏断裂和西侧的岷江断裂处引起的同震位移为厘米级.同震位移的反演结果显示：断层面上滑动量主要集中在7 km深度,最大量值约为0.4 m,平均滑动角为-15°,利用滑动分布计算的相应矩震级为M_W6.4,与地震波反演结果相当.结合同震滑动分布、同震主应变分布、余震分布和震源机制解等特征,推测此次地震破裂极值区累积的能量得到较充分释放,进一步分析得出此次地震在塔藏断裂、岷江断裂和虎牙断裂处产生了一定的应力变化,值得持续关注.     

Numerical Simulation on Coseismic Effect of the November 14,2001 Great Kunlun Earthquake,Northern Tibet,China

The November 14, 2001 Ms8.1 Kunlun Mountains earthquake in northern Tibet is the largest earthquake occurring on the Chinese mainland since 1950. We apply a three-dimensional （3-D） finite element numerical procedure to model the coseismic displacement and stress fields of the earthquake based on field investigations. We then further investigate the stress interaction between the Ms8.1 earthquake and the intensive aftershocks. Our primary calculation shows that the coseismic displacement field is centralized around the east Kunlun fault zone. And the attenuation of coseismic displacements on the south side of Kunlun fault zone is larger than that on the north side. The calculated coseismic stress field also indicates that the calculated maximal shear stress field is centralized around the east Kunlun fault zone; the directions of the coseismic major principal stress are opposite to that of the background crustal stress field of the Qinghai-Xizang （Tibet） Plateau. It indicates that the earthquake relaxes the crustal stress state in the Qinghai-Xizang （Tibet） Plateau. Finally, we study the stress interaction between Ms8.1 earthquake and its intensive aftershocks. The calculated Coulomb stress changes of the Ms8.1 great earthquake are in favor of triggering 4 aftershocks.     

日本9.0级大震体应变同震分析

以江苏溧阳地震台体应变仪记录到的2011年3月11日日本本州东海岸附近海域9.0级地震为基础,分析记录的应变波特征,计算出震中距和震级等参数。通过初步了解体应变记录远震、大震的特征,对进一步分析、研究体应变映震能力和震源机制有一定的借鉴作用。     

重庆荣昌井水温同震-震后响应特征及机理研究

利用位于华蓥山断裂带的重庆荣昌井水温数字化分钟值观测资料,统计分析了该井水温对2008年1月—2021年9月全球M_S≥7.0、川滇地区M_S≥6.0、重庆及周边地区M_S≥4.0共273次地震的同震-震后响应动态特征,对井水温同震优势方向成因和机理进行了深入研究,获得以下认识:①荣昌井水温同震-震后响应能力较好,对近震和远震均可记录到;该井水温同震响应由深及浅的顺序发生,响应持续时间随观测深度的增加而增加,响应幅度随观测深度的增加而减小,且该井水温同震-震后响应持续时间较同井水位的长;②自观测以来,荣昌井多层水温同震响应方向均为上升,说明单个井水温对不同地震的同震响应存在优势响应方向,水温的同震特征更依赖于井孔自身观测条件的影响;荣昌井水温同震响应优势方向上升可能是地震波的扰动造成井下深部气体释放,并沿裂隙上升进入井含水层系统而引起;③荣昌井水位-水温对中、远场地震的同震为同向上升正相关关系或振荡—上升,对近场地震的同震为水位下降—水温上升的反相关关系,可能是近场地震和中、远场地震引起的水位同震响应变化机制不一致所致;引起荣昌井水温同震变化的地震波能量密度e(r)＞10^-5J/m³,而引起荣昌井水温和水位同震反向变化的地震波能量密度e(r)＞1J/m³。     

用多种数据构建2008年汶川特大地震同震位移场

本文主要以GPS、精密水准观测和卫星SAR遥感图像分析2008年汶川特大地震同震位移特征.GPS数据包括:(1)四川盆地和川西高原地区各类国家等级GPS网点复测;(2)沿破裂带国家天文大地网GPS复测.前者推算的同震位移测定精度优于2 cm,后者6~8 cm.SAR遥感资料包括:(1)ALOS 卫星升轨相位干涉图像,精度优于8 cm;(2)ALOS和ENVISAT卫星影像合成的三维位移图,精度优于0.5 m.同震位移场显示,断层下盘(四川盆地)变形总体呈扇形集中指向震中,断层上盘(龙门山)变形总体上呈逆时针旋转态势,最大的实测水平位移5.5 m.汶川、理县、茂县等地测站位移指向破裂带方向,而平武、青川等地测站逐渐转变为平行,乃至远离破裂带方向,与汶川地震逆冲兼走滑的破裂特征一致.断层上盘大幅隆升,下盘靠近断层的区域以下沉为主,远场表现为幅度很小的隆升,垂直升降区域间,有一条与龙泉山断裂带平行的升降过渡带,调节龙泉断层的应力状态.用实测变形场检验多个地震波破裂模型表明,近场(距离断层50 km) 模型形变准确度可达40~50 cm, 远场精度优于5 cm.     

强震变形的地震大地测量监测与研究

地震大地测量是将大地测量（特别是空间大地测量）与地震学及构造地质学进行融合的新兴交叉学科,其可用于监测地震孕育的地球物理背景场及动态变化过程,对相关形变实现了10²a～10^-2s的宽频带监测,基本弥补了地震学与构造地质学间的频率空白。以多频带的地震大地测量技术（GNSS、InSAR、高频GNSS）为支撑的陆态网络工程,不仅获得了中国大陆长期的地壳运动图像,而且在强震形变监测中发挥了重要作用。汶川、芦山、尼泊尔廓尔喀及九寨沟等地震的研究成果表明,高频/静态GNSS、InSAR、精密水准相融合的多频大地测量,极大地拓展了地震形变监测的时空频域,促进了大陆型地震的相关研究,为地震预测预警研究奠定了基础。然而,目前使用的地震大地测量资料有限,同时,我国地震大地测量监测网络也有待不断加密和优化。     

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

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

基于子带干涉测量技术的巴基斯坦地震形变获取研究

2013年9月24日发生在巴基斯坦俾路支省(Balochistan)境内阿瓦兰县(Awaran)的MW7.7级地震,在地表产生了最大达10m的滑动量.利用TerraSAR-X短波雷达数据获取的InSAR同震形变场产生了密集且大范围的干涉条纹,给后续的相位解缠带来困难.而子带干涉法是一种无需或只需进行少量相位解缠,即可获得绝对相位差的新方法.其主要思路是通过缩减带宽以增长波长,从而减少干涉条纹数,降低解缠难度或不需解缠直接得到绝对相位差.但由于带宽的缩减,导致噪声的增大和旁瓣带来的额外干扰,使干涉图质量下降,因此在子带干涉参数选取、噪声滤波以及处理流程等方面需要特殊处理,特别是子带的中心频率和带宽的选取会很大程度地影响测量精度.首先选取典型DEM实验区,以干涉图相干性和误差为评价指标,利用逐步参数选取法,研究相关参数对子带干涉测量的影响,制定最优的参数方案,认识参数选取的原则和方法.在此基础上,将子带干涉应用于巴基斯坦地震的同震形变场获取.最后,将子带干涉、Landsat 8光学影像的交叉频谱相关法、offset-tracking、常规DInSAR获取的同震形变场进行比较,并与模型拟合的形变场进行对比分析.结果表明,子带干涉虽然会受失相干的影响,其提取的形变场范围相较于Landsat 8和offset-tracking有所缺失,但在共同覆盖的区域其精度和噪声水平更优,相比较于常规DInSAR,更适用于条纹密集和形变量大的地区.     

山西离石中心地震台洞体形变观测异常分析

通过对离石中心地震台洞体形变近6年观测资料中典型异常形态及与其相关地震的对应关系进行分析,结果表明,异常形态主要表现为短期趋势转折,突跳、毛刺、加粗等的数据波动,固体潮畸变伴掉格,且大多数异常是3种形态的组合;多数异常在震前发生,震后迅速恢复,固体潮分钟值曲线上有同步变化;水平摆和伸缩仪对应国内外6级以上的地震出现不同程度的异常反映,水管仪表现不明显。