基于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内可得到稳定的震级。     

尼泊尔MW7.9地震震后断层滑动时空特征分析

利用2015年尼泊尔MW7.9地震震后16个站531 d的连续观测数据,基于弹性位错理论,使用主成分反演方法（PCAIM）分析该地震震后断层滑动的时空变化特性。结果表明：1）震后断层滑动服从对数衰减模式;2）余滑主要分布在主震破裂带北部区域,最大滑动量达到20.6 cm,位于地表以下26.7 km处;3）余滑释放的地震矩为3.36×1020 N·m,对应矩震级约为M_W7.6。     

基于GNSS观测研究2015年尼泊尔MW7.8地震震后地壳形变特征及其机制

收集及处理尼泊尔境内的GPS连续观测站和中国藏南地区的GPS基准站数据,获得2015年尼泊尔MW7.8地震震后3 a的GPS水平形变场。结果显示,尼泊尔地震的震后形变主要分布于尼泊尔北部及中尼边境区域,且东西方向形变较小,南北方向形变较大,整体继续向南运动,最大震后位移约为10.93 cm。采用孔隙弹性回弹模型计算的理论地表位移远小于GPS观测值,无法解释GPS观测到的震后形变。采用震后余滑模型反演的结果表明,震后余滑主要集中在断层的下倾延伸部分,且空间分布较广,余滑释放的地震矩为1.09×1020 Nm。采用PSGRN/PSCMP程序计算粘弹性引起的理论地表形变结果显示,粘弹性松弛模型不能解释近场GPS观测值,但在远场区域的运动方向与GPS观测值一致。采用粘弹性松弛和震后余滑组合机制模型进行反演,余滑释放的地震矩降为1.08×1020 Nm,且空间分布更加集中。研究结果表明,组合机制模型在保证了模型拟合精度的基础上,反演结果与应力驱动模型反演结果更接近。     

2021年青海玛多7.4级地震前应变演化特征分析

计算1999～2016年、2017～2020年青海玛多及其周边地区GNSS速度场,同时采用连续滑移回归方法获取该地区震前面膨胀率、主应变率及最大剪切应变率,分析其动态演变特征。结果表明：1)震前该地区构造运动并未出现明显调整;2)震前震中附近区域面膨胀率下降,说明震中附近区域应变能累积降低;3)震前震中附近区域最大剪切应变率和主应变率有所增强,说明震中附近出现应力调整现象,剪切作用增强。     

汶川地震震后余滑和震后粘弹性松弛数值模拟分析

采用2010～2015年汶川地区GNSS震后形变资料,利用有限元法建立三维震后粘弹性松弛模型,通过二维格网搜索获得龙门山断裂带上盘最佳弹性层厚度和中下地壳最佳粘滞系数,并分析汶川地震震后2～7 a粘弹性松弛影响下的震后形变特征;然后采用2008～2009年GNSS震后形变资料,根据最佳参数建立粘弹性松弛与余滑组合模型,...     

联合光学和雷达遥感数据的2021年玛多MW7.3地震同震形变提取

利用Sentinel-1升降轨数据,基于D-InSAR技术获取2021-05-22青海玛多M_W7.3地震LOS向形变场。由于地震破裂达到地表,沿断层剖面处出现了干涉条纹不连续、破碎重叠等现象,为准确提取地震发震断层的几何展布,联合地震前后的Landsat-8光学影像数据,基于频率域互相关算法提取地震水平向形变场及断层的几何展布特征。结果显示,玛多地震的地表破裂轨迹长达155.6 km,在首端及末端存在分支破裂,地表主体破裂带可分为3段,段间走向差异较大,整体呈NWW向展布。发震断层主要以左旋走滑为主,推测本次地震由昆仑山口-江错断裂控制。     

基于GPS观测值研究日本M_W9.0地震震后形变机制

基于GPS观测的震后水平位移对2011年日本MW9.0地震的震后形变特征进行研究。震后近5a(截至2015-12),震后水平位移累积达到东向60~165cm,南向20~65cm,距离震中较远处的G104、G105及J192站点观测到的震后位移累积变化已超过同震位移,且震后形变还在持续。联合震后余滑和粘弹性位错理论模拟震后形变,利用这2种作用机制对震后GPS水平形变进行解释。研究表明,震后余滑在震后形变最初阶段起主要作用,但随着时间增长而逐渐衰减,粘滞性松弛作用的贡献随着时间增长而变大,GPS观测到的震后地表形变可由这2种机制结合得到较合理的解释。利用震后GPS形变模拟估算地震区域的地幔粘滞系数在1.5×1019 Pa·s量级。     

玛多MW7.3地震同震滑动分布的三维有限元模拟分析

