Generally,tsunami waves become hazardous only when approaching the coast.Studying the runup and inundation of tsunami waves is important for understanding the tsunami evolution and for tsunami hazard assessment.Here,we simulated the 1993 Hokkaido-Nansei-Oki(HNO) tsunami using a finite-difference model based on nonlinear shallow-water equations.We focused on the runup and inundation of tsunami wave propagation onto coastal area of Okushiri Island near Hokkaido,Japan,and investigate the relationship of different runup heights with the morphology and bathymetry of the seashore.In the simulation,a nested 4-layer grid system and moving boundary technique are adopted to study runup and inundation.The calculated tsunami heights and inundations in the region agreed well with field measurements.The local bathymetric and topographic characteristics had a first-order effect on the runup.Numerical experiments show that the focusing of certain local bathymetric features would amplify both wave height and current velocity remarkably.The results show that computation on dense grids is necessary to reproduce the observed runup heights,and inundation velocity is an important factor preventing tsunami devastation.In addition,we discussed the potential capability of sediment transport to illustrate the impact of tsunami waves on coastal geomorphology.