Open AccessA three‐dimensional semianalytic viscoelastic model for time‐dependent analyses of the earthquake cycle
Author(s) -
Smith Bridget,
Sandwell David
Publication year - 2004
Publication title -
journal of geophysical research: solid earth
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.67
H-Index - 298
eISSN - 2156-2202
pISSN - 0148-0227
DOI - 10.1029/2004jb003185
Subject(s) - viscoelasticity , geology , convolution (computer science) , fault (geology) , asthenosphere , computation , fourier transform , time domain , computer science , seismology , tectonics , lithosphere , mathematical analysis , algorithm , physics , mathematics , machine learning , artificial neural network , computer vision , thermodynamics
Exploring the earthquake cycle for large, complex tectonic boundaries that deform over thousands of years requires the development of sophisticated and efficient models. In this paper we introduce a semianalytic three‐dimensional (3‐D) linear viscoelastic Maxwell model that is developed in the Fourier domain to exploit the computational advantages of the convolution theorem. A new aspect of this model is an analytic solution for the surface loading of an elastic plate overlying a viscoelastic half‐space. When fully implemented, the model simulates (1) interseismic stress accumulation on the upper locked portion of faults, (2) repeated earthquakes on prescribed fault segments, and (3) the viscoelastic response of the asthenosphere beneath the plate following episodic ruptures. We verify both the analytic solution and computer code through a variety of 2‐D and 3‐D tests and examples. On the basis of the methodology presented here, it is now possible to explore thousands of years of the earthquake cycle along geometrically complex 3‐D fault systems.