首次构建2021年玛多M_W7.3地震三维有限元模型,分析同震滑动分布特征。首先解算高精度GNSS同震形变观测数据;然后建立玛多地震三维有限元模型,并用弹性半空间Okada模型验证其准确性;最后以40个近场和远场GNSS同震形变观测数据为约束,利用最小二乘法反演玛多地震同震滑动分布模型,从而模拟同震位移。结果表明,玛多地震引起的破裂主要分布在野马滩和黄河乡附近,最大滑动值约为3.4 m。该结果与现场考察结果一致,能够很好地解释GNSS同震形变观测数据。     

基于H-C方法的地震断层面快速识别方法研究——以2021-05-22青海玛多MS7.4地震为例

基于震后各机构给出的玛多M_S7.4地震震源机制解,利用震源-矩心(H-C)方法对地震发震断层进行快速判断,在不需要参考余震情况下,快速识别地震破裂方向。18组测定结果均显示,NWW走向的节面Ⅰ为发震断层面。综合利用震源区地质构造特征、震后余震序列分布、区域构造应力场、烈度等震线以及野外地质考察结果,验证断层测算结果的准确性,最终推断玛多M_S7.4地震断层面为NWW走向的节面Ⅰ、昆仑山口-江错断裂为发震断层,震源机制解显示是一次以左旋走滑为主的地震事件。     

2021年阿拉斯加MW8.2地震震源破裂过程反演

基于GPS近场形变和远场P波,利用有限断层方法反演2021-07-29阿拉斯加M_W8.2地震的震源破裂过程。结果显示,此次地震持续时间约100 s,释放的地震矩约为1.59×10²¹ Nm。断层破裂滑移自震中沿断层向NEE延伸,震中东北侧的断层滑动量较大。基于GPS数据和位错理论模型分析地震的地表形变,同震形变整体上符合海沟特大地震的逆冲断层弹性回跳理论模式。受到断层闭锁的影响,太平洋板块俯冲北美板块,使得上盘北美板块受到长时间持续挤压,地震发生瞬间北美板块转变为拉张松弛状态。     

2012年苏门答腊MW8.6级地震震后效应模拟

基于弹性-粘弹性分层半空间中矩形位错理论和地震波反演的断层模型,结合研究区地壳-上地幔平均波速分层结构,利用PSGRN/PSCMP软件模拟计算2012年苏门答腊MW8.6地震产生的同震效应以及考虑2004年MW9.1地震效应的震后地表形变和重力变化,同时给出震后形变和重力的年变化率。结果表明,同震形变和重力变化显示发震断层左旋走滑特征,其变化主要集中于断层附近;海水质量的重新分布对同震垂直位移和重力变化影响显著;震后50 a年均变化率揭示了地表形变和重力变化的过程, 其图像明显受2004年MW9.1地震震后效应影响,呈现复杂的正负6区域分布;震后50 a近场年均形变量达10 mm/a,年均重力变化量达1.2 μGal/a,而远场年均形变量低于1 mm/a,年均重力变化量低于0.1 μGal /a;水平形变在震后400 a内变化明显,变化率逐步减小,800 a后基本稳定;垂直位移和重力变化在震后100 a内变化尤为显著,变化率快速衰减,400 a后基本稳定不变。     

基于2015年尼泊尔MW7.8地震震后形变探测青藏高原南缘的岩石圈流变结构

基于有限元方法,采用2015年尼泊尔MW7.8地震震后5 a的GPS观测资料约束青藏高原南缘的岩石圈流变结构,利用通过主前缘逆冲断裂带(MFT)间的距离来约束边界位置的垂直边界结构和印度弹性俯冲板片结构探测印度板块和青藏高原板块的边界结构。结果表明,2种结构均能产生与观测数据一致的南-西南水平运动模式,并能较好地解释远场地面位移,但垂直边界结构模型得到的垂直形变偏大,而俯冲板片结构模型能很好地解释中尼边境及以北地区的隆升现象。青藏高原下地壳稳态和瞬态粘滞系数的最优结果分别为1×1018 Pa·s和1×1017 Pa·s。     

2015年尼泊尔MW7.8地震近场GNSS时间序列分析

收集尼泊尔地震近场密集的GNSS测站数据,采用GAMIT/GLOBK软件包解算2011～2017年尼泊尔近场58个GNSS连续观测站和41个震后GNSS流动观测站数据,对获得的数据进行周期性、共模误差以及噪声的时间序列分析。结果显示：1)GRACE改正GNSS时间序列周期项,可使E、N、U三个分量的WRMS平均降低15.52%、26.41%和45.06%;2)针对时间序列中存在大量缺失数据,改进传统的PCA方法,利用第一主成分计算共模误差,3个分量有效率分别达到61.91%、53.91%和53.88%;3)噪声分析表明,尼泊尔地震近场GNSS站点最佳噪声模型为WN+FN,相比单纯的WN模型,WN+FN模型估计震间速度、同震和震后位移误差分别约为其6～9倍、4～7倍和2～4倍。该研究结果对区域坐标框架维护、尼泊尔地震的同震破裂和震后余滑的地球物理模型反演具有重要意义。     

青海玛多MS7.4地震安西站VP宽频带倾斜仪数据特征分析

