skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

Journal Article · · Geophysical Journal International
DOI:https://doi.org/10.1093/gji/ggac050· OSTI ID:1862887

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

Research Organization:
Univ. of Texas, Austin, TX (United States)
Sponsoring Organization:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); National Science Foundation (NSF)
Grant/Contract Number:
SC0019303; EAR-2121666; EAR-2045292; 19214743; 1927216; ACI-1550593
OSTI ID:
1862887
Alternate ID(s):
OSTI ID: 1979491
Journal Information:
Geophysical Journal International, Journal Name: Geophysical Journal International Vol. 230 Journal Issue: 2; ISSN 0956-540X
Publisher:
Oxford University PressCopyright Statement
Country of Publication:
United Kingdom
Language:
English

References (115)

Coseismic Rupture Process of the Large 2019 Ridgecrest Earthquakes From Joint Inversion of Geodetic and Seismological Observations journal November 2019
Inverse Problem Theory and Methods for Model Parameter Estimation book January 2005
Towards adjoint-based inversion for rheological parameters in nonlinear viscous mantle flow journal September 2014
Bayesian Dynamic Finite‐Fault Inversion: 2. Application to the 2016 M w 6.2 Amatrice, Italy, Earthquake journal July 2019
Slow slip near the trench at the Hikurangi subduction zone, New Zealand journal May 2016
A family of higher order mixed finite element methods for plane elasticity journal February 1984
Inferring the basal sliding coefficient field for the Stokes ice sheet model under rheological uncertainty journal January 2021
Dynamically Triggered Changes of Plate Interface Coupling in Southern Cascadia journal November 2019
A scalable algorithm for MAP estimators in Bayesian inverse problems with Besov priors journal January 2015
Pre- and postseismic slow slip surrounding the 2011 Tohoku-oki earthquake rupture journal July 2013
The geophysics, geology and mechanics of slow fault slip journal August 2018
OpenMP: an industry standard API for shared-memory programming journal January 1998
Mixed finite element methods for elliptic problems journal September 1990
Internal deformation due to shear and tensile faults in a half-space journal April 1992
Domain decomposition and multiscale mortar mixed finite element methods for linear elasticity with weak stress symmetry journal November 2019
Variations in earthquake rupture properties along the Gofar transform fault, East Pacific Rise journal April 2012
deal.II---A general-purpose object-oriented finite element library journal August 2007
hIPPYlib: An Extensible Software Framework for Large-Scale Inverse Problems journal October 2018
Frictional Afterslip Following the 2005 Nias-Simeulue Earthquake, Sumatra journal June 2006
Hessian-based model reduction for large-scale systems with initial-condition inputs journal February 2008
Least Squares Methods for Ill-Posed Problems with a Prescribed Bound journal February 1970
An overview of the Trilinos project journal September 2005
Code Verification by the Method of Manufactured Solutions journal November 2001
Performance and Scalability of the Block Low-Rank Multifrontal Factorization on Multicore Architectures
  • Amestoy, Patrick R.; Buttari, Alfredo; L'Excellent, Jean-Yves
  • ACM Transactions on Mathematical Software, Vol. 45, Issue 1 https://doi.org/10.1145/3242094
