基于GPS的2021年玛多MW7.4地震同震静、动态形变场

利用青海省CORS网GPS连续观测数据获取2021年玛多M_W7.4地震的同震形变场以及震时地表运动状态。结果表明,断层南北两侧的近场测站分别展现出南东向和北西向运动,符合左旋走滑机制。同震形变集中在震中距300 km范围内,震中距150 km以内的站点均能反映出cm级的位移,最大水平向位移为28.3 cm(JDUO站)。高频GPS动态形变提取出的永久位移与静态解算结果相当,其动态波形最大峰值为49.9 cm(KANQ站)。依据震中距和波形初动时间估计得到地震波速为2.8~4.9 km/s,断层东端的站点估算速度大于其他站点,可能与此次破裂的方向性或者破裂传播速度有关。依据震级经验公式,动态波形估算震级在M6.8~7.6区间,拟合平均震级为M7.35±0.15。若在实时条件下,震后70 s内可得到稳定的震级。     

高频GNSS观测快速估计地震震级及其精度评估

以2016-11-13新西兰地震为例,利用PANDA软件的精密单点定位(PPP)模块,基于高频GNSS数据分析地表形变并计算震级,在此基础上定量评估GNSS地震波测量精度对震级估计的影响。结果表明,当GNSS测站的峰值地面位移（PGD）精度在2 cm范围内时,对震级估计的误差影响在0.1以内。另外,由于震级估计与测站震中距和PGD有关,震中距越大,估计的震级对PDG越敏感。
     

高频GPS测定的尼泊尔MW7.8和MW7.3地震震时地表动态变形过程

以陆态网络和尼泊尔境内高频（1 Hz）GPS观测数据为基础,采用动态双差相对定位方法,获取2015年尼泊尔M_W7.8和M_W7.3地震震时近场地表动态变形过程。结果表明,M_W7.8地震震中东侧高频GPS站动态位移幅度明显大于震中西侧;各高频GPS站动态位移幅度不仅与测站震中距有关,而且与地震破裂传播方向有关;M_W7.3地震引起的水平动态位移相对较小。将高频GPS与邻近强震仪动态位移时序进行对比发现,二者在振幅和相位上具有较好的一致性。     

利用高频GPS数据研究汶川地震引起的近场地表运动学特征

利用单历元定位方法处理2008年汶川Ms8.0级地震震中附近GPS连续站1 Hz数据,获取近场地壳运动学特征。观测结果表明,测站运动的方向和幅度不仅与测站与震中距离相关,而且与破裂带的距离相关,近场运动学的复杂性反映了震源破裂过程具有逆冲和走滑的特征。距震中最近的郫县站最大运动学位移EW方向超过100 cm,NS方向达到70 cm,其GPS位移和强震仪得到的位移一致。对12个GPS测站的运动学位移与永久静态位移进行对比,发现6个测站前者明显大于后者,5个测站二者大小相当,成都站前者小于后者。     

日本Mw9.0强震的高频GPS形变波特征初探

利用高频GPS数据长基线双差精密定位技术,获取2011年3月日本Mw9.0强震近场和远场GPS连续参考站的瞬时动态形变。分析表明：距震源较近的IGS站（USUD）的同震动态形变水平方向幅度超过60 cm,震后水平形变达到30 cm;位于远场的GPS连续参考站（JLYJ）的动态形变,清晰地给出了地震面波引起的地表形变波形。对JLYJ站的谱分析表明,地震面波引起的动态形变周期为10～17 s,并出现了两次峰值,东西向和垂向上的动态位移量远大于北南向,说明JLYJ站的形变主要由Rayleigh波及其谐波所主导,Rayleigh波的传播速度为3.5 km/s,波长约50 km。     

球体位错理论在2011年日本强震中的应用研究

利用球体位错理论计算了2011年日本强震产生的远场同震位移与应变,并利用GPS远场数据修正了该强震的总地震矩。结果表明：1）利用球体位错理论计算得到的理论水平位移场显示,垂直于发震断层的广大区域同震位移较大,位移矢量总体都指向震中地区,震中距约5 000千米的地方亦产生了3 mm以上的同震水平位移。理论位移与远场GPS观测结果具有良好的一致性;2）比较两个独立断层模型对应的理论同震位移场发现,震源西部地区远场位移总体上只有1%～4%的微小差异,而东部广大海域的差异则达到同震信号的6%～15%,震中周围差异更大。该差异表明,相对于震源仅局域覆盖的日本本土GPS观测数据对2011年日本强震的断层滑动分布模型的约束能力有限;3）依据中国及邻区的远场GPS同震观测数据修正2011年日本强震的总地震矩,把该地震释放的总能量约束在(3.24～4.96)×10 ²² Nm,相应的矩震级为Mw8.97～9.10;4）2011年日本强震在华北地区产生的同震应变与该区的长期应力变化背景场大体相反,表明该强震使华北地区的地壳产生了松弛效应。     

