2015年尼泊尔Mw7.9地震震前断层闭锁与同震位移特征分析

对尼泊尔震前的GPS水平速度场进行融合处理,获得跨喜马拉雅中东段沿N15°E方向的GPS水平速度场剖面。采用弹性半空间矩形位错模拟反演,结果表明,尼泊尔地震前喜马拉雅主前缘断裂带存在浅部闭锁,深部无震蠕滑,闭锁深度为22.5km,蠕滑段滑动速率为19.0mm/a,断层倾角为10°。GPS观测的同震形变场揭示了尼泊尔地震引起的地壳形变特征,最大同震位移位于尼泊尔境内的KKN4,该站向南移动1.89m;地震还在我国藏南地区造成最大0.54m的永久形变。距离震中400km以外GPS观测到的同震形变微弱,在误差范围之内。依据GPS同震水平位移在N15°E方向的位移剖面,采用矩形位错模型简单模拟了尼泊尔地震的同震破裂。结果显示,GPS同震位移剖面可以用喜马拉雅冲断带前缘主断裂(MFT)以北的基底低角度逆断层引起的弹性位错模拟,浅部的主前缘逆冲断裂、主边界逆冲断裂和主中央逆冲断裂等分支并没有破裂。     

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

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

汶川地震断层滑动模型反演及地震区域构造的大地测量约束

以汶川地震前后GPS观测得到的同震形变场为约束,利用带约束的最小二乘反演方法,分别反演考虑与不考虑滑脱层的最优同震断层滑动模型,并在此基础上正演了同震重力变化。结果表明,含滑脱层的断层滑动模型模拟的GPS台站的同震水平形变和重力台站的同震重力变化都更接近实际的观测结果,证实龙门山的中央断裂带南段存在大规模的近水平的滑脱层参与了汶川地震的同震滑动。     

基于非均匀地壳模型的格林函数法反演2001年昆仑山口西Ms8.1地震的同震滑移

建立基于地壳不均匀的三维有限元模型,将昆仑山口西Ms8.1地震产生的长约400 km的破裂带在走向上分为16段,生成各个破裂段滑移相对于地表位移的格林函数,并以GPS观测的水平位移为约束,用阻尼最小二乘法反演破裂滑移分布,结果表明反演结果与地质考察的地表破裂数据接近.用反演结果重构的地表位移场与观测值一致,并体现了断裂...     

Sentinel-1 SAR数据约束的2018年谢通门MW5.8地震同震破裂模型

通过Sentinel-1卫星升降轨数据获取谢通门地震的同震形变场,并基于均匀弹性半无限位错模型反演地震的同震滑动分布模型。InSAR同震形变场表明,升降轨视线向最大形变量分别为0.049 m和0.051 m,形变场长轴大致呈南北方向,位于甲岗-定结断裂西侧。通过对倾角和倾向进行格网搜索发现,西倾节面更可能为该地震的发震节面。反演结果表明,滑动分布主要位于2~10 km深度范围内,平均滑动量为0.02 m,最大滑动量为0.10 m,发震断层倾角为47°,平均滑动角为-81.60°,显示该地震以正倾滑动为主。大地测量数据约束的该地震震中为30.27°N、87.75°E,震源深度为6.58 km,释放地震矩为5.056×10¹⁷ Nm,对应矩震级为M_W5.7,与GCMT、USGS公布的震级基本一致。综合分析震中位置和滑动机制认为,甲岗-定结断裂的分支断层为本次谢通门地震的发震断层。     

利用InSAR技术研究2016年西藏定结Mw5.3地震发震构造

利用Sentinel-1A卫星数据和D-InSAR技术，获得2016-05-22西藏日喀则市定结县M_W5.3地震LOS方向的同震形变场图像，分析InSAR形变结果的变化特征，并以此为约束反演该地震的断层几何参数和同震滑动分布特征。结果表明，定结M_W5.3地震发震断层走向近SN，断面倾向E，倾角约43°，破裂长度约9 km，同震滑动主要集中在2~5 km深度范围内，以正断倾滑为主，最大滑动量为0.24 m，矩震级为5.3级。2016年定结地震发震构造是定结-申扎伸展断裂系中的一条新生盲断层。     

九寨沟Ms7.0地震的InSAR观测及断层破裂研究

为提高九寨沟2017-08-08 MS7.0地震的InSAR同震形变场精度，厘清发震断层的构造形态及形成机制，对形变区相干系数进行统计分析并确定相干性阈值，利用GACOS对形变场进行大气校正，再根据余震分布和地质背景确定发震断层的基本形状，最后基于Okada弹性半空间位错模型反演发震断层的滑动分布。改正后的InSAR同震形变场显示，视线向最大下沉为25 cm，最大抬升为10 cm，分别位于震中西北和东南，形变长轴为北西向,主要形变区位于发震断层西部。改正后的InSAR形变场残差均方根较改正前小，最大滑动量为0.9 m，平均滑动角为-0.5°，破裂主要集中在地下1~20 km范围内，矩震级为6.5，与USGS和GCMT等机构的结果一致。研究表明，利用GACOS改正九寨沟InSAR同震形变场对提高形变场精度具有一定的作用，反演断层滑动分布的结果较改正前差别不明显，发震断层的属性与虎牙断裂北段的性质基本一致。结合余震重定位结果可以推断，发震断层为虎牙断裂的北向延伸部分，此次地震事件为巴颜喀拉地块南东向扭转与华南地块碰撞的结果，中下地壳粘性流体的差异分布是导致虎牙断裂倾角变化的主要原因。     

