Rapid and simultaneous estimation of fault slip and heterogeneous lithospheric viscosity from post-seismic deformation

Abstract

Post-seismic deformation is commonly attributed to viscoelastic relaxation and/or afterslip, although discerning between the two driving mechanisms can be difficult. A major complication in modeling post-seismic deformation is that forward models can be computationally expensive, making it difficult to adequately search model space to find the optimal fault slip distribution and lithospheric viscosity structure that can explain observable post-seismic deformation. We propose an inverse method which uses coseismic and early post-seismic deformation to rapidly and simultaneously estimate a fault slip history and an arbitrarily discretized viscosity structure of the lithosphere. Our method is based on an approximation which is applicable to the early post-seismic period and expresses surface deformation resulting from viscoelastic relaxation as a linearized function with respect to lithospheric fluidity. We demonstrate this approximation using two-dimensional earthquake models. We validate the approximation and our inverse method using two three-dimensional synthetic tests. The success of our synthetic tests suggests that our method is capable of distinguishing the mechanisms driving early post-seismic deformation and recovering an effective viscosity structure of the lithosphere.