Regularized inversion for coseismic slip distribution with active constraint balancing

Abstract

Estimating the spatial distribution of coseismic slip is an ill-posed inverse problem, and solutions may be extremely oscillatory due to measurement errors without any constraints on the coseismic slip distribution. In order to obtain stable solution for coseismic slip inversion, regularization method with smoothness-constrained was imposed. Trade-off parameter in regularized inversion, which balances the minimization of the data misfit and model roughness, should be a critical procedure to achieve both resolution and stability. Then, the active constraint balancing approach is adopted, in which the trade-off parameter is regarded as a spatial variable at each model parameter and automatically determined via the model resolution matrix and the spread function. Numerical experiments for a synthetical model indicate that regularized inversion using active constraint balancing approach can provides stable inversion results and have low sensitivity to the knowledge of the exact character of the Gaussian noise. Regularized inversion combined with active constraint balancing approach is conducted on the 2005 Nias earthquake. The released moment based on the estimated coseismic slip distribution is 9.91×10²¹ N·m, which is equivalent to a moment magnitude of 8.6 and almost identical to the value determined by USGS. The inversion results for synthetic coseismic uniform-slip model and the 2005 earthquake show that smoothness-constrained regularized inversion method combined with active constraint balancing approach is effective, and can be reasonable to reconstruct coseismic slip distribution on fault.

Keywords

coseismic slip inversion / trade-off parameter / active constraint balancing / model resolution matrix / spread function

References

[1] MasterlarkT. Finite element model predictions of static deformation from dislocation source in a subduction zone: Sensitivities to homogeneous, isotropic, Poisson-solid, and half-space [J]. Journal of Geophysical Research, 2003, 108(B11): 2540-2556

[2] MaertenF, ResorP, PollardD, MaertenL. Inverting for slip on three-dimensional fault surfaces using angular dislocations [J]. Bulletin of the Seismological Society of America, 2005, 95(5): 1654-1665

[3] WangY-z, ZhuJ-j, OuZ-q, LiZ-w, XingX-m. Coseismic slip distribution of 2009 L’Aquila earthquake derived from InSAR and GPS data [J]. Journal of Central South University, 2012, 19(1): 244-251

[4] FarquharsonC G, OldenburgD W. A comparison of automatic techniques for estimating the regularization parameter in non-linear inverse problems [J]. Geophysical Journal International, 2004, 156(3): 411-425

[5] BarnhartW D, LohmanR B. Automated fault model discretization for inversions for coseismic slip distributions [J]. Journal of Geophysical Research, 2010, 115(10): B10419

[6] XuC J, DengC Y, ZhouL X. Coseismic slip distribution inversion method based on the variance component estimation [J]. Geomatics and Information Science of Wuhan University, 2016, 41137-44

[7] YabukiT, Matsu’UraM. Geodetic data inversion using a Bayesian information criterion for spatial distribution of fault slip [J]. Geophysical Journal International, 1992, 109(2): 363-375

[8] JiangG Y, XuC J, WenY M. Inversion for coseismic slip distribution of the 2010 Mw6.9 Yushu Earthquake from InSAR data using angular dislocations [J]. Geophysical Journal International, 2013, 181(3): 1441-1458

[9] HreinsdóttirS, FreymuellerJ T, FletcherH J, LarsenC F. Coseismic slip distribution of the 2002 Mw7.9 Denali fault earthquake, Alaska, determined from GPS measurements [J]. Geophysical Research Letters, 2003, 30(13): 1670

[10] ChenS-h, ZhuZ-q, LuG-y, CaoS-j. Inversion of gravity gradient tensor based on preconditioned conjugate gradient [J]. Journal of Central South University: Science and Technology, 2013, 442619-625

[11] OkadaY. Surface deformation due to shear and ensile faults in a half-space [J]. Bulletin of the Seismological Society of America, 1985, 75(4): 1135-1154

[12] OkadaY. Internal deformation due to shear and ensile faults in a half-space [J]. Bulletin of the Seismological Society of America, 1992, 82(2): 1018-1040

[13] SunJ J, LiY G. Adaptive Lp inversion for simultaneous recovery of both blocky and smooth features in a geophysical model [J]. Geophysical Journal International, 2014, 197(2): 882-899

[14] TongX-z, LiuJ-x, GuoR-w, LiuH-fei. Three-dimensional magnetotelluric regularized inversion based on smoothness-constrained model [J]. Transactions of Nonferrous Metals Society of China, 2014, 24(2): 509-513

[15] YiM J, KimJ H, ChungS H. Enhancing the resolving power of least-squares inversion with active constraint balancing [J]. Geophysics, 2003, 68(3): 931-941

[16] MenkeWGeophysical data analysis: discrete inverse theory [M], 2012, New York, Academic Press, Inc.

[17] AnM. A simple method for determining the spatial resolution of a general inverse problem [J]. Geophysical Journal International, 2012, 191(2): 849-864

[18] JooY, SeolS J, ByunJ. Acoustic full-waveform inversion of surface seismic data using the Gauss-Newton method with active constraint balancing [J]. Geophysical Prospecting, 2013, 61(S1): 166-182

[19] LeeS K, KimH J, SongY, LeeC K. MT2DinvMatlab-A program in MATLAB and FORTRAN for two-dimensional magnetotelluric inversion [J]. Computers and Geosciences, 2009, 35(8): 1722-1734

[20] ChaveA D, LezaetaP. The statistical distribution of magnetotelluric apparent resistivity and phase [J]. Geophysical Journal International, 2007, 17(1): 127-132

[21] LiD-z, HuangQ-h, ChenX-b. Error effects on magnetotelluric inversion [J]. Chinese Journal of Geophysics, 2009, 52(1): 268-274

[22] KreemerC, BlewittG, MaertenF. Co- and postseismic deformation of the 28 March 2005 Nias Mw8.7 earthquake from continuous GPS data [J]. Geophysical Research Letters, 2006, 33(7): L07307

[23] KoncaA O, HiorleifsdottirV, SongT A. Rupture kinematics of the 2005, Mw 8.6, Nias-Simeulue earthquake from the joint inversion of seismic and geodetic data [J]. Bulletin of the Seismological Society of America, 2007, 971307-322

[24] HsuY J, SimonsM, WilliamsC, CasarottiE. Threedimensional FEM derived elastic Greens’s functions for the coseismic deformation of the 2005 Mw8.7 Nias-Simeulue, Sumatra earthquake [J]. Geochemistry Geophysics Geosystems, 2011, 127Q07013