2016年新疆阿克陶地震Sentinel-1 InSAR同震形变特征

地震 ›› 2018, Vol. 38 ›› Issue (1): 17-25.

2016年新疆阿克陶地震Sentinel-1 InSAR同震形变特征

宋瑞庆¹, 孟国杰¹, 张奎², 董彦芳¹, 龚发明²

1.中国地震局地震预测重点实验室地震预测研究所, 北京 100036;
2.重庆大学通信工程学院, 重庆 400044

收稿日期:2017-05-08 出版日期:2018-01-31 发布日期:2019-08-14
通讯作者: 孟国杰,研究员。E-mail:mgj@cea-ies.ac.cn
作者简介:宋瑞庆(1992-),男,甘肃武威人,在读硕士研究生,主要从事现代地壳运动与形变研究
基金资助:
国家自然科学基金(41461164004, 41404027, 41304002);国家国际科技合作专项项目
(2015DFR21100)

InSAR Co-seismic Deformation Characteristics of the 2016 Aketao Earthquake using Sentinel-1 Data

SONG Rui-qing¹, MENG Guo-jie¹, ZHANG Kui², DONG Yan-fang¹, GONG Fa-ming²

1.Key Laboratory of Earthquake Prediction, Institute of Earthquake Science, CEA, Beijing 100036, China;
2.College of Communication Engineering, Chongqing University, Chongqing 400044, China

Received:2017-05-08 Online:2018-01-31 Published:2019-08-14

摘要/Abstract

摘要： 2016年11月25日新疆克孜勒苏州阿克陶县发生M_W6.6地震。本文利用合成孔径雷达差分干涉测量技术, 对Sentinel-1卫星获取的升、降轨雷达数据进行了处理, 提取了该次地震的同震形变场, 并结合形变场特征与震源机制解, 采用梯度下降法反演发震断层的滑动分布。结果表明, 升、降轨LOS向同震形变场在发震断层两侧具有明显不同的形变特征, 主要形变区域分布在断层两侧, 升轨LOS向形变量可达-8.2 cm与11.2 cm, 降轨LOS向形变量可达-21.4 cm与13.1 cm; 反演的升、降轨干涉形变场与InSAR测量值之间的残差得到有效控制, 大部分的残差介于±5 cm之间; 断层滑动分布主要集中于沿断层面深约2~18 km处, 最大滑动量位于沿断层面深约7 km处可达0.96 m; 平均滑动角约182.29°, 最大滑动处的滑动角约197.13°, 两个滑动分布中心的滑动角均接近180°, 表明阿克陶地震为一典型的右旋走滑破裂性事件; 当剪切模量取32 Gpa时, 反演的发震断层地震矩M₀可达9.75×10¹⁸, 相当于矩震级M_W6.60, 与地震波形反演结果一致。

关键词: InSAR, 同震形变场, Sentinel-1卫星, 2016年阿克陶地震

Abstract: An M_W6.6 earthquake occurred on November 25, 2016 in Aketao, Xinjang region. By utilizing differential interferometric synthetic aperture radar technique, we retrieved co-seismic deformation field of the Aketao earthquake through processing the ascending and descending data obtained from Sentinel-1. Besides, with the combination of deformation field features and focal mechanism solution, we inverted effectively the seismic slip distribution using steepest descent method. The results demonstrate that the ascending and the descending LOS co-seismic deformation fields have distinct characteristics on both sides of the seismogenic fault. The ascending LOS displacement is -8.2 cm and 11.2 cm, and the descending LOS displacement is -21.4 cm and 13.1 cm. The residuals between the inversed ascending/descending deformation and the InSAR measurements are effectively controlled. Most residuals are within -5 cm to 5 cm. Slips are mainly distributed 2~18 km beneath the fault plane. The maximum slip is about 0.96 m, which is located 7 km beneath the fault plane. The average slip angle is about 182.29°. The maximum slip angle is about 197.13°. The slip angle of two sliding distribution centers is close to 180°, indicating that the Aketao earthquake is a typical dextral strike-slip rupture event. When the shear modulus is set to 32 Gpa, the seismic moment M₀ of the inversed seismogenic fault is up to 9.75×10¹⁸, which is equivalent to moment magnitude M_W6.6 and is mostly identical to the focal mechanism solution.

Key words: InSAR, Co-seismic deformation field, Sentinel-1, The Aketao earthquake

