Mapping three-dimensional co-seismic surface deformations associated with the 2015 M_W7.2 Murghab earthquake based on InSAR and characteristics of crustal strain

Three-dimensional(3 D) co-seismic surface deformations are of great importance to interpret the characteristics of coseismic deformations and to understand the geometries and dynamics of seismogenic faults. In this paper, we propose a method for mapping 3 D co-seismic deformations based on InSAR observations and crustal strain characteristics. In addition, the search strategy of correlation points is optimized by adaptive correlation distance, which greatly improves the applicability of the proposed method in restoring deformations in decorrelation areas. Results of the simulation experiment reveal that the proposed method is superior to conventional methods in both the accuracy and completeness. The proposed method is then applied to map the 3 D co-seismic surface deformations associated with the 2015 MW7.2 Murghab earthquake using ascending and descending ALOS-2 PALSAR-2 images. The results show that the seismogenic fault is the Sarez-Karakul fault(SKF), which is dominated by NE-SW strike slips with an almost vertical dip angle. The north section and the south segment near the epicentre have obvious subsidence along with a southwestward motion in the northwest wall, and the southeast wall has northeast movement and surface uplift trend along the fault zone. The strain field of the earthquake is also obtained by the proposed method. It is found that the crustal block of the seismic area is obviously affected by dilatation and shear forces, which is in good agreement with the movement character of the sinistral slip.     

Preliminary analysis on the source properties and seismogenic structure of the 2017 <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>7.0 Jiuzhaigou earthquake

At GMT time 13:19, August 8, 2017, an M_s7.0 earthquake struck the Jiuzhaigou region in Sichuan Province, China, causing severe damages and casualties. To investigate the source properties, seismogenic structures, and seismic hazards, we systematically analyzed the tectonic environment, crustal velocity structure in the source region, source parameters and rupture process, Coulomb failure stress changes, and 3-D features of the rupture plane of the Jiuzhaigou earthquake. Our results indicate the following: (1) The Jiuzhaigou earthquake occurred on an unmarked fault belonging to the transition zone of the east Kunlun fault system and is located northwest of the Huya fault. (2) Both the mainshock and aftershock rupture zones are located in a region where crustal seismic velocity changes dramatically. Southeast to the source region, shear wave velocity at the middle to lower crust is significantly low, but it rapidly increases northeastward and lies close to the background velocity across the rupture fault. (3) The aftershock zone is narrow and distributes along the northwest-southeast trend, and most aftershocks occur within a depth range of 5–20 km. (4) The focal mechanism of the Jiuzhaigou earthquake indicates a left-lateral strike-slip fault, with strike, dip, and rake angles of 152°, 74° and 8°, respectively. The hypocenter depth measures 20 km, whereas the centroid depth is about 6 km. The co-seismic rupture mainly concentrates at depths of 3–13 km, with a moment magnitude (M_w) of 6.5. (5) The co-seismic rupture also strengthens the Coulomb failure stress at the two ends of the rupture fault and the east segment of the Tazang fault. Aftershocks relocation results together with geological surveys indicate that the causative fault is a near vertical fault with notable spatial variations: dip angle varies within 66°–89° from northwest to southeast and the average dip angle measures ~84°. The results of this work are of fundamental importance for further studies on the source characteristics, tectonic environment, and seismic hazard evaluation of the Jiuzhaigou earthquake.     

Rupture behavior and deformation localization of the Kunlunshan earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>7.8) and their tectonic implications

Earthquake surface rupture is the result of transformation from crustal elastic strain accumulation to permanent tectonic deformation. The surface rupture zone produced by the 2001 Kunlunshan earthquake (M _w 7.8) on the Kusaihu segment of the Kunlun fault extends over 426 km. It consists of three relatively independent surface rupture sections: the western strike-slip section, the middle transtensional section and the eastern strike-slip section. Hence this implies that the Kunlunshan earthquake is composed of three earthquake rupturing events, i.e. the M _w =6.8, M _w =6.2 and M _w ⩽=7.8 events, respectively. The M _w =7.8 earthquake, along the eastern section, is the main shock of the Kunlunshan earthquake, further decomposed into four rupturing subevents. Field measurements indicate that the width of a single surface break on different sections ranges from several meters to 15 m, with a maximum value of less than 30 m. The width of the surface rupture zone that consists of en echelon breaks depends on its geometric structures, especially the stepover width of the secondary surface rupture zones in en echelon, displaying a basic feature of deformation localization. Consistency between the Quaternary geologic slip rate, the GPS-monitored strain rate and the localization of the surface ruptures of the 2001 Kunlunshan earthquake may indicate that the tectonic deformation between the Bayan Har block and Qilian-Qaidam block in the northern Tibetan Plateau is characterized by strike-slip faulting along the limited width of the Kunlun fault, while the blocks themselves on both sides of the Kunlun fault are characterized by block motion. The localization of earthquake surface rupture zone is of great significance to determine the width of the fault-surface-rupture hazard zone, along which direct destruction will be caused by co-seismic surface rupturing along a strike-slip fault, that should be considered before the major engineering project, residental buildings and life line construction. Supported by the National Natural Science Foundation of China (Grant No. 40474037) and the National Basic Research Program of China (Grant No. 2004CB418401)     

