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

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

芦山地震和九寨沟地震台站实测震级差异性分析

本文选取2013年芦山地震和2017年九寨沟地震波形,重新量取垂直向振幅,计算宽频带面波震级M_S（BB）,分析各台站实测震级出现方向性差异的原因。其中,通过572个宽频带台站实测芦山地震震级M_S（BB）7.1,通过603个宽频带台站实测九寨沟地震震级M_S（BB）6.9。芦山地震实测震级大于M_S（BB）7.3的台站呈现WN-ES向分布,与断层倾向一致;实测震级小于M_S（BB）7.0的台站呈现NE-WS向分布,与其所在断层走向一致。九寨沟地震实测震级大于M_S（BB）7.0的台站分布呈现NE向分布,与断层倾向一致;实测震级小于M_S（BB）6.8的台站总体分布较为离散,大体呈现NW-SE向分布,与树正断裂走向一致。实测震级偏大的台站方向性分布与多普勒效应和P波辐射花样联系不明显。对比分析芦山地震和九寨沟地震,去除场地响应和仪器自身影响,台站实测震级差异性仍然存在,因此,台站实测震级差异性是由于受到了多普勒效应、辐射花样、仪器和场地响应之外的因素影响。综合考虑地震震级涉及的影响因素,芦山地震和九寨沟地震的台站实测震级差异性可能与地震波的传播路径有关。     

2017年8月8日九寨沟MS7.0地震序列分析

2017年8月8日四川九寨沟发生M_S7.0地震,该地震发生在巴颜喀喇块体的东北边界,震中区域构造条件复杂,是巴颜喀喇块体北侧左旋走滑环境向东侧逆冲挤压环境过渡的位置,附近地区历史强震较多。九寨沟地震是一次主-余型地震,余震活动水平较弱,主震发生后短时间内M_L≥4.0余震的“等待时间”存在异常,震后较长时间余震活动恢复到正常状态,序列h值、余震视应力等符合主-余型序列特征。序列b值为0.84,G-R关系推测序列最大余震的震级约为M_L5.4（M_S5.0）,8月9日发生的M_S4.8地震是目前该序列的最强余震。通过与1970年以来附近地区7级左右地震序列的对比认为,九寨沟地震与1976年松潘-平武2次7.2级地震序列在余震空间位置、发震构造和震源机制等方面存在较大差异,因此,不具备发育为震群型序列的条件。九寨沟地震主震视应力为0.36～0.38MPa,属于应力下调模型,序列余震的平均视应力水平接近龙门山断裂带附近中小地震的平均背景水平。     

基于地震预警参数阈值评估内蒙古地区地震潜在破坏区

利用内蒙古地区测震台网记录到的地震波形数据,基于地震预警参数阈值（特征周期τ_c和位移幅值P_d）开展评估内蒙古地区中强地震潜在破坏区范围的研究。按照一定的筛选条件进行筛选后,拟合出适合内蒙古地区的τ_c与震级M_L、P_d与τ_c和震源距R的统计关系式;根据现有地震事件的最大震级,得到仪器烈度为Ⅳ度、震级M_L=5.0时对应的参数阈值为P_d=0.1075cm和τ_c=0.686s。最后对近年来内蒙古境内2次中强地震开展线下模拟,结果表明该方法需要台网具有一定密度,如能均匀分布,效果更好,否则会出现结果仅由经验关系式决定、而实际观测值没有发挥作用的问题。     

2017年8月8日四川九寨沟MS7.0地震前甘肃南部地区视应力变化

利用甘肃省区域台网提供的地震波形资料,计算甘肃南部地区2010年1月～2017年8月M_L≥2.0地震的视应力,分析九寨沟地震前甘南地区视应力的时空演化特征,所得结果如下：①M_L2.6～3.1、M_L2.0～2.5震级段地震视应力的空间分布具有较好的一致性,九寨沟M_S7.0地震前,距震中较近的川甘交界地区的舟曲、文县等地呈现出较为集中的视应力高值异常;②M_L2.0～2.5、M_L2.6～3.1地震的视应力在大区域范围内随时间的变化趋势有一定差异,M_L2.6～3.1地震视应力变化更为明显,在较大区域范围内表现为震前显著升高,而M_L2.0～2.5地震视应力则在震前一两年内开始下降;③随着区域划分范围向震中靠近,M_L2.0～2.5、M_L2.6～3.1地震视应力的变化逐渐趋于一致,至震中附近时,2个震级段视应力均表现为“震前几年长时间升高—临震前几个月下降”的同步变化。     

2017年8月8日九寨沟MS7.0地震前成都台地磁谐波振幅比异常分析

