青藏高原东缘地震序列的流变与构造启示：基于2021-2022年三次中强震后早期余震序列与滑移的动力学解析

本研究针对青藏高原东缘2021-2022年三次中强震（Mw 6.1-7.4）的早期余震序列，通过InSAR反演早期震后滑移（afterslip）与高精度地震重定位分析，揭示了余震时空演化受控于同震应力驱动的震后滑移机制，首次量化了走滑断层流变参数（a-b）与复发周期的深度依赖性差异，为理解脆-韧性转换带动态响应及区域地震危险性评估提供了新范式。

　　

青藏高原作为全球最活跃的构造带之一，其东缘复杂的地震活动一直是地球科学界的焦点。2021-2022年间，该区域接连发生三次中强走滑地震：2021年漾濞Mw 6.1、玛多Mw 7.4和2022年门源Mw 6.7地震，这些事件不仅造成显著灾害，更暴露出对余震触发机制、断层流变特性认知的不足。传统研究多依赖地质或大地测量数据估算地震复发周期，但难以捕捉震后短期动力学过程；而余震深度分布的时空演变如何反映深部流变属性，仍是未解之谜。

针对这些问题，成都理工大学地质灾害防治与地质环境保护国家重点实验室的研究团队开展了系统性研究。通过整合卫星雷达干涉测量（InSAR）与高精度地震重定位技术，首次揭示了青藏高原东缘走滑断层系统的震后滑移-余震耦合机制，相关成果发表于《International Journal of Applied Earth Observation and Geoinformation》。

研究采用三大关键技术：1）利用Sentinel-1卫星数据反演早期震后形变场，构建震后滑移模型；2）应用双差定位法（hypoDD）处理台网数据，获得米级精度的余震重定位目录；3）基于速率-状态摩擦定律，建立余震迁移速度与流变参数（a-b）的定量关系。

早期震后形变与滑移模型

InSAR数据显示玛多和门源地震后分别出现最大3 cm和1 cm的雷达视线向位移，反演获得等效矩震级Mw 5.95和Mw 5.7的震后滑移。滑移主要分布在同震破裂下盘（12-26 km和8-14 km深度），与余震形成空间互补格局（图2,5），表明速度强化（VS）区主导震后松弛。

余震序列时空演化

重定位目录揭示三类特征：1）沿断层走向的差异性迁移速度（漾濞和门源SE向5 km/day快于NW向1 km/day，玛达达100 km/day）；2）深度依赖性迁移，浅部速度显著高于深部（图8）；3）脆-韧性过渡带深度随时间快速浅化（漾濞0.36 km/day、玛多0.73 km/day、门源0.09 km/day），反映应变率衰减诱导的流变响应（图6,7）。

流变与构造启示

通过建立余震累积数与震后滑移的物理模型（公式1-3），首次量化了青藏高原走滑断层的深度依赖性参数：1）流变参数A′（0.1-1 MPa）与实验室测量一致；2）摩擦参数（a-b）随深度递减（10^-4-10^-3），符合速率-状态摩擦理论；3）余震推导的复发周期（漾濞124年、玛多383年、门源207年）显著短于地质/大地测量估算值，揭示地震周期中加载速率的非稳态特性（表S4,S5）。

这项研究的意义在于：1）提出余震时空模式可作为深部流变属性的"动态探针"，尤其适用于缺乏震后形变数据的区域；2）证实震后早期滑移主导余震触发，为地震危险性评估提供物理基础；3）建立的跨尺度分析方法（从米级余震定位到千米级滑移反演）可推广至全球类似构造背景区。该成果不仅深化了对青藏高原东缘孕震机制的理解，更为活动断层长期行为预测提供了新思路。