中图分类号:

P315.7

宋瑞庆, 孟国杰, 张奎, 董彦芳, 龚发明. 2016年新疆阿克陶地震Sentinel-1 InSAR同震形变特征[J]. 地震, 2018, 38(1): 17-25.

SONG Rui-qing, MENG Guo-jie, ZHANG Kui, DONG Yan-fang, GONG Fa-ming. InSAR Co-seismic Deformation Characteristics of the 2016 Aketao Earthquake using Sentinel-1 Data[J]. EARTHQUAKE, 2018, 38(1): 17-25.

参考文献

[1] 陈杰, 李涛, 孙建宝, 等. 2016年11月25日新疆阿克陶M_W6.6地震发震构造与地表破裂[J]. 地震地质, 2016, 38(4): 1160-1174.
[2] Robinson A C, Yin A, Manning C E, et al. Tectonic evolution of the northeastern Pamir: Constraints from the northern portion of the Cenozoic Kongur Shan extensional system, western China[J]. Geological Society of America Bulletin, 2004, 116(7): 953.
[3] Strecker M R, Frisch W, Hamburger M W, et al. Quaternary deformation in the Eastern Pamirs, Tadzhikistan and Kyrgyzstan[J]. Tectonics, 1995, 14(5): 1061-1079.
[4] Chevalier M, Li H, Pan J, et al. Fast slip-rate along the northern end of the Karakorum fault system, western Tibet[J]. Geophysical Research Letters, 2011, 38(22): 116-120.
[5] Massonnet D, Rossi M, Carmona C, et al. The displacement field of the Landers earthquake mapped by radar interferometry [J]. Nature, 1993, 364(6433): 138-142.
[6] Peltzer G, Crampé F, King G. Evidence of nonlinear elasticity of the crust from the M_W7.6 Manyi (Tibet) earthquake[J]. Science, 1999, 286(5438): 272-276.
[7] Simons M, Fialko Y, Rivera L. Coseismic deformation from the 1999 M_W7.1 Hector Mine, California, earthquake as inferred from InSAR and GPS observations[J]. Bulletin of the Seismological Society of America, 2002, 92(4): 1390-1402.
[8] 单新建, 宋小刚, 韩宇飞, 等. 汶川M_S8.0地震前InSAR垂直形变场变化特征研究[J]. 地球物理学报, 2009, 52(11): 2739-2745.
[9] 李奇, 孟国杰, 张奎, 等. 基于分区解缠方式研究2010年玉树M_W6.9地震同震形变[J]. 地震, 2015, 35(3): 22-30.
[10] 洪顺英, 刘智荣, 申旭辉, 等. 多视线向InSAR形变场约束的改则地震滑动分布反演[J]. 遥感学报, 2015, 19(2): 288-301.
[11] Wang R, Diao F, Hoechner A. SDM-A geodetic inversion code incorporating with layered crust structure and curved fault geometry[C]. EGU General Assembly Conference, EGU General Assembly Conference Abstracts, 2013.
[12] Rocca F. Perspectives of Sentinel-1 for InSAR applications[C]. Geoscience and Remote Sensing Symposium, IEEE, 2012, 53: 1738-1740.
[13] Prats-Iraola P, Scheiber R, Marotti L, et al. TOPS interferometry with TerraSAR-X[J]. IEEE Transactions on Geoscience and Remote Sensing, 2012, 50(8): 3179-3188.
[14] Scheiber R, Moreira A. Coregistration of Interferometric SAR Images using spectral diversity[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(5): 2179-2191.
[15] Sansosti E, Berardino P, Manunta M, et al. Geometrical SAR image registration[J]. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(10): 2861-2870.
[16] Zhang K, Ge L, Hu Z, et al. Phase unwrapping for very large interferometric data sets[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(10): 4048-4061.
[17] VPU用户手册. http:∥insar.cqu.edu.cn/vpu/user_guide_cn.html.
[18] Okada Y. Surface deformation due to shear and tensile faults in a half-space[J]. Bulletin of the Seismological Society of America, 1985, 75(4): 1135-1154.
[19] Jonsson S, Zebker H, Segall P, et al. Fault slip distribution of the 1999 M_W7.1 Hector Mine, California, earthquake, estimated from satellite radar and GPS measurements[J]. Bulletin of the Seismological Society of America, 2002, 92(4): 1377-1389.