GPS time-series and its response to <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 Kunlunshan earthquake

In this paper, observation data in 25 GPS reference stations of China have been analyzed by calculating GPS position coordinate time-series with GIPSY. Result shows there is an obvious trend variation in such time-series. The trend variations of time series along the longitude and latitude coordinate reflect the motion of each position in the global-plate, in which the trend variation in the vertical direction reveals some large-scale construction information or reflects the local movement around the positions. The analysis also shows that such time-series have a variation cycle of nearly 1.02 a, but the reason still remains to be further studied. At the end of this paper, response of the time-series of M _S=8.1 Kunlunshan earthquake was analyzed, and the seismogenic process of M _S=8.1 Kunlunshan earthquake, according to the time proceeding and the feature of anomaly, was divided into 3 phases—changes in blocks with forces, strain accumulation, quick accumulation and slow release of energy. At the initial stage of seismogenic process of M _S=8.1 earthquake and at the imminent earthquake, coseismic process as well as during the post earthquake recovery, anomaly in vertical direction is always in a majority. The anomalous movement in vertical direction at the initial stage resulted in a blocking between faults, while at the middle stage of seismogenic process, the differential movement between blocks are in a majority, which is the major reason causing energy accumulating at the blocking stage of faults.     

Preliminary analysis on characteristics of coseismic deformation associated with <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 western Kunlunshan Pass earthquake in 2001

Based on the analysis of coseismic deformation in the macroscopic epicentral region extracted by Differential Interferometric Synthetic Aperture Radar (D-InSAR), and combined with the seismic activity, focal mechanism solutions of the earthquake and field investigation, the characteristic of coseismic deformation of M _S=8.1 western Kunlunshan Pass earthquake in 2001 was researched. The study shows that its epicenter lies in the northeast side of Hoh Sai Hu; and the seismogenic fault in the macroscopic epicentral region can be divided into two central deformation fields: the west and east segments with the lengths of 42 km and 48 km, respectively. The whole fault extends about 90 km. From the distribution of interferometry fringes, the characteristic of sinistral strike slip of seismogenic fault can be identified clearly. The deformations on both sides of the fault are different with an obviously higher value on the south side. In the vicinity of macroscopic epicenter, the maximum displacement in look direction is about 288.4 cm and the minimum is 224.0 cm; the maximum sinistral horizontal dislocation of seismogenic fault near the macroscopic epicenter is 738.1 cm and the minimum is 551.8 cm.     

2022年青海门源6.9级地震引起的地表同震位移场研究

中国青海省门源县于2022年1月8日发生6.9级地震。依据该地震震源断层信息设置4种不同的位错分布模式,基于Okada提出的地表位移解析解分别计算4种模式下地表同震位移场,结合现场观测数据,探讨发震断层的滑动形式及其对周边地表产生的影响。结果表明,此次地震发震断裂初步判断主要为冷龙岭断裂西侧延伸至托莱山断裂,以左旋走滑断层为主,断层面上最大位错量达到4 m左右;震中西南侧向NE方向运动,东南侧向SE方向运动,西北侧和东北侧分别向NW以及SW方向运动;震中附近小范围区域产生了超过1.5 m的地表水平位移,破裂带上存在竖向地表位移超过0.5 m的区域;现场监测到局部产生最大约2.1～2.3 m的水平位错,以及部分区段垂直位错量最大达到0.7 m;以震中位置为中心,断层引起的地表位移影响范围达到约30 km×36 km,此范围内产生的地表位移大于0.1m。研究为此次地震的震后恢复工作以及此区域后续的工程设防等提供参考。     

基于InSAR的青海大柴旦地震三维同震形变场获取与震源特征分析

