基于GPS观测的2015年尼泊尔MS8.1地震多尺度震后形变特征

摘要/Abstract

摘要： 2015年4月25日青藏高原南缘的喜马拉雅地震带发生尼泊尔M_S8.1大地震, 中国藏南区域和尼泊尔境内的GPS连续站均观测到了明显的同震和震后形变。收集中国大陆构造环境监测网络藏南区域GPS基准站和尼泊尔境内的GPS连续观测站数据, 采用统一策略精密解算后获得了测站坐标时序。利用震前两年以上的数据拟合测站长期构造运动和周期性非构造运动, 并从震后的时间序列扣除之后获得了GPS测站的震后形变。采用震后余滑的对数型模型、粘弹性松弛的指数型模型以及二者组合的综合模型等3种震后形变模型拟合了GPS测站观测到的震后形变。采用基于贝叶斯估计的马尔科夫链蒙特卡洛(Markov Chain Monte Carlo, MCMC)算法反演了模型参数。结果显示, 对数型模型加指数型模型的综合模型能够更好地拟合GPS观测到的震后形变, 表明GPS观测到的震后形变至少包含两种不同的变形机制, 具有多尺度特征。指数型模型的特征时间为876.623年, 对应的等效粘滞系数为8.30×10²⁰ Pa·s。对比不同测站两年内所积累的震后形变与同震比值, 相比震中附近的测站, 中国藏南区域GPS测站观测到的震后形变所占比例更高, 可能指示喜马拉雅地震带主要断裂的持续性加载, 其地震危险性值得进一步关注。

关键词: 尼泊尔地震, GPS, 贝叶斯估计, 对数型模型, 指数型模型, 多尺度特征

Abstract: On April 25th 2015, M_S8.1 Nepal earthquake occurred in the Himala earthquake belt that locates in the south margin of the Tibetan Plateau. We analyzed the continuous GPS data from Crustal Movement Observation Network of China and Nepal continuous GPS network. We firstly adopted the multi-scales wavelet method to investigate the pre-seismic effect of tectonic activities on the new set GPS stations, and corrected the coordinate time series of GPS sites. Then,retrieved the postseismic deformations in these time series, and studied the effects of the afterslip based on logarithmic model and visco-ealstic relaxation based on exponential model induced by this event. The parameters of three poseismic deformation models including afterslip, viscoelastic relaxion and the combination of both were inverted by MCMC algorithm of Bayesion frame, showing that observed postseismic deformation involved at least two different mechanisms. We found that the postseismic deformation was featured with variable temporal scales, and controlled by effects of both the afterslip and visco-ealstic relaxation. The charasteristic time of exponential model is 876.623 a, with a viscosity of 8.30×10²⁰ Pa·s. The inverted cumulative deformation shows that this effect is ongoing along the main fault, while the seismic risk in the upper part of unruptured area still deserved further attention.

Key words: 2015 Nepal earthquake, GPS, Bayesian framewor, Logarithmic model, Exponential model, Multi-scale characteristic

中图分类号:

P315.7

苏小宁, 孟国杰. 基于GPS观测的2015年尼泊尔M_S8.1地震多尺度震后形变特征[J]. 地震, 2018, 38(4): 1-13.

SU Xiao-ning, MENG Guo-jie. Multiscale Postseismic Deformation based on cGPS Time Series Following the 2015 Gorkha Nepal M_S8.1 Earthquake[J]. EARTHQUAKE, 2018, 38(4): 1-13.

参考文献

