Based on 3-D finite difference method simulations of seismic wave propagation, we examined the processes by which the complex, scattered high-frequency (f > 1 Hz) seismic wavefield during crustal earthquakes is developed due to heterogeneous structure, which includes smallscale velocity inhomogeneity in subsurface structure and irregular surface topography on the surface, and compared with obse…