In the present study, the local tsunami amplification observed in Ryori Bay, located on the Sanriku coast of Japan, was investigated using numerical simulations. Large-scale tsunami propagation simulations and tsunami inundation simulations for the bay were systematically conducted to estimate and model the 2011, 1933, and 1896 tsunamis that occurred off the Sanriku coast and which resulted in large run-ups. The simulation results, which are moderately consistent with observations, presented larger run-up heights and inundations for the 1933 and 1896 tsunamis (which followed relatively small earthquakes) compared to those of the 2011 tsunami (which followed a larger earthquake). Furthermore, the frequency analysis indicated that the former two tsunamis comprised higher predominant components. A tsunami inundation simulation using parametrized synthetic waveforms was conducted to identify the contributing factors associated with the large amplification and run-ups. The results indicated that the predominant components are significantly amplified in the bay and the initial decrease in the water surface elevation prior to the primary waves of the two tsunamis leads to an increase in their run-up heights. Furthermore, the simulated waveforms of the tsunamis revealed that the 1933 and 1896 tsunamis had their wavefronts changed into a steep wavefront, i.e., a bore-like wave, during their wave developments in the bay, attributed to shoaling, narrowing bay width, and the nonlinear effect of the wave. These results, therefore, indicate that bores which are known to generate large run-up heights were generated in the bay during the two tsunamis.