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. 87). 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.