A non-linear geodetic data inversion using ABIC for slip distribution on a fault with an unknown dip angle

Abstract

S U M M A R YWe have developed a method of geodetic data inversion for slip distribution on a fault with an unknown dip angle. A common strategy for obtaining slip distribution in previous studies is to first determine the fault geometry by minimizing the square misfit under the assumption of a uniform slip on a rectangular fault, and then apply the usual linear inversion technique to estimate a slip distribution on the determined fault. It is not guaranteed, however, that the fault determined under the assumption of a uniform slip gives the best fault geometry for a spatially variable slip distribution. The inverse problem is non-linear for cases with unknown fault geometries, but the non-linearity of the problems is actually weak, when we can assume the fault surface to be a flat plane. In particular, when a clear trace of coseismic faults is observed on the Earth's surface, only the dip angle is an unknown parameter to determine the fault geometry. Then, we regarded the dip angle as an hyperparameter that prescribed the structure of parametric models, and obtained the best estimate of the dip angle using Akaike's Bayesian Information Criterion (ABIC). With the best estimate of the dip angle, we can obtain the slip distribution on the fault based on the maximum-likelihood principle. We applied the method to the InSAR data of the 1995 Dinar, Turkey earthquake and obtained a much lower dip angle than the previous analyses.

A New Strategy to Compare Inverted Rupture Models Exploiting the Eigenstructure of the Inverse Problem

“…This example highlights the abstract nature of the coimage subspace: a truncated model, although robust, may display source patterns that are not directly amenable for physical interpretation. The enhanced sensitivity of the near-field observations is particularly visible when considering static displacement data, such as Interferometric Synthetic Aperture Radar data (e.g., Fukahata and Wright, 2008).…”

Deformation and Source Parameters of the 2015 Mw 6.5 Earthquake in Pishan, Western China, from Sentinel-1A and ALOS-2 Data

“…Previous studies (e.g., [44,47,48]) indicate that the fault geometry (especially the dip angle) derived from the uniform slip model is not optimal for the distribution slip model. In this study, the grid search method proposed by Feng et al [44] is used to determine the optimal dip and smoothing factor simultaneously.…”

InSAR-derived crustal deformation and fault models of normal faulting earthquake (M j 7.0) in the Fukushima-Hamadori area

“…We incorporate the cross-terms of the covariance matrix in the inversion scheme, following the equation presented by Fukahata and Wright (2008), and now take the characteristic correlation distance of errors to be 10 km (Wright et al, 2003;Fukahata and Wright, 2008). Figure 8 shows the calculated LOS displacement predicted from our model and the residual between the observations and the calculations.…”

“…We cannot neglect the contribution of the cross-terms of the covariance matrix for the InSAR data, because they have a strong spatial correlation in general, which largely results from the variations of atmospheric water vapour (Lohman and Simons, 2005;Fukahata and Wright, 2008). We incorporate the cross-terms of the covariance matrix in the inversion scheme, following the equation presented by Fukahata and Wright (2008), and now take the characteristic correlation distance of errors to be 10 km (Wright et al, 2003;Fukahata and Wright, 2008).…”