Preliminary Finite Fault Results for the Jan 22, 2017 Mw 7.9 -6.2669,155.1189 Earthquake (Version 1)


DATA Process and Inversion

We used GSN broadband waveforms downloaded from the NEIC waveform server. We analyzed 37 teleseismic broadband P waveforms and 9 broadband SH waveforms selected based on data quality and azimuthal distribution. Waveforms are first converted to displacement by removing the instrument response and are then used to constrain the slip history using a finite fault inverse algorithm (Ji et al., 2002). We begin modeling using a hypocenter matching or adjusted slightly from the initial NEIC solution (Lon. = 155.1 deg.; Lat. = -6.3 deg., Dep. = 136.0 km), and a fault plane defined using either the rapid W-Phase moment tensor (for near-real time solutions), or the gCMT moment tensor (for historic solutions).


Result

The solutions of two nodal planes explain the data equally well. Both solutions are presented.

Result of Fault Plane 1

Seismic moment Mo=8.79e+27 dyne.cm
Moment Magnitude Mw=7.9
Strike= 135.0 Dip= 55.0 Rake= 100.0




Surface projection of the slip distribution superimposed on GEBCO bathymetry. Red lines indicate major plate boundaries [Bird, 2003]. Gray circles, if present, are aftershock locations, sized by magnitude.


Cross-section of slip distribution



Cross-section of slip distribution. The strike direction is indicated above each fault plane and the hypocenter location is denoted by a star. Slip amplitude is shown in color and the motion direction of the hanging wall relative to the footwall (rake angle) is indicated with arrows. Contours show the rupture initiation time in seconds.


Moment Rate Function




Source time function, describing the rate of moment release with time after earthquake origin, relative to the peak moment rate (listed in the top right corner of the plot).


Comparison of data and synthetic seismograms




Comparison of teleseismic body waves. Data are shown in black and 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-meters. The number above the beginning of each trace is the source azimuth; below is the epicentral distance. Shading describes relative weighting of the waveforms.




Comparison of teleseismic body waves. Data are shown in black and 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-meters. The number above the beginning of each trace is the source azimuth; below is the epicentral distance. Shading describes relative weighting of the waveforms.


Result of Fault Plane 2

Seismic moment Mo=8.78e+27 dyne.cm
Moment Magnitude Mw=7.9
Strike= 298.0 Dip= 37.0 Rake= 76.0




Surface projection of the slip distribution superimposed on GEBCO bathymetry. Red lines indicate major plate boundaries [Bird, 2003]. Gray circles, if present, are aftershock locations (up to 7 days), sized by magnitude.


Cross-section of slip distribution



Cross-section of slip distribution. The strike direction of the fault plane is indicated by the black arrow and the hypocenter location is denoted by the red star. The slip amplitude are showed in color and motion direction of the hanging wall relative to the footwall is indicated by black arrows. Contours show the rupture initiation time in seconds.


Moment Rate Function




Source time function, describing the rate of moment release with time after earthquake origin.


Comparison of data and synthetic seismograms




Comparison of teleseismic body waves. Data are shown in black and 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-meters. The number above the beginning of each trace is the source azimuth; below is the epicentral distance. Shading describes relative weighting of the waveforms.




Comparison of teleseismic body waves. Data are shown in black and 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-meters. The number above the beginning of each trace is the source azimuth; below is the epicentral distance. Shading describes relative weighting of the waveforms.


Scientific Analysis:

Not available yet.

Slip Distribution:

The plots above and a variety of data files for the finite fault solution in different formats can be obtained by clicking on the Downloads tab below.

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.

Ji, C., D. V. Helmberger, D. J. Wald, and K. F. Ma (2003), Slip history and dynamic implications of the 1999 Chi-Chi, Taiwan, earthquake, J Geophys Res-Sol Ea, 108(B9).

Shao, G. F., X. Y. Li, C. Ji, and T. Maeda (2011), Focal mechanism and slip history of the 2011 M-w 9.1 off the Pacific coast of Tohoku Earthquake, constrained with teleseismic body and surface waves, Earth Planets Space, 63(7), 559-564.


Acknowledgement and Contact Information

This work is supported by the National Earthquake Information Center (NEIC) of United States Geological Survey. This web page is built and maintained by Dr. G. Hayes at the NEIC.