青海玛多M_S7.4地震前甘肃省安西县地震监测站VP宽频带倾斜仪出现明显异常变化且记录到该地震的同震波形。结合CTS-1地震计进行融合分析可知,2套仪器同震波形的震相到时基本同步,表明2套仪器具有可靠性和一致性,但震相记录特征存在差异。2021-04-17 VP宽频带倾斜仪数据打破原有倾斜趋势,年变向西北方向倾斜,05-13恢复为向西南方向倾斜,震前记录到倾角急剧增大的方位与震中至台站的方位一致,倾斜旋转变化方向与区域构造应力方向一致。在同时段时均值数据时频谱中,异常时段存在（1.085~4.340）×10^-6 Hz(周期为64~128 h)的低频信号。此次地震中2套仪器的优势频率和时频变化特征存在显著差异,这种差异可能与2套仪器不同的机械结构、频段宽度、采样频率有关。     

新疆于田MS7.3地震同震与震后形变机制研究

利用InSAR技术获取2008-03-21新疆于田MS7.3地震的同震和震后形变场。同震分布式滑动反演结果表明,同震断层最大滑动量达5.4 m,主要分布在南部断层的0~5 km深度附近,地震以正断错动为主,兼有左旋走滑分量。震后形变结果表明,发震断层北段两侧存在差异性运动,最大累积差异形变在震后782 d达15 cm。进一步分析表明,震后断层余滑可能是震后形变的主要机制。余滑反演结果表明,震后2 a断层余滑量相对较小,滑移区范围明显减小且均位于浅部区域,北部断层能量释放较彻底,南部仍存有少量能量,整体能量基本释放完。     

利用长期GPS数据研究2008年汶川地震震后形变

采用2008年汶川M_W7.9地震震后7a的GPS连续站和流动站的长期观测资料,基于粘弹性松弛模型,利用PSGRN/PSCMP程序计算汶川震后形变,反演巴颜喀拉块体下的粘弹性参数。其中分层模型考虑巴颜喀拉块体下存在低速带,最优拟合的粘弹性松弛模型显示该低速带粘弹性参数为2.51×10¹⁸ Pa·s,而中下地壳粘弹性参数为3.98×10¹⁸ Pa·s。     

顾及黏弹性松弛效应的2015年尼泊尔MW7.9地震震后早期余滑模型

处理尼泊尔境内及中国藏南的连续GPS数据,获得尼泊尔地震震后1 a的三维形变场。结果表明,震后形变主要发生在尼泊尔北部及中尼边境区域,水平形变以向南运动为主,垂直形变以隆升为主。采用有限断层模型反演的震后余滑分布范围十分广泛,主要集中在同震破裂下倾部分,向下延伸至30 km,平均余滑深度为20 km。震后1 a余滑释放的地震矩为8.8×1019 Nm,等价矩震级为MW7.3,占同震释放总能量的12%。扣除黏弹性松弛效应之后,反演的余滑分布更加集中于同震破裂下倾区域,断层模型底部无余滑,平均余滑深度减小为16 km。经黏弹性松弛效应改正后的余滑模型释放的地震矩为 5.7 ×1019Nm,仅为同震释放能量的8%左右。该余滑模型不但提高了数据的拟合精度,且其分布特征更接近应力驱动的模型。由于震后余滑发生在深部而非浅部未破裂区域,说明断层浅部仍处在闭锁状态,其未来的地震危险性仍需密切关注。     

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

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

InSAR约束下2001年昆仑山MS 8.1地震震后粘弹性松弛效应

根据2001年昆仑山MS8.1地震震后2~6 a的InSAR形变场,采用震后粘弹性松弛模型进行模拟,使用半解析平面分层模型研究青藏高原昆仑山地震震区流变结构。Maxwell半空间粘弹性应力松弛模型模拟结果表明,最优的弹性上地壳厚度为15~20 km,莫霍面深度为70 km,该地区下地壳有效粘滞系数在（1~2.5）×1019 Pas之间,上地幔有效粘滞系数在（1~6.3）×1019 Pas之间。该结果与青藏高原玛尼震区的震后粘弹性松弛反演结果基本一致。     

2008年汶川MW7.9地震库仑应力变化的可靠性研究

基于2008年汶川M_W7.9地震精细的同震滑动分布模型,计算此次地震引起的周边地区库仑应力变化;利用变化系数(CV)评价库仑应力变化的可靠性,并评估接收断层走向、倾角、滑动角、摩擦系数、Skempton系数对库仑应力变化可靠性的影响;根据CV可靠阈值剔除不可靠库仑应力变化,分析可靠库仑应力变化与余震分布之间的关系。结果表明,库仑应力变化在正值区或负值区内部的可靠性较高,在接近正负值交界处库仑应力变化的可靠性急剧下降;接收断层倾角对库仑应力变化可靠性的影响最大;79.3%的余震位于可靠库仑应力变化分布区域,剔除不可靠库仑应力变化及不可靠库仑应力变化分布区域内余震后,余震触发率下降3.1%。