A Caltech Library Service

Gamra: Simple meshing for complex earthquakes

Landry, Walter and Barbot, Sylvain (2016) Gamra: Simple meshing for complex earthquakes. Computers and Geosciences, 90 . pp. 49-63. ISSN 0098-3004.

[img] PDF - Published Version
Creative Commons Attribution.

[img] PDF - Erratum
Creative Commons Attribution.

[img] Archive (ZIP) (Application 1) - Supplemental Material
Creative Commons Attribution.

[img] PDF - Submitted Version
See Usage Policy.


Use this Persistent URL to link to this item:


The static offsets caused by earthquakes are well described by elastostatic models with a discontinuity in the displacement along the fault. A traditional approach to model this discontinuity is to align the numerical mesh with the fault and solve the equations using finite elements. However, this distorted mesh can be difficult to generate and update. We present a new numerical method, inspired by the Immersed Interface Method (Leveque and Li, 1994), for solving the elastostatic equations with embedded discontinuities. This method has been carefully designed so that it can be used on parallel machines on an adapted finite difference grid. We have implemented this method in Gamra, a new code for earth modeling. We demonstrate the correctness of the method with analytic tests, and we demonstrate its practical performance by solving a realistic earthquake model to extremely high precision.

Item Type:Article
Related URLs:
URLURL TypeDescription Paper
Barbot, Sylvain0000-0003-4257-7409
Alternate Title:Gamra: Simple Meshes for Complex Earthquakes
Additional Information:© 2016 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license( Received 5 October 2015, Revised 18 February 2016, Accepted 19 February 2016, Available online 23 February 2016. We thank two anonymous reviewers for their comments that improved this paper. This research used resources of the Argonne Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract DE-AC02-06CH11357. This research was supported by the Gordon and Betty Moore Foundation, the National Research Foundation of Singapore under the NRF Fellowship scheme (National Research Fellow Award No. NRF-NRFF2013-04) and by the Earth Observatory of Singapore and the National Research Foundation and the Singapore Ministry of Education under the Research Centres of Excellence initiative. This is EOS publication 113.
Errata:The authors would like to acknowledge two more important sources of support for this research. This research was supported by NSF Award EAR-0949446 and the Tectonics Observatory. The authors would like to apologise for any inconvenience caused.
Group:Infrared Processing and Analysis Center (IPAC)
Funding AgencyGrant Number
Department of Energy (DOE)DE-AC02-06CH11357
Gordon and Betty Moore FoundationUNSPECIFIED
National Research Foundation (Singapore)NRF-NRFF2013-04
Earth Observatory of SingaporeUNSPECIFIED
Ministry of Education (Singapore)UNSPECIFIED
Caltech Tectonics ObservatoryUNSPECIFIED
Subject Keywords:Immersed Interface Method; Earthquakes; Fault systems; Elastostatics; Adaptive mesh refinement; Parallel computing
Other Numbering System:
Other Numbering System NameOther Numbering System ID
Earth Observatory of Singapore113
Record Number:CaltechAUTHORS:20160531-100116966
Persistent URL:
Official Citation:Walter Landry, Sylvain Barbot, Gamra: Simple meshing for complex earthquakes, Computers & Geosciences, Volume 90, Part A, May 2016, Pages 49-63, ISSN 0098-3004, (
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:67489
Deposited By: Tony Diaz
Deposited On:31 May 2016 17:08
Last Modified:18 Nov 2020 18:55

Repository Staff Only: item control page