Smooth surface reconstruction of earthquake faults from distributed moment-potency-tensor solutions

TL;DR

This paper presents a novel method to reconstruct smooth 3D earthquake fault surfaces from potency tensor data, addressing overdeterminacy issues and improving accuracy over previous methods.

Contribution

It introduces an analytical solution for 3D fault surface reconstruction from potency tensor data, incorporating an a priori constraint to handle overdeterminacy.

Findings

Successfully applied to 2013 Balochistan earthquake data

Reproduces fault surfaces more accurately than previous methods

Demonstrates robustness with noisy synthetic data

Abstract

Earthquake faults as observed by seismic motions primarily manifest as displacement discontinuities within elastic continua. The displacement discontinuity and the surface normal vector (n-vector) of such an idealized earthquake source are measured by the tensor of potency, which is seismic moment normalized by stiffness. This study formulates an inverse problem to reconstruct a smooth 3D fault surface from an areal density field of the potency tensor. Here, the surface is represented by an elevation field, while nodal planes of the potency density represent the surface normal (n-vector) field, reducing the problem to an n-vector-to-elevation transform. Although this transform is a one-to-one mapping in 2D, it becomes overdetermined in 3D because the n-vector has two degrees of freedom while the scalar elevation has only one, admitting no solution in general. This overdeterminacy originates from modeling the potency density, the inelastic strain with six degrees of freedom, as a displacement discontinuity of five degrees of freedom. Whereas this overdeterminacy appears as the violation of the determinant-free constraint in point potency sources, it raises a conflict with the global consistency of the n-vector field in areal potency densities. Recognizing this capacity of the potency density to describe inelastic strain incompatible with displacement discontinuity, we introduce an a priori constraint to define the fault as the smooth surface that best approximates inelastic strain as displacement discontinuity. We derive an analytical solution for this formulation and demonstrate its ability to reproduce 3D surfaces from noisy synthetic n-vectors. We integrate this formula into potency density tensor inversion and apply it to the 2013 Balochistan earthquake. The estimated 3D geometry shows better agreement with observed fault traces than previous quasi-2D methods, validating our proposal.