Determination of Smoothing Factor for the Co-seismic Slip Distribution Inversion

doi:10.11947/j.AGCS.2018.20170724

Abstract

Abstract: For the determination of the smoothing factor (also known as the regularization parameter) in the co-seismic slip distribution inversion, the compromise curve between the model roughness and the data fitting residual is generally used to determine (in order to distinguish the method proposed in this paper, the method is called "L curve" according to its shape). Based on the L-curve, the eclectic intersection curve as a new method are proposed to determine the smoothing factor in this paper. The results of the simulated experiment show that the inversion accuracy of the parameters of the seismic slip distribution with the smoothing factor determined by the eclectic intersection curve method is better than that of the L curve method. Moreover, the eclectic intersection curve method and the L curve method are used to determine the smoothing factor of L'Aquila and Taiwan Meinong earthquake slip distribution inversion respectively, and the inversion results are compared and analyzed. The analysis results show that the L'Aquila and Meinong of Taiwan actual earthquake slip distribution results are in the range of other scholars at home and abroad, and compared with the L curve method, the eclectic intersection curve method has advantages of high computation efficiency, no need to depend on data fitting degree and more appropriate of smoothing factor and so on.

Key words: co-seismic slip distribution, smoothing factor, eclectic intersection curve method, L curve method, the L'Aquila earthquake, the Meinong earthquake

CLC Number:

P228

WANG Leyang, ZHAO Xiong. Determination of Smoothing Factor for the Co-seismic Slip Distribution Inversion[J]. Acta Geodaetica et Cartographica Sinica, 2018, 47(12): 1571-1580.

References

