巴颜喀拉块体东北缘GPS应变率空间分布特征及其与2017年九寨沟MS7.0地震的关系

巴颜喀拉块体东北缘是构造变形和地震活动较强的区域, 2017年九寨沟M_S7.0地震就发生在该区域内。 利用多尺度球面小波方法解算GPS应变率场, 分析巴颜喀拉块体东北缘2009年至2017年的应变率场分布特征, 该方法的优点是可以将GPS应变率场按照不同的空间尺度进行小波分解, 给出不同空间尺度的应变率场。 结果表明在2017年九寨沟地震之前, 震中附近应变积累显著, 虎牙断裂北延断裂的左旋走滑速率为3.0 mm/a, 拉张速率为3.1 mm/a, 表明该条断裂以左旋走滑为主兼有拉张特征, 与九寨沟地震的震源机制解一致。 除九寨沟震中附近外, 在岷县与漳县交界处、 理县和汶川、 青川等地区主应变率、 面应变率、 最大剪应变率也较大, 这可能与2013岷县漳县(M_S6.6)、 2008年汶川(M_S8.0)、 2014年理县(M_S4.8)以及2014青川县(M_S4.8)地震有关。     

2020年土耳其MW5.7地震InSAR同震形变特征与破裂滑动分布

为准确认识2020-06-14土耳其MW5.7地震的发震位置、构造特点以及地震危险性,利用D-InSAR技术对Sentinel-1A数据进行处理,基于GACOS进行大气校正,获得视线向（line of sight,LOS）同震形变场。降轨LOS向同震形变场显示,断层北侧抬升,最大抬升形变量约8.87 cm;南侧沉降,最大沉降量约-7.75 cm。以LOS向同震形变为约束,先采用贝叶斯自举优化法反演发震断层几何参数,然后使用有限断层方法反演地震破裂滑动分布。结果显示,断层走向约257.48°±0.65°、倾角约79.69°±0.98°、滑动角约154.2°±3.8°,震中位置为（40.754°E,39.389°N）,破裂区域长度约8 km、宽度约6 km,破裂的最浅埋深约0.8 km、最大埋深约8.9 km,最大滑动量约0.57 m,对应深度约4.278 km。地震释放地震矩约4.54×1017 Nm,对应矩震级MW5.7,与土耳其灾害和紧急情况管理局公布的结果一致。此次土耳其地震受近东西向的潜伏断层控制,以右旋走滑为主兼具少量的逆冲特性。地震造成部分地区的库伦破裂应力增量超过0.1 bar,这些区域未来的地震危险性值得关注。     

2020年西藏尼玛Mw6.3地震InSAR同震形变特征与破裂滑动分布

2020年7月22日(UTC)西藏尼玛发生Mw6.3地震.为准确确定地震的发震断层、反演破裂滑动分布以及评估地震危险性,本文利用合成孔径雷达差分干涉(DInSAR)技术对Sentinel-1A/B数据进行处理,获取视线向(line of sight,LOS)同震形变场,形变场显示以沉降为主,长轴呈NE向,升、降轨最大沉降量分别为～30 cm和～26 cm.采用基于序列蒙特卡罗采样的贝叶斯方法反演发震断层的位置和几何参数,并在此基础上通过最速下降法反演同震破裂滑动分布.反演结果显示:发震断层走向～30.41°±1.25°,倾角～49.52°±1.06°,属东倾断层.震中位置为(86.89°E,33.18°N),对应深度7.23 km.同震破裂只存在一个滑动中心,主要集中在3.4～11.8 km深度,最大滑动量～0.98 m,平均滑动角-74.54°,表明本次地震是以正断为主兼具少量左旋走滑,释放地震矩约为3.86×1018 N·m,对应矩震级Mw6.36,略大于USGS(United States Geological Survey,美国地质调查局)提供的震级.对于深度5 km和10 km位于破裂区南北端,以及深度15 km和20 km位于破裂区东西侧,其库仑破裂应力变化≥0.01 MPa,未来危险性值得更多的关注.     

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

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

2022年门源6.9级地震的GPS形变特征

