Source parameters of the 2009 L'Aquila earthquake,Italy from Envisat and ALOS satellite SAR images

2009年4月6日意大利L’Aquila地区发生了Mw6.3级地震,该地震造成了300余人的人员死亡.本文联合不同波长、不同入射倾角的升降轨Envisat和ALOS卫星的差分干涉数据对该地震进行震源机制解的反演研究.研究首先对卫星雷达影像进行二通差分干涉处理,获取了覆盖L' Aquila地震震区的完整InSAR同震形变场,然后结合四叉树和均匀采样方法对原始观测数据进行降采样.在此基础上,联合GPS形变观测数据,利用弹性半空间矩形和三角位错模型,以及断层自动剖分技术对断层面进行最优离散剖分,反演获取了发震断层的精确几何参数和最优断层滑动分布,结果显示分布式三角位错滑动模型能够很好地解释观测到的地表形变场.反演结果表明发震断层是一个以正倾滑为主兼有少量右旋走滑的盲断层;基于观测数据最优确定的断层单元的最短边长为0.4 km,最长边长为6.3 km;此次地震的滑动分布主要发生在5～14 km深度的范围内,最大滑移量为1.07 m,释放的能量为3.43×1018N·m(Mw6.32),与地震学的研究结果非常一致.     

利用D-InSAR技术研究西藏改则地震同震形变场

针对2008年1月9日MW6.4西藏改则地震和2008年1月16日的MW5.9余震,通过两通(2-pass)加外部DEM差分干涉处理技术(D-InSAR),提取了地震区域2次地震累积的视线向(LOS)同震形变场。结果表明:发震断层均为正断层,位于依布茶卡-日干配错断裂端点附近。主震发震断层走向为N30°E,余震发震断层走向为N21°E,两断层距离约7km;在影像上主震发震断层有造成地表破裂的痕迹,余震未见地表破裂的痕迹;这次地震造成的同震形变场长约30km,宽约20km,主震断层上盘和下盘视线向最大形变量分别为39.2cm和11.2cm,两盘相对位错达50.4cm,余震造成的视线向形变量为9.4cm     

基于方差分量估计的2015年尼泊尔MW7.8地震同震滑动分布

基于三角位错的断层自动剖分方法,联合ALOS-2 InSAR和GPS同震位移数据反演了2015年尼泊尔M_W7.8地震的发震断层参数及滑动分布.反演中基于赫尔默特方差分量估计方法确定InSAR、GPS水平向与垂直向的3种观测位移数据的相对权比.反演结果显示,方差分量估计方法定权得到的最大滑动量5.5 m,大于验前方差定权给出的5.1 m的结果,同时前者显示同震破裂沿东南方向传播至接近M_W7.3最大余震区域,然后分别继续向北和东南传播约30 km,并形成一个滑动空区（slip gap）,而后者的结果没有表现出这种显著的滑动空区特征.方差分量估计定权得到的滑动模型能更好地解释观测到的InSAR同震形变场,其拟合残差标准差较验前方差定权的结果减小了2.5 cm,GPS的东和北分量的残差标准差分别减小了0.4 cm和0.9 cm.与两种不同的基于矩形位错的方法相比,利用方差分量估计的三角位错断层自动剖分方法的滑动误差均值和标准差相对最小,其中均值为0.037 m,标准差为0.028 m,表明该方法得到的滑动模型具有相对更好的无偏性和有效性.

     

基于多期DInSAR数据的2016年定结地震序列发震断层特征研究

2016年5月22日,在西藏定结县发生四次M_W4~5地震,研究本序列地震的发震断层几何和运动特征对于认识周边活动断裂性质具有重要意义.由于发震地区偏远,且观测台网分布稀疏,本研究采用星载雷达干涉测量（DInSAR）技术进行了同震形变场重建,但是定结地震震级较小,单干涉像对获取的形变场受相位噪声影响较大.为了解决这一问题,本研究基于时间序列Sentinel-1A干涉数据生成多期同震与非同震干涉图,并利用叠加平均法对本次定结地震同震形变场进行重建,提取了定结2016年5月22日多次地震产生的同震累计整体形变场.基于InSAR同震形变场和区域地质特征,研究进行了滑动分布反演,确定其主要贡献的发震断层几何参数及滑动分布：断层走向为188°,倾角为43°,平均滑动角为78°,发震断层的运动性质以正断为主兼具少量左旋走滑分量,滑动主要集中在断层垂直深度0~9 km处,最大滑动量约为25 cm,位于断层倾向深度3 km处,反演得到的矩震级为M_W5.58.本研究结果表明采用星载InSAR叠加平均技术可以较好地压制相位噪声,有效提取此类中小型浅源地震同震微弱形变场.最后,我们认为本次定结地震与藏南拆离断层与申扎-定结断层的活动密切相关.

     

