基于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。综合分析认为,本次地震的震源特征为右旋走滑兼有少量正断分量,是发生在公格尔拉张系内拉张环境下的构造地震事件。     

2016年新疆阿克陶MS6.7地震同震形变特征和断层滑动分布反演

本文针对2016年11月25日在新疆阿克陶发生的地震, 用差分干涉测量技术(D-InSAR)对3种不同观测模式的升、 降轨数据进行处理, 提取了多视线向的同震形变场; 根据不同模式的LOS向形变量, 构建形变分解模型, 将其分解为垂直向形变量和沿断层走向形变量; 结合同震形变场特征与震源机制解, 采用单断层模型, 利用梯度下降法(SDM), 以Multi-LOS向形变进行约束, 反演了阿克陶地震的同震滑动分布特征。 研究结果表明, 升、 降轨LOS向同震形变场在发震断层两侧具有不同的形变特征, 发震断层走向近EW向; LOS向形变量分解表明, 此次地震破裂以右旋走滑为主; 滑动分布反演的形变残差介于0~5 cm之间, 发震断层的滑动量主要位于2~16 km深部, 最大滑动量可达1.02 m, 位于断层面深部5.83 km处, 最大滑动量处的滑动角为185.24°; 平均滑动角为181.32°, 平均滑动量为0.12 m; 滑动分布反演也证明该地震为右旋走滑破裂事件, 与LOS向形变分解结果一致; 当剪切模量μ=3.2 GPa时, 反演得到的地震矩震级约M_W6.6。     

利用InSAR技术研究2016年青海门源MW5.9地震同震形变场及断层滑动分布

2016年1月21日, 青海省门源县冷龙岭断裂带附近发生了M_W5.9地震。 基于Sentinel-1A影像, 采用差分干涉雷达测量技术研究了此次地震产生的同震形变场, 结果表明, 门源地震的形变影响范围约20~30 km, 形变态势在升降轨道形变场均显示为隆升, 基本沿冷龙岭断裂呈近似同心圆展布, 推测可能是冷龙岭断裂与民乐—大马营断裂之间的一条逆断层, 沿雷达视线方向最大形变量级约为6 cm。 均匀滑动反演显示门源发震断层长7.3 km, 宽6.2 km, 走向298.6°, 倾角34.5°, 倾向宽度9.5 km, 沿走向滑动量为170 mm, 沿倾向滑动量为460 mm, 矩震级为M_W5.97; 分布式滑动反演显示门源地震以逆冲为主, 兼具少量右旋走滑分量, 滑动量主要集中在沿断层倾向方向, 距离地表5～15 km处, 最大滑动量约0.3 m, 位于断层倾向深度10 km处, 矩震级为M_W5.93。     

尼泊尔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数据反演的结果一致.     

Sentinel-1A及ALOS2数据约束下的阿克陶地震震源参数反演及发震构造研究

2016年11月25日,在新疆阿克陶县发生了M_W6.6地震.本文利用Sentinel-1A宽幅数据和ALOS2精细数据获得了同震形变场,干涉形变场沿木吉断裂展布,显示本次地震破裂长度可达70 km,在地表形成两个形变中心,且震中东部形变场条纹密集而西部稀疏、影响范围广,跨断层剖面显示视线向最大形变量可达12 cm.利用一种结合先验知识的多视角最小二乘迭代分解法求解出地表的三维形变场,结果显示震中东部形变中心垂向向下运动最大可达20 cm,木吉断裂南侧西向运动可达10 cm,断层表现为右旋走滑兼具正断作用.采用前向模拟的方法确定二段式分段断层模型能够较好地恢复观测形变场,进而以InSAR观测数据为约束,基于弹性半空间形变模型采用两段非均匀断层滑动模型来反演断层面上的精细滑动值.反演结果显示本次地震可能为两次地震事件,发生在断层西段的第一次事件以右旋走滑为主,走向103°,倾角76°,滑动角-167°,震源深度10.1 km,累计地震矩为7.2×10¹⁸ N·M;东段的第二次事件为右旋走滑兼具正断作用,走向109°,倾角略缓约55°,滑动角-160°,震源深度5.3 km,累计地震矩7.76×10¹⁸ N·M.本次地震是一次发生在公格尔拉张系的拉张环境下的构造地震.

     

2008年MS7.1于田地震InSAR同震形变场及其震源滑动反演

