In this study, stochastic earthquake catalog of the Erzincan region in Turkey is generated based on synthetic ground motions. Monte Carlo simulation method is used to identify the spatial and temporal distribution of events. Ground motion time histories are generated using stochastic simulation methodology. Annual exceedance rate of each ground motion amplitude is calculated through statistical distribution of the complete set of ground motions. The results are compared with classical probabilistic seismic hazard analysis (PSHA). Classical PSHA generally produces larger spectral amplitudes than the proposed study due to wide range of aleatory variability. The effects of near-field forward directivity and detailed site response are also investigated on the results.