本文利用Envisat ASAR的升、降轨和宽幅数据,通过基于先验知识的最小二乘迭代逼近获取大柴旦2次地震的地表三维同震形变.结果表明,2008年MW6.3地震垂直向形变主要发生在断层南盘,以隆升形变为主,最大隆升量约10cm,北盘沉降量小于等于-1cm.东西向形变在南盘呈向东运动的特征,最大运动量约4cm,北盘向西运动,最大运动量约为-2cm.2009年MW5.8地震垂直向形变显示断层南盘抬升的特征,最大抬升量约27cm,北盘最大沉降量约-3cm.东西向形变表现为南盘向东运动,最大约10cm,北盘向西运动,约为-4cm.可以看出这两次地震均表现为逆冲为主,兼少量左旋走滑的震源特征.视线向结果无法判定同震形变的少量走滑特征,而地表三维分量可以有效地识别出少量左旋还是右旋走滑的震源特性.本文以视线向、垂直向、东西向形变量作为约束条件,利用Okada模型正演了2008年地震同震三维形变场.结果显示,采用逆冲兼少量左旋走滑的发震断层参数,视线向、垂直向、东西向正演结果与观测结果吻合.这也表明采用分解后的地表三维同震形变场可以有效地识别出发震断层的少量左旋走滑特征.     

2021年甘肃阿克塞5.5级地震前区域地壳变形异常特征分析

为研究2021年甘肃阿克塞5.5级地震前震中周边的地壳变形运动特征,根据CMONOC提供的GPS观测资料,通过GPS速度场、面应变率场和站间基线时间序列分析,探讨了震前的地壳变形特征。研究结果表明,震中处于发震断层西南侧远场速度大、近场速度小且断裂两侧速度场方向明显不同的位置,同时位于与发震构造类型一致的面应变率场的压缩高值区,即挤压应变累积区。跨发震断裂的基线长度在震前3年多时间出现的转折变化,表明震前发震断裂两侧左旋走滑速率明显减缓;而断裂两侧的挤压速率有所加快则反映出本次地震前震中区域存在一定程度的应变能积累。研究结果有助于认识该地震的孕震过程,并可为相似构造类型区域的强震预测提供参考依据。     

Constraints on rupture speed of the 2001 <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript> 8.1 West Kunlun Mountain Pass earthquake by co-seismic surface rupture slip displacements based on the slip-weakening mechanism with frictional undershoot involved

With co-seismic surface rupture slip displacements provided by the field observation for the 2001 M_S8.1 West Kunlun Mountain Pass earthquake, this paper estimates the rupture speed on the main faulting segment with a long straight fault trace on the surface based on a simple slip-weakening rupture model, in which the frictional overshoot or undershoot are involved in consideration of energy partition during the earthquake faulting. In contrast to the study of Bouchon and Vallée, in which the rupture propagation along the main fault could exceed the local shear-wave speed, perhaps reach the P-wave speed on a certain section of fault, our results show that, under a slip-weakening assumption combined with a frictional undershoot (partial stress drop model), average rupture speed should be equal to or less than the Rayleigh wave speed with a high seismic radiation efficiency, which is consistent with the result derived by waveform inversion and the result estimated from source stress field. Associated with the surface rupture mechanism, such as partial stress drop (frictional undershoot) associated with the apparent stress, an alternative rupture mechanism based on the slip-weakening model has also been discussed.     

2016年11月13日新西兰凯库拉(Kaikoura)MW7.8地震同震位移和应力数值分析

2016年11月13日新西兰南岛北端凯库拉（Kaikoura）发生了M_W7.8大地震,造成了强烈的地表变形并引发大面积滑坡和海啸的发生.基于美国地质调查局（USGS）断层滑动模型,建立全球同震横向不均匀并行椭球型地球模型,计算了此次新西兰凯库拉大地震产生的同震形变和应力及库仑应力变化.初步计算结果表明：新西兰凯库拉M_W7.8地震造成断层上盘东北向抬升,下盘西南俯冲;引起发震区域同震位移较大,从凯库拉到坎贝拉（Campbell）以及首都惠灵顿（Wellington）整体上东北向抬升,最大同震水平位移1.2 m,垂直位移1.1 m.此次大地震释放了发震断层上积累的压应力,但增加了发震断层两端的挤压力;同时,同震剪应力变化增加了NE-SW向断层发生右旋滑动的危险性;采用此次地震发震断层参数计算得出的最大库仑应力变化增加区域集中在发震断层两端,可达到MPa量级.当分别采用新西兰北岛Awatere断裂系和南岛Wellington断裂系参数计算库仑应力变化时,发现新西兰北岛和南岛震中以南区域的库仑应力均增加,可触发部分余震的发生.     

