青海门源6.9级地震同震及震前GNSS变形特征分析

北京时间2022年1月8日1点45分,在我国青海省门源县发生了6.9级地震.通过震中附近陆态网络GNSS连续观测数据得到的同震位移场显示,距离震中最近的QHME站同震位移最大,EW向达到20.31 mm,SW向达到-35.45 mm,震中附近五个站的同震位移反映出此次地震的左旋同震破裂特征;GNSS站间基线时间序列结果...     

GNSS约束的2022年泸定M6.8地震滑动分布及同震应力变化

2022年9月5日四川泸定县发生M6.8级强震,震中位于鲜水河断裂东南末端,野外地质调查初步结果显示本次地震并未发现明显的地表破裂迹象.本文基于震后科考GNSS流动观测和震中周边连续站观测资料,解算并提取震中90 km范围内31个测站的静态水平向同震位移.结果显示：GNSS同震形变场空间分布呈现明显的左旋走滑特征,观测到的最大水平向同震位移达23 cm,震中距50 km范围内同震位移量普遍大于1 cm.基于GNSS观测资料反演得到的同震滑动分布显示泸定地震地表破裂主要集中在磨西至田湾之间,主破裂深度2～8 km,最大滑动量～1.96 m,地震矩9.25×10¹⁸N·m,对应矩震级M_W6.6.静态库仑应力结果显示本次地震增强了震源破裂区周边活动断层的库仑应力,并触发了大量余震,余震主要发生在库仑应力增强区域.结合震间闭锁分布、历史地震及库仑应力变化,我们认为未来需要密切关注与磨西断裂交接的安宁河、大凉山断裂以及康定—磨西段的地震危险性.     

GNSS观测的2021年5月22日玛多MS7.4地震同震位移及其约束反演的滑动破裂分布

2021年5月22日青海省果洛州玛多县发生M_W7.4地震,此次地震产生的地表破裂在空间上表现出明显的分段特征.本文基于不同来源的GNSS连续观测网数据获取了此次地震的精细三维同震形变场,结果显示:观测到的最大水平位移量达到280 mm,最大垂直形变量仅为25 mm,暗示此次地震的逆冲分量较小;此次地震具有较为明显的左旋走滑特征,同震形变基本对称,在NW-SE向的影响范围更广,该方向上水平同震形变大于3 mm的震中距范围超过500 km.进而,本文以余震精定位结果和GNSS观测到的三维同震形变场为约束,构建了地表破裂线为折线、倾角为85°、倾向西南的断层模型,反演了滑动破裂分布.结果显示:滑动破裂分布在震中两侧不均匀,均破裂到地表,破裂深度达到15 km左右,最大滑移量为4.73 m,计算的矩震级为M_W7.37.该结果与余震精定位结果具有很好的一致性,破裂的极值区正好位于早期余震空区,推测该余震空区未来的发震风险性较低.最后基于反演结果模拟计算了震中区域形变和应变场,结合应变值在断层地表迹线东南侧呈现挤压特征和已有的研究成果,推测此次地震增强了巴颜喀拉块体在东部地区挤压应力的积累特征,导致东部地区发震危险性增强,值得后续跟踪研究.     

2022年1月8日青海门源MS6.9地震GNSS同震形变场及其断层滑动分布

2022年1月8日在青海省海北州门源县发生M_S6.9地震,本次地震是继2016年门源M_S6.4地震后冷龙岭断裂周边发生的又一次强震。确定本次地震的破裂分布对分析该地区震害风险具有重要意义。通过收集震中及周边12个GNSS连续站点和震后加密观测的17个流动站点观测资料,获取了震中100 km范围内29个测站的GNSS静态同震形变场,并以此为约束反演了本次地震同震滑动分布。结果显示,近场GNSS观测到的最大形变量可达1.3 m。反演的最优破裂模型显示该地震主破裂区深度位于0～10 km,滑动破裂出露地表,最大滑动量为4.07 m,地震矩释放能量约1.1×10¹⁹ N·m,对应矩震级M_W6.7。门源地震破裂至地表是造成该地区基础设施破坏的直接原因。     

GNSS测得的大地震前兆水平形变

