Crustal Deformation Modeling
|
Elastostatic Deformation in Heterogeneous Media
The deformation due to faulting in the presence of arbitrary spatial variations in elastic properties can in general be evaluated with fully numerical techniques like finite elements or in some simple cases with semi-analytic techniques, like the propagator matrix. Several investigators have proposed to tackle the general case that includes gradual variations as well as lateral and vertical changes in elastic moduli. Here we propose an iterative method that uses equivalent body forces and equivalent surface tractions to represent the effect of inclusions. The formulation for 2-D antiplane deformation can be found in (Barbot et al. 2008) and approach for the general 3-D case in (Barbot et al. 2009). The method has been benchmarked in numerous cases for strike-slip faults with lateral and vertical variations in elastic moduli. The iterative solution can be obtained using a perturbation approach or a successive over-relaxation method. The SOR technique allows one to 1) accelerate convergence in cases of small variations in elastic properties and 2) evaluate deformation in media with large elastic contrasts. However, the convergence of SOR is not always garantied and depends upon the degree of elastic contrast and the size of the inclusions and its proximity from the source.
- Barbot S., Y. Fialko and D. Sandwell, "Three-Dimensional Models of Elasto-Static Deformation in Heterogeneous Media, with Applications to the East California Shear Zone", Geophysical Journal International, in press, 2009
- Barbot S., Y. Fialko and D. Sandwell, "Effect of a compliant fault zone on the inferred earthquake slip distribution", Journal of Geophysical Research, Volume 113, Issue B6, CiteID B06404
Heterogeneous Half Space Deformation Code: STATICI (last released 02/2010)
The source code of STATICI, used to evaluate static deformation due to fault in a heterogeneous half space with spatially-varying elastic properties can be downloaded here
Some programs for post processing and visualization with gmt can be downloaded here
Some compilation options relate to the choice of FFT (fourt, FFTW, SGI FFT), output formats (text and/or GMT) and parallelization method (MPI or OpenMP). These options are explained in the readme file contained in the tar ball. See also some related comments about the homogeneous version STATICU. Instructions and run examples can be found in the readme file.
The latest release includes an additional text output format compatible with GMT, fixes bugs that were asssociated with thrust faults, and provides some post-processing tools for visualization and data conversion.
Some programs for post processing and visualization with gmt can be downloaded here
Some compilation options relate to the choice of FFT (fourt, FFTW, SGI FFT), output formats (text and/or GMT) and parallelization method (MPI or OpenMP). These options are explained in the readme file contained in the tar ball. See also some related comments about the homogeneous version STATICU. Instructions and run examples can be found in the readme file.
The latest release includes an additional text output format compatible with GMT, fixes bugs that were asssociated with thrust faults, and provides some post-processing tools for visualization and data conversion.