利用1999～2021年的GPS观测资料研究2022年门源6.9级地震震中附近区域的地壳形变特征。速度场在祁连山-海原地震带两侧呈现出明显的“南大北小”特性。发震断层速度剖面结果显示,断层两侧挤压速率为3.1 mm/a,走滑速率为3.9 mm/a,这与走滑为主的震源机制解一致。门源地震位于主应变率场和面膨胀率场高值区向低值区过渡的梯度带上,发生在最大剪应变率场和应变率第二张量不变量高值区边缘地区。同震库仑应力触发了大部分的余震,对冷龙岭断裂到海原断裂的地区有一定的应力加载作用。综合地壳形变特征和库仑应力等分析结果认为,冷龙岭断裂到海原断裂地区的地震危险性有所提高。     

2021年双湖MS5.8地震InSAR同震形变场及断层滑动分布反演

基于Sentinel-1A数据,利用合成孔径雷达差分干涉(DInSAR)技术提取视线向(LOS)同震形变场,利用基于序列蒙特卡罗采样的贝叶斯方法反演发震断层几何参数与断层滑动分布。结果表明,发震断层走向约47.810°±2.218°,倾角约36.664°±2.499°;震中位置为87.715°E、34.365°N,震源深度为5.673 km。同震滑动分布主要集中在3.9~7.5 km深度,最大滑动量约0.42 m,平均滑动角约-66.598°±3.258°。本次地震以正断为主,兼具少量左旋走滑性质,反演得到的矩震级为M_W5.5,略小于USGS提供的震级。     

基于GPS观测研究中国东北地区现今地壳形变特征

基于中国东北和俄罗斯远东东南部2012—2017年的GPS观测数据, 利用包含年周期、 半年周期、 线性项和阶跃项的函数模型拟合GPS站坐标时间序列, 得到ITRF2014下的速度场, 并进一步转换到欧亚参考框架下得到相对欧亚板块的速度场。 基于多尺度球面小波方法解算应变率场, 并分析了其空间分布特征, 同时研究了各GPS站对2011年日本东北M_W9.0大地震的震后松弛响应特征和背景形变场特征。 结果表明: ① 若不扣除日本东北大地震的松弛效应, 相对欧亚板块中国东北主体上表现为东南方向运动, 在依兰—伊通断裂和嫩江断裂带之间, 地壳表现为逆时针旋转, 其他区域向东南方向运动, 方向一致性较好, 在敦化—密山断裂东侧速度大小明显增加。 敦化—密山断裂和依兰—伊通断裂两侧拉张量分别为3.96±0.04 mm/a和0.71±0.05 mm/a, 两条断裂的剪切运动不明显。 总体上, 面应变率显示出NW—SE向的拉张和NE—SW向的挤压, 面应变率显示出依兰—伊通断裂南端、 嫩江断裂带北端和俄罗斯远东东南部呈挤压状态。 在依兰—伊通断裂、 敦化—密山断裂南侧以及俄罗斯远东东南部最大剪应变率相对较大。 ② 各GPS测站对2011年日本东北M_W9.0大地震震后松弛的响应整体上表现为东南向运动, 松弛形变量随震中距增加而减小。 松弛效应的面应变率总体上表现为NW—SE向的拉张和NE—SW向的挤压, 面应变率显示出依兰—伊通、 敦化—密山断裂南端、 嫩江断裂带北端以及俄罗斯远东地区具有挤压特征, 其他地区表现为拉张特征。 中国与俄罗斯远东边界南端存在一个明显的最大剪应变率高值区。 ③ 扣除日本东北M_W9.0大地震引起的松弛变形后, 总体上面应变率仍然表现为NW—SE向的拉张和NE—SW向的挤压, 面应变率最大值仍然位于依兰—伊通断裂和敦化—密山断裂南端、 第二松花江断裂带以及俄罗斯远东和中国边界最南段。 在依兰—伊通断裂、 敦化—密山断裂南端, 中国与俄罗斯远东边界南端的最大剪应变率高值区仍然存在, 表明这些地区应变积累较快, 并且一直在持续。