A three‐step maximum a posteriori probability method for InSAR data inversion of coseismic rupture with application to the 14 April 2010 Mw 6.9 Yushu, China, earthquake

• DOI: 10.1002/jgrb.50244
• Publication Date: 2013-06-11T07:12:54Z
• Authors: Jianbao Sun; Zheng‐Kang Shen; Roland Bürgmann; Min Wang; Lichun Chen; Xiwei Xu
• Keywords: Synthetic data; Magnetic dip

Abstract We develop a three‐step maximum posteriori probability method for coseismic rupture inversion, which aims at maximizing the posterior density function (PDF) of elastic deformation solutions earthquake rupture. The originates from fully Bayesian inversion and mixed linear‐nonlinear methods shares same PDF with them, while overcoming difficulties convergence when large numbers low‐quality data are used greatly improving rate using optimization procedures. A highly efficient global algorithm, adaptive simulated annealing, is to search (“ mode ” in statistics) first step. second step approaches “true” solution further Monte Carlo technique positivity constraints, all parameters obtained as initial solution. Then slip artifacts eliminated models third procedure step, fixed fault geometry parameters. design model 45° dip angle oblique slip, produce corresponding synthetic interferometric aperture radar (InSAR) sets validate reliability efficiency new method. then apply this InSAR distribution 14 April 2010 M w 6.9 Yushu, China earthquake. Our preferred composed three segments most occurring within 15 km depth reaches 1.38 m surface. seismic moment released estimated be 2.32e+19 Nm, consistent estimate 2.50e+19 Nm.