2017年8月8日四川省九寨沟县发生M_S7.0地震,成都地震基准台距此次地震震中约255km,震前该台地磁谐波振幅比出现趋势性异常。对成都台GM4磁通门磁力仪秒数据进行谐波振幅比计算,结果显示,谐波振幅比在九寨沟M_S7.0地震前表现出下降—转折—恢复上升的异常变化形态,地震发生在异常恢复期,同时,异常表现出由长周期向短周期迁移、SN向与EW向变化不同步的特征,芦山M_S7.0地震前成都台地磁谐波振幅比亦呈现了相似的异常变化特征。     

2018年9月4日新疆伽师MS5.5地震序列及发震构造讨论

2018年9月4日新疆伽师发生M_S5.5地震,震中处于塔里木地块西北缘,位于1997~1998年伽师强震群震区内。此次伽师地震前发生了M_S4.7前震,截至9月30日最大余震震级为M_S4.6（M_L5.0）,初步判定为前-主-余型地震序列。序列精定位结果显示,余震沿近NE向展布,主震震源深度与1997~1998年伽师强震主震基本一致,发震断层陡立。本文从区域的构造环境、地震震源机制解和余震分布特征等方面分析认为,地震发生在伽师隐伏断裂东南端部,为1997~1998年伽师强震群震区的一次新的构造活动。序列参数、视应力等计算结果显示,伽师M_S5.5地震的预测最大余震震级与最大余震震级M_S4.6接近,表明序列最大余震已经发生。     

2015年7月3日新疆皮山MS6.5地震烈度与地震地质灾害特点

2015年7月3日新疆和田地区皮山县发生的M_S6.5地震极震区烈度达到Ⅷ度,砂土液化、地裂缝和崩塌等地质灾害发育。极震区内砂土液化发生范围广泛,地表发生喷砂冒水和边坡的侧向移动,对道路、水利等设施造成了严重的破坏;地裂缝断错多处路面,宽度最大达30 cm,裂缝多呈NWW走向;崩塌主要发生在高阶地边坡地带,人烟稀少,未造成人员、财产损失。本次地震地质灾害是新疆近十年来M_S6以上较为严重的一次。     

利用反投影技术快速测定中强震矩震级及破裂过程——以2023年9月9日摩洛哥MW6.9地震为例

利用远场台阵反投影技术可以较少地依赖先验经验,于震后快速获得震源破裂过程时空分布特征,并据此估算矩震级。本文使用美国阿拉斯加宽频带台阵记录到的远场直达P波数据,基于反投影方法,对2023年9月9日摩洛哥M_W6.9地震震源破裂过程及震级进行分析。结果显示,地震破裂走向以NE-SW向为主,地震能量在10～20s集中在震中距10km范围内释放,与震源机制解和震中区亚特拉斯山脉走向基本一致。震源破裂持续时间结合在全球分布均匀的GSN台网记录到的P波最大位移数据,在震后约15min即可快速估算出本次地震矩震级为M_W7.0,与中国地震台网中心和美国地质调查局发布震级基本一致。     

2022年1月8日青海门源MS6.9地震同震形变场及地下位错反演

针对2022年1月8日青海门源M_S6.9地震,基于Sentinel-1A卫星的SAR影像获取了本次地震的同震形变场,并进行计算分析,运用SDM程序反演地下静态位错。初步结果如下：①本次地震震中位于蝴蝶状分布的形变图中心,抬升盘与下降盘之间的破裂迹线明显,破裂迹线长约24km,形变中心位于托莱山与冷龙岭断裂交汇处。②降轨数据显示,在断层区域视线向（LOS）形变量最小约-0.55m（远离卫星）,最大约0.68m（靠近卫星）;升轨数据显示,在断层区域视线向（LOS）形变量最小约-0.49m（远离卫星）,最大约0.42m（靠近卫星）。③对中国地震局发布的门源6.9级地震烈度图和本研究获取的地震形变图进行叠加分析,可以看出二者吻合度较高,地表有较明显形变所涉及区域约0.4万km²。④运用降轨InSAR数据反演得出地下最大位错量约3.29m,深度4.75km,在近地表位错破裂长度可达30km,地表有明显破裂的区域长约25km,与中国地震局科考工作每日情况的结果基本一致。     

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地震的发震断裂为冷龙岭断裂     

InSAR数据约束的2021年5月21日云南漾濞MS6.4地震发震构造研究

