Vertical Cross-Section of Subduction Zone (2D)#
The files are in the directory examples/subduction-2d
.
The files and directories for this set of examples includes:
README.md
README file containing a brief description of the various examples.
*.cfg
PyLith parameter files.
generate_gmsh.py
Python script to generate mesh using Gmsh.
*.msh
Gmsh mesh files generated by Gmsh.
*.jou
Files used to construct the finite-element mesh using CUBIT/Trelis.
*.exo
Exodus II mesh files generated by Cubit.
*.spatialdb
Spatial database filesFiles associated with the spatial databases.
viz
Directory containing ParaView Python scripts and other files for visualizing results.
output
Directory containing simulation output. It is created automatically when running the simulations.
Overview#
This example examines quasistatic interseismic and coseismic deformation in 2D for a subduction zone (see Fig. 85). It is based on the 2011 M9.0 Tohoku-oki earthquake off the east coast of Japan. We build on the previous examples and add complexity through a series of steps:
- Step 1
Static coseismic slip on the subduction interface.
- Step 2
Quasistatic interseismic deformation with creep on the top and bottom of the slab, except in the zone of coseismic slip.
- Step 3
Quasistatic earthquake cycle with prescribed coseismic slip and creep.
Warning
Steps 4-6 have not yet been updated for PyLith v3.
Important
We decribe how to generate the finite-element mesh using both Gmsh and Cubit. The files for both methods are included. We use the Gmsh files in the PyLith parameter files. See examples/strikeslip-2d/step01_slip_cubit.cfg for a description of how to modify the parameter files to switch from using mesh files from Gmsh to mesh files from Cubit.