利用InSAR研究乌恰地震同震形变

2008年10月5日在新疆乌恰地区39.50°N、73.64°E处发生MW6.7地震,用ROI_PAC软件处理4景日本ALOS卫星数据得到乌恰地震的三维同震形变场。距离向,断层北侧最大位移达39 cm,南侧最大位移达36 cm;方位向,北侧和南侧的最大位移分别为1.5 m和2 m。采用Okada弹性位错模型分析该地震的滑动分布,模拟得出的断层走滑角为48°,倾滑角为53°,深度为10.5 km,最大形变量2.5 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.     

基于线弹性位错模型及干涉雷达同震形变场反演1997年玛尼Mw75级地震参数－Ⅰ.均匀滑动反演

1997年11月8日西藏玛尼Mw75级地震是干涉雷达技术应用于地震观测以来的一次重要事件.Peltzer等认为该形变场体现了地壳的非线性弹性形变响应，挤压象限和拉张象限的杨氏模量比值为2.相同的情形并没有在其他地震同震形变场中发现，通常用线弹性理论就能够很好地解释观测数据.考虑到他们所用断层模型的简化程度和纯走滑约束等，本研究认为非线性弹性解释是牵强的.应用广泛使用的Okada线弹性位错模型，采用卫星观测得到的地表断层位置，去除倾滑为零的约束，基于非线性优化反演方法寻求拟合雷达观测的最佳断层几何参数和均匀滑动参数.结果表明，线弹性模型能够满意解释观测数据.断层在朝阳湖以东的段落最深达到165 km，随着断层向两侧延伸，深度逐渐减小；反演得到的断层倾     

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...     

基于InSAR形变观测反演2015年皮山MS6.5地震滑动分布

2015年7月3日在新疆皮山县发生了M_S6.5地震, 该地震使当地遭受了巨大的经济损失。 本文利用欧空局提供的Sentinel-1A卫星差分干涉数据对该地震的震源机制情况进行了反演研究, 首先运用两轨法对卫星雷达影像进行差分干涉处理, 获取了覆盖皮山地震震区的同震形变场, 然后利用弹性半空间的均匀滑动模型反演获取了发震断层的几何参数, 并对原始观测数据进行降采样处理, 在此基础上运用分布式滑动模型反演获取了更为精细的断层滑动分布, 结果显示分布式滑动模型与观测结果有很高的拟合度。 反演结果表明发震断层是一个以逆冲为主兼有极少量左旋走滑的盲断层, 此次地震断层面的同震活动分布主要集中在7~15 km深度范围内, 同震的地震矩为6.28×10¹⁸N·m, 矩震级为M_W6.46, 与前人的研究结果非常一致。     

2010年青海玉树MS7.1级地震地表破裂带和形变特征分析

通过对SPOT卫星影像上地表破裂的目视解译,以及对ALOS PALSAR卫星数据进行D-InSAR形变提取和分析,结合地震活动性、震源机制、活动构造等资料,确定了发震断层空间分布、断层性质和同震形变场分布特征.结果显示,玉树地震发生在甘孜—玉树断裂带上,总体走向约为300°,断层近乎直立.根据相干性强弱将Ⅰ区地表破裂划分为三段:北段长22 km,中段长5 km,南段长6 km,破裂带总长度约33 km.Ⅱ区内非相干带长约10 km.同震形变场分布在78 km×55 km范围内,主震所在的形变Ⅰ区断层两侧视线向相对位错约为0.78 m,转换成水平位错约为1.5 m;余震所在的Ⅱ区形变相对较小.     

2010年青海玉树MS7.1级地震地表破裂带和形变特征分析

通过对SPOT卫星影像上地表破裂的目视解译,以及对ALOS PALSAR卫星数据进行D-InSAR形变提取和分析,结合地震活动性、震源机制、活动构造等资料,确定了发震断层空间分布、断层性质和同震形变场分布特征.结果显示,玉树地震发生在甘孜—玉树断裂带上,总体走向约为300°,断层近乎直立.根据相干性强弱将Ⅰ区地表破裂划分为三段:北段长22 km,中段长5 km,南段长6 km,破裂带总长度约33 km.Ⅱ区内非相干带长约10 km.同震形变场分布在78 km×55 km范围内,主震所在的形变Ⅰ区断层两侧视线向相对位错约为0.78 m,转换成水平位错约为1.5 m;余震所在的Ⅱ区形变相对较小.     

基于线弹性位错模型反演1997年西藏玛尼Mw7.5级地震的干涉雷达同震形变场——Ⅱ滑动分布反演

