昆仑山口西MS8.1地震INSAR斜距向同震位错分解

通过INSAR获得的地震同震形变场是地面像元相对卫星雷达波照射方向的视线向(或斜距向)变化量，INSAR视线向变化量分解为水平分量及垂直分量是多解的。为解决INSAR斜距向分解为水平位错及垂直位移的不唯一性，利用GPS定点的现场实测位移观测值，在昆仑山口西8．1级地震破裂带主断面上，建立起一种理论与实测相结合的INSAR斜距向位移分解方程，获得唯一解析解。该方法既保持了现场观测的真实性与精度优势，又利用了INSAR实时获取连续形变场的技术特点，通过一系列近似计算最终获得主破裂带上连续变化的水平及垂直同震形变曲线。     

昆仑山MS81地震震源参数的多破裂段模拟研究

昆仑山MS81地震的已有研究结果在破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数方面存在较大差异.本文采用D-InSAR技术首次获得昆仑山MS81地震干涉同震形变场，结合野外科学考察的实测值，进行了主破裂带InSAR视线向变化量的分解，通过对InSAR分解结果、野外科学考察、遥感解译等多源数据综合分析，重新划分了昆仑山地震的次级破裂段.进而通过对地震南北盘同震应变的分析，发现了昆仑山地震的南北两盘分别受挤压和拉张两种应力作用，研究表明多种岩石在拉张和压力作用下其最小主应力下的杨氏模量表现出非线性弹性特征，从而提出对昆仑山地震地表位移及震源特征参数分析时应考虑非线弹性介质导致的非线性弹性位移分布特征.基于上述原因，本文对Okada线弹性位错模型的算法进行了改进，提出了“多震源、非均一位错分量、多破裂段叠加”的线弹性模型，该模型模拟出的形变场干涉纹图较好地体现了地震形变场的分布特征，并由此获得了一套较为完整的地震发震断层的几何学特征参数，为破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数研究提供了较好的解释.     

昆仑山MS81地震震源参数的多破裂段模拟研究

昆仑山MS81地震的已有研究结果在破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数方面存在较大差异.本文采用D-InSAR技术首次获得昆仑山MS81地震干涉同震形变场，结合野外科学考察的实测值，进行了主破裂带InSAR视线向变化量的分解，通过对InSAR分解结果、野外科学考察、遥感解译等多源数据综合分析，重新划分了昆仑山地震的次级破裂段.进而通过对地震南北盘同震应变的分析，发现了昆仑山地震的南北两盘分别受挤压和拉张两种应力作用，研究表明多种岩石在拉张和压力作用下其最小主应力下的杨氏模量表现出非线性弹性特征，从而提出对昆仑山地震地表位移及震源特征参数分析时应考虑非线弹性介质导致的非线性弹性位移分布特征.基于上述原因，本文对Okada线弹性位错模型的算法进行了改进，提出了“多震源、非均一位错分量、多破裂段叠加”的线弹性模型，该模型模拟出的形变场干涉纹图较好地体现了地震形变场的分布特征，并由此获得了一套较为完整的地震发震断层的几何学特征参数，为破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数研究提供了较好的解释.     

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

2022年青海门源MW6.6地震发震构造——来自InSAR和高分影像约束

2022年1月8日,在青海省门源县海原断裂带冷龙岭断裂西段和托莱山断裂东端发生了M_W6.6地震,发育了明显地震地表破裂带.本文基于欧空局哨兵1号雷达影像,利用InSAR技术获取了门源地震的同震形变场,升降轨InSAR同震视线向位移表现出相反的形变特征,量级达到～60 cm左右,同震形变存在非对称性分布特征,结合高分7号观测解译了近断层地表形变.利用InSAR观测反演获得了发震断层参数及详细滑动分布,计算了同震库仑应力变化,并对发震构造及震中区域未来地震危险性进行了分析讨论.结果表明：门源地震至少有两条断裂发生了破裂,主断层对应冷龙岭断裂西段,InSAR确定的最优断层模型显示主断层东段存在沿走向变化特征,西段则在地质解译断层基础上向西延伸,次断层对应地质解译的托莱山断裂东端,两个断裂组成一个平躺的Y型分布.断层最大滑动量约为～3.7 m,断层浅部存在滑动,表明该地震破裂到了地表,地表破裂长度约19 km.主断层滑动主要集中在0～9 km深度范围;次断层滑动主要集中在0～4 km深度范围,InSAR确定的矩震级为M_W6.6.库仑应力变化结果显示民...     

