A generalization of the stochastic summation scheme of small earthquakes following a general form of the source spectrum (ω−n) model is developed to estimate the proper probability density functions for rupture propagation over the fault and to simulate site-specific strong ground motions in two independent stages. The probability density functions for rupture propagation are first developed for a single-corner-frequency source model, and then extended to multiplicative and additive double-corner-frequency source models to partially account for the effects of finite-fault rupture such as fault geometry, rise time, and rupture time on ground motions, and to better match observed Fourier amplitude spectra (FAS). The generalized two-stage stochastic summation scheme developed in this study can be used to simulate realistic strong-ground-motion time histories from small events, represented as empirical Green’s functions (EGFs), by specifying only two source parameters of the seismic moment and stress drop for a given site. The proposed two-stage stochastic summation schemes are tested using numerical examples following the theoretical source models as well as observed ground motions of the 2011 Mw 9.0 Tohoku earthquake in Japan and the 1994 Mw 6.7 Northridge earthquake in California to verify the applicability of the approach. Comparisons between the observed and simulated time histories for these two earthquakes, their response spectra, and FAS show satisfactory performance of the proposed approach for engineering purposes. The results of these comparisons indicate that the proposed stochastic summation of small events as EGFs can adequately replicate observed ground motions and evaluate uncertainties in ground motions with the generation of many time history realizations corresponding to a multitude of possible rupture processes . The proposed approach can also be used to generate multiple realistic acceleration time histories matched on average to the target response spectra for use in earthquake performance design, and to develop ground-motion models for low-seismicity regions, where there is a lack of recording of strong ground motions.