 |
|
|
|
 |
 |
 |
|
 |
|
 |
|
 |
|
 |
|
 |
|
 |
|
 |
|
 |
|
 |
|
 |
|
>Return to Research
Geodynamic modeling codes
15 Coseismic Deformation from the 1999 Mw 7.1 Hector Mine, California, Earthquake as Inferred from InSAR and GPS Observations, M. Simons, Y. Fialko, and L. Rivera, Bull. Seismol. Soc. Am., 92, 1390-1402, 2002. [PDF]
We use a variety of codes to model lithospheric deformation. Generally, the quick and dirty calculations adopt Cartesian elastic half space models with simple sources constructed by superposition of dislocations. These codes are all from the literature and usually implemented in Matlab. For static elastic problems, with Luis Rivera, we developed a 1D layered elastic code (15) which we refer to as EDKS (elementary displacement Kernel Structure). This code is a Fortran 90 code with an available Matlab front end. This software was designed for computational efficiency since we use it for inversions with InSAR data. Generally, Green’s functions a pre-calculated for a given layered structure. Then, when they are actually used, we simply interpolate for the exact source and receiver locations from a table. These codes are available upon request.
Recently, with Charles Williams (RPI) we are focusing on finite element based approaches for handling the effects of 3D elastic structure (such as in the case of the Chi Chi earthquake in Taiwan or most subduction zones) as well as post-seismic inelastic deformation. This effort is an integral part of the larger Computational Infrastructure for Geodynamics (CIG) effort. For more information on this effort, please consult the CIG – short term deformation group web pages.
Mark Simons' Paper Collection: Entire paper including figures are all made available online (within the bounds of copyright restrictions).