基于InSAR与GPS观测的汶川同震垂直形变场的获取

发震断层的形变是断层活动的重要参数之一,对认识断层性质、震源机理有重要作用。文中以逆冲性质为主的汶川地震为例,采用符合地表水平形变特征的Biharmonic样条插值对GPS水平形变矢量插值,然后再分解为EW和SN向分量。利用可靠的GPS观测值对InSAR参考点进行校正,统一两者的坐标系。通过对汶川地震视线向形变场剖面与GPS对比分析发现,断层上盘GPS与InSAR观测参考点相差9.93cm,而下盘则为-11.49cm。在此研究基础上,通过GPS水平形变场与InSAR视线向形变场联合解算,获取了汶川地震垂直连续形变场。结果表明,断层两侧垂直形变衰减较快,横跨断裂带形变量30cm的宽度不超过50km;沿发震断层附近垂直形变高值区分布不均匀,主要集中分布在发震断裂的汶川县城至都江堰段、茶坪—北川—南坝段和青川段。这3段各有特色,南段断层两侧垂直形变极不对称,主要以上盘剧烈抬升为主,最大抬升区域在映秀镇至连山坪一带,抬升量达到5.5m。中段表现为较强的反对称性,断层一侧抬升另一侧沉降。该段上盘最大抬升区域在茶坪东侧,抬升量为255cm,下盘最大沉降量在永庆,沉降量为-215cm。北端垂直形变相对较小,主要分布在青川北侧,呈对称分布,在发震断层最北端,最大抬升量为120cm。     

亚像素相位相关法在获取汶川地震近场形变中的应用

震后地表实际破裂带的分布及其近场的形变特征,是理解块体运动学特性、断层破裂特征、地震发生机制等科学问题的十分重要的约束条件。基于InSAR获取的汶川地震同震形变场,由于发震断层附近同震形变梯度巨大,沿断层带出现了非相干条带,以致于无法获得断层附近的形变量。而基于亚像素级的光学影像偏移量法为获取断层附近大形变分布提供了可能。文中以SPOT卫星影像为数据源,采用光学影像偏移量法获得了什邡及茂县地区的水平位移形变场。结果显示龙门山断裂带上至少2条断裂同时发生破裂,形成了主要地表破裂带(龙门山镇-高川破裂带)和次级地表破裂带(汉旺破裂带),沿龙门山镇-高川破裂带平均位移量为4～6m,在高川附近伴随的平均右旋水平位移为1～3m; 汉旺破裂带因逆冲导致水平缩短,平均位移量一般为1～2m。汶川-茂县断裂带没有明显的地表破裂带。研究表明,利用光学影像相位相关法能够获得近断层位错量,可以成为InSAR手段的重要补充。     

基于BFGS法融合InSAR和GPS技术监测地表三维形变

虽然InSAR技术具有高精度、大范围和高空间分辨率等优点,但只能监测雷达视线方向上的一维地表形变;而GPS技术虽可以监测地表的三维形变,但其空间分辨率很低.本文针对融合InSAR和GPS技术监测地表高空间分辨率三维形变展开研究.首先证明了简单的局部最优化迭代算法就能求得综合InSAR和GPS监测地表形变速率的能量函数模型的全局最优估值.随后提出了利用BFGS局部最优算法反演最优的地表三维形变速率.该方法既能避免全局最优化算法计算复杂且难以收敛的问题,又能克服传统的解析法中数值计算不稳定的缺点.最后,通过模拟实验和美国南加州真实数据实验表明,该方法能够得到高精度的地表三维形变速率场.而且当观测或插值误差导致解析法误差较大时,BFGS方法仍能得到高精度、稳定的全局最优解.     

2021年青海玛多7 4级地震的同震变形分析

2021年5月22日,青海省玛多县发生了 7.4级地震,该地震发生在巴颜喀拉地块北部边界东昆仑断裂带以南约70 km,属于块体内部断裂带地震.根据中国大陆构造环境监测网络提供资料,距离震源30多公里的玛多台站记录到东西向永久位移约25 cm.同时,InSAR也观测到明显的形变场,升轨和降轨的最大相对形变量分别约1.87 m和2.32 m.为了解释这些大地测量观测数据,本文利用该地震的三个断层滑动模型,基于不同地球模型的地震位错理论,计算同震变形场,并分别与GNSS观测数据和InSAR视线向形变量对比分析,结果显示基于InSAR数据反演的断层滑动模型产生的位移场与球形地球模型的理论计算结果最为吻合.进一步,利用较优断层模型计算2021年青海玛多7.4级地震的理论同震位移、大地水准面、重力和应变等变化,该结果为玛多地震的GNSS和重力观测的解释提供理论参考依据.     