本研究利用InSAR技术与ALOS PALSAR雷达数据,获取了2008年3月20日于田MS7.1地震视线向同震形变场,并基于该数据集和限制性最小二乘算法反演了此次地震的断层滑动分布;通过构造四大类反演方案,详细分析了InSAR观测系统中的入射角与方位角对反演结果的影响.结果表明:入射角随点位变化对反演结果有较大影响,使用其平均值将对破裂细节产生一定影响;而方位角对反演结果的影响不大,使用其平均值是一种较为理想的选择;引入入射角与方位角变化后,反演获得了较佳的于田地震同震滑动,主要集中分布于0~14 km深度附近,最大滑动量达3.2 m,矩张量为3.3×1019 N·m,相当于矩震级MW7.0.     

2022年青海门源MW6.6地震InSAR同震形变场与震源特征

本文利用合成孔径雷达差分干涉(Differential Interferometric Synthetic Aperture Radar, D-InSAR)技术处理哨兵卫星升降轨影像,获取了2022年青海门源M_S6.9地震的同震变形场,变形结果显示此次地震导致冷龙岭断裂西段两侧约30 km×20 km范围明显变形,卫星视线向最大同震位移达50 cm.基于有限断层模型,本文构造断层滑动分布模型,反演同震变形场数据.模型显示,冷龙岭断裂的破裂以左旋走滑为主,破裂长度约23 km,主要分布在10 km以上深度,整体形态表现为不规则椭圆盘,其中约3.9 m的最大破裂位于5～6 km深度的椭圆中心.基于反演模型估算的地震平均应力降约5.9 MPa,地震标量矩为1.0×10¹⁹ N·m,对应矩震级为M_W6.6.本文计算了1927年古浪地震、1986年及2016年门源地震对本次地震破裂面的同震静态库仑应力扰动,结果表明冷龙岭断裂西段受到历次地震的库仑应力加载作用,库仑应力增加累积达0.12 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地震的破裂空区完全破裂.考虑到塔藏断裂和虎牙断裂的运动性质,可初步判定发震断层为虎牙断裂北侧延伸分支.基于三种不同接收断层模型的同震库仑应力变化计算结果反映出该区域以应力释放为主,进一步触发较大走滑型余震的可能性不大.

     

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

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

2016年1月21日门源Ms6.4级地震InSAR同震形变及发震构造研究

2016年1月21日青海省门源县发生M_S 6.4级地震,距离震中位置最近的已知断裂为位于其北侧7 km的冷龙岭断裂,该断裂晚第四纪以来主要表现为左旋走滑运动,局部兼具倾滑分量.但本次地震的震源机制解显示为逆冲型地震,与人们认知的冷龙岭断裂走滑运动性质有所差异.本研究利用Sentinel-1A数据获取了升降轨方向的形变场,形变场特征显示本次地震以逆冲为主,最大形变量在6 cm左右.以此数据为约束,采用两步法反演了发震断层参数,结果显示仅以InSAR数据为约束并不能唯一确定断层的倾向,因此本文对比了NE倾向和SW倾向等不同断层模型的反演结果,综合分析该地区的地质构造背景,认为SW倾向的断层模型更加合理,本次地震由冷龙岭北侧的弧形次生断层引起,发震断层面走向141°,倾向SW,倾角43°,平均滑动角为72.7°,最大滑动量为0.44 m,反演矩震级为M_W 6.0,震源深度10 km.该次生断层与冷龙岭断裂一起构成正花状构造,冷龙岭主干断层近直立插入基底,夹持部分形成隆起断块,推测本次地震是青藏高原向NW推挤生长,在压扭性作用下隆起断块的一种表现.     

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形变研究与模拟

采用欧洲空间局ERS-2的星载干涉雷达数据,选取1997年11月8日MW7.6级玛尼地震作为研究对象,采用了差分干涉方法,在通过对覆盖同一地区的SAR数据进行差分干涉处理,得到了玛尼地震的视线向同震形变场。经研究发现:该地震形变场呈长轴近北东东向不规则椭圆形分布,地表破裂带长度约为130km,发震断层走向约为78°,断裂为左行走滑特征。断层以南为隆起区,在发震断层附近最大视线向隆起位移量为113.6cm,断层以北为沉降区,最大视线向沉降位移量为170.4cm。基于Okada模型实现了具有复杂结构的4段断层段参数的InSAR形变场数据模拟,获得断层的最大走滑为6m,估计出玛尼地震的标量地震矩M0为2.69×10^20Nm,计算得到的矩震级MW为7.6。证明了研究方法的正确性和研究结论的可靠性。     

基于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。     

基于GIS的当雄同震InSAR形变场分析