journal March 2019
Surface deformation due to shear and tensile faults in a half-space journal August 1985
Firedrake: Automating the Finite Element Method by Composing Abstractions journal December 2016
Stress change before and after the 2011 M9 Tohoku-oki earthquake journal December 2018
Tectonic tremor and deep slow slip on the Alpine Fault: ALPINE FAULT TREMOR AND SLIP journal May 2012
Efficient D-Optimal Design of Experiments for Infinite-Dimensional Bayesian Linear Inverse Problems journal January 2018
A domain decomposition approach to implementing fault slip in finite-element models of quasi-static and dynamic crustal deformation: FAULT SLIP VIA DOMAIN DECOMPOSITION journal June 2013
Rupture Process of the 2019 Ridgecrest, California Mw 6.4 Foreshock and Mw 7.1 Earthquake Constrained by Seismic and Geodetic Data journal July 2020
Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities journal September 2009
A decadal-scale deformation transient prior to the 2011 M w 9.0 Tohoku-oki earthquake: MAVROMMATIS ET AL. journal July 2014
A Computational Framework for Infinite-Dimensional Bayesian Inverse Problems Part I: The Linearized Case, with Application to Global Seismic Inversion journal January 2013
Fast Algorithms for Bayesian Uncertainty Quantification in Large-Scale Linear Inverse Problems Based on Low-Rank Partial Hessian Approximations journal January 2011
MUQ: The MIT Uncertainty Quantification Library journal December 2021
On the existence, uniqueness and approximation of saddle-point problems arising from lagrangian multipliers journal January 1974
Decade‐scale decrease in b value prior to the M 9‐class 2011 Tohoku and 2004 Sumatra quakes journal October 2012
Seismic evidence of fluid migration in northeastern Japan after the 2011 Tohoku-Oki earthquake journal June 2021
Three-dimensional viscoelastic finite element model for postseismic deformation of the great 1960 Chile earthquake journal January 2004
A compiler for variational forms journal September 2006
A preconditioner for the equations of linear elasticity discretized by the PEERS element journal June 2004
hIPPYlib: An Extensible Software Framework for Large-Scale Inverse Problems Governed by PDEs: Part I: Deterministic Inversion and Linearized Bayesian Inference journal April 2021
Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet journal September 2015
Unified form language: A domain-specific language for weak formulations of partial differential equations journal February 2014
Geometric MCMC for infinite-dimensional inverse problems journal April 2017
The Mechanics of Earthquakes and Faulting book December 2018
A family of mixed finite elements for the elasticity problem journal August 1988
DOLFIN: Automated finite element computing journal April 2010
Theoretical Displacements and Stresses near Fractures in Rock: with Applications to Faults, Joints, Veins, Dikes, and Solution Surfaces book January 1987
Slow Earthquakes and Nonvolcanic Tremor journal May 2011
Tremor-tide correlations and near-lithostatic pore pressure on the deep San Andreas fault journal December 2009
Episodic slow slip events in the Japan subduction zone before the 2011 Tohoku-Oki earthquake journal July 2013
The use of a priori data to resolve non-uniqueness in linear inversion journal April 1979
Preceding, coseismic, and postseismic slips of the 2011 Tohoku earthquake, Japan: AFTERSLIP OF THE TOHOKU EARTHQUAKE journal July 2012
Reduced symmetry elements in linear elasticity journal January 2009
A Silent Slip Event on the Deeper Cascadia Subduction Interface journal April 2001
A community benchmark for subduction zone modeling journal December 2008
A Stochastic Newton MCMC Method for Large-Scale Statistical Inverse Problems with Application to Seismic Inversion journal January 2012
A new elasticity element made for enforcing weak stress symmetry journal March 2010
Intense foreshocks and a slow slip event preceded the 2014 Iquique Mw 8.1 earthquake journal July 2014
A new family of mixed finite elements in ?3 journal January 1986
Geodynamic inversion to constrain the non-linear rheology of the lithosphere journal July 2015
Geodetically Inferred Locking State of the Cascadia Megathrust Based on a Viscoelastic Earth Model journal September 2018
Mixed finite element methods for linear elasticity with weakly imposed symmetry journal October 2007
FEM-based discretization-invariant MCMC methods for PDE-constrained Bayesian inverse problems journal October 2016
Fluid flow in subduction zones: The role of solid rheology and compaction pressure journal September 2014
A very long-term transient event preceding the 2011 Tohoku earthquake journal January 2015
Automated Derivation of the Adjoint of High-Level Transient Finite Element Programs journal January 2013
A new nonlinear finite fault inversion with three-dimensional Green's functions: Application to the 1989 Loma Prieta, California, earthquake: FAULT INVERSION WITH 3-D GREEN'S FUNCTION journal February 2004
A Randomized Maximum A Posteriori Method for Posterior Sampling of High Dimensional Nonlinear Bayesian Inverse Problems journal January 2018
Two families of mixed finite elements for second order elliptic problems journal June 1985
Terra FERMA : The T ransparent F inite E lement R apid M odel A ssembler for multiphysics problems in E arth sciences journal February 2017
Resolving depth-dependent subduction zone viscosity and afterslip from postseismic displacements following the 2011 Tohoku-oki, Japan earthquake journal February 2017
Geodetic determination of relative plate motion in central California journal February 1973
Mixed Finite Element Methods for Linear Viscoelasticity Using weak Symmetry journal June 2010
A Technique for the Numerical Solution of Certain Integral Equations of the First Kind journal January 1962
PEERS: A new mixed finite element for plane elasticity journal December 1984
Finding Structure with Randomness: Probabilistic Algorithms for Constructing Approximate Matrix Decompositions journal January 2011
Goal-oriented optimal design of experiments for large-scale Bayesian linear inverse problems journal July 2018
Dual hybrid methods for the elasticity and the Stokes problems: a unified approach journal June 1997
Three-Field Block Preconditioners for Models of Coupled Magma/Mantle Dynamics journal January 2015
Benchmarking FEniCS for mantle convection simulations journal January 2013
Spud 1.0: generalising and automating the user interfaces of scientific computer models journal January 2009
jInv--a Flexible Julia Package for PDE Parameter Estimation journal January 2017
Viscoelastic Flow in the Lower Crust after the 1992 Landers, California, Earthquake journal November 1998
A community benchmark for viscoplastic thermal convection in a 2-D square box: BENCHMARK FOR VISCOPLASTIC CONVECTION journal July 2015
A Computational Framework for Infinite-Dimensional Bayesian Inverse Problems, Part II: Stochastic Newton MCMC with Application to Ice Sheet Flow Inverse Problems journal January 2014
A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling journal January 2001
Modeling the poroelastic response to megathrust earthquakes: A look at the 2012 M 7.6 Costa Rican event journal April 2018
Mixed and nonconforming finite element methods : implementation, postprocessing and error estimates journal January 1985
Weakly Imposed Symmetry and Robust Preconditioners for Biot’s Consolidation Model journal July 2017
A Decade of Lessons Learned from the 2011 Tohoku‐Oki Earthquake journal May 2021
A Multipoint Stress Mixed Finite Element Method for Elasticity on Simplicial Grids journal January 2020
A generalized poroelastic model using FEniCS with insights into the Noordbergum effect journal February 2020
A Mixed Finite Element Method for Plasticity Problems with Hardening journal September 1977
How do “ghost transients” from past earthquakes affect GPS slip rate estimates on southern California faults?: “GHOST TRANSIENTS” AND FAULT SLIP RATES journal April 2013
Optimal experimental design for infinite-dimensional Bayesian inverse problems governed by PDEs: a review journal March 2021
Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries book January 1997
Stress before and after the 2011 great Tohoku-oki earthquake and induced earthquakes in inland areas of eastern Japan: STRESS CHANGE BY 2011 TOHOKU EARTHQUAKE journal February 2012
Mechanics of graben formation in crustal rocks: A finite element analysis journal October 1989
Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate journal October 2020
Constraining effective rheology through parallel joint geodynamic inversion journal September 2014
Effects of material property variations on slip estimates for subduction interface slow-slip events: Elastic heterogeneity and SSE inversions journal February 2015
A Fast and Scalable Method for A-Optimal Design of Experiments for Infinite-dimensional Bayesian Nonlinear Inverse Problems journal January 2016
Viscoelastic relaxation following subduction earthquakes and its effects on afterslip determination journal February 2015
Learning physics-based models from data: perspectives from inverse problems and model reduction journal May 2021
A simple and efficient method for introducing faults into finite element computations journal October 1981
On the solvability of incompressible Stokes with viscoplastic rheologies in geodynamics: ISSUES IN VISCOPLASTICITY journal June 2016
Computational Methods for Inverse Problems book January 2002
Three-dimensional FEM derived elastic Green's functions for the coseismic deformation of the 2005 M w 8.7 Nias-Simeulue, Sumatra earthquake: THREE-DIMENSIONAL GREEN'S FUNCTIONS FOR THE EARTHQUAKE journal July 2011
Finite element study of uplift and strain across Vancouver Island journal October 1994
Coseismic deformation due to the 2011 Tohoku-oki earthquake: influence of 3-D elastic structure around Japan journal September 2016
Randomization and Reweighted $\ell_1$-Minimization for A-Optimal Design of Linear Inverse Problems journal January 2020
An Element-Based Preconditioner for Mixed Finite Element Problems journal January 2021

Similar Records

Dynamics, Radiation, and Overall Energy Budget of Earthquake Rupture With Coseismic Off-Fault Damage
Journal Article · Tue Aug 06 00:00:00 EDT 2019 · Journal of Geophysical Research. Solid Earth · OSTI ID:1862887

Effects of Fault Roughness on Coseismic Slip and Earthquake Locations
Journal Article · Thu Jul 18 00:00:00 EDT 2019 · Journal of Geophysical Research. Solid Earth · OSTI ID:1862887

An improved geodetic source model for the 1999 MW 6.3 Chamoli earthquake, India
Journal Article · Tue Feb 16 00:00:00 EST 2016 · Geophysical Journal International · OSTI ID:1862887