GPS和InSAR同震形变约束的尼泊尔M_W7.9和M_W7.3地震破裂滑动分布

2015年4月25日尼泊尔爆发MW7.9地震,继而引发5月12日MW7.3级余震,GPS、InSAR监测到震源区及周边大范围同震形变.本文以国内外的GPS和InSAR同震形变为约束,考虑喜马拉雅断裂带岩石圈垂向分层和横向差异的影响,反演主喜马拉雅逆冲断裂在这次主震和余震中破裂面形状和滑动分布.结果显示,主震从USGS确定的震中位置向东偏南延伸100km以上,破裂地面迹线与主前缘逆冲断裂迹线基本一致.破裂面倾角约7°~11°,大部分破裂集中在深度8~20km,同余震分布深度一致.主震最大滑动量约6.0~6.6m,位于14km深处.余震破裂集中在震中附近30km范围内,填补了主震东部破裂空区,最大滑动约3.6~4.6 m,位于13km深.深度20km以下基本没有破裂.地壳介质不均匀性对破裂滑动分布的影响较大,介质不均匀模型的观测值不符值比各向同性弹性半空间模型降低10%以上.本文地震破裂模型特征与地震反射剖面、以及根据震间期大地测量数据反演的喜马拉雅深部蠕滑剖面极其相似.跨喜马拉雅断裂剖面的震间形变量与地震破裂滑移量直接相关.以此推算,尼泊尔中部大震原地复发周期在300年以上.     

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.     

2021年青海玛多MS7.4地震同震应力场与形变台站同震阶变分布关系探讨

利用基于升、降轨InSAR形变场及余震精定位结果反演得到的同震滑动模型,通过PSGRN/PSCMP程序获得同震水平形变场及应力场分布特征,结合玛多M_S7.4地震周边形变同震阶变台站分布特征,探讨同震应力场变化与同震阶变台站分布间的关系。模拟得到的水平形变场结果显示,此次玛多地震为左旋走滑运动特征,水平形变量主要集中在巴颜喀拉块体内,其次是北部的柴达木块体;羌塘块体以及祁连块体同震水平位移量较小;昆仑山口-江错断裂作为一条NE倾向的走滑型断裂,断层上盘区域滑动量明显大于下盘,模拟得到的最大水平形变量达1380mm;形变同震阶变的台站主要集中分布在祁连山断裂带中东段以及西秦岭等地区,祁连山断裂带中东段位于此次玛多地震同震正应力变化正值区域,而西秦岭等地区则处于玛多地震同震剪切应力变化的正值区域,即出现同震阶变的台站与同震应力场变化的正值区域具有较好的一致性。     

2008年汶川M_S8.0地震地表破裂变形定量分析——北川-映秀断裂地表破裂带

2008年汶川MS8.0地震在北川-映秀断裂产生了长达240km的同震地表破裂。通过详细的测量、基于测量标志与断裂变形的几何关系对数据的分析,给出了观测点的断裂同震地表变形的垂直位移、倾向水平缩短、走向滑动、断层上盘水平运动方向等参数。结果显示,断裂同震变形分布的空间变化很大,目前获得的最大水平位移位于虹口乡深溪沟,为4.98m,同时也是最大右旋走滑位移点,走滑量4.5m,而目前获得的最大垂直位移在其东北的支沟,为5.7~6.7m。NE向断裂水平位移多为1~2m,垂直位移多为3m左右,而小鱼洞-草坝分支断裂水平位移和垂直位移都更小,只有0.5~1.5m。擂鼓镇附近的数据则反映与断裂相关的巨型滑坡可能将重力变形叠加到构造变形中。由断层水平缩短和垂直位移计算的断层倾角表明,北川-映秀断裂是浅部陡倾的具有走滑分量的逆断层     

2016年日本熊本地震破裂时空过程联合反演

