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. Issue 8 (27th August 2013)

• Record Type: Journal Article
• Authors: Sun, Jianbao; Shen, Zheng‐Kang; Bürgmann, Roland; Wang, Min; Chen, Lichun; Xu, Xiwei
• Is Part Of: Journal of geophysical research. Volume 118:Issue 8(2013:Aug.)
• Page Start: 4599
• Page End: 4627
• Publication Date: 2013-08-27
• Subjects: three‐step maximum a posteriori probability inversion -- fully Bayesian inversion -- mixed linear‐nonlinear Bayesian inversion -- Ganzi‐Yushu fault -- Yushu earthquake -- coseismic slip; Geomagnetism -- Periodicals; Geochemistry -- Periodicals; Geophysics -- Periodicals; Earth sciences -- Periodicals; 551.1
• DOI: 10.1002/jgrb.50244
• Languages: English
• ISSNs: 2169-9313

Abstract

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