以21世纪全球各具特色且最典型的6次大地震,尤其是巨大地震前后GNSS观测得到的地壳水平形变为重点,回顾探寻前兆形变,探索地震预测预报。6次地震为:2008年汶川8.0级、2010年智利8.8级、2011年东日本9.0级、2015年尼泊尔8.1级、2016年日本九州岛7.3级地震及2023年土耳其7.8级双震。地壳形变是物理问题,区域参考框架是利用GNSS连续观测研究前兆位移物理问题的基本前提,而同震水平位移是探寻前兆形变的关键现象。大震前震中及其附近既无明显的垂直位移积累,也无明显的水平剪切位移积累。大地震震中及其附近前兆水平位移主要有两种形态,水平位移达到峰值或为闭锁;但近震中震前介质已非弹性,仅远场介质为弹性。GNSS观测结果表明震前震中主要为水平挤压,震时产生剪切破裂,这与岩石破裂实验结果一致,因此“压-剪”(弹性)回跳模型符合GNSS观测结果,但至今仅获得少数地震短临前兆形变观测结果。GNSS依然是探索地震预测的主要观测手段,但由于与地震断层破裂临界状态密切相关,利用各种观测技术探索地震断层破裂的临界现象,即地震短临前兆,是突破地震短临预报的关键,仍需要利用多学科各种观测技术...     

基于GNSS和InSAR约束的2021年玛多MS7.4地震同震滑动分布及应用

2021年5月22日青海省果洛藏族自治州玛多县发生M_S7.4地震,震中位于青藏高原中部的巴颜喀拉块体,这是近20多年来在巴颜喀拉块体周边发生8次M≥7级强震后,块体内部的一次强震,也是汶川地震以来中国大陆发生的最大一次地震,因此该地震的成因及周边地区未来的地震危险性值得重点关注.本文利用震后及时获取的39个近场流动GNSS观测,联合61个GNSS连续观测、Sentinel-1和ALOS-2 InSAR观测获取了本次地震精细的同震形变场,以此为约束,基于均匀弹性半无限位错模型,反演了发震断层的滑动分布,并计算了同震库仑应力变化.GNSS水平同震形变十分显著,断层南北两侧的GNSS点位,最大水平形变分别达0.7 m和-1.2 m,距震中200 km的测点仍有1 cm左右的同震形变.Sentinel-1和ALOS-2的升降轨InSAR同震形变场显示此次地震造成了约160 km长的地表破裂,最大视线向形变分别达0.9 m和1.2 m.同震滑动分布模型显示,发震断层由主段和次段组成,长度分别为170 km和20 km,主段倾向北,倾角85°,平均滑动角为-4.36°,表明玛多地震是一次典型的走滑型地震.次段倾向南,倾角68°,平均滑动角为-11.84°.地震破裂主要集中在0~15 km深度范围,最大滑动量为4.4 m,对应深度6.97 km.反演给出的矩震量为1.61×10²⁰N·m,对应矩震级M_W7.4.主发震断层上存在4个凹凸体,玛多地震是一次不对称双侧破裂事件.结合余震精定位、野外调查及地质资料,我们认为主发震断裂为昆仑山口—江错断裂,东部的次级破裂与主破裂机制不同.同震库仑应力结果显示,东昆仑断裂玛沁段应力有所增加(>0.01 MPa),处于应力加载状态,未来发生强震的危险性较高.     

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和重力观测的解释提供理论参考依据.     

2021年青海玛多MW7.4地震GNSS同震形变场及其断层滑动分布

2021年5月22日青海玛多M_W7.4地震作为发生在巴颜喀拉块体内部的一次强震,再次引起了人们对该地区地震活动性的强烈关注.本文基于震后GNSS流动观测和区域连续GNSS站资料,解算了106个站点的同震形变及其中17个站点的高频形变波形.同震形变场显示玛多地震具有典型的左旋走滑特征,GNSS观测到的最大同震位移达到1.2 m.GNSS与InSAR数据相符度较高,GNSS提供了准确的近场形变信息.基于GNSS同震形变场,本文反演了断层滑动分布,并计算了发震断层上产生的库仑应力变化.结果表明,发震断层的滑动破裂存在多个凹凸体,破裂分段特征明显且出露地表,与野外地表破裂考察和余震分布吻合,主体破裂位于断层面0~10 km的浅部区域,最大滑动量达到4.6 m,地震矩1.63×10²⁰N·m,矩震级为M_W7.4;发震断层上静态库仑应力增加区域与余震分布具有一致性,说明余震主要是由静态库仑应力加载而触发的.     