为了深入认识2016年4月15日日本熊本地震破裂的复杂性,利用远场体波资料和同震InSAR资料联合反演了此次地震的震源破裂时空过程. 联合反演结果表明:熊本地震的震源破裂持续时间约为25 s,整个破裂过程释放的总标量矩为6.03×1019 N·m,对应于矩震级M_W7.1;同震滑动主要集中分布于浅部,破裂以右旋走滑为主,但在沿倾向0—5 km范围内,破裂呈较强的正断特征;此次地震破裂的最大同震滑动量约为4.9 m,且最大同震位错区位于背离断层走向上、距离起始破裂点约5—10 km的区域;破裂前期（0—7 s）,在倾向上向浅表发生破裂,在走向上向东北和西南两侧扩展;大约7 s后,破裂背离断层走向主要向东北方向扩展. 根据有限断层联合反演结果推测,此次熊本地震破裂可能出露至地表.      

提取地震同震形变场特征的最佳雷达波段问题研究

利用弹性形变模型, 并以走滑型断层为例, 对地震同震形变场进行了模拟研究。 结合雷达成像几何关系, 得到干涉雷达的LOS向形变场以及干涉相位条纹图像。 对比分析了不同雷达波段(X波、 C波、 L波)以及不同雷达入射角(30°和50°)的同震形变场和干涉相位图。 结果表明, 同一雷达波段不同雷达入射角的LOS向形变场会有不同的效果, 入射角逐渐增大时, LOS向形变量随之增大, 体现了发震断层走滑特征; 采用不同雷达入射角, 其探测到的地表形变垂向分量和水平分量信息量不同, 表明入射角较大或较小时能够分别较好地对地表形变的水平变化和垂向变化进行监测; 同一雷达入射角, 不同雷达波段的干涉相位条纹有明显区别, 波长越短, 干涉条纹越密集, 容易出现去相关现象, 长波长的干涉相位条纹比较清晰。 从某种意义上讲, 尽管长波长的干涉相位反映地表形变信息的细节会减少, 但受噪声影响较小, 抗干扰性能较短波长的强, 干涉像对相干性高。 地表形变表现为长期缓慢的小幅度渐变性形态, 也表现为短期迅速的大面积突发性形态。 不同雷达波长探测到的地表形变有效信息不同, 为了更好地利用干涉雷达对各种尺度和不同变化速率的地表形变进行监测, 需要在满足精度的条件下尽可能地选择波长较长及多样化的雷达数据。     

THE CORRELATION BETWEEN GEOMETRIC FEATURE OF CO-SEISMIC RUPTURE AND CO-SEISMIC DISPLACEMENT

The existence of asperity has been confirmed by heterogeneously distributed seismic activities along the slipping surface associated with recent huge earthquakes, such as the M8.0 2008 Wenchuan earthquake and M9.0 2011 Tohoku-Oki earthquake. The location of asperity embedded in the seismogenic depth always corresponds to the area of high value of the co-seismic displacement and stress drop where the elastic energy is accumulated during the inter-seismic periods. Fault segmentation is an essential step for seismic hazard assessment. So far, the fault trace is dominantly segmented by considering its geometric features, such as bends and steps. But the connection between the asperity and geometric feature of the slipping surface is under dispute. Research on correlation between geometric feature of surface rupture and co-seismic displacement is of great significance to understand the relationship of seismicity distribution to geometric morphology of sliding surface. To scrutinize the correlation between the geometric feature and co-seismic displacement, we compiled 28 earthquake cases among which there are 19 strike-slip events and 9 dip-slip events. These cases are mainly collected from the published investigation reports and research papers after the earthquake occurred. All the earthquakes' magnitude is between M_W5.4~8.1 except for the M_W5.4 Ernablla earthquake. The range of the rupture length lies between 4.5~426km. Each case contains surface rupture trace mapped in detail with corresponding distribution of co-seismic displacement, but the rupture maps vary in projected coordinate system. So, in order to obtain uniform vector graphics for the following data processing, firstly, vectorization of the surface rupture traces associated with each case should be conducted, and secondly, the vector graphics are transformed into identical geographic coordinate system, i.e. WGS1984-UTM projected coordinate system, and detrended to adjust its fitted trend line into horizontal orientation. The geometric features of surface rupture trace are characterized from three aspects, i.e. strike change, step and roughness. Previous studies about the rupture geometry always describe the characteristics from the whole trace length, consequently, the interior change of the geometric characteristics of the rupture is overlooked. In order to solve this problem, a technique of moving window with a specified window size and moving step is performed to quantify the change of feature values along the fault strike. The selected window size would directly affect the quantified result of the geometric feature. There are two contrary effects, large window size would neglect the detail characteristics of the trace, and small window size would split the continuity of the target object and increase the noise component. So we tested a set of sizes on the Gobi-Altay case to select a proper value and choose 1/25 of the whole rupture length as a proper scaling. Here, we utilize the included angle value of the fitted line in the adjoining windows, Coefficient of variation and the intercept value of the PSD(Power Spectra Density)for characterizing the change of strike, step size and roughness. The rupture trace is extracted within every moving window to calculate the aforementioned feature values. Then we can obtain three sets of data from every rupture trace. The co-seismic displacement is averaged in piecewise with uniform interval and moving step along the fault strike. Then, the correlations between three kinds of feature value and the co-seismic displacement are calculated respectively, as well as the P-value of correlation coefficient significant test. We divided cases into two groups according to the slip mode, i.e. strike-slip group and dip-slip group, and contrast their results. In the correlation result list, there is an apparent discrepancy in correlation values between the two groups. The values of the strike-slip group mostly show negative, which indicates that geometric feature of the rupture trace is in inverse proportion to the displacement. In dip-slip group, the values distribute around zero, which suggests the geometric features is irrelevant to the displacement. Through the analysis of the correlation between the surface rupture and co-seismic displacement, the following conclusions can be reached:1)In comparison with the dip-slip earthquake type, the characteristics of surface rupture of strike-slip earthquakes have a higher-level of correlation with the distribution of the co-seismic displacement, which suggests that the geometric features of strike-slip active faults may have a higher reference value in the fault-segmentation research than the dip-slip type; 2)In most strike-slip events, there is a negative correlation between the geometric features and the co-seismic displacement, which implicates that the higher the feature values of the steps, strike change and roughness, the lower the corresponding co-seismic displacement is; 3)Among the three quantified features of the surface rupture trace, the ranking of relevancy between them and the co-seismic displacement is:step size > strike change > roughness.     