1997年11月8日西藏M_w7.5级玛尼地震是干涉雷达技术应用于地震观测以来的一次重要事件.在第一部分中，我们应用广泛使用的Okada线弹性位错模型，假设断层的各个分段滑动量均匀，反演得到断层各个分段的几何参数和均匀滑动量.本部分的反演进一步去除滑动均匀假设，并利用更能反映断层真实状态的角形元位错模型（线弹性），在第一部分反演得到断层几何的基础上，反演断层面的静态位错分布.反演结果表明，线弹性滑动分布模型能够更好地解释观测数据，进一步提高反演的数据拟合程度.最终得到了断层面上的走滑和倾滑位错分布.首次得到的断层面滑动分布显示断层面滑动在浅部（0～12 km）比较集中，地震破裂长度约170 km，最大左旋走滑位移达4.8 m；反演结果还表明局部段落存在较大倾滑位移，量值达到1.9 m，这在断层模型中是不能忽略的，它可能是断层两侧形变不对称的原因之一；反演得到的标量地震矩为2.18×10²⁰ N·m，相当于矩震级M_w7.5，与Velasco等利用地震波形反演得到的结果一致.     

2021年青海玛多MW7.3地震InSAR的同震形变场、断层滑动分布及其对周边区域的应力扰动

文中基于D-InSAR技术,利用欧空局C波段升降轨哨兵SAR数据,获取了2021年5月22日青海玛多MW7.3地震的InSAR同震形变场,并对同震形变的空间特征、量级和断层破裂的分段性进行了分析.哨兵卫星的高质量观测数据清晰地描绘了玛多地震的地表破裂迹线,地表破裂长度约210km.为了进一步认识玛多地震断层深部的同震滑...     

基于线弹性位错模型反演1997年西藏玛尼Mw7.5级地震的干涉雷达同震形变场——Ⅱ滑动分布反演

1997年11月8日西藏M_w7.5级玛尼地震是干涉雷达技术应用于地震观测以来的一次重要事件.在第一部分中，我们应用广泛使用的Okada线弹性位错模型，假设断层的各个分段滑动量均匀，反演得到断层各个分段的几何参数和均匀滑动量.本部分的反演进一步去除滑动均匀假设，并利用更能反映断层真实状态的角形元位错模型（线弹性），在第一部分反演得到断层几何的基础上，反演断层面的静态位错分布.反演结果表明，线弹性滑动分布模型能够更好地解释观测数据，进一步提高反演的数据拟合程度.最终得到了断层面上的走滑和倾滑位错分布.首次得到的断层面滑动分布显示断层面滑动在浅部（0～12 km）比较集中，地震破裂长度约170 km，最大左旋走滑位移达4.8 m；反演结果还表明局部段落存在较大倾滑位移，量值达到1.9 m，这在断层模型中是不能忽略的，它可能是断层两侧形变不对称的原因之一；反演得到的标量地震矩为2.18×10²⁰ N·m，相当于矩震级M_w7.5，与Velasco等利用地震波形反演得到的结果一致.     

近场位移数据约束的2013年芦山地震破裂模型及其构造意义

以往的研究显示了2013年芦山MS7.0级地震发震断层的隐伏逆冲断层基本特征,但是破裂深部细节差异较大.本文以近场密集的同震形变数据约束芦山地震破裂面几何形状及滑动分布, 结果显示芦山地震破裂面具有铲状结构,上部16 km为43°~50°高角度断层,深部16~25 km为小于27°的低角度断层,破裂深度与重定位的余震分布深度一致.破裂分布模型清楚显示上下两个断层上各有一个滑动幅度大于0.5 m的峰值破裂区,最大滑动量1.5 m位于13 km深处.重定位的余震分布基本都落在最大滑动量等值线外部库仑应力增加的区域.芦山地震破裂面几何形状和滑动分布特征与2008年汶川MS8.0级地震映秀—北川破裂相似,支持龙门山冲断带发育大规模的近水平滑脱层, 是青藏高原东缘地壳缩短增厚、龙门山挤压隆升的重要证据.     

基于线弹性位错模型及干涉雷达同震形变场反演1997年玛尼Mw75级地震参数－Ⅰ.均匀滑动反演

1997年11月8日西藏玛尼Mw75级地震是干涉雷达技术应用于地震观测以来的一次重要事件.Peltzer等认为该形变场体现了地壳的非线性弹性形变响应，挤压象限和拉张象限的杨氏模量比值为2.相同的情形并没有在其他地震同震形变场中发现，通常用线弹性理论就能够很好地解释观测数据.考虑到他们所用断层模型的简化程度和纯走滑约束等，本研究认为非线性弹性解释是牵强的.应用广泛使用的Okada线弹性位错模型，采用卫星观测得到的地表断层位置，去除倾滑为零的约束，基于非线性优化反演方法寻求拟合雷达观测的最佳断层几何参数和均匀滑动参数.结果表明，线弹性模型能够满意解释观测数据.断层在朝阳湖以东的段落最深达到165 km，随着断层向两侧延伸，深度逐渐减小；反演得到的断层倾     

