Evaluation of Strong-Motion Parameters in Qom Province in Iran Using Generalized Inversion Method

Document Type : Articles

Authors

1 Nuclear Science and Technology Research Institute, Tehran, Iran

2 International Institute of Earthquake Engineering and Seismology, Tehran, Iran

Abstract

Being located on the seismic belt, Iran has experienced frequent and sometimes intense seismic activities. Many researchers discussed about seismotectonic characteristics of Iran, and based on similarities in seismotectonic parameters, they divided Iran into different tectonic provinces. Most acceptable seismotectonic divisions for Iran is provided by Mirzaie which introduces five tectonic provinces; Alborz-Azerbaijan (North and Northwest of Iran), Kopet Dagh (Northeast of Iran), Zagros (Southwest of Iran), Makran (Southeast of Iran) and Central-East of Iran. In this study, Qom region located in Central-East province is discussed, which is a seismic active region and faced large earthquakes such as Avaj (22/6/2002) with magnitude 6.1 epicentered in 35.71˚N and 49.02˚E, Silakhor(30/3/2006) with magnitude 6.2 epicentered in 33.69˚N and 48.95˚E and Kahak (18/6/2007) with magnitude 5.7 epicentered in 34.54˚N and 50.94˚E, hence it is very important to have estimates of earthquake parameters such as anelastic and geometrical attenuation, average stress parameter, and site amplification for this region, which has high seismicity rates and especially has populated areas with high seismic risk. In this study, generalized inversion method is used which is common approach for deriving the modeling parameters from the empirically observed data. Generalized inversion method introduced by Andrews. Andrews recast the spectral ratios approach into an inverse problem for earthquakes and stations with multiple records. This method has been later improved by other researchers. Ground-motion amplitudes are affected by source, path and site effects. Source term includes the effect of magnitude, stress parameter and the type of faulting on the observed ground-motion amplitude. Path term includes both anelastic and geometrical attenuation. Site term includes the effect of the near surface features like geology and topography on the observe data that consists of site amplification and high frequency attenuation. High frequency attenuation is calculated by estimating kappa using Hough-Anderson method. Convolution of these effects is what stations are recording as earthquake records. Generalized inversion method introduces a way to divide recorded spectrum into mentioned effects. At first, by transforming time-domain data into frequency-domain and then with the aid of natural logarithm operator, a linear equation is achieved. Afterward, we can derive modeling parameters for each of the effective terms.
In this study, for source term, the Brune source model (omega square) is used. To remove the trade-off between anelastic attenuation and geometrical spreading in path effect, geometrical spreading model from previous studies is used. In this formulation, all the site spectra may be multiplied by an arbitrary function of frequency and all the source spectra divided by the same function without changing the fit to the data that adds a degree of freedom to our system, which in this study is resolved by using predefined site term for reference site (Arak station) as it is located on hard rock (Vs30 >700 m/s) and can be evaluated using Boore-Joyner generic rock site model. Database is constructed from 39 number of strong-motion records, recorded in 14 stations from 9 earthquakes with magnitude greater than 4 within target region. At first, horizontal and vertical components are extracted (horizontal component is considered as geometric mean of two horizontal components) and corrected using wavelet de-noising method. Shear wave extracted from all records using Husid method for beginning of S-wave and Kinoshita method for end of S-wave and Fourier spectrum is sampled in frequencies between 0.4 and 15 Hz. A system of over-determined equations is resulted.
Stress drop, corner frequency, moment magnitude and seismic radius is calculated for each earthquake using Grid-search method. The relation between S-wave energy and moment magnitude is expressed by , regional stress drop is 81 bar and frequency dependent quality factor is  that aligns well with previous studies in this region. Calculated quality factor indicates that the region of study has high seismicity rates. Also site amplification is calculated for all the stations in the selected frequency range of 0.4-15 Hz. Results of this study can further be used for a better evaluation of seismic hazard in a region of interest.

Keywords