地震预警中的单台综合定位方法

基于日本基盘强震观测网中某个地震事件首个地震台记录的震级在4级以上、震中距为30~100 km的强震数据,运用B-Δ法和Voronoi多边形法相结合的方法计算震中距。结果表明,B-Δ法和Voronoi多边形法相结合计算出的震中距与单用B-Δ法的计算结果相比,单台定位所用时间相同,但其定位精度明显提高。     

应用GPS观测资料反演潜在震源参数

根据点源孕震力学模型 ,应用福建 台湾地区GPS位移速度资料反演 1999年台湾南投Mw7.6大震前的孕震参数。反演所得的震中与实际发生的地震震中相差 6 7km ,主压应力场方向为 319.9°,倾角 - 2 8.6° ,即仰角2 8.6° ,与用震源机制解得到的P轴方向和倾角分别相差 17°和 8.6° ,与菲律宾板块向亚洲板块挤压的方向是一致的。用土耳其地区GPS位移速率资料反演 1999年土耳其Mw7.6地震前的孕震参数 ,反演得到的震中与实际震中相差 70km ,主应力方位角和倾角分别为 16 2 .5°、11.5°,与震源机制解结果分别差 2 4 .5°和 1.5°。     

加速度记录测定近震震级在江苏地区的计算研究

基于江苏强震台网和江苏测震台网的地震记录,采用时间域内实时仿真WA记录方法,对2007~2018年记录到的M2.3以上地震进行近震震级测定。为了定量分析基于加速度测定近震震级的可靠性,首先采用国标规定的新震级测定方法和量规函数测定了基于加速度记录的震级,并将结果与江苏测震台网基于速度记录测定的结果进行对比分析,得到线性拟合关系式：ML(ACC)=0.987ML＋1.371;另外计算了15 324组基于加速度测定的震级与人工编目震级之间的绝对偏差,结果显示,偏差在±0.5级以内的结果占87%;通过4阶多项式拟合后发现,在小于200 km范围内测定的震级略低于人工编目的结果,震中距超过200 km后无系统偏差。总体而言,两者结果的一致性较好,表明利用加速度记录仿真WA位移记录来计算近震震级是可行的,得到的结果也是可靠的。     

库尔勒台JCZ-1T型地震计近震震相特征分析

分析库尔勒台超宽频带JCZ-1T型地震计记录的不同震中距的近震波形特征,发现新疆部分区域记录的地震波形有其特殊的特征：1）震中距小于230 km时,记录到Pg、Sg波;震中距大于230 km时,可以记录到首波Pn;2）除昆仑山特殊区域外,震中距大约为500 km时,面波发育不完全;3）震中距约大于700 km的喀什-乌恰交汇区范围,Pg、Sn波发育不清晰。     

利用GPS高频信号分析玛多地震产生的动态位移特征

为研究玛多地震对周围地表动态位移的影响,利用超长基线模式解算震中距400 km范围内的11个高频GPS观测数据,并对各站水平方向的位移时间序列进行分析.结果表明,玛多地震对震中距最近的青海玛多站的东西向产生约25 cm的位移,并且没有恢复到震前状态;对青海玛沁站和甘肃玛曲站的南北向产生约27 cm的位移,但很快恢复到震...     

高频GNSS资料的强震震相识别与震中确定

利用2011日本东北Mw9.0地震和2008汶川Mw7.9地震若干台站的高频（1 Hz）GNSS观测资料,通过精密单点定位方法解算出地震发生时的位移形变波形。对位移波形一阶差分得到震时高频GNSS速度时间序列,采用S变换谱的方法对该速度时序进行地震波到时的拾取,即震相识别。将地震波拾取结果用于震中位置和发震时刻反演,并与美国地质调查局USGS公布的资料对比。结果发现,两者日本地震震中位置相差约16 km,发震时刻相差约0.7 s;汶川地震震中相差约4 km,发震时刻相差约0.4 s。     

重力仪与地震计记录的地震波信号频谱特征比较——以芦山7.0地震为例

以芦山MS7.0地震的地震波记录为例,比较分析我国12个地震台的同址gPhone重力仪与CTS-1地震计记录的地震波信号的频谱特征。观测结果表明,震中距大于2 000 km时,两种仪器记录的P波波形基本一致,只在幅度上稍有区别。P波的频谱分析结果表明,震中距较小时gPhone重力仪会遗失一些高频信号,震中距大于2 000 km时,两种仪器记录的P波优势频带均为0~0.3 Hz。     

