<p>Earthquake-induced landslides often have unique regional distribution characteristics and cause significant damage to the natural environment and severe loss of life and property. The aim of this study is to comprehensively investigate the contribution of ground motion time histories to the unique landslide distribution of the 2018 Hokkaido Eastern Iburi earthquake: although the geological and geomorphological conditions are similar on the northern and southern sides of the epicenter, the landslide density is obviously greater on the northern side. Overall, the synthetic waveforms in the near-source area are generated using a distance-weighted spatial interpolation method and a spectrum-matching technique, and dynamic finite difference simulations are then employed to further reproduce the actual landslide distribution pattern. The workflow contains two steps. First, the seismic ground motions in the near-source area are spatially evaluated by combining the ground motion model and the phase characteristics extracted from the downhole records. Then, dynamic analysis of slope stability using synthetic waveforms is conducted to form landslide susceptibility maps. As a result, the actual spatial pattern in which landslide density is greater on the northern side of the epicenter is reproduced and explained successfully. Differences in phase components are found to be responsible for the different deformation processes of slopes under seismic loading. The proposed method in this study can be applied to risk analysis of regional landslide hazards in regard to future earthquakes.</p>