[1] Ader T, Avouac J P, Jing L Z, et al. Convergence rate across the Nepal Himalaya and interseismic coupling on the Main Himalayan Thrust: Implications for seismic hazard[J]. J. Geophys. Res, 2012, 117, B04403.
[2] 邓起东, 程绍平, 马冀, 等. 青藏高原地震活动特征及当前地震活动形式[J]. 地球物理学报, 2014, 57(7): 2025-2042.
[3] 薛艳, 刘杰, 刘双庆, 等. 2015年尼泊尔M_S8.1与M_S7.5级地震活动特征[J]. 科学通报, 2015, 60: 3617-3626.
[4] 赵文津. 尼泊尔大地震发生的构造背景及发展趋势[J]. 科学通报, 2015, 60: 1953-1957.
[5] Galetzka J, Melgar D, Genrich J F. et al. Slip pulse and resonance of the Kathmandu basin during the 2015 Gorkha earthquake, Nepal[J]. Science, 2015, 349(6252): 1091-1095.
[6] Tung S, Masterlark T. Coseismic slip distribution of the 2015 M_W7.8 Gorkha, Nepal, earthquake from joint inversion of GPS and InSAR data for slip within a 3-D heterogeneous Domain[J]. J. Geohpys. Res., 2016, 121(5): 3479-3503.
[7] Zhang H, Van der lee S, Ge Z X. Multi-array rupture imaging of the devastating 2015 gorkha, nepal earthquake sequence[J]. Geophys. Res. Let., 2016, 43(2): 584-591.
[8] 刘刚, 杨少敏, 师宏波, 等. 2015年尼泊尔地震破裂过程的统一模型[J]. 地球物理学报, 2017, 60(7): 2663-2679.
[9] 单新建, 张国宏, 汪驰升, 等. 基于InSAR和GPS观测数据的尼泊尔地震发震断层特征参数联合反演研究[J]. 地球物理学报, 2015, 58(11): 4266-4276.
[10] 熊维, 谭凯, 刘刚, 等. 2015年尼泊尔M_W7.9地震对青藏高原活动断裂同震、震后应力影响[J]. 地球物理学报, 2015, 58(11): 4305-4316.
[11] 苏小宁, 王振, 孟国杰, 等. GPS观测的2015年尼泊尔M_S8.1级地震震前应变积累及同震变形特征[J]. 科学通报, 2015, 60(22): 2115-2123.
[12] 张勇, 许力生, 陈运泰. 2015年尼泊尔M_W7.9地震破裂过程快速反演与初步联合反演[J]. 地球物理学报, 2015, 58(5): 1804-1811.
[13] Jouanne F, Awan A, Madji A, et al. Postseismic deformation in Pakistan after the 8 October 2005 earthquake: Evidence of afterslip along a flat north of the Balakot-Bagh Thrust[J]. J. Geophys. Res., 2011, 116(B7): 86-94.
[14] 孟国杰, 苏小宁, 徐婉桢, 等. 基于GPS观测研究2010年青海玉树M_S7.1地震震后地壳形变特征及其机制[J]. 地球物理学报, 2016, 59(12): 4570-4583.
[15] Azúa B M, Demets C, Masterlark T. Strong interseismic coupling, fault afterslip, and viscoelastic flow before and after the Oct. 9, 1995 Colima-Jalisco earthquake: Continuous GPS measurements from Colima, Mexico[J]. Geophys. Res. Lett., 2002, 29(8): 122-121.
[16] Barbot S, Fialko Y. A unified continuum representation of post-seismic relaxation mechanisms: semi-analytic models of afterslip, poroelastic rebound and viscoelastic flow[J]. Geophys. J. Int., 2010, 182(3): 1124-1140.
[17] Bedford J, Moreno M, Li S, et al. Separating rapid relocking, afterslip, and viscoelastic relaxation: An application of the postseismic straightening method to the Maule 2010 cGPS[J]. J. Geophys. Res., 2016, 121(10): 7618-7638.
[18] Suito H, Freymueller J T. A viscoelastic and afterslip postseismic deformation model for the 1964 Alaska earthquake[J]. J. Geophys. Res., 2009, 114, B11404.
[19] Tobita M. Combined logarithmic and exponential function model for fitting postseismic GNSS time series after 2011 Tohoku-Oki earthquake[J]. Earth Planets Space, 2016, 68: 41.
[20] Sreejith K M, Sunil P S, Agrawal R, et al. Coseismic and early postseismic deformation due to the 25 April 2015 M_W7.8 Gorkha, Nepal, earthquake from InSAR and GPS measurements[J]. Geophys.Res. Lett., 2016, 43: 3160-3168.
[21] Gualandi A, Avouac J P, Galetzka J, et al. Pre- and post-seismic deformation related to the 2015, M_W7.8 Gorkha earthquake, Nepal[J]. Tectonophysics, 2017, 714-715: 90-106.
[22] Dong D N, Bock Y. Global Positioning System Network analysis with phase ambiguity resolution applied to crustal deformation studies in California[J]. J. Geophys. Res., 1989, 94(B4): 3949-3966.
[23] Dieterich J H. Kilgore B D. Direct observation of frictional contects: new insights for state-dependent propertities[J]. Pure and Applied Geophys., 1994, 143(1-3): 283-302.
[24] Fukuda J. Johuson K M. Larson K M. Fault friction parameters inferred from the early stages of afterslip following the 2003 Tokachi-oki earthquake[J]. J. Geophys. Res., 2009, 114, B04412.
[25] Marone C J, Scholz C H, Bilham R. On the mechanics of earthquake afterslip[J]. J. Geophys. Res., 1991, 96(B5): 8441-8452.
[26] Perfettini H, Avouac J P. Modeling afterslip and aftershocks following the 1992 Landers earthquake[J]. J. Geophys. Res., 2007, 112: B07409.
[27] Calais E, Vergnolle M, Deverchere J, et al. Are post-seismic effects of the M8.4 Bolnay earthquake(1905 July 23) still influencing GPS velocities in the Mongolia-Baikai area? [J]. Geophys. J. Int., 2002, 149(1): 157-169.
[28] 石耀霖, 曹建玲. 中国大陆岩石圈等效粘滞系数的计算和讨论[J]. 地学前沿(中国地质大学(北京); 北京大学), 2008, 25(3): 82-95.