利用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, 与地震学反演的结果比较吻合.      

漾濞MS6.4地震同震形变特征及发震构造探讨

2021年5月21日云南省漾濞县发生MS6.4地震.文中基于升、降轨Sentinel-1 SAR影像,利用InSAR技术获取了此次地震的同震形变场,反演获取了发震断层的精细滑动分布,计算了区域应变分配及同震位错引起的周边各断裂上的库仑应力变化,对发震构造及周边断裂的地震危险性进行了讨论.结果表明:InSAR同震形变场显...     

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

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

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

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

2020年1月19日新疆伽师县MS6.4级地震InSAR同震形变场与断层滑动分布反演

利用欧空局Sentinel-1A SAR数据,重建了2020年1月19日新疆伽师县MS6.4地震InSAR同震形变场.以升降轨InSAR形变场联合约束,反演了发震断层参数与同震滑动分布.结果 表明,地震同震形变场(长轴走向近EW向)发生在柯坪塔格褶皱带与奥兹格尔它乌褶皱带之间的区域内,升降轨观测的视线向(LOS)形变量...     

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

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

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

综合InSAR和应变张量估计2016年MW7.0熊本地震同震三维形变场

利用合成孔径雷达(Synthetic Aperture Radar,SAR)影像提取与地质活动相关的三维地表形变场,对深入理解地质灾害的形成机制及其潜在灾害风险评估非常重要.目前,利用SAR影像的同震三维形变场提取主要利用单个像素点的多次观测构建观测方程,然后基于加权最小二乘(Weighted Least Squares,WLS)方法分解从而获得同震三维形变场,因此该方法缺乏对相邻像素点空间相关性的约束.考虑相邻同震位移点的应力连续性,研究学者提出了顾及大地测量应变张量和卫星形变观测的SAR同震三维形变场方法(Extended Simultaneous and Integrated Strain Tensor Estimation from geodetic and satellite deformation Measurements,ESISTEM).本文以2016年MW7.0熊本地震为例,收集了覆盖此次地震的ALOS-2卫星升降轨影像,利用传统差分InSAR(DInSAR)方法和子孔径雷达干涉测量(Multiple Aperture InSAR,MAI)方法分别对升降轨SAR影像对进行处理,得到视线向(LOS)形变和方位向形变,最后利用ESISTEM方法获取此次地震的三维同震形变场.此外,利用GPS和野外考察观测对本文的三维形变场结果进行结果精度分析.研究结果表明,与传统WLS方法相比,ESISTEM方法不仅能有效抑制奇异像素点对形变结果的干扰,同时对近断层的失相干信号能进行较好的恢复,更有助于解释地表破裂区的地震形变特征和掌握地震发生机制.本文确定的三维同震形变场结果显示主形变区发生在Futagawa断层中部和Hinagu断层最北端,最大水平位移为2m,抬升为0.55m.断层破裂以NE-SW走向的右旋走滑为主兼有部分正断成分.应变张量分析表明发震断层处受到了明显的收缩力和剪切力的作用.     

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

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

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年尼泊尔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同震形变场局部形变细节, 结合震区地质背景、 断裂分布及断层运动特征, 获得了同震破裂拟出露地表迹线。     

基于D-InSAR技术的伊朗巴姆地震地表形变监测

以伊朗巴姆地区为例,对伊朗巴姆地震造成的地表形变进行了差分干涉测量,得到了垂直向的同震三维形变场,并运用GIS三维分析技术对形变场进行了分析。实验结果表明,地震在巴姆城市的东侧造成了较大形变,在西侧也产生了微量形变。巴姆城市北部地块沉降,南部地块隆起。同时在巴姆城市南部可明显看到地震造成的断层。实验结果验证了基于C波段的SAR数据的D-InSAR技术在干燥地区监测地表形变方面的可行性。本文对产生去相关效应的原因进行了解释,认为对于干燥少植被的地区干涉效果较好。并指出,如果能够通过技术进步提高雷达干涉测量的精度并降低观测成本,同时将该技术与GPS、GIS等技术相结合,从而更好地研究形变机理,这将对地质灾害的研究产生重大意义。     

2014年美国加州纳帕M_W6.1地震断层参数的Sentinel-1A InSAR反演

2014年8月24日,在美国加州旧金山海湾北部的纳帕地区发生了MW6.1地震.发震断层是西纳帕断裂系统中的一部分,但是该断层之前并未被足够重视.本文利用欧洲空间局最近发射成功并刚刚投入使用的Sentinel-1A卫星获取的第一对同震干涉像对(20140807-20140831),得到了该地震的地表同震形变场,结合震后24h内区域GPS同震形变资料作为约束条件,反演了纳帕地震的断层几何参数以及滑动分布.Sentinel-1A干涉结果表明,此次地震造成了明显的地面形变,视线向最大抬升和最大沉降量均达到了10cm.联合反演结果表明,该发震断层的走向为344°,倾角为80°.主要破裂以右旋走滑为主,平均倾滑角为-146.5°,最大倾滑量达到了1.1m,位于地表下约4km,存在明显的滑动亏损现象.此次地震,累计释放地震矩达1.5×1018 N·m,约合矩震级MW6.1.该结果略小于InSAR单独约束结果,可能与Sentinel-1A像对中包含的快速震后形变分量有关.     

InSAR analysis of a blind thrust rupture and related active folding: the 1999 Ain Temouchent earthquake (M w 5.7, Algeria) case study