2022年泸定6.8级地震GNSS同震形变场及其约束反演的破裂滑动分布

2022年9月5日四川省甘孜州泸定县发生6.8级地震,地震发生在川滇菱形块体东边界鲜水河断裂磨西段西侧附近.本文基于震中350 km范围内中国大陆构造环境网络和中国地震科学实验场118个GNSS连续站数据,观测获得了精细的同震形变场,结果显示：同震形变跨发震断层呈空间对称分布,表明本次地震具有显著的左旋走滑特征;记录到的最大同震形变发生在震中距40 km的SYD5(石棉安顺场)站,东西向和南北向形变量分别达到-22.0±1.2 mm和11.6±0.9 mm,震中距100 km以外测站的水平同震形变均小于5 mm;垂向同震形变不显著.结合走向163°和倾角77°的发震断层模型,本文对发震断层面的同震破裂滑动分布进行了反演,结果显示：同震滑动主要集中在主震东南侧余震空区0～15 km深度范围内,且破裂达到了地表,同震释放的地震矩与矩震级为M_W6.57地震相当.结合理论同震形变场、主应变场和邻近区域主要活动断裂库仑应力变化的空间分布特征,本次可能导致处于强闭锁和地震空区的安宁河断裂石棉—冕宁段未来发震风险性增强.     

利用空间观测技术研究青海玛多7.4级地震孕育发生变形时空特征

为揭示2021年青海玛多7.4级地震震前变形与同震破裂特征,本文系统收集与处理了青海省CORS站和中国地壳运动观测网络等GNSS资料、以及Sentinel-1和ALOS-2卫星SAR数据,通过多源数据综合处理获取了玛多地震形变场动态结果.震前GNSS形变场结果显示,玛多震源区处于其北侧剪切形变高值区边缘,震前存在左旋剪...     

2010年青海玉树地震同震-震后形变场特征及演化过程

基于星载合成孔径雷达差分干涉测量技术(DInSAR)和4期ENVISAT/ASAR雷达数据,获得了不同时间基线的三个同震干涉形变场和两个震后干涉形变场,并对这五个在时段上互有重叠的形变场进行了综合分析.结果表明,玉树地震同震形变场为围绕发震断层NW展布的椭圆形干涉条纹,覆盖范围约89 km×59 km.断层运动性质为左旋走滑.两盘最大视线向相对形变量至少达45 cm,最大形变出现在结古镇附近.时间基线不同的同震形变场总体上基本一致,但两盘最大相对形变量和局部形变存在差异.震后时间较长的干涉对反映的最大形变量反而减小;在震后时间较短的干涉对上于结古镇西南侧观测到的局部形变,在震后时间较长的干涉对上却没有出现.分析认为在形变量最大的结古镇附近可能出现了震后快速弹性回弹,导致随震后时间延续,形变量反而减小的现象.玉树地震震后形变主要出现在断层附近、震后不久的时段内,形变量在8 cm以下,具有与同震方向一致和相反的两种震后形变方式.在结古镇西南观测到一个与同震形变相反的局部沉降,应为震后弹性回弹.在微观震中处的断层附近观测到与同震方向一致的震后形变,可能是震后余滑.通过对地震前后不同时间基线的多个干涉对的联合对比分析,可以在一定程度上区分同震形变与震后形变,更好地研究地震引起的变形过程,特别是地震断层附近短期震后形变场的演化过程,为进一步研究断层带的岩性特征、物理力学及运动特性提供约束.     

基于GNSS与InSAR约束的九寨沟MS7.0地震滑动模型及其库仑应力研究

