In order to extrapolate the frequency of severe tsunami run-up events, we assumed a generalized Pareto distribution and estimated its parameters from historical records on tsunamis spanning 1,000 years by using either the maximum likelihood method or the probability-weighted moments-matching method.