[1] 许才军, 邓长勇, 周力璇. 利用方差分量估计的地震同震滑动分布反演[J]. 武汉大学学报(信息科学版), 2016, 41(1):37-44. XU Caijun, DENG Changyong, ZHOU Lixuan. Coseismic Slip Distribution Inversion Method Based on the Variance Component Estimation[J]. Geomatics and Information Science of Wuhan University, 2016, 41(1):37-44.
[2] FUNNING G J, PARSONS B, WRIGHT T J, et al. Surface Displacements and Source Parameters of the 2003 Bam (Iran) Earthquake from Envisat Advanced Synthetic Aperture Radar Imagery[J]. Journal of Geophysical Research:Solid Earth, 2005, 110(B9):B09406.
[3] JIANG Zaisen, WANG Min, WANG Yanzhao, et al. GPS Constrained Coseismic Source and Slip Distribution of the 2013M_w6.6 Lushan, China, Earthquake and Its Tectonic Implications[J]. Geophysical Research Letters, 2014, 41(2):407-413.
[4] 许才军, 刘洋, 温扬茂. 利用GPS资料反演汶川M_w7.9级地震滑动分布[J]. 测绘学报, 2009, 38(3):195-201, 215. XU Caijun, LIU Yang, WEN Yangmao. M_w7.9 Wenchuan Earthquake Slip Distribution Inversion from GPS Measurements[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(3):195-201, 215.
[5] 张国宏, 屈春燕, 宋小刚, 等. 基于InSAR同震形变场反演汶川M_w7.9地震断层滑动分布[J]. 地球物理学报, 2010, 53(2):269-279. ZHANG Guohong, QU Chunyan, SONG Xiaogang, et al. Slip Distribution and Source Parameters Inverted from Co-seismic Deformation Derived by InSAR Technology of Wenchuan M_w7.9 Earthquake[J]. Chinese Journal of Geophysics, 2010, 53(2):269-279.
[6] 温扬茂, 何平, 许才军, 等. 联合Envisat和ALOS卫星影像确定L'Aquila地震震源机制[J]. 地球物理学报, 2012, 55(1):53-65. WEN Yangmao, HE Ping, XU Caijun, et al. Source Parameters of the 2009 L'Aquila Earthquake, Italy from Envisat and ALOS Satellite SAR Images[J]. Chinese Journal of Geophysics, 2012, 55(1):53-65.
[7] 李志才, 张鹏, 金双根, 等. 基于GPS观测数据的汶川地震断层形变反演分析[J]. 测绘学报, 2009, 38(2):108-113, 119. LI Zhicai, ZHANG Peng, JIN Shuanggen, et al. Wenchuan Earthquake Deformation Fault Inversion and Analysis Based on GPS Observations[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(2):108-113, 119.
[8] WANG Leyang, ZHAO Xiong, GAO Hua. A Method for Determining the Regularization Parameter and the Relative Weight Ratio of the Seismic Slip Distribution with Multi-source Data[J]. Journal of Geodynamics, 2018(118):1-10.
[9] 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.
[10] OKADA Y. Internal Deformation due to Shear and Tensile Faults in a Half-space[J]. Bulletin of the Seismological Society of America, 1992, 92(2):1018-1040.
[11] SUN Wenke, OKUBO S, VANÍČEK P. Global Displacements Caused by Point Dislocations in a Realistic Earth Model[J]. Journal of Geophysical Research:Solid Earth, 1996, 101(B4):8561-8577.
[12] 孙建宝, 徐锡伟, 沈正康, 等. 基于线弹性位错模型及干涉雷达同震形变场反演1997年玛尼M_w7.5级地震参数-Ⅰ. 均匀滑动反演[J]. 地球物理学报, 2007, 50(4):1097-1110. SUN Jianbao, XU Xiwei, SHEN Zhengkang, et al. Parameter Inversion of the 1997 Mani Earthquake from InSAR Co-seismic Deformation Field Based on Linear Elastic Dislocation Model-I. Uniform Slip Inversion[J]. Chinese Journal of Geophysics, 2007, 50(4):1097-1110.
[13] 刘洋, 许才军, 温扬茂, 等. 2008年大柴旦M_w6.3级地震的InSAR同震形变观测及断层参数反演[J]. 测绘学报, 2015, 44(11):1202-1209. DOI:10.11947/j.AGCS.2015.20140628. LIU Yang, XU Caijun, WEN Yangmao, et al. The InSAR Coseismic Deformation Observation and Fault Parameter Inversion of the 2008 Dachaidan M_w6.3 Earthquake[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(11):1202-1209. DOI:10.11947/j.AGCS.2015.20140628.
[14] TIKHONOV A N. Regularization of Incorrectly Posed Problems[J]. Soviet Mathematics Doklady, 1963, 4(1):1624-1627.
[15] GOLUB G H, HEATH M, WAHBA G. Generalized Cross-Validation as a Method for Choosing a Good Ridge Parameter[J]. Technometrics, 1979, 21(2):215-223.
[16] 王振杰. 大地测量中不适定问题的正则化解法研究[D]. 武汉:中国科学院测量与地球物理研究所, 2003. WANG Zhenjie. Research on the Regularization Solutions of Ill-posed Problems in Geodesy[D]. Wuhan:Institute of measurement and Geophysics of the Chinese Academy of Sciences, 2003.
[17] 王振杰, 欧吉坤. 用L-曲线法确定岭估计中的岭参数[J]. 武汉大学学报(信息科学版), 2004, 29(3):235-238. WANG Zhenjie, OU Jikun. Determining the Ridge Parameter in a Ridge Estimation Using L-curve Method[J]. Geomatics and Information Science of Wuhan University, 2004, 29(3):235-238.
[18] HANSEN P C. Analysis of Discrete Ill-posed Problems by Means of the L-curve[J]. SIAM Review, 1992, 34(4):561-580.
[19] HANSEN P C, O'LEARY D P. The Use of the L-curve in the Regularization of Discrete Ill-posed Problems[J]. SIAM Journal on Scientific Computing, 1993, 14(6):1487-1503.
[20] FAN Qingbiao, XU Caijun, YI Lei, et al. Implication of Adaptive Smoothness Constraint and Helmert Variance Component Estimation in Seismic Slip Inversion[J]. Journal of Geodesy, 2017, 91(10):1163-1177.
[21] 鲁洋为, 王振杰. 用U曲线法确定岭估计中的岭参数[J]. 导航定位学报, 2015(3):132-134, 138. LU Yangwei, WANG Zhenjie. Determining the Ridge Parameter in Ridge Estimation Using U-curve Method[J]. Journal of Navigation and Positioning, 2015(3):132-134, 138.
[22] 袁强强, 沈焕锋, 李平湘, 等. 自适应正则化多幅影像超分辨率重建[J]. 中国图象图形学报, 2010, 15(12):1720-1727. YUAN Qiangqiang, SHEN Huanfeng, LI Pingxiang, et al. Adaptively Regularized Muti-frame Image Super-resolution Reconstruction[J]. Journal of Image and Graphics, 2010, 15(12):1720-1727.
[23] XU Peiliang. Truncated SVD Methods for Discrete Linear Ill-posed Problems[J]. Geophysical Journal International, 1998, 135(2):505-514.
[24] JONSSON, ZEBKER H, SEGALL P, et al. Fault Slip Distribution of the 1999M_w 7.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.
[25] 王乐洋, 许才军, 张朝玉. 一种确定联合反演中相对权比的两步法[J]. 测绘学报, 2012, 41(1):19-24. WANG Leyang, XU Caijun, ZHANG Chaoyu. A Two-step Method to Determine Relative Weight Ratio Factors in Joint Inversion[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(1):19-24.
[26] WANG Leyang, XU Guangyu. Variance Component Estimation for Partial Errors-in-variables Models[J]. Studia Geophysica et Geodaetica, 2016, 60(1):35-55.
[27] WALTERS R J, ELLIOTT J R, D'AGOSTINO N, et al. The 2009 L'Aquila Earthquake (Central Italy):A Source Mechanism and Implications for Seismic Hazard[J]. Geophysical Research Letters, 2009, 36(17):L17312.
[28] 王乐洋, 李海燕, 温扬茂, 等. 地震同震滑动分布反演的总体最小二乘方法[J]. 测绘学报, 2017, 46(3):307-315. DOI:10.11947/j.AGCS.2017.20160212. WANG Leyang, LI Haiyan, WEN Yangmao, et al. Total Least Squares Method Inversion for Coseismic Slip Distribution[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(3):307-315. DOI:10.11947/j.AGCS.2017.20160212.
[29] 李海燕. 震源位错模型参数反演方法研究[D]. 南昌:东华理工大学, 2016. LI Haiyan. Research on Inversion Method for Hypocenter Parameter Base on Dislocation Model[D]. Nanchang:East China University of Technology, 2016.
[30] ANZIDEI M, BOSCHI E, CANNELLI V, et al. Coseismic Deformation of the Destructive April 6, 2009 L'Aquila Earthquake (Central Italy) from GPS Data[J]. Geophysical Research Letters, 2009, 36(17):L17307.
[31] CHELONI D, D'AGOSTINO N, D'ANASTASIO E, et al. Coseismic and Initial Post-seismic Slip of the 2009M_w6.3 L'Aquila Earthquake, Italy, from GPS Measurements[J]. Geophysical Journal International, 2010, 181(3):1539-1546.
[32] 王乐洋, 高华, 冯光财. 2016年台湾美浓M_w 6.4地震震源参数的InSAR和GPS反演[J]. 地球物理学报, 2017, 60(7):2578-2588. WANG Leyang, GAO Hua, FENG Guangcai. InSAR and GPS Inversion for Source Parameters of the 2016M_w6.4 Meinong, Taiwan Earthquake[J]. Chinese Journal of Geophysics, 2017, 60(7):2578-2588.
[33] FENG Guangcai, LI Zhiwei, SHAN Xinjian, et al. Geodetic Model of the 2015 April 25M_w 7.8 Gorkha Nepal Earthquake and M_w7.3 Aftershock Estimated from InSAR and GPS Data[J]. Geophysical Journal International, 2015, 203(2):896-900.
[34] HUANG M H, TUNG H, FIELDING E J, et al. Multiple Fault Slip Triggered Above the 2016 M_w 6.4 Meinong Earthquake in Taiwan[J]. Geophysical Research Letters, 2016, 43(14):7459-7467.
[35] LEE S J, YEH T Y, LIN Y Y. Anomalously Large Ground Motion in the 2016 ML 6.6 Meinong, Taiwan, Earthquake:A Synergy Effect of Source Rupture and Site Amplification[J]. Seismological Research Letters, 2016, 87(6):1319-1326.