GPS约束下的2022年台湾MW6.9地震破裂滑动分布

利用9个GPS站水平和垂直向形变数据反演2022年台湾M_W6.9地震破裂滑动分布。结果表明,此次地震破裂至地表,滑移以走滑为主兼逆冲分量,主要沿NNE向延伸,有2个破裂集中区,共释放地震矩约5.73×10¹⁹ Nm。基于GPS观测和位错理论模型对此次地震的地表形变进行分析,认为同震形变整体上符合台湾岛东海岸区域构造运动特征,菲律宾海板块俯冲欧亚板块是中央山脉断裂活动的主要动力来源。     

改则Ms6.9地震InSAR同震形变场模拟

以InSAR同震形变场为约束，采用各向同性弹性半空间位错模型，利用试错法对2008年西藏改则县扎西错\%M\%s6.9地震断层几何形状进行格网搜索寻优模拟。观测值离散化得到1 010个数据总集，并重采样80次得到80个数据子集，每个子集之间存在微小的噪声差异，采用Bootstrap方法进行80次搜索寻优。首先用单断裂矩形断层模型进行模拟，最大残差为30 cm。然后用双断裂矩形断层模型进行模拟，最大残差为15 cm。双断裂模型反映了观测形变场的基本特征，并与地震资料符合。根据80次搜索寻优得到主断裂倾滑量95%置信区间为1.6~2.4 m。模拟表明，地表变形是主、余震两次破裂的结果。实际观测和数值模拟结果表明，InSAR可以作为监测青藏高原大陆动力学的有效技术。     

InSAR约束下的2021年西藏比如M_(W)5.7地震断层几何参数及滑动分布北大核心CSCD

利用Sentinal-1A卫星SAR数据获取覆盖2021-03-19西藏比如M_(W)5.7地震的降轨同震视线向形变场,反演此次地震的发震断层几何参数及同震滑动分布。结果表明,此次地震破裂在一个之前未被探测到的隐伏断层上;破裂以正断为主,伴有少量左旋走滑分量;滑动主要分布在2.4~11.2 km深度处,且在6.4 km深度处滑动量达到最大,约为0.17 m;获取的地震矩约为3.6×10^(17) Nm,相当于M_(W)5.7。这些精确的断层参数及震源模型可为评估区域地震危险性及构造演化提供重要参考。     

GPS�ⶨ���Ĵ�«ɽMw6.6����ͬ���α�

??????????????Χ??GPS??????????????з?????????t???????λ?????????α?????50 km??Χ????GPS????α???????????????????70 mm????????????30 s???????λ??????????????????α????????λ???????????У???С?????о???????????????29 km????????λ???21 mm????????150 km???λ???1??2 mm????GPS???????α????????????????????????????????????????????λ?????и???20 km????????廬?????0.6 m??GPS????????α??о????????????????????????????????ε?????????????????     

����GPS��InSAR�����봨Mw7.9����ϲ㻬���ֲ�

??????8.0????????????????????????????????????λ??????????????η????????ν??????????????????8.0??????GPS??InSAR????α??????????????????????????????????????????1?????????????????GPS??InSAR??????????α????????????97.27%;2??????????????????????????????????????0??20 km??????????????λ???????????????????????????????????????????0??20 km???????????????????10 m????????10??15 km?????????????6.5 m??3??????????????????????????????????????????????????????????103??????????????????????121????????????????????138??????????????????????????105??4?????????????????8.7??10 ²⁰ Nm???????Mw7.9??     

远场GPS同震位移的变形特征分析——以2011年日本东北大地震为例

利用2011日本东北大地震的亚洲东部GPS远场水平位移数据,通过计算及比较分析4个独立震源解和倾/滑角效应,研究了远场地震变形的空间分布特征和规律。结果表明,对于低角度逆冲型地震来说,水平位移格林函数在远场位移场中占主导作用。基于此确定了远场GPS水平位移场能在反演中有效利用的空间范围,为联合远场GPS数据和近场GPS数据反演断层滑动分布作理论准备,并为远场数据反演位错Love数（h、l）提供理论依据。     

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

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

基于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,且空间分布更加集中。研究结果表明,组合机制模型在保证了模型拟合精度的基础上,反演结果与应力驱动模型反演结果更接近。     

昆仑山口西Ms8.1地震的地壳变形特征

利用全球定位系统1991－2001年及震后4期GPS观测数据。获得了2001年11月14日昆仑山口西地震（Ms8.1)的同震和震后形变运动图像，跨破裂带GPS网最近两点测得的最大同震形变为1.9m左右，而震后4个月断层蠕滑引起的变形约为80mm，破裂带两侧震后变形幅度具有非对称性。南侧震后变形基本是北侧的2－3倍。研究结果显示破裂带南盘在震后向偏东方向有明显移动，预示本次地震后能量的重新分配与积累，根据近几十年以来东昆仑断裂带的大地震由西向东扩展的特点，未来地震有向东迁移的可能。     

汶川地震前后川滇块体应变与断裂变形特征研究

利用1999—2007和2007—2009年GPS资料,采用块体变形模型和GPS速度剖面研究了川滇地区的分阶段变形特征,结果显示川西北-滇中-滇西南块体主压应变率方向由北向南呈顺时针方向旋转;安宁河断裂、则木河断裂存在剪切应变积累,小江断裂南段以走滑为主、北段以应变积累为主;汶川地震对丽江-小金河断裂中南段、安宁河断裂和金沙江断裂有一定影响,对则木河断裂、小江断裂和红河断裂的影响较小。