2017年8月8日的九寨沟M_S7.0地震发生在岷江断裂、塔藏断裂及虎牙断裂交汇地区,地处青藏高原东北部的川甘交界地区,位于巴颜喀拉地块的东缘,地质构造复杂,对于九寨沟地震震中位置和发震断层的确定,存在不同意见.本文利用GNSS及升降轨InSAR观测,在获取九寨沟地震同震形变场的基础上,基于均匀弹性半无限位错模型,联合反演了发震断层的滑动分布模型,并计算了同震库仑应力变化.InSAR同震形变场显示,视线向最大沉降量和抬升量分别为0.21 m和0.16 m,形变场长轴为NW向,形变主要集中在断层西侧.距震中40 km和65 km的九寨和松潘两县,水平向的GNSS同震位移分别达14.31 mm和8.22 mm.联合GNSS和InSAR同震形变场反演得到的滑动分布主要集中在沿走向5~33 km,倾向2~20 km的范围内,平均滑动量为0.18 m,最大滑动量为0.91 m.发震断层长40 km,宽30 km,走向155°,倾角81°,滑动角-9.56°.同震位移场及滑移分布模型表明此次地震为一次左旋走滑为主的地震事件,地震破裂并未完全到达地表,与虎牙断裂北段的几何产状和运动学性质更为接近,结合精定位余震的分布,我们确定虎牙断裂北段为此次地震的发震断层,震中位于北纬33.25°,东经103.82°,震源深度10.86 km,矩震量为7.754×10¹⁸ Nm,相应的矩震级为M_W6.5,与美国地调局和哈佛大学给出的震源机制解基本一致.同震库仑应力导致了虎牙断裂北段延长线的东北和西南两端应力增强,其中塔藏断裂的罗叉段和马磨段未来强震的危险性值得关注.     

联合GNSS和InSAR观测位移反演2008年汶川大地震断层位错模型参数

利用现代空间大地测量技术,尤其是卫星合成孔径雷达干涉测量,能够获取高精度、高空间分辨率的同震和孕震形变,为地震断层形变和破裂机制研究提供了前所未有的机遇。本文介绍了利用大地测量观测数据反演地震断层位错模型参数的贝叶斯反演方法。联合运用2008汶川大地震前后GNSS和InSAR技术观测获得的同震位移,反演了地震断层的几何参数和滑动位错分布。研究结果表明,汶川地震的断层滑动主要集中在倾角较陡的浅部,同时包含逆冲和右旋走滑,其中最大逆冲6.1m,最大右旋6.5m。根据断层滑动分布正演计算得到的上盘同震位移明显小于下盘,预示该断层两侧孕震形变可能存在较大的不对称性。     

利用GPS观测资料分析2008年于田Ms7.3地震的同震位移及震后形变

2008年3月21日新疆于田发生Ms7.3级地震.本文通过处理、分析GPS数据,得到破裂断层北侧100 km附近的同震位移及震后形变信息.在观测区域GPS点监测到10 mm左右的同震位移,其中最大为南向14 mm,东向5 mm.同震位移呈现一致性的东南向运动特征,证实于田地震存在显著的左旋走滑分量.震后台站向西南方向运...     

Dynamic pattern characteristics of fault deformation and gravity field in the development process of Yongdeng M S=5.8 earthquake

In this paper, the fault deformation abnormality, dynamic evolution features of gravity and vertical deformation field in the seismogenic process of the Yongdeng, Gansu Province earthquake on July 22, 1995 are studied primarily. There appeared α-β-γ tri-stage anomaly at three sites near the epicenter, and there appeared anomalies of step and sudden jump at more than 10 sites in outer region since 1993. The high value area before shock, coseismic effect and process of recovery aftershock were monitored by portable gravity data. Data reflects the changing process of fault movement from the quasi-linear to the nonlinear in the near source region during seismogenic development of the Yongdeng earthquake and evolution of gravity field from heterogeneity of seismogenic term to quasi-homogeneity of postseismic term. There exists close relationship between strong earthquake and dynamic evolution of regional stress-strain field. Considering all above, the experience and lessons in this medium-short term prediction test are summarized.     

2014年新疆于田MS7.3地震同震位移及位错反演研究

