Relative moment tensor inversion method using synthetic data and its application to the aftershocks of the Faryab Earthquake on February 28, 2006 (Mw 6.0)

Document Type : Research Article

Authors

Abstract

Relative Moment Tensor Inversion, RMTI, method is a suitable method to determine the focal mechanism of some close-lying earthquakes (a cluster of earthquakes) recorded by common seismic stations. With enough data on the propagation path, a Green's function is needed in classical source studies. However, the effect of the propagation path is minimized in RMTI applications using relative data. The observed data in an RMTI method consists of the relative amplitudes of the direct P phase and direct S phase on vertical and rotated horizontal components of some close-lying earthquakes recorded by common seismic stations. The effect of the propagation paths from the source-region to a given seismic station is minimized using the relative amplitudes of the corresponding phases (e.g. P-phase) at a given station. The focal mechanisms of the earthquakes are then determined using a linear weighted least-squares approach for the six components of the seismic moment tensor without knowing the complete Green's functions. The calculation is done using only a simple velocity model at the source region. In this research, we have implemented the RMTI method for four groups of synthetic data. Each group consists of 14 events that their corresponding locations and focal-mechanisms are taken according to the events of March 25, 2006 occurred in the Fin region, South of Iran. These events were recorded by Iranian National Broad-Band Seismic Network (INSN) in local distances and their magnitudes (ML) vary between 3.5 and 5.6. In this synthetic data-set, we tried to analyze the effect of a velocity-model and also non-double couple components (CLVD) in the moment tensor inversion using different scenarios. The synthetic data were firstly inverted using a classical moment tensor inversion in a time domain and then were inverted using the RMTI method in both time and frequency domains. The effects of the velocity model and CLVD components in the data were then analyzed using the synthetic data at different frequency bands used in the inverse algorithm. Our results indicate that the RMTI method can be easily implemented to retrieve the focal mechanism of close-lying earthquakes in local distances. Our results also indicate that an unrealistic increase of non-double-couple components causes wrong results, but the crustal model does not cause any significant effect and the effect of the propagation path in a relative inversion is minimized. The aftershock focal mechanisms of the February 28, 2006 earthquake were determined using an RMTI method. The focal mechanisms of the aftershocks were consistent with main shock.
 
 

Keywords


Bouchon, M., 1981, A simple method to calculate Green's functions for elastic layered media: Bull. Seism. Soc. Am., 71)4(,959-971.
Dahm, T., 1996, Relative moment tensor inversion based on ray-theory: theory and synthetic tests: Geophys. J. Int. 124, 245–257.
Dziewonski, A. M., Chou, T. A., and Woodhouse, J. H., 1981, Determination of earthquake source parameters from waveform data for studies of global and regional seismicity: J. Geophys. Res., 86, 2825-52.
Gilbert, F., and Dziewonski, A. M., 1975, An application of normal mode theory to the retrieval of structural parameters and source mechanisms from seismic spectra: Phil. Trans. R. Soc. A 274, 369- 71.
Global CMT Catalog Search, http://www.globalcmt.org/CMTsearch.html
Hatzfeld, D., Tatar, M., Priestley, K., and Ghafory-Ashtyany, M., 2003, Seismological constraints on the crustal structure beneath the Zagros mountain belt (Iran): Geophys. J. Int., 155, 403-410.
Herrmann, R. B., 1996, Computer Programs in Seismology, (last accessed August 2008), Saint Louis University, http://www.eas.slu.edu/People/RBHerrmann/CPS330.html
Incorporated Research Institutions for Seismology (IRIS), http://www.iris.edu.
International Institute of Seismology and Earthquake Engineering (IIEES), http://www.iiees.ac.ir.
ISC Bulletin, http://www.isc.ac.uk/search/bulletin/index.html
Kanamori, H., and Given, J. W., 1981, Use of long period surface waves for rapid determination of earthquake source parameters: Phys. Earth Planet. In. 27, 8-31.
Kennett, B. L. N., and Engdahl, E. R., 1991, Traveltimes for global earthquake location and phase identification: Geophysical Journal International, 122, 429–465.
Klein, F., 2002, Hypoinverse earthquake location program: ftp://ehzftp.wr.usgs.gov/klein/hyp2000
Langston, C. A., and HeImberger, D. V., 1975, A procedure for modelling shallow dislocation sources: Geophys. J. R. Astr. Soc. 42, 117- 30.
Lay, T., and Wallace, T. C., 1995, Modern Global Seismology: Academic Press, New York.
Mendiguren, J. A., 1977, Inversion of surface wave data in source mechanism studies: J. Geophys. Res., 82, 889-94.
Shomali, Z. H., and Roberts, R. G., 2009, An application of Relative Moment Tensor Inversion to the 26 December 2003 Mw 6.6 Iran-Bam Earthquake: Bull. Seism. Soc. Am., 99(1), 159–171
Sokos, E., Zahradník, J., 2008, ISOLA a Fortran code and a Matlab GUI to perform multiple-point source inversion of seismic data: Computers and Geosciences, 34, 967-977.
Strelitz, R. A., 1978, Moment tensor inversions and source models: Geophys. J. Roy. Astron. Soc., 52, 359-64.
Stump, B. W., and Johnson, L. R., 1977, The determination of source properties by the linear inversion of seismograms, Bull. Seism. Soc. Am., 67, 1489-502.
Zahradnik, J., Jansky, J., and Papatsimpa, N., 2001, Focal mechanisms of weak earthquakes from amplitude spectra and polarities: Pure Appl. Geophys., 158, 647-665.