Publications

Export 2 results:
Sort by: Author Title Type [ Year  (Desc)]
2010
Tong, XP, Sandwell DT, Fialko Y.  2010.  Coseismic slip model of the 2008 Wenchuan earthquake derived from joint inversion of interferometric synthetic aperture radar, GPS, and field data. Journal of Geophysical Research-Solid Earth. 115   10.1029/2009jb006625   AbstractWebsite

We derived a coseismic slip model for the M(w) 7.9 2008 Wenchuan earthquake on the basis of radar line-of-sight displacements from ALOS interferograms, GPS vectors, and geological field data. Available interferometric synthetic aperture radar (InSAR) data provided a nearly complete coverage of the surface deformation along both ascending (fine beam mode) and descending orbits (ScanSAR to ScanSAR mode). The earthquake was modeled using four subfaults with variable geometry and dip to capture the simultaneous rupture of both the Beichuan fault and the Pengguan fault. Our model misfits show that the InSAR and GPS data are highly compatible; the combined inversion yields a 93% variance reduction. The best fit model has fault planes that rotate from shallow dip in the south (35 degrees) to nearly vertical dip toward the north (70 degrees). Our rupture model is complex with variations in both depth and rake along two major fault strands. In the southern segment of the Beichuan fault, the slip is mostly thrust (<13 m) and occurred principally in the upper 10 km of the crust; the rupture progressively transformed to right-lateral strike slip as it propagated northeast (with maximum offsets of 7 m). Our model suggests that most of the moment release was limited to the shallow part of the crust (depth less than 10 km). We did not find any "shallow slip deficit" in the slip depth distribution of this mixed mechanism earthquake. Aftershocks were primarily distributed below the section of the fault that ruptured coseismically.

2004
Smith, B, Sandwell D.  2004.  A three-dimensional semianalytic viscoelastic model for time-dependent analyses of the earthquake cycle. Journal of Geophysical Research-Solid Earth. 109   10.1029/2004jb003185   AbstractWebsite

[ 1] 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.