运用Sentinel-1A卫星数据和D-InSAR技术,获取2021-05-21云南漾濞M_S6.4地震的同震形变场。结果显示,漾濞地震同震形变场长轴近NW展布升降轨形变场符号相反,视线向最大沉降量和抬升量为0.1 m。InSAR同震形变场反演的滑动分布主要集中在沿走向2～12 km,倾向1～9 km的范围内,最大滑动量0.35 m,发震断层长9.8 km、宽4 km,滑动量主要集中在地下3～6 km范围内,滑动角-146.7°。同震位移场及滑动分布模型反映本次地震为发震断层的右旋走滑事件,地震破裂未达到地表。断层模型反演结果显示,矩震级为M_W6.1,发震断层以北西走向右旋走滑运动为主,初步认为本次M_W6.1地震发震断裂可能是一条NW向的维西—乔后断裂西侧的隐伏次生断裂。     

Preliminary analysis on characteristics of coseismic deformation associated with M S=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. Foundation item: National Natural Science Foundation of China (40374013) and “Researching on the Disaster Earthquake” (2003) of Public Welfare Research Item, Ministry of Science and Technology of China.     

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.     

新10井水位对九寨沟M_S7.0、精河M_S6.6地震同震响应

基于新10井水位对九寨沟M_S7.0、精河M_S6.6地震的同震响应特征,分别利用水震波和潮汐分析法,反演得到了新10井水位在地震波作用期间和地震波作用前后含水层参数特征,探讨了2次地震引起的新10井同震响应变化机理。结果表明,新10井水位对九寨沟地震和精河地震具有不同的同震响应形态,利用水平流模型反演所得的新10井含水层渗透系数,在九寨沟地震时为61m/d,而在精河地震时为147m/d,表明地震波作用导致水平向的渗透性增强;利用垂向流模型反演得到的新10井含水层渗透系数,在2次地震之前约为49×10~(-5)m/d,而在地震之后约为18×10~(-5)m/d,表明地震波作用导致垂直向的渗透性减弱。因此,新10井水位的同震响应变化与地震波引起的含水层渗透率的改变有关。     

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

九寨沟M_S7.0地震对巴颜喀拉块体东北缘活动断裂影响的有限元模拟

在详细调研地震地质资料的基础上,构建了巴颜喀拉地块东北缘三维有限元模型。以九寨沟M_S7.0地震同震位错为荷载,模拟计算了九寨沟地震的发生对巴颜喀拉块体东北缘主要活动断裂加卸载效应的影响。模拟结果显示,九寨沟地震的发生对龙日坝断裂、虎牙断裂、青川-平武断裂西段、迭部-白龙江断裂西段和东段、临潭-宕昌断裂东段,以及处于甘青川交界危险区内的东昆仑断裂东段、塔藏断裂西段,处于六盘山南-西秦岭东危险区的西秦岭北缘断裂东段表现为库仑应力加载;对岷江断裂、塔藏断裂东段库仑应力卸载效应显著。     

九寨沟MS7.0地震前地磁异常

分析距离九寨沟M_S7.0地震震中分别为65 km、254 km和238 km的松潘地磁台、成都地磁台和天水地磁台2017年连续观测结果,发现九寨沟地震前,松潘地磁台出现异常变化,主要表现为2017年5月开始,松潘地磁台长时间垂直分量累积变化量达40 nT。短时间内也存在地磁场突变,最高可达10 nT,而成都地磁台和天水地磁台地磁场无明显异常。通过对观测数据的详细检查发现松潘地磁台很多观测日的最后一个观测值与后一观测日的第一个观测值相差超过1 nT,因此认为九寨沟地震前观测数据的异常是由观测系统异常导致。对三个台站地磁场垂直分量日变化的统计分析发现垂直分量在不同时段相关性保持一致,说明地震前短时间内不存在地磁垂直分量的明显异常变化,三个台站地磁场日变化的差异是由地下介质性质不同导致。岩石圈磁场在震中西侧区域出现正负异常交替现象,这可能与九寨沟地震的孕育有关。     

九寨沟M_S7.0地震震源机制解及构造应力场研究

2017年8月8日四川九寨沟县发生M_S7.0地震。根据中国地震台网固定台站的记录波形,利用gCAP方法和P波初动符号反演方法求解了主震的震源机制解,其结果与哈佛CMT震源机制解、美国地质勘探局发布的震源机制以及中国地震局发布的各震源机制解有很好的一致性。此外,还求解了九寨沟地震震中附近2010—2016年间28次小震的震源机制解,利用得到的小震震源机制解反演了该区域的构造应力场。结果表明,这一区域构造应力场的最大应力轴与最小应力轴均为水平方向,其中最大应力轴方向为NWW方向,最小应力轴方向为SSW方向。