We study the surface deformation associated with the 22 December 1999 earthquake, a moderate sized but damaging event at Ain Temouchent (northwestern Algeria) using Interferometric Satellite Aperture Radar images (InSAR). The mainshock focal mechanism solution indicates reverse faulting with a NE–SW trending rupture comparable to other major seismic events of this section of the Africa–Eurasia plate boundary. Previously, the earthquake fault parameters were, however, poorly known because no aftershocks were precisely determined and no coseismic surface ruptures were observed in the field. Using a pair of ERS data with small baseline and short temporal separation in the ascending orbit we obtained an interferogram that shows the coseismic surface displacement field despite poor coherence. The interferogram measures four fringes and displays an ellipse-shaped lobe with ∼11 cm peak line-of-sight displacement. The elastic modeling using a boundary element method (Poly3Dinv) indicate coseismic slip reaching up to 1 m at 5 km depth on the N 57° E trending, dipping 32° NW Tafna thrust fault. The geodetic estimate of seismic moment is 4.7 × 10¹⁷ N m. (Mw 5.7) in is good agreement with seismological results. The elliptical shape of the surface displacement field coincides with the NE–SW trending Berdani fault-related fold. The consistency between the geological observations and InSAR solution shed light on the precise earthquake location and related Tafna fault parameters.     

联合强震记录和InSAR/GPS结果的四川九寨沟7.0级地震震源滑动分布反演及其地震学应用

2017年8月8日我国四川九寨沟发生里氏7.0级地震.本研究利用基线校正方法获得距震中100km范围内9个强震台站同震位移,基于Sentinel-1卫星干涉SAR影像对获取了InSAR同震形变场.结合GPS形变数据,本研究进行了震源滑动模型联合反演,结果显示此次地震整体以走滑运动为主,释放地震矩约为7.60×1018 N·m(～MW6.52).通过对比模拟形变场和观测值显示,联合反演结果优于单独基于InSAR形变场的反演结果.静态应力变化计算结果显示断层平均静态应力降为1.07MPa.反演滑动模型沿走向和倾角方向拐角波数值分别为0.99×10-4和1.10×10-4.同震静态库仑应力变化计算结果显示共有83.6%的余震位于库仑应力增加的区域,被主震所触发的余震占总数的77.9%,主震对后续余震具有显著触发作用.强地面运动模拟结果显示模拟结果在烈度分布范围和等级方面与调查烈度符合度很高,模拟结果能够很好地反映断层破裂的方向性效应等特征.本研究计算结果显示九寨沟地震无论是平均静态应力降还是拐角波数均低于同类型地震的平均水平,这可能是造成本次地震强地震动水平相对不高的原因.     

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

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

基于升降轨Sentinel-1 SAR数据研究2020年新疆于田MS6.4地震震源机制

2020年6月26日新疆于田县发生M_S6.4地震。 本文先利用Sentinel-1A/B卫星升、 降轨SAR数据, 结合InSAR技术提取此次地震的同震形变场; 并以同震形变场为约束, 利用贝叶斯方法的均匀滑动模型反演发震断层的几何参数; 最后基于梯度下降法(Steepest Decent Method, SDM)来确定发震断层更为精细的滑动分布。 研究结果表明: InSAR形变场的分布主要沿北北东方向, 东西方向跨度约40 km, 南北方向跨度约30 km, 形变特征为西侧沉降与东侧隆升, 西侧最大沉降约200 mm, 东侧最大隆升约70 mm。 于田地震发震断层为一条走向187.68°、 倾角59.78°、 滑动角77.76°的隐伏断层; 于田地震的宏观震源深度为6.29 km, 距震级M_W6.19。 断层滑动分布主要集中在沿走向7~21 km与倾向4~11 km的区域内, 平均滑动量约0.2 m, 在沿倾向7~8 km处的最大滑动量约0.97 m, 同震位错主要表现为正断运动。     

中小地震三维形变场重构方法研究与同震滑动分布反演——以2016年5月22日定日MW5.3地震为例

地震三维形变场对于研究地震发震机制等具有重要意义。 已有的InSAR三维形变场重构研究中, 只有中强地震的实例。 由于形变量级小、 InSAR方位向形变的误差较大, 中小地震的三维形变场重构易受噪声等影响。 本文以2016年5月22日西藏定日M_W5.3地震为例, 开展中小地震三维形变场重构的尝试。 首先基于InSAR技术获取了Sentinel-1升轨和降轨观测模式下的同震形变场, 再结合同震形变场的特点、 区域构造特征等, 添加限定方程(走滑运动为0), 重构了同震三维形变场。 结果显示, 震中附近以下降为主, 幅度达7 cm, 南北方向形变较小(约2 mm), 此区域还伴有2 cm的西向水平运动; 形变中心区域东西两侧部分区域均出现少许东向运动(1.5 cm)。 由同震形变场特征判断此次地震以正断破裂为主。 本文提出了基于连续性分层采样选取样本点方法, 以适应本地震形变场的实际情况。 对所得的LOS向位移场和重构的三维形变场进行降采样, 反演得到了断层面上的滑动分布, 两种数据得到的结果相似, 最优发震断层的走向约181°, 倾角约45°, 断层错动平均滑动角约-87.1°, 平均滑动量约为3.6 cm, 最大滑动量位于深度6.5 km处, 相当于一次M_W5.4的地震。