An M_S6.2 earthquake struck Jishishan County, Gansu Province, China, on December 18, 2023, leading to a tremendous sand boil and causing more than 150 deaths. As the trigger of the sand boils, the ground shaking during the earthquake attracts extensive concerns. In this study, ground motions of the Jishishan earthquake were reproduced by using a stochastic finite-fault modeling approach and compared with observations. To maintain the far-field received energy independent of subfault size, two improvements were made to the source spectral model. The comparison shows that: Simulated ground motions agree well with observations in terms of the waveform, peak, and duration of the acceleration time histories at reported stations. The response spectra of the observed time series are in general well reproduced by the simulation except for the remarkably large amplitude at certain periods that may result from a site response. Both the simulated and observed PGA are higher than the predictions of the empirical model but remain close to +1 standard deviation. The contour map of the simulated PGA and PGV shows a similar pattern to the observation except for the slight underestimation near the northeast and northwest corners of the study area. The maximum intensity derived from the simulated PGA is degree Ⅺ, which is also consistent with the reported shaking intensity.