利用考虑地形起伏的反演方法获取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等机构给出的结果一致。     

2010年玉树MS7.1地震同震破裂、余震分布特征及其与构造的关系

针对2010年4月14日玉树发生M_S7.1地震,本文利用InSAR数据给出同震视线向位移确定出的发震断层空间展布,并以该同震位移为约束反演得到主震和最大余震的同震位错分布.结果表明,主震同震位错发生在东玉树断裂,最大余震同震位错发生在西玉树断裂东端;基于位错分布计算了同震库仑应力变化与西部余震集中区地震活动之间的关系,结果反映玉树主震后最为活跃的余震活动可能受控于近东西向的次级断层(走向约为85°),而非玉树主干断裂;玉树断裂带整体呈现为左旋走滑运动,但其具体断层运动形式表现出主干断裂典型走滑运动、走滑断裂间的拉张和逆冲性质的次级运动、次级断裂与主干断裂相互作用下更为复杂的多方向次级断层活动等等不同变形特征,而主震同震破裂与余震空间分布均与这些不同断层变形方式有着密切关系.     

Analysis of Anomalous Characteristics of Regional Crustal Deformation before the September 16th, 2019 Zhangye MS5.0 Earthquake, Gansu, China

To study the crustal movement in the vicinity of the epicenter before the Zhangye M_S5.0 earthquake in 2019, the characteristics of crustal deformation before the earthquake are discussed through the GPS velocity field analysis based on the CMONOC data observed from GPS. The baseline time series between two continuous GPS stations and the strain time series of an area among several stations are analyzed in the epicenter area. The resulting time series of baseline azimuth around the epicenter reflects that the energy of the fault in the northern margin of Qilian Mountain is accumulated continuously before 2017. Besides, the movement trend of azimuth slows down after 2017, indicating the stress accumulation on both sides of the seismogenic fault zone has reached a certain degree. The first shear strain and EW-direction linear strain in the epicentral area of the Zhangye M_S5.0 earthquake remain steady after 2017, and the surface strain rate decreases gradually after 2016. It is illustrated that there is an obvious deformation loss at the epicentral region three years before the earthquake, indicating that a certain degree of strain energy is accumulated in this area before the earthquake.     

Regional tectonic deformation setting before the M<Subscript>s</Subscript>8.1 earthquake in the west of the Kunlun Mountains Pass

This paper gives a preliminarily study of the regional tectonic deformation setting before the M_s8.1 earthquake that occurred in the west of the Kunlun Mountains Pass; in the study, the data of the velocity field of crustal horizontal movement during 1991-2000 observed by GPS in and around the Qinghai-Tibet block and those of gravity reiteration in 1998 and 2000 were used. Analysis shows that the preparation and occurrence of this large earthquake are related to the horizontal movement and deformation setting in a large region and might be attributed to the block activity on a relatively large scale. Within the Qinghai-Tibet block, the region of left-lateral shear deformation is of a very large extent. This large earthquake occurred right in such a place where the left-lateral shear strain along the fault strike had the highest rate and the planar dilatation strain was tensile, which was on the margin of negative value region of abnormal gravity variation. The regional tectonic deformation setting can help the huge left-lateral strike-slip rupture to develop.     

关于芦山7.0级地震在龙门一带是否存在同震地表破裂的讨论

4.20芦山地震后,有学者在芦山县龙门乡发现一系列的线性裂缝和砖块的旋转变形等"地震地表破裂迹象",由此推测芦山—龙门一线存在隐伏逆断裂,并认为该断裂有可能是此次地震的发震断裂。因此,进一步探讨芦山—龙门一线是否存在潜在的发震断裂,无论是对研究芦山7.0级地震的发震断裂,还是对灾区的重建指导都十分重要。在龙门乡开展了地质灾害调查、跨谷地的地质剖面实测,槽探和人工地震勘探等工作。结果显示:至少在800m深度范围内,不存在芦山-龙门隐伏断裂。此带上的地裂缝等现象不是由断层位错引起,而更可能是地震动在阶地陡坎附近造成的地基或边坡效应所致。     

甘孜-玉树断裂当江段地表破裂遗迹与1738年玉树西北地震的讨论

