2017年九寨沟MS7.0地震同震地表三维形变场解算研究

利用覆盖九寨沟M_S7.0地震的Sentinel-1升、降轨和Radarsat-2升轨数据,分别提取了3个轨道沿雷达视线向形变,运用多平台联合观测方法解算了九寨沟地震沿地表真实的垂直向、SN向和EW向形变信息。结果表明,3个轨道InSAR数据均监测到了LOS向同震形变,范围约55km×45km,呈“果仁状”,靠近卫星飞行方向最大形变量为12.9cm (降轨),远离卫星飞行方向最大形变量为19.5cm (升轨)。三维形变结果显示,垂直向上位移达23.4cm,垂直向下位移达17.6cm;北向位移达141.8cm,南向位移达100.2cm;东向位移达22.0cm,西向位移达10.7cm。树正断裂两端地块呈非对称水平相对运动,上盘单侧向东形变较剧烈,符合左旋走滑型地震事件的运动特征。     

2016年新疆阿克陶地震Sentinel-1 InSAR同震形变特征

2016年11月25日新疆克孜勒苏州阿克陶县发生M_W6.6地震。 本文利用合成孔径雷达差分干涉测量技术, 对Sentinel-1卫星获取的升、 降轨雷达数据进行了处理, 提取了该次地震的同震形变场, 并结合形变场特征与震源机制解, 采用梯度下降法反演发震断层的滑动分布。 结果表明, 升、 降轨LOS向同震形变场在发震断层两侧具有明显不同的形变特征, 主要形变区域分布在断层两侧, 升轨LOS向形变量可达-8.2 cm与11.2 cm, 降轨LOS向形变量可达-21.4 cm与13.1 cm; 反演的升、 降轨干涉形变场与InSAR测量值之间的残差得到有效控制, 大部分的残差介于±5 cm之间; 断层滑动分布主要集中于沿断层面深约2~18 km处, 最大滑动量位于沿断层面深约7 km处可达0.96 m; 平均滑动角约182.29°, 最大滑动处的滑动角约197.13°, 两个滑动分布中心的滑动角均接近180°, 表明阿克陶地震为一典型的右旋走滑破裂性事件; 当剪切模量取32 Gpa时, 反演的发震断层地震矩M₀可达9.75×10¹⁸, 相当于矩震级M_W6.60, 与地震波形反演结果一致。     

Using Envisat InSAR time-series to investigate the surface kinematics of an active salt extrusion near Qum,Iran

Lower Oligocene-Miocene rock salt is exposed in several diapirs in the central plateau of Iran along the northern margin of the Great Kavir basin. These include a small but mature salt extrusion known as Qum Kuh located near the city of Qum. We use small-baseline interferograms to study the surface displacements of Qum Kuh. The interferometric dataset consisted of 149 interferometric displacement maps derived from 35 C-band ASAR images collected by Envisat satellite in both descending and ascending orbits from 2003 through 2009.The results of the displacement time series analyses show significant long-term trend of horizontal motion toward gentle slopes of Qum Kuh, punctuated by seasonal variations during dry and wet seasons. Interferometric Synthetic Aperture Radar (InSAR) time series results suggest that the salt near the spreading center, where the salt first surfaces from depth, is extruded to the more gentle flanks of Qum Kuh and that down-slope gravitational spreading of the extruded salt accelerates when it is wet.     

Surface deformation analysis in the Ischia Island (Italy) based on spaceborne radar interferometry

The island of Ischia, located at the northwestern end of the Gulf of Napoli (Italy), is a volcanic area that is historically active (the Arso eruption, in 1302 and the Casamicciola earthquake, in 1883) and has diffuse hydrothermal phenomena. We present in this work a study of the surface deformation occurring in the island, which is based on applying the Differential Synthetic Aperture Radar Interferometry (DInSAR) algorithm referred to as Small BAseline Subset (SBAS) technique. This study is focused on the 1992–2003 time interval and SAR data acquired by the European Remote Sensing (ERS) satellites from ascending and descending orbits have been used, thus allowing us to discriminate the vertical and east–west components of the displacements. A validation of the DInSAR results has been carried out first by comparing the vertical deformations estimated from the SAR data with those measured from the spirit leveling network that is present in the area. In particular, we computed the difference between the mean vertical deformation velocities estimated from the SAR and the corresponding geodetic measurements along three main leveling lines; the maximum value of the root mean square difference is of about 1 mm/yr. The final discussion is dedicated to the interpretation of the detected displacements, benefiting from the overall information extracted from the ascending and descending DInSAR measurements. In particular, DInSAR data relative to the vertical deformation component show that the present-day subsidence of Ischia mainly develops in areas characterized by active landsliding and along faults; moreover, the deflation of the island, which is recorded by the horizontal displacement component, is probably related to the de-pressurization of the hydrothermal system.     

