Hybrid broadband simulation of strong-motion records from the September 16, 1978, Tabas, Iran, earthquake (M w 7.4)

This paper presents a simulation of three components of near-field ground shaking recorded during the main shock at three stations of the September 16, 1978, Tabas (M w = 7.4), Iran, earthquake, close to the causative fault. A hybrid method composed of a discrete wavenumber method developed by Bouchon (Bouchon in Bull Seismol Soc Am 71:959–971, 1981; Cotton and Coutant in Geophys J Int 128:676–688, 1997) and a stochastic finite-fault modeling based on a dynamic corner frequency proposed by Motazedian and Atkinson (Bull Seismol Soc Am 95:995–1010, 2005), modified by Assatourians and Atkinson (Bull Seismol Soc Am 97:935–1949, 2007), is used for generating the seismograms at low (0.1–1.0 Hz) and high frequencies (1.0–20.0 Hz), respectively. The results are validated by comparing the simulated peak acceleration, peak velocity, peak displacement, Arias intensity, the integral of velocity squared, Fourier spectrum and acceleration response spectrum on a frequency-by-frequency basis, the shape of the normalized integrals of acceleration and velocity squared, and the cross-correlation with the observed time-series data. Each characteristic is compared on a scale from 0 to 10, with 10 being perfect agreement. Also, the results are validated by comparing the simulated ground motions with the modified Mercalli intensity observations reported by reconnaissance teams and showed reasonable agreement. The results of the present study imply that the damage distribution pattern of the 1978 Tabas earthquake can be explained by the source directivity effect.