对历史记载的公元1738年玉树西北地震的震级及其发震构造目前仍存有争议。卫星影像解译和野外调查发现沿甘孜-玉树断裂当江段分布一条长约75km的左旋走滑地震地表破裂带,其最大同震水平位移约2.1m。综合分析该地表破裂带特征、探槽揭露信息、测年结果以及历史文献记载等资料,认为当江段应为1738年玉树西北地震的发震断层,基于震例类比和经验公式估算该次地震的震级为71/2级。沿甘孜-玉树断裂的历史地震破裂分布显示,玉树段在隆宝镇以西存在近50km长的破裂空段;当江段距1738年地震的离逝时间也可能已经接近其地震复发周期,上述两个段落未来均存在大震危险。     

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,发震构造为具有少量走滑性质的逆冲断裂。从发震构造特征、同震滑动分布推测,伽师地震发震构造是柯坪塔格褶皱带滑脱面以上沉积盖层内的逆冲断裂,支持了柯坪推覆体的薄皮构造模型观点。     

Focal depths for moderate-sized aftershocks of the Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>8.0 earthquake and their implications

Sliding-window cross-correlation method is firstly adopted to identify sPn phase, and to constrain focal depth from regional seismograms, by measuring the time separation between sPn and Pn phases. We present the focal depths of the 17 moderate-sized aftershocks (M _S⩾5.0) of the Wenchuan M _S8.0 earthquake, using the data recorded by the regional seismic broadband networks of Shaanxi, Qinghai, Gansu, Yunnan and Sichuan. Our results show focal depths of aftershocks range from 8 to 20 km, and tend to cluster at two average depths, separate at 32.5°N, i.e., 11 km to the south and 17 km to the north, indicating that these aftershocks are origin of upper-to-middle crust. Combined with other results, we suggest that the Longmenshan fault is not a through-going crustal fault and the Pingwu-Qingchuan fault may be not the northward extension of the Longmenshan thrust fault. Supported by the National Natural Science Foundation of China (Grant Nos. 40604009 and 40574040) and Special Project for the Fundamental R & D of Institute of Geophysics, China Earthquake Administration (Grant No.DQJB08B20)     

Numerical simulation of strong ground motion for the <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>8.0 Wenchuan earthquake of 12 May 2008

The Wenchuan earthquake of 12 May 2008 is the most destructive earthquake in China in the past 30 years in terms of property damage and human losses. In order to understand the earthquake process and the geo-morphological factors affecting the seismic hazard, we simulated the strong ground motion caused by the earthquake, incorporating three-dimensional (3D) earth structure, finite-fault rupture, and realistic surface topography. The simulated ground motions reveal that the fault rupture and basin structure control the overall pattern of the peak ground shaking. Large peak ground velocity (PGV) is distributed in two narrow areas: one with the largest PGV values is above the hanging wall of the fault and attributed to the locations of fault asperities and rupture directivity; the other is along the northwestern margin of the Sichuan Basin and caused by both the directivity of fault rupture and the amplification in the thick sediment basin. Rough topography above the rupture fault causes wave scattering, resulting in significantly larger peak ground motion on the apex of topographic relief than in the valley. Topography and scattering also reduce the wave energy in the forward direction of fault rupture but increase the PGV in other parts of the basin. These results suggest the need for a localized hazard assessment in places of rough topography that takes the topographic effects into account. Finally, had the earthquake started at the northeast end of the fault zone and ruptured to the southwest, Chengdu would have suffered a much stronger shaking than it experienced on 12 May, 2008. Supported by the U.S. National Science Foundation (Grant Nos. EAR 0738779 and OCE 0727919), the National Basic Research Program of China (Grant No. 2004CB418404), and partially by the National Nature Science Foundation of China (Grant No. 40521002)     

1679年三河—平谷8级地震地表破裂端部特征及其地质意义

地震地表破裂端部的几何结构与运动学特征研究有助于科学认识断裂的破裂传播与终止过程。夏垫断裂是华北平原区最为重要的隐伏强震构造之一,于1679年发生了三河—平谷M8历史大地震,但其同震地表破裂长度及端部变形特征仍存争议。基于前人研究结果,在野外地质调查的基础上,跨1679年三河—平谷8级地震地表破裂端部布设了2条浅层地震勘探剖面,研究断裂端部的新活动特征。结果显示,断裂端部的最新活动时代为全新世,运动方式以走滑为主兼正断,且呈现出明显的滑动亏损特征。结合同震垂直位移分布等数据,分析认为该次地震的地表破裂仅长10余公里,与8级地震不匹配,其发震构造和机制仍需深入研究。