基于GPS、InSAR和强震数据联合反演2017年九寨沟MS7.0地震同震滑动分布

以多源观测为约束,利用近场GPS、InSAR同震位移和强震观测数据联合反演2017年九寨沟M_S7.0地震发震断层的几何形状和破裂模型,并初步分析其地震成因。结果表明,九寨沟地震的发震断层走向154°,倾角80°,以左旋走滑为主,兼有少量正断分量;地震破裂主要集中在2~12 km深度范围,最大滑移量为1.6 m,位于地下6 km深处,断层在近地表没有滑移;地震释放的地震矩为5.59×10¹⁸ Nm,矩震级为M_W6.44。结合余震分布认为,该地震的发震断层为隐伏的虎牙断裂北段。     

利用PPP反演海潮负荷位移参数

选取高纬地区的3个CORS站以及中纬地区的SHAO测站,利用PPP浮点解计算测站未经海潮负荷位移改正的坐标时间序列,采用傅里叶变换方法提取分潮的周期信息,进而反演其中8个主要分潮的海潮负荷位移参数,将其与5个全球海潮模型DTU10、EOT11a、GOT4.8、TPXO.7.2、FES2004计算得到的海洋负荷位移参数比较。结果表明,PPP反演得到的8个主要分潮的海潮负荷位移参数与海潮模型计算结果具有较好的一致性,两者均方根差异为mm级,说明利用PPP反演海潮负荷位移参数是可行的。     

用trackRT进行实时GPS数据处理

利用trackRT实时定位模块对GPS连续站数据进行处理,得到测站的实时动态位移,精度水平方向约为1 cm、垂向约为2 cm;利用其RINEX模拟模块trackRTr对2011年日本Mw9.0地震的1 Hz GPS数据进行模拟实时处理,与事后处理结果较为一致;分析卫星轨道星历和基线长度对解算精度的影响,结果表明基线长度对解算精度影响较大。     

2016年杂多MS6.2地震的InSAR形变及断层滑动模型

利用欧空局Sentinel-1A卫星的升轨观测,基于InSAR技术获取2016-10-17青海杂多地区MS6.2地震的同震形变场,并分别采用最速降法和三角元法反演地震的位错分布模型。结果表明,本次地震发生在杂多上拉秀断裂带,走向北东东,地震破裂以正断倾滑为主,兼具少许左旋走滑运动分量,最大滑移量0.17 m,震源深度为9 km,平均滑动角约为-76°,最大形变中心位于32.85°N、94.86°E附近,地震位错主要分布在深度8~12 km,与余震的深度分布较为一致,反演给出的矩震级为MW5.9,与USGS和CENC测定的矩震级和面波震级基本一致。     

联合高频GNSS和强震仪观测获取同震位移

采用动态精密单点定位技术处理高频GNSS数据获取GNSS位移,利用空间叠加滤波法消除其共模误差的影响,基于Kalman滤波将高频GNSS位移和强震观测加速度融合,利用GNSS位移校正强震仪加速度的基线漂移,得到融合同震位移和速度,建立顾及基线偏差改正的高频GNSS和强震仪组合模型。以2016年意大利北部6.0级地震为例,探讨本文模型的有效性。实验结果表明,本文模型具有高精度、宽频带特点,能准确描述地震引起的永久阶跃和震时地表位移变化。     

GPS实时单站速度位移解算方法及其精度分析

基于载波相位变率测速模型,采用相邻历元差分计算载波相位变率,并利用消电离层组合削弱电离层的影响,结合广播星历,实现单站GPS速度的实时解算,最后将速度在短时间内进行积分,得到测站在该时间内的实时位移。通过静态GPS实测数据和日本2011-03-11大地震期间IGS站观测数据对测站实时速度和位移进行动态解算,并分析解算结果的精度。结果表明,测速精度可以达到mm级,积分得到的测站位移精度可以达到亚dm级甚至cm级。     

2017年四川九寨沟MS7.0地震震前地壳密度变化特征

基于紧凑约束最小二乘优化反演算法,利用2014~2017年地表重力变化观测成果,反演得到2017年九寨沟MS7.0地震前震中及周边地区的地壳三维密度变化。结果表明：1）在震中及周边地区出现大范围的重力变化,并沿张掖-西宁-九寨沟一带形成最大值达9 μGal的重力变化梯度带;2） 地壳密度变化量随深度增加而增大,其中在45 km附近达到最大,为0.2 g/cm3;3） 震中及邻区范围内形成密度变化的正负四象限分布（2014-04~2017-04）,其表现的物质运移方向呈左旋走滑式分布,与震源机制解基本一致。据此推断,在震前产生了与地震孕育相关的深部构造运动,而强震易发生在物质运移汇聚地方。