  We modify the new method for the statistical estimation of the tail
distribution of earthquake seismic moments introduced by Pisarenko et al.
[2009] and apply it to the earthquake catalog of Japan (1923-2007). The method
is based on the two main limit theorems of the theory of extreme values and on
the derived duality between the Generalized Pareto Distribution (GPD) and
Generalized Extreme Value distribution (GEV). We obtain the distribution of
maximum earthquake magnitudes in future time intervals of arbitrary duration
tau. This distribution can be characterized by its quantile Qq(tau) at any
desirable statistical level q. The quantile Qq(tau) provides a much more stable
and robust characteristic than the traditional absolute maximum magnitude Mmax
(Mmax can be obtained as the limit of Qq(tau) as q tends to 1, and tau tends to
infinity). The best estimates of the parameters governing the distribution of
Qq(tay) for Japan (1923-2007) are the following: Form parameter for GEV =
-0.1901 +- 0.0717; position parameter GEV(tau=200)= 6.3387 +- 0.0380; spread
parameter for GEV(tau=200)= 0.5995 +- 0.0223; Q_0.90,GEV(tau=10)= 8.34 +- 0.32.
We also estimate Qq(tau) for a set of q-values and future time periods in the
range for tau between 1 and 50 years from 2007. For comparison, the absolute
maximum estimate Mmax from GEV, which is equal to 9.57 +- 0.86, has a scatter
more than twice that of the 90 percent quantile Q_{0.90,GEV}(tau=10) of the
maximum magnitude over the next 10 years counted from 2007.
