The empirical Green’s function method is applied to the foreshock and the mainshock of the 2016 Gyeongju earthquake to simulate strong ground motions of the mainshock and scenario earthquake at seismic stations of seven metropolises in South Korea, respectively. To identify the applicability of the method in advance, the mainshock is simulated, assuming the foreshock as the empirical Green’s function. As a result of the simulation, the overall shape, the amplitude of PGA, and the duration and response spectra of the simulated seismic waveforms are similar with those of the observed seismic waveforms. Based on this result, a scenario earthquake on the causative fault of Gyeongju earthquake with a moment magnitude 6.5 is simulated, assuming that the mainshock serves as the empirical Green’s function. As a result, the amplitude of PGA and the duration of simulated seismic waveforms are significantly increased and extended, and the spectral amplitude of the low frequency band is relatively increased compared with that of the high frequency band. If the empirical Green’s function method is applied to several recent well-recorded moderate earthquakes, the simulated seismic waveforms can be used as not only input data for developing ground motion prediction equations, but also input data for creating the design response spectra of major facilities in South Korea.