汶川Ms8．0级地震InSAR同震形变场观测与研究

利用InSAR技术,采用地震前后日本ALOS/PALSAR数据,提取了2008年5月12日四川汶川地震7个条带的地表同震形变场.每个形变条带南北向500 km,东西向70 km,7个形变场覆盖了映秀镇、都江堰、茂县、北川、平武和青川.结果显示,此次地震地表破裂带在北川-映秀断裂带上.地表破裂带从汶川县映秀镇西南震中附近一直到青川县苏河北侧,全长约为230 km.发震断层西北盘为抬升盘,南东盘断层附近,仍然表现为隆起区,显示出以逆冲为主的断层性质.在汶川县映秀镇西侧震中区,最大相对卫星视线向形变量达260 cm,如果全部换算成垂直形变,则两个区域的垂直相对形变达3.3 m.从北川至平通一带,卫星视线向形变范围在120～180 cm的隆起带.其中,擂鼓镇隆起形变范围170～180 cm,换算成垂直形变约在2.2～2.3 m之间.在青川苏河北附近,有70～80 cm范围的隆起形变.在雅安、峨眉山一带,以及射洪至重庆北侧一带有大范围沉降区.在重庆及其南侧区域有幅度在20～30 cm小范用隆起.由青川向东至广元、宁强一带,有形变幅度在60～70 cm的隆起区.整个形变场影响范围较大,四川盆地出现了不同程度的地表形变.     

InSAR约束下的2008年汶川地震同震和震后形变分析

2008年5月12日,青藏高原东缘的龙门山断裂带上发生了Mw7.9级汶川地震.本文通过分析覆盖汶川地震震中区域的ALOS/PALSAR像对的方位向偏移量来选择无明显电离层扰动影响的像对进行干涉处理,获取了高精度、连续的InSAR地表形变场.在此基础上,结合高精度GPS同震形变数据,采用同震、黏弹性松弛震后形变联合反演模型同时确定了汶川地震的同震滑动分布和龙门山地区的流变结构参数.研究结果表明,汶川地震是一个断层破裂非常复杂的地震事件,其中,北川段、岳家山段、虹口段和汉旺段的滑动以逆冲为主,而青川段以右旋走滑为主.滑动主要发生在10 km深度以上的区域,最大滑动量位于虹口段的东北端,达10.7 m.地震释放的总能量为9.28×1020 N·m（Mw7.91）,与地震学的结果一致.联合反演模型确定的龙门山地区中下地壳的黏性系数为2×1018 Pa·s,为青藏高原东部地区的黏性系数提供了一个可靠的下限值.如果有更长时间的震后形变观测时间序列,将为该区域提供更为可靠的流变结构.     

2010年1月12日海地 Mw7.0地震InSAR同震形变观测及同震滑动分布反演

2010年1月12日GMT时间21时53分,在海地境内(72.57°W,18.44°N)发生了Mw7.0地震.文中利用干涉合成孔径雷达(InSAR)方法获得了覆盖整个震区的高精度形变观测资料,用以研究该地震的发震机理.采用ALOs PALSAR数据,分析了轨道、大气等误差源对干涉信号的影响,最终获得了雷达视线向(LOS...     

汶川MS8.0级地震InSAR同震形变场观测与研究

利用InSAR技术,采用地震前后日本ALOS/ PALSAR数据,提取了2008年5月12日四川汶川地震7个条带的地表同震形变场.每个形变条带南北向500 km,东西向70 km,7个形变场覆盖了映秀镇、都江堰、茂县、北川、平武和青川.结果显示,此次地震地表破裂带在北川—映秀断裂带上.地表破裂带从汶川县映秀镇西南震中附近一直到青川县苏河北侧,全长约为230 km.发震断层西北盘为抬升盘,南东盘断层附近,仍然表现为隆起区,显示出以逆冲为主的断层性质.在汶川县映秀镇西侧震中区,最大相对卫星视线向形变量达260 cm,如果全部换算成垂直形变,则两个区域的垂直相对形变达3.3 m.从北川至平通一带,卫星视线向形变范围在120～180 cm的隆起带.其中,擂鼓镇隆起形变范围170～180 cm,换算成垂直形变约在2.2～2.3 m之间.在青川苏河北附近,有70～80 cm范围的隆起形变.在雅安、峨眉山一带,以及射洪至重庆北侧一带有大范围沉降区.在重庆及其南侧区域有幅度在20～30 cm小范围隆起.由青川向东至广元、宁强一带,有形变幅度在60～70 cm的隆起区.整个形变场影响范围较大,四川盆地出现了不同程度的地表形变.