InSAR技术对同震形变的量测达到了厘米级的精度,但在数据后处理和结果分析方面仍然存在明显的不足。采用可对空间信息进行存储、管理及分析等功能的GIS,对InSAR数据获取的2008年10月6日16时30分西藏自治区当雄县MW6.3地震的同震形变场进行分析。结果表明:(1)GIS可对多源数据进行有效管理;(2)可以确定当雄地震的震中位置;(3)可获取沉降区的最大形变量;(4)可确定主要的形变区间;(5)可将形变结果进行三维展示。GIS可有效地弥补InSAR数据后处理、数据分析及成果展示方面的不足。     

2008年于田MW7.2地震滑动分布反演及三维同震形变场解算

基于ENVISAT ASAR升降轨数据, 利用InSAR获取2008年于田M_W7.2地震同震形变场; 采用SDM反演本次地震断层滑动分布; 使用PSCMP正演获取于田M_W7.2地震南北向模拟形变量, 并结合升降轨同震形变场, 解算三维同震形变场。 同震形变场分析表明, 2008年于田M_W7.2地震以正断为主, 且带有走滑运动特征, 破裂带走向为NNE向。 同震滑动分布反演结果显示, 断层沿走向被分为4段F1、 F2、 F3、 F4, 其滑动分布集中在0~14 km区间, 以F2、 F3段为主, 最大滑动量约5.31 m, 位于F2段深部2.76 km处; 沿破裂带走向, 左旋走滑位移与垂直位移比值有增大的趋势; 反演获得的地震矩M₀=5.58×10¹⁹N·m, 相当于矩震级M_W7.1。 三维同震形变场解算结果显示, 断层上盘整体表现为沉降, 断层下盘整体表现为隆升, 且沉降量明显大于隆升量, 表明地震以正断破裂为主; 除靠近断裂带中上部表现为向东南运动外, 上盘整体上表现为向西南运动; 断层下盘则整体表现为向东北运动, 证明破裂兼有左旋走滑运动。 滑动分布反演、 正演与三维同震形变场解算结果皆表明, 于田M_W7.2地震破裂以正断为主, 且带有一定的左旋走滑。     

基于InSAR同震形变观测反演2010年新西兰南岛Mw7.1 Darfield地震同震破裂分布

2010年9月4日新西兰南岛Canterbury平原发生了M_w7.1地震,震源深度约为10 km.本次地震发生在一条震前不为人所知的断层上.我们利用覆盖整个震区的合成孔径雷达(SAR)观测资料,通过干涉处理分析获得雷达视线向(LOS)同震形变场;以此资料为约束反演了断层的几何参数以及同震破裂分布.结果显示,该地震造成四条相对独立断层的破裂.大部分的地震矩释放发生在Greendale断层(编号1-4),其错动以右旋走滑为主,最大破裂约为8.5 m.其它三条断层中,经过震源的逆冲断层最大破裂为5.1 m (编号6),位于Greendale断层以西的逆冲断层最大破裂为3.5 m (编号5),位于Greendale断层北面的走滑断层最大破裂为1.9 m(编号7).反演的Greendale断层地表滑动与地质调查得到的地表破裂在形态和数值上均吻合较好.本次地震释放的地震矩为5.0×10¹⁹N·m,矩震级为7.1.板块边界带形变场分析表明,Darfield地震的发生受边界带应变分配在该地区残留构造应力场控制,其复杂性体现了区域构造应力场的特点.地震对其周围地区的应力场影响较大,库仑应力增加区与余震分布有一定对应关系,并在2011年Christchurch 6.3级地震发震断层区域造成约0.1bar的库仑应力增加,对此地震有一定的触发作用.     

青海门源MS6.9地震震前OLR异常与InSAR同震形变关联性分析

对比分析利用涡度距平法提出的2022年1月8日青海门源 M_S6.9 地震震前射出长波辐射(OLR)短期异常分布和震后InSAR技术提取的门源地震同震形变空间分布,结果显示,震前红外辐射增强区与InSAR同震破裂形变区的空间位置基本吻合,扩展形式基本相似(同震破裂形变区分布在红外辐射异常区内部)。在震前的全国范围OLR空间分布上,仅青海德令哈—西宁—甘肃武威一带出现了呈“哑铃”状近WE向展布的OLR热辐射增强区,空间可辨识度高,OLR异常时空演化过程遵循了岩石应力加载破裂过程中的热异常规律,显示热异常变化与应力变化存在关联; InSAR技术提取的同震形变同样位于肃南—祁连断裂(俄堡段)、托莱山断裂和冷龙岭断裂的交汇区。InSAR同震形变结果揭示了地表形变以水平方向为主,断层运动具有典型的走滑变形特征。InSAR同震形变结果为红外遥感反映地震形变提供可检验的地质实体监测证据,验证了门源地震前辐射增强异常是地震构造地应力强度变化的遥感物理参量反映。