This study sought to simulate strong broadband seismic motions beyond the corner frequencies used, for the same events, in previous studies. To correct discrepancies among the corner frequencies of events, a scaling law based on the ω À2 model was assumed and the spectral amplitude decay beyond the corner frequency was compensated. The observations were also corrected for location, focal mechanism, and time of occurrence. After estimating the empirical Green's tensor spatial derivative (EGTD) from 11 aftershock events, using 0.2-10 Hz band-pass-filtered waveforms, the strong motion records for the mainshock and aftershocks were simulated. In the simulation of each event, the EGTD elements were multiplied by the moment tensor elements followed by summation and then corrected in their spectral amplitude taking the corner frequency of each event into account. At the closest epicentral distance, and for most events, an acceptable agreement was found between calculated and observed waveforms. The results were also compared with the outcomes of simulations using the empirical Green's function method. As EGTD elements are determined by local underground structures, they could prove useful for future structural studies.