References

  1. Andrews, D.J. (1986) Objective determination of source parameters and similarity of earthquakes of different size. Earthquake Source Mechanics. American Geophysical Union, 37, 259-267.
  2. Motazedian, D. (2006) Region-specific key seismic parameters for earthquakes in Northern Iran. Bull. Seismol. Soc. Am., 96(4),1383-1395.
  3. Hartzell, S. (1992) Site response estimation from earthquake data. Bull. Seism. Soc. Am., 82(6), 2308-2327.
  4. Berberian, M. (2006) Contribution to the seismotectonics of Iran (part 2). Geological Survey of Iran, 39, 429-517.
  5. Nowroozi, A. (1976) Seismotectonic provinces of Iran. Bull. Seism. Soc. Am., 66(4), 1249-1276.
  6. Khademi, M.R., Nayeri, A. (1976) Seismotectonic provinces of Iran. Iranian Committee on Large Dams; Bulletin 6.
  7. Mirzaei, N., Mengtan, G., Yuntai, C. (1998) Seismic source regionalization for seismic zoning of Iran: Major seismo tectonic provinces. J. Earthquake Pred. Res., 7(4), 465-95.
  8. Babaahmadia, A., Safaei, H., Yassaghi, A., Vafa, H., Naeimi, A., Madanipour, S., Ahmadi, M. (2010) study of Quaternary structures in the Qom region, West Central Iran. Journal of Geodynamics, 50, 355-367.
  9. Ansari, A., Noorzad, A., Zare, M., (2007) Application of wavelet multi-resolution analysis for correction of seismic acceleration records, J. Geophys. Eng., 4(4), 1-16.
  10. Husid, P. (1967) Gravity Effects on the Earthquake Response of Yielding Structures. Report of Earthquake Engineering Research Laboratory, California Institute of Technology, Pasadena, California.
  11. Kinoshita, SH. (1994) Frequency-Dependent Attenuation of Shear Waves in the Crust of the Southern Kanto Area, Japan. Bull. Seism. Soc. Am., 84(5), 1387-1396.
  12. Boore, D.M. (1983) Stochastic simulation of high-frequency ground motions based on seismological models of the radiated spectra. Bull. Seism. Soc. Am., 73(6), 1865-1894.
  13. Anderson, J.G. and Hough, S.E. (1984) A Model for the Shape of the Fourier Amplitude Spectrum of Acceleration at High Frequencies. Bull. Seismol. Soc. Am., 74(5), 1969-1993.
  14. Aki, K. (1967) Scaling law of seismic spectrum. J. Geophys. Res., 72(4), 1217-1231.
  15. Yoshimoto, K., Sato, H., and Ohtake, M. (1993) Frequency-dependent attenuation of P and S waves in the Kanto area, Japan, based on the coda normalization method. Geophys. J. Int., 114, 165-174.
  16. Zafarani, H., Hassani, B. and Ansari, A. (2012) Estimation of earthquake parameters in the Alborz seismic zone, Iran using generalized inversion method. Soil Dynam. Earthq. Eng., 42, 197-218.
  17. Ghasemi, H., Kamalian, N., Hamzehloo, H., Beitollahi, A. (2006) Determination of quality factor of direct shear waves, Qs, in Alborz region using near-field strong motion records of Kojor earthquake in frequency range of 1~32 Hz. Physics of the Earth and Space, 31,103-112.
  18. Mousavi, M., Zafarani, H., Noorzad, A., Ansari, A., Bargi, KH. (2007) Analysis of Iranian strong motion data using the specific barrier model. J. Geophys. Eng., 4(4), 1-14.
  19. Raoof, M., Herrmann, R.B., and Malagnini, L. (1999) Attenuation and excitation of three-component ground motion in southern California. Bull. Seismol. Soc. Am., 89(4), 888-902.
  20. Brune, J.N. (1970) Tectonic stress and the spectra of seismic shear waves from earthquakes. J. Geophys. Res., 75(26), 4997-5009.
  21. Boore, D. (2003) Simulation of ground motion using the stochastic method. Pure Appl. Geophys., 160, 635-76.
  22. Engdahl, E.R., Jackson, J.A., Myers, S.C., Bergman, E.A., Priestley, K. (2006) Relocation and assessment of seismicity in the Iran region. Geophys. J. Int., 167(2), 761-778.

Files

History

  • Receive Date: 18 July 2017
  • Revise Date: 27 February 2021
  • Accept Date: 26 December 2020

Share

How to cite

Statistics