GPS测定的四川九寨沟MW6.6地震同震形变

2017年8月8日四川阿坝州九寨沟发生M_W6.6地震,震源机制解显示该地震为左旋走滑型地震。对震中周围的GPS连续站观测资料进行处理,获得高频GPS动态形变和静态同震水平位移。震中100km范围内四川松潘和甘肃武都站观测到1 Hz动态形变。距离震中约69km的松潘站观测的同震水平位移为7.4mm。根据少量的GPS静态同震位移反演的同震破裂模型显示本次地震的最大滑动量为376mm,地震矩为7.25×1018 N·m,等效矩震级为M_W6.6。正演计算的同震三维形变场显示本次地震的最大水平位移可达4~5cm,垂直位移呈四象限分布,最大可达1.5cm,区域内10个流动GPS站可观测到同震形变。     

Co-seismic ground deformation and source parameters of Mani M7.9 earthquake inferred from spaceborne D-InSAR observation data

In recent ten years, Differential Interferometric Synthetic Aperture Radar (D-InSAR) has become a major technique of space-based geodesy together with GPS, VLBI and SLR. Interferometric Synthetic Aper-ture Radar (InSAR) has many advantages, such as all-weather, all-time, strong stability and dynamic survey property, and no requirements for ground sta-tions. In particular, the surveying results by InSAR can cover a large range of the ground deformation field in succession and has gr…     

Modification of fault slip models of the Mw 9.0 Tohoku Earthquake by far field GPS observations

GPS data from Crustal Movement Observation Network of China (CMONOC) are used to derive far-field co-seismic displacements induced by the M_w 9.0 Tohoku Earthquake. Significant horizontal displacements about 30 mm, 10 mm, and 20 mm were caused by this large event in northeast China, north China, and on the Korean peninsula respectively. Vectors of relatively large horizontal displacements with dominant east components pointed to the epicenter of this earthquake. The east components show an exponential decay with the longitude, which is characteristic of the decay of the co-seismic horizontal displacements associated with earthquakes of thrust rupture. The exponential fit of the east components shows that the influence of the co-seismic displacements can be detected by GPS at a distance of about 3200 km from the epicenter of the earthquake. By considering the capability of the far field displacements for constraining the inversion of the fault slip model of the earthquake, we use spherically stratified Earth models to simulate the co-seismic displacements induced by this event. Using computations and comparisons, we discuss the effects of parameters of layered Earth models on the results of dislocation modeling. Comparisons of the modeled and observed displacements show that far field GPS observations are effective for constraining the fault slip model. The far field horizontal displacements observed by GPS are used to modify the slips and seismic moments of fault slip models. The result of this work is applicable as a reference for other researchers to study seismic source rupture and crustal deformation.