The southernmost portion of the Ryukyu Trench near the island of Taiwan potentially generates tsunamigenic earthquakes with magnitudes from 7.5 to 8.7 through shallow rupture. The fault model for this potential region dips 10° northward with a rupture length of 120 km and a width of 70 km. An earthquake magnitude of Mw 8.15 is estimated by the fault geometry with an average slip of 8.25 m as a constraint on the earthquake scenario. Heterogeneous slip distributions over the rupture surface are generated by a stochastic slip model, which represents the decaying slip spectrum according to k -2 in the wave number domain. These synthetic slip distributions are consistent with the abovementioned identical seismic conditions. The results from tsunami simulations illustrate that the propagation of tsunami waves and the peak wave heights largely vary in response to the slip distribution. Changes in the wave phase are possible as the waves propagate, even under the same seismic conditions. The tsunami energy path not only follows the bathymetry but also depends on the slip distribution. The probabilistic distributions of the peak tsunami amplitude calculated by 100 different slip patterns from 30 recording stations reveal that the uncertainty decreases with increasing distance from the tsunami source. The highest wave amplitude for 30 recording points is 7.32 m at Hualien for 100 different slips. Compared with the stochastic-slip distributions, the uniform slip distribution will be highly underestimated, especially in the near field. In general, the uniform slip assumption only represents the average phenomenon and will consequently ignore the possibility of tsunami waves. These results indicate that considering the effects of heterogeneous slip distributions is necessary for assessing tsunami hazards to provide additional information about tsunami uncertainties and facilitate a more comprehensive estimation.