An experimental quantification of the strength and volume of real, heterogeneous, fractured rock masses is crucial when assessing rock slope stability. In order to quantitatively characterize the internal structure of fractured rock slopes, we present three‐dimensional numerical simulations of seismic wave propagation and compare with observations. We introduce a simple, effective model for fractured rock mass, which can easily be applied to simulate weak‐motion seismic wave propagation. The macroscopic compliant fractures cutting the rock mass are modeled as finite‐width zones of reduced elastic parameters characterized by shear and normal stiffness. The widths of such zones are not fixed and can be adjusted to fit the grid step in the numerical method. The proposed rock mass model is applied and tested for the Walkerschmatt site in southwest Switzerland. Synthetic ambient vibrations are generated using a finite‐difference method for the fractured rock mass, shaped by the real terrain geometry, and compared with the measurements. The observed seismic response is satisfactorily reproduced in a broad frequency range (0.5–10 Hz). The synthetized response is primarily controlled by the stiffness, depth, number of fractures, and inertial mass of the fractured rock. The simulated amplification and ground‐motion directionality correspond with the observed levels, unless (1) the simplified cracks reach depths of 200–300 m; and (2) the fracture network is larger with respect to the mapped network. This illustrates the potential of ambient vibration methods in combination with numerical simulations to infer depth, volume, and mechanical characteristics of slope instability.