The goal of this study is to provide a stochastic method to investigate the effects of the randomness of soil properties due to their natural spatial variability on the response spectra spatial variation at sites with varying conditions.For this purpose,Monte Carlo Simulations are used to include the variability of both incident ground motion and soil parameters in the response spectra by mean of an appropriate coherency loss function and a site-dependent transfer function,respectively.The approach is built on the assumption of vertical propagation of SH type waves in soil strata with uncertain parameters.The response spectra are obtained by numerical integration of the goving equation of a single-degree-of-freedom (SDOF) system under non-stationary site-dependent and spatially varying ground motion accelerations simulated with non-uniform spectral densities and coherency loss functions.Numerical examples showed that randomness of soil properties significantly affects the amplitudes of the response spectra,indicating that as the heterogeneity induced by the randomness of the parameters of the medium increases,the spectral ordinates attenuate.