利用于田震中300 km范围内的1个GPS连续站和12个GPS流动站数据,解算得到了2014年新疆于田M_S7.3地震地表同震位移,并反演了发震断层滑动分布,探讨此次地震对周边断裂的影响.地表同震位移结果显示,GPS观测到的同震位移范围在平行发震断裂带的北东-南西向约210 km,垂直发震断裂带的北西-南东方向约为120 km,同震位移量大于10 mm的测站位于震中距约120 km以内;同震位移特征整体表现为北东-南西方向的左旋走滑和北西-南东方向的拉张特征,其中在北东-南西方向,I069测站位移最大,约为32.1 mm,在北西-南东方向,XJYT测站位移最大,约为28.1 mm;位错反演结果表明,最大滑动位于北纬36.05°,东经82.60°,位于深部约16.6 km,最大错动量为2.75 m,反演震级为M_W7.0,同震错动呈椭圆形分布,以左旋走滑为主并具有正倾滑分量,两者最大比值约为2.5：1,同震错动延伸至地表,并向北东方向延伸,总破裂长度约50 km,地表最大错动约1.0 m;同震水平位移场模拟结果显示贡嘎错断裂、康西瓦断裂和普鲁断裂等不同位置主应变特征具有差异性,这种差异特征是否影响断裂带以及周围区域的应力构造特征,值得关注.     

Coseismic and postseismic subsurface displacements and strains for a vertical strike-slip fault in a three-layer elastic medium

A three-layer elastic-gravitational fault displacement model using dislocation theory has been developed and used to examine the effect of layering of earth elastic moduli on surface and subsurface displacement fields for a vertical strike-slip fault. The model has been used to examine the effect of depth variation of elastic properties at coseismic and postseismic time scales. For pure strike-slip motion the effect of gravity on coseismic and postseismic horizontal deformation is negligible. For coseismic deformation the model predicts that (for constant Poisson's ratio) an increase in elastic moduli with depth attenuates the displacements within the upper layers with respect to displacement distribution for a uniform half-space, while an inclusion of a soft layer between the top layer and lower half-space amplifies upper layer displacements. The effect of variation in Poisson's ratio on surface and subsurface displacements has also been examined.The effect of postseismic stress relaxation on surface and subsurface displacements for a three-layer model has been calculated and compared with that of a uniformly relaxed half-space model. Layer 1 is assumed to correspond to the upper crust, layer 2 the lower crust and layer 3 the upper mantle. The effect of postseismic stress relaxation within a uniform half-space and within just the lower crust and upper mantle has been examined. Stress relaxation within the whole half-space decreases the amplitude and shortens the wavelength of displacements, while stress relaxation within the lower two layers increases the amplitude and broadens the wavelength of displacements. The difference between uniform and layered postseismic relaxation is particularly pronounced at the base of the crust.Coseismic and postseismic normal and volumetric strains for a vertical strike-slip fault have also been examined. For a uniformly relaxed half-space model, an increase in normal strains is shown with respect to the coseismic elastic solution, whereas the postseismic volumetric strain is effectively zero. For a three-layer model with stress relaxation in the lower layers only, the normal and volumetric strains within the top elastic layer resemble coseismic strains, while in the lower layers which suffer a rigidity decrease, the postseismic volumetric strain is effectively zero.     

2001年昆仑山口西8.1级地震同震形变场特征的初步分析a

利用差分干涉雷达测量技术获取的宏观震中区的同震形变场，结合对地震活动性、震源机制、野外考察等资料分析，对昆仑山口西8.1级地震同震形变场特征进行了研究. 结果表明:宏观震中位于库赛湖东北侧，宏观震中区发震断层可分为两个形变中心区域，其中西段长约42 km，东段长约48 km，整个发震断层主破裂段长90 km；由干涉形变条纹分布格局可清楚地判断出发震断层的左旋走滑特征；断层两盘变形特征不同，南盘变形程度明显大于北盘；宏观震中附近最大斜距向位移量为288.4 cm，最小斜距向位移量为224.0 cm，宏观震中发震断层最大左旋水平位错为738.1 cm，最小地面左旋水平位错为551.8 cm.      

Preliminary analysis on characteristics of coseismic deformation associated with M_S=8.1 western Kunlunshan Pass earthquake in 2001

IntroductionOnNovember14,2001,aMS=8.1earthquakeoccurredonthewestofKunlunshanPassintheborderareaofQinghaiandXinjiang,whichwasthestrongestearthquakeinChinesemainlandsincetheMS=8.0earthquakeoccurredinDangxiongdistrictofXizangAutonomousRegiononNovember18,1951.TheearthquakeoccurredontheEasternKunlunTectonicZone,whichwasapalaeoplatejunctionzoneinsideTibetanPlateau.ItdividedTibetanPlateauintothesouthandnorthparts.ThezoneplayedaveryimportantroleinTibetanPlateausdeformationprocessanddynamicev…