Comparison of geomorphological field mapping and 2D‐InSAR mapping of periglacial landscape activity at Nordnesfjellet,northern Norway

The ability to continuously monitor the dynamic response of periglacial landforms in a climate change context is of increasing scientific interest. Satellite radar interferometry provides information on surface displacement that can be related to periglacial processes. Here we present a comparison of two‐dimensional (2D) surface displacement rates and geomorphological mapping at periglacial landform and sediment scale from the mountain Nordnesfjellet in northern Norway. Hence, 2D Interferometric Synthetic Aperture Radar (InSAR) results stem from a 2009–2014 TerraSAR‐X dataset from ascending and descending orbits, decomposed into horizontal displacement vectors along an east–west plane, vertical displacement vectors and combined displacement velocity. Geomorphological mapping was carried out on aerial imagery and validated in the field. This detailed landform and sediment type mapping revealed an altitudinal distribution dominated by, weathered bedrock blockfields, surrounded primarily by slightly, to non‐vegetated solifluction landforms at the mountain tops. Below, an active rockslide and associated rockfall deposits are located on the steep east‐facing side of the study area, whereas glacial sediments dominate on the gentler western side. We show that 2D InSAR correctly depicts displacement rates that can be associated with typical deformation patterns for flat‐lying or inclined landforms, within and below the regional permafrost limit, for both wet and dry areas. A net lowering of the entire landscape caused by general denudation of the periglacial landforms and sediments is here quantified for the first time using radar remote sensing. Copyright © 2018 John Wiley & Sons, Ltd.     

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和应变张量估计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走向的右旋走滑为主兼有部分正断成分.应变张量分析表明发震断层处受到了明显的收缩力和剪切力的作用.     

基于Sentinel-1A数据反演漾濞MS6.4地震的同震形变场及断层几何参数

2021年5月21日晚21时48分,云南省大理州漾濞县(震中:25.67°N,99.87°E)发生M_S6.4地震,震源深度8 km。为快速获得此次地震同震形变场及断层几何参数,研究该次地震的发震构造等,文章基于震前、震后的sentinel-1A卫星升降轨SAR数据进行二轨法差分雷达干涉测量(DInSAR),并基于Okada弹性半空间位错模型反演断层几何参数。研究结果如下:(1)此次地震造成的同震形变场长约19 km,宽约20 km;(2)升轨雷达视线向最大形变约为8.2 cm,降轨雷达视线向最大形变约为8.7 cm;(3)地震断层走向为313.7°,倾角为87°,滑动角为175°,为右旋走滑型断层,最大滑动量为0.79 m,反演得出的地震矩为1.48×10~(18) N·m,矩震级为M_W6.1。在川滇块体向南挤出的构造背景下,块体西边界的维西—乔后断裂、红河断裂发生右旋走滑,本次地震便是维西—乔后断裂南段分支断裂右旋走滑活动的体现。     

利用考虑地形起伏的反演方法获取2017年伊拉克MW7.3地震的断层参数及滑动分布

2017年11月13日伊拉克北部地区苏莱曼尼亚省发生了M_W7.3地震, 造成了重大的人员和经济损失。 本文利用升降轨的Sentinel-1和降轨的ALOS-2卫星的SAR数据, 通过差分干涉测量技术获取了该地震的同震形变场, 联合DInSAR和MAI技术, 采用抗差最小二乘法求解该地震的同震三维形变场。 基于改进的考虑地形起伏的均匀位错模型反演确定了发震断层的断层参数, 最后基于非均匀位错模型得到了发震断层的分布式滑动分布模型。 结果显示: ALOS-2卫星降轨轨道观测到的伊拉克地震引起的LOS向地表形变最大为55.8 cm抬升和47.9 cm下沉; Sentinel-1卫星观测到的伊拉克地震引起的LOS向地表形变为: 升轨轨道最大为87.9 cm抬升和17.1 cm下沉; 降轨轨道最大为55.6 cm抬升和38 cm下沉; 相对于前人的研究成果, 本文利用改进的考虑地形起伏的反演方法得到的发震断层几何参数表明发震断层为NNW走向, 倾向角为352°, 同震破裂以逆冲为主, 同时兼有一定的左旋走滑分量。 基于均匀位错模型反演得到的断层滑动分布结果表明, 同震破裂未延伸至地表, 主要滑动量集中在12~18 km, 最大滑动量位于15 km深度, 达到4.3m, 反演得到的矩震级为M_W7.35, 与UGSG、 GCMT等机构给出的结果一致。     

Some problems of landslide monitoring using satellite radar imagery with different wavelengths: Case study of two landslides in the region of Greater Sochi

The problems of processing and interpreting the data provided by radar satellite interferometry for the conditions of landslides covered by vegetation are analyzed in two case studies of landslides in the Northern Caucasus in the region of Kepsha and Mamaika villages in the vicinity of the railway tunnels. The estimates of the displacement fields are obtained by the method of persistent scatterers using the StaMPS program package. The five-year experience of landslide monitoring shows that in the unfavorable conditions of satellite radar interferometry, proper selection of the strategy of satellite image processing is vital. In the present paper, we discuss, in particular, the crop selection, the selection of the master image, reference area, and digital elevation model. For the landslide located in the sparsely populated region near Kepsha village, we used the data from the ascending and descending tracks of the long-wavelength ALOS and shorter-wavelength ENVISAT satellites. For the landslide in the region of Mamaika village with a large number of different buildings serving as good scatterers for radar signals, we used the images from the ENVISAT and from TerraSAR satellite, which transmits even shorter waves. The average line-of-sight (LOS) displacement velocities V _LOS for the landslide near Kepsha village measure at most 10 mm per annum, which means that this landslide has remained stable at least since 2004. The landslide in Mamaika village is significantly more active. The average LOS displacement velocities in the active part of this landslide attain 60 mm per annum. The artificial corner reflector installed on the segment of the landslide where natural scatterers of radar signal are absent made it possible to estimate the LOS displacement velocity on this segment of the slope at 49 mm per annum.     

Decomposing InSAR LOS displacement into co-seismic dislocation with a linear interpolation model： A case study of the Kunlun Mountain Ms=8.1 earthquake

It has always been a difficult problem to extract horizontal and vertical displacement components from the InSAR LOS （Line of Sight） displacement since the advent of monitoring ground surface deformation with InSAR technique. Having tried to fit the firsthand field investigation data with a least squares model and obtained a preliminary result, this paper, based on the previous field data and the InSAR data, presents a linear cubic interpolation model which well fits the feature of earthquake fracture zone. This model inherits the precision of investigation data; moreover make use of some advantages of the InSAR technique, such as quasi-real time observation, continuous recording and all-weather measurement. Accordingly, by means of the model this paper presents a method to decompose the InSAR slant range co-seismic displacement （i.e. LOS change） into horizontal and vertical displacement components. Approaching the real motion step by step, finally a serial of curves representing the co-seismic horizontal and vertical displacement component along the main earthquake fracture zone are approximately obtained.     

2020年1月19日新疆伽师MS6.4地震InSAR同震形变场特征及发震构造初步探讨

基于Sentinel-1 SAR升、降影像,利用D-InSAR技术获取新疆伽师M S6.4地震的同震形变场,结果表明,本次地震引起的同震形变场整体呈近椭圆状分布,形变区东西长约66 km,南北宽约40 km,整个形变场由南部隆升区和北部沉降区组成,南部最大隆升量约7 cm,北部最大沉降量约3 cm。本次地震发生在块体俯冲界面处的低倾角逆冲推覆构造带上,隆升和沉降两个中心均位于逆冲推覆体的上盘,形变主要以隆升形变为主,符合低倾角逆断层中强震的变形特征。在沉降区与隆升区之间干涉条纹连续分布,未出现表征地表破裂位置的空间失相关带,表明地震未引起明显的地表破裂。结合震源机制、余震精定位及区域构造特征,初步推断认为伽师地震的发震构造可能为柯坪塔格推覆构造前缘的N倾的柯坪断裂。     

InSAR数据约束的2020年西藏定日Mw5.7地震源模型及构造意义

2020年3月20日,处在北喜马拉雅断裂与申扎一定结断裂交会区的西藏定日发生Mw5.7地震,此次地震发震断层源模型的研究对于认识该地区的复杂地质构造具有重要意义.本文利用升、降轨Sentinel-1A星载SAR数据,通过合成孔径雷达干涉测量(InSAR)处理获得了地震引起的地表位移场,其中,雷达视线方向(LOS)最大地...     

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数据自适应融合形变测量

近年来,合成孔径雷达干涉测量(Interferometric Synthetic Aperture Radar,InSAR)技术在地面沉降监测方面展现了巨大的应用潜力,但受其重访周期和一维形变测量能力的限制,仅利用单一轨道卫星观测数据很难揭示真实的地表形变特征及其演化规律.随着在轨运行的SAR卫星系统不断增加,使得融合相同时间段内覆盖同一区域的多源多轨道InSAR数据成为可能.然而目前普遍采用的多源InSAR数据融合方法均为针对大尺度形变监测设计,或者忽略南北向形变甚至水平形变,容易造成误判.为此,本文对经典小基线集(Small Baseline Subset, SBAS)时序InSAR分析方法进行改进,在其形变反演模型中加入东西向和南北向形变参数,采用方差分量估计方法解算多源观测数据验后方差,通过迭代精化确定权重矩阵,从而获得形变参数的最优估值.使用美国南加州地区的ALOS PALSAR和ENVISAT ASAR数据开展实验,利用南加州综合GPS网(SCIGN)位于研究区域内的9个站点观测数据进行验证,结果表明本文方法得到的融合形变测量结果在垂直向上能够准确反映地表形变波动,周期性与GPS观测比较一致;同时,融合得到的三维形变场显示南加州洛杉矶地区存在不可忽略的水平形变,东西向形变测量精度略高于南北向.因此,基于方差分量估计的多源InSAR融合方法在提高形变测量时间序列连续性的同时,能够更准确地反演研究区域三维形变特征.     

InSAR干涉形变场远程模拟系统实现与地震地表形变场模拟算例分析

利用InSAR(Interferometric SAR)干涉测量技术可以获得地表形变场(视线向, LOS), 将其与弹性半空间中一定断层模型模拟计算出的地表形变场进行正演分析, 以便获取发震断层的几何学和运动学参数, 是一种典型的前向模拟模式。 该计算模式由于模拟计算程序调试、 安装过程的复杂性, 使得在每台微机上安装此类程序不仅费时费力, 且浪费大量的计算资源。 文中首先介绍了实现InSAR干涉形变场模拟模型由单机推广到Internet/Intranet上进行远程计算最终建立“InSAR干涉形变场远程模拟系统”的过程, 并基于此系统平台对1997年玛尼M_S7.9地震的InSAR地表形变场进行了模拟试算、分析等。     

InSAR数据约束的2022年1月8日青海门源MS6.9地震发震构造研究

利用Sentinel-1A卫星升降轨道数据和D-InSAR技术获得青海门源2022年1月8日M_S6.9地震的同震形变场,并基于弹性半空间位错模型反演其震源参数,利用分布滑动模型确定断层面上的滑动分布。结果表明,2022年1月8日青海门源地震的同震形变场沿NWW-SEE方向分布;断裂带南缘升轨影像和降轨影像最大视距分别为61 cm和62 cm,断裂带北缘升轨影像和降轨影像最大视距地表形变量分别为43 cm和56 cm。InSAR同震形变场断裂尺度模型断层长30 km,宽18 km,最大滑移量3.5 m;断层滑动分布模型表明该地震为左旋走滑地震。结合冷龙岭断裂的运动特征和几何特征,初步确定此次M_S6.9地震的发震断裂为冷龙岭断裂     

STUDY ON THE SEISMOGENIC FAULT CHARACTERSTICS OF 2016 MW5.9 MENYUAN EARTHQUAKE BASED ON Sentinel -1A DATA

In this paper, we processed and analyzed the Sentinel-1A data by "two-pass" method and acquired the surface deformation fields of Menyuan earthquake. The results show the deformation occurred mainly in the south wall of fault, where uplift deformation is dominant. The uplift deformation is significantly larger than the subsidence and the maximum uplift of ascending and descending in the LOS is 6cm, 8cm respectively. Meanwhile, based on the Okada model, we use the ascending and descending passes data as constraints to invert jointly the fault distribution and source parameters through constructing fault model of different dip directions. The optimum fault parameters are:The dip is 43°, the strike is 128°with the mean rake of 85°. The maximum slip is about 0.27m. The inverted seismic moment M₀ is 1.13×10¹⁸N·m, and the moment magnitude M_W is 5.9. The SW-dipping Minyue-Damaying Fault is possibly the seismogenic fault, based on the comprehensive analysis of the focal mechanisms, aftershocks relocation results and the regional tectonic background. The focus property is dominated by thrust movement with a small amount of dextral strike-slip component. The earthquake is the result of local stress adjustment nearby the Lenglongling Fault under the background of northeastward push and growth of Tibet Plateau.     

利用PSInSAR时间序列研究拉奎拉地震位移场变化特征

对2009年4月6日意大利拉奎拉(LAquila)Mw6.3级地震9景ENVISAT重轨单视复数据采用双轨模式进行永久散射体(Persistent Scatterer,PS)处理,获得了此次地震时序位移场,结合PS位移场Delaunay算法数值分析,结果表明:1)Envisat干涉雷达完整而清晰地探测到了地震前后位移场变化过程及其在不同阶段与震源断裂相关的不同形变特征:震前蠕动位移-明显变形-震期快速突变-震后量级明显减缓的持续变形;2)在此次分析的约54×59km2区域内,震中以西为大面积上升区,视线向最大上升量为130mm;震中以东为下降区域,但下降集中在破裂区,主要在发震过程及震后形成,视线向最大下沉量为210mm,与GPS监测结果一致;3)地震产生的破裂主要分布在约22×14km2范围内,沿NW走向、SW倾斜的Paganica-S.Demetrio正断层展布,方向约135°;4)本文为探索利用PSInSAR时间序列进行地震趋势预测提供了一个较好的案例.     

Analysis of Ground Deformation Detected Using the SBAS-DInSAR Technique in Umbria, Central Italy

Ground deformation affecting the Umbria region (central Italy) in the 9-year period from 1992 to 2000 was investigated through multi-temporal Differential Synthetic Aperture Radar Interferometry (DInSAR). For the purpose, the Small BAseline Subset (SBAS) technique was adopted, which allows studying the temporal evolution of the detected deformation at two spatial scales: a low-resolution (regional) scale, and a full-resolution (local) scale. For the analysis, SAR data acquired by the European Remote Sensing (ERS-1/2) satellites along ascending and descending orbits were used. The detected deformation was analysed to investigate its relevance to geophysical, geomorphologic, and human-induced processes that may result in hazardous conditions to the population of Umbria. Low-resolution deformation data were used to: (i) determine the amount of displacement caused by the Umbria-Marche earthquake sequence from September 1997 to April 1998 in the Foligno area, (ii) determine the number and percentage of the known landslides that can be monitored by the DInSAR technology in the investigated area, and (iii) identify and measure subsidence induced by exploitation of a confined aquifer in the Valle Umbra. Results indicate that earthquakes moved through the Foligno area westwards up to 3.9 cm and with an uplift reaching 1.7 cm. Intersection in a GIS of the low-resolution deformation maps with a detailed landslide inventory map allowed the determination that the portion of landslides that can be monitored by the SBAS-DInSAR technique in Umbria ranges from 2.7% to 3.4%, and the percentage of the total landslide area ranges from 10.4% to 12.8%. In the Valle Umbra, a dependency was found between the time and the amount of detected ground deformation, and the record of water withdrawal. The full-resolution deformation data were used to investigate the movement of the Ivancich landslide, in the Assisi Municipality. Joint analysis of the spatial and the temporal characteristics of the ground displacement allowed the formulation of a hypothesis on the landslide geometry and deformation pattern.