Link to USGS home page.
Earthquake Hazards Program

Finite Fault Model

Magnitude 8.6 NORTHERN SUMATRA, INDONESIA
Monday, March 28, 2005 at 16:09:36 UTC

Preliminary Earthquake Report
U.S. Geological Survey, National Earthquake Information Center
World Data Center for Seismology, Denver

Preliminary Rupture Model
Contributed by Chen Ji, Caltech

DATA Process and Inversion

We used the GSN broadband data downloaded from the IRIS DMC. We analyzed 13 teleseismic P waveforms and 8 SH waveforms selected based upon data quality and azimuthal distribution. Waveforms are first converted to displacement by removing the instrument response and then used to constrain the slip history based on a finite fault inverse algorithm (Ji et al, 2002). We use the hypocenter of the USGS (Lon.=97.01 deg.; Lat.=2.08 deg.). The fault planes are defined by slightly modifying the quick moment tensor solution of the Dr. Polet.

Result

After comparing the waveform fits based on two planes, we find that the nodal plane (strike=325 deg., dip=15 deg.) fits the data better. The seismic moment release based upon this plane is 1.43x10**29 dyne.cm using the PREM earth. The centriod location (Lat.=1.48 deg.; Lon.=97.01 deg. depth =31 km) is consistent with that of Harvard CMT.

Cross-section of slip distribution Cross-section of slip distribution

Figure: The big black arrow shows the fault's strike. The colors show the slip amplitude and white arrows indicate the direction of motion of the hanging wall relative to the footwall. Contours show the rupture initiation time and the red star indicates the hypocenter location.

Comparison of data and synthetic seismograms

Comparison of data and synthetic seismograms

Figure: The Data are shown in black and the synthetic seismograms are plotted in red. Both data and synthetic seismograms are aligned on the P or SH arrivals. The number at the end of each trace is the peak amplitude of the observation in micro-meter. The number above the beginning of each trace is the source azimuth and below it is the epicentral distance.

Surface projection of the slip distribution supperimposed on topography and bathymetry map ETOPO2

Figure: Surface projection of the slip distribution supperimposed on topography and bathymetry map ETOPO2. The ocean is plotted in blue and land is plotted in green. The black line indicates the plate boundary.

CJ's Comments:

The static displacement field on the surface (depth=0.0) predicted by this model is up to 3 m for uplift and 4 m for horizonal motion. The sea floor should suffer larger peak motion considering smaller distances to the fault plane (where the depth is negative). The stars from south to north indicate the epicenters of this event, 2002 Mw=7.2 earthquake, and 2004 Sumatra-Andanman island earthquake, respectively.




Download (Slip Distribution, Not available yet)

References

Ji, C., D.J. Wald, and D.V. Helmberger, Source description of the 1999 Hector Mine, California earthquake; Part I: Wavelet domain inversion theory and resolution analysis, Bull. Seism. Soc. Am., Vol 92, No. 4. pp. 1192-1207, 2002.

Bassin, C., Laske, G. and Masters, G., The Current Limits of Resolution for Surface Wave Tomography in North America, EOS Trans AGU, 81, F897, 2000.


SUBFAULT FORMAT CMTSOLUTION FORMAT
FirstGov button  Take Pride in America button