In this article, stochastic finite-fault simulation combined with site response analysis is used to understand the spatial distribution of ground motion in Guwahati city due to three damaging earthquakes. The rock level ground motion for the scenario earthquakes is generated based on the stochastic finite-fault methodology. These simulated motions are further amplified up to the surface by equivalent linear site response analyses using the available borelog data at 100 different locations in Guwahati city. A set of twenty simulated rock level time histories for each event, are used to compute the surface level ground motion. Response spectra are computed and the results are presented in the form of contour maps, at selected natural periods. The mean amplification due to local soil deposit is as high as 2.2 at most of the sites in Guwahati city. Based on these simulated motions, an average site correction factor is obtained for soil sites in Guwahati city. The standard error in the simulated response spectra is also reported. The contour maps obtained will be useful in identifying vulnerable places in Guwahati city.