Earthquakes occur because of abrupt slips on faults due to accumulated stress in the Earth’s crust. Because most of these faults and their mechanisms are not readily apparent, deterministic earthquake prediction is difficult. For effective prediction, complex conditions and uncertain elements must be considered, which necessitates stochastic prediction. In particular, a large amount of uncertainty lies in identifying whether abnormal phenomena are precursors to large earthquakes, as well as in assigning urgency to the earthquake. Any discovery of potentially useful information for earthquake prediction is incomplete unless quantitative modeling of risk is considered. Therefore, this manuscript describes the prospect of earthquake predictability research to realize practical operational forecasting in the near future.

## References

*IEEE Trans. Automat. Control.*

**AC-19**716–723. MR423716 10.1109/TAC.1974.1100705Akaike, H. (1974). A new look at the statistical model identification.

*IEEE Trans. Automat. Control.*

**AC-19**716–723. MR423716 10.1109/TAC.1974.1100705

*Bayesian Statistics*(

*Valencia*, 1979) (J. Bernard, M. De Groot, D. Lindley and A. Smith, eds.) 143–166. Univ. Press, Valencia. MR638876Akaike, H. (1980a). Likelihood and the Bayes procedure. In

*Bayesian Statistics*(

*Valencia*, 1979) (J. Bernard, M. De Groot, D. Lindley and A. Smith, eds.) 143–166. Univ. Press, Valencia. MR638876

*Ann. Inst. Statist. Math.*

**32**311–324. MR609025 10.1007/BF02480336Akaike, H. (1980b). On the use of the predictive likelihood of a Gaussian model.

*Ann. Inst. Statist. Math.*

**32**311–324. MR609025 10.1007/BF02480336

*A Celebration of Statistics*(A. Atkinson and S. Fienberg, eds.) 1–24. Springer, New York. MR816143 10.1007/978-1-4613-8560-8_1Akaike, H. (1985). Prediction and entropy. In

*A Celebration of Statistics*(A. Atkinson and S. Fienberg, eds.) 1–24. Springer, New York. MR816143 10.1007/978-1-4613-8560-8_1

*Bull. Earthq. Res. Inst.*

**43**237–238.Aki, K. (1965). Maximum likelihood estimate of $b$ in the formula $\log\mathrm{N}=a-bM$ and its confidence limits.

*Bull. Earthq. Res. Inst.*

**43**237–238.

*Earthquake Prediction*(

*Maurice Ewing Series*, 4) (D. Simpson and P. Richards, eds.) 566–574. American Geophysical Union, Washington, DC.Aki, K. (1981). A probabilistic synthesis of precursory phenomena. In

*Earthquake Prediction*(

*Maurice Ewing Series*, 4) (D. Simpson and P. Richards, eds.) 566–574. American Geophysical Union, Washington, DC.

*J. Geophys. Res.*

**118**1–14.Bansal, A. and Ogata, Y. (2013). Non-stationary epidemic type aftershock sequence model for seismicity prior to the 26 December 2004, M9.1 Sumatra-Andaman Islands mega-earthquake.

*J. Geophys. Res.*

**118**1–14.

*Biometrika*

**62**269–276. MR400509 10.1093/biomet/62.2.269Cox, D. R. (1975). Partial likelihood.

*Biometrika*

**62**269–276. MR400509 10.1093/biomet/62.2.269

*An Introduction to the Theory of Point Processes. Vol. I*:

*Elementary Theory and Methods*, 2nd ed. Springer, New York. MR1950431Daley, D. J. and Vere-Jones, D. (2003).

*An Introduction to the Theory of Point Processes. Vol. I*:

*Elementary Theory and Methods*, 2nd ed. Springer, New York. MR1950431

*J. Geophys. Res.*

**99**2601–2618.Dieterich, J. (1994). A constitutive law for rate of earthquake production and its application to earthquake clustering.

*J. Geophys. Res.*

**99**2601–2618.

*Geophys. J. Int.*

**100**19–32.Frohlich, C. and Davis, S. D. (1990). Single-link cluster analysis as a method to evaluate spatial and temporal properties of earthquake catalogues.

*Geophys. J. Int.*

**100**19–32.

*Biometrika*

**58**255–277. MR319314 10.2307/2334515Good, I. J. and Gaskins, R. A. (1971). Nonparametric roughness penalties for probability densities.

*Biometrika*

**58**255–277. MR319314 10.2307/2334515

*J. Geophys. Res.*

**103**24347–24358.Harris, R. A. (1998). Introduction to special section: Stress triggers, stress shadows, and implications for seismic hazard.

*J. Geophys. Res.*

**103**24347–24358.

*Nature Geoscience*

**2**14–144.Hashimoto, C., Noda, A., Sagiya, T. and Matsu’ura, M. (2009). Interplate seismogenic zones along the Kuril–Japan trench inferred from GPS data inversion.

*Nature Geoscience*

**2**14–144.

*Biometrika*

**58**83–90. MR278410 10.1093/biomet/58.1.83Hawkes, A. G. (1971). Spectra of some self-exciting and mutually exciting point processes.

*Biometrika*

**58**83–90. MR278410 10.1093/biomet/58.1.83

*Geophysical Research Letters*

**33**L07302.Iwata, T. and Katao, H. (2006). Correlation between the phase of the moon and the occurrences of microearthquakes in the Tamba region through point-process modeling.

*Geophysical Research Letters*

**33**L07302.

*Rep. Coord. Comm. Earthq. Predict.*

**80**80–99.JMA (2008). Japan Meteorological Agency, Seismic activity in and around Kanto and Chubu Districts (May 2007–April 2008).

*Rep. Coord. Comm. Earthq. Predict.*

**80**80–99.

*Ann. Geophys.*DOI:10.4401/ag-5350. Available at http://www.annalsofgeophysics.eu/index.php/annals/article/view/5350/5371.Jordan, T. H., Chen, Y. T., Gasparini, P., Madariaga, R., Main, I., Marzocchi, W., Papadopoulos, G., Sobolev, G., Yamaoka, K. and Zschau, J. (2011). Operational earthquake forecasting. State of knowledge and guidelines for utilization.

*Ann. Geophys.*DOI:10.4401/ag-5350. Available at http://www.annalsofgeophysics.eu/index.php/annals/article/view/5350/5371.

*Nature*

**335**690–694.Keilis-Borok, V., Knopoff, L., Rotwain, I. and Allen, C. (1988). Intermediate-term prediction of occurrence times of strong earthquakes.

*Nature*

**335**690–694.

*J. Geophys. Res.*

**115**B10312. MR2683216Kumazawa, T., Ogata, Y. and Toda, S. (2010). Precursory seismic anomalies and transient crustal deformation prior to the 2008 Mw $=6.9$ Iwate-Miyagi Nairiku, Japan, earthquake.

*J. Geophys. Res.*

**115**B10312. MR2683216

*J. Geophys. Res.*

**116**B08303.Lei, X., Xie, C. and Fu, B. (2011). Remotely triggered seismicity in Yunnan, southwestern China, following the 2004 Mw9.3 Sumatra earthquake.

*J. Geophys. Res.*

**116**B08303.

*J. Geophys. Res.*

**113**B10310.Lei, X., Yu, G., Ma, S., Wen, X. and Wang, Q. (2008). Earthquakes induced by water injection at $\sim 3$ km depth within the Rongchang gas field, Chongqing, China.

*J. Geophys. Res.*

**113**B10310.

*New Zealand Journal of Geology and Geophysics*

**40**77–89.Ma, L. and Vere-Jones, D. (1997). Application of M8 and Lin–Lin algorithms to New Zealand earthquake data.

*New Zealand Journal of Geology and Geophysics*

**40**77–89.

*Bull. Disas. Prey. Res. Inst.*,

*Kyoto Univ.*

**36**43–98.Matsumura, K. (1986). On regional characteristics of seasonal variation of shallow earthquake activities in the World.

*Bull. Disas. Prey. Res. Inst.*,

*Kyoto Univ.*

**36**43–98.

*Earth Planets Space*

**63**159–169.Nanjo, K., Tsuruoka, H., Hirata, N. and Jordan, T. (2011). Overview of the first earthquake forecast testing experiment in Japan.

*Earth Planets Space*

**63**159–169.

*Geophysical Journal International*

**191**653–658.Nanjo, K., Tsuruoka, H., Yokoi, S., Ogata, Y., Falcone, G., Hirata, N., Ishigaki, Y., Jordan, T., Kasahara, K., Obara, K., Schorlemmer, D., Shiomi, K. and Zhuang, J. (2012). Predictability study on the aftershock sequence following the 2011 Tohoku-Oki, Japan, Earthquake: First results.

*Geophysical Journal International*

**191**653–658.

*J. Geophys. Res.*

**116**B04315.Nomura, S., Ogata, Y., Komaki, F. and Toda, S. (2011). Bayesian forecasting of the recurrent earthquakes and its predictive performance for a small sample size.

*J. Geophys. Res.*

**116**B04315.

*J. Phys. Earth.*

**31**115–124.Ogata, Y. (1983). Estimation of parameters in the modified Omori formula for aftershock frequencies by the maximum likelihood procedure.

*J. Phys. Earth.*

**31**115–124.

*Int. J. Forecasting*

**11**539–541.Ogata, Y. (1995). Evaluation of probability forecasts of events; invited discussion as a commentary on “Forecasting Earthquakes and Earthquake Risk” by Prof. D. Vere-Jones.

*Int. J. Forecasting*

**11**539–541.

*J. Geophys. Res*

**109**B04305.Ogata, Y. (2004b). Seismicity quiescence and activation in western Japan associated with the 1944 and 1946 great earthquakes near the Nankai trough.

*J. Geophys. Res*

**109**B04305.

*J. Geophys. Res.*

**110**B08305.Ogata, Y. (2005b). Synchronous seismicity changes in and around the northern Japan preceding the 2003 Tokachi-oki earthquake of M8.0.

*J. Geophys. Res.*

**110**B08305.

*Rep. Coord. Comm. Earthq. Predict.*

**73**666–669.Ogata, Y. (2005c). Simultaneous estimation of b-values and detection rates of earthquakes for the application to aftershock probability forecasting (in Japanese).

*Rep. Coord. Comm. Earthq. Predict.*

**73**666–669.

*Rep. Coord. Comm. Earthq. Predict.*

**74**523–528.Ogata, Y. (2005d). Seismicity changes in and around Kyushu District before the 2005 earthquake of M7.0 in the western offshore of Fukuoka Prefecture (in Japanese).

*Rep. Coord. Comm. Earthq. Predict.*

**74**523–528.

*Geophys. Res. Let.*

**33**L01303.Ogata, Y. (2006a). Monitoring of anomaly in the aftershock sequence of the 2005 earthquake of M7.0 off coast of the western Fukuoka, Japan, by the ETAS model.

*Geophys. Res. Let.*

**33**L01303.

*Tectonophysics*

**424**291–306.Ogata, Y. (2006b). Seismicity anomaly scenario prior to the major recurrent earthquakes off the east coast of Miyagi Prefecture, northern Japan.

*Tectonophysics*

**424**291–306.

*J. Geophys. Res.*

**112**B10301.Ogata, Y. (2007). Seismicity and geodetic anomalies in a wide preceding the Niigata-Ken-Chuetsu earthquake of 23 October 2004, central Japan.

*J. Geophys. Res.*

**112**B10301.

*Pure and Applied Geophysics*

**167**1115–1127.Ogata, Y. (2010a). Anomalies of seismic activity and transient crustal deformations preceding the 2005 M7.0 earthquake west of Fukuoka.

*Pure and Applied Geophysics*

**167**1115–1127.

*Rep. Coord. Comm. Earthq. Predict.*

**86**123–125.Ogata, Y. (2011a). Operational probability foreshock forecasts up until Tohoku-Oki earthquake (in Japanese).

*Rep. Coord. Comm. Earthq. Predict.*

**86**123–125.

*Earth*,

*Planets and Space*

**63**217–229.Ogata, Y. (2011b). Significant improvements of the space–time ETAS model for forecasting of accurate baseline seismicity.

*Earth*,

*Planets and Space*

**63**217–229.

*Geophys. J. Int.*

**186**331–348.Ogata, Y. (2011c). Pre-seismic anomalies in seismicity and crustal deformation: Case studies of the 2007 Noto Hanto earthquake of M6.9 and the 2007 Chuetsu-oki earthquake of M6.8 after the 2004 Chuetsu earthquake of M6.8.

*Geophys. J. Int.*

**186**331–348.

*J. R. Stat. Soc. Ser. B Stat. Methodol.*

**44**102–107. MR655379Ogata, Y. and Akaike, H. (1982). On linear intensity models for mixed doubly stochastic Poisson and self-exciting point processes.

*J. R. Stat. Soc. Ser. B Stat. Methodol.*

**44**102–107. MR655379

*Ann. Inst. Statist. Math.*

**34**373–387.Ogata, Y., Akaike, H. and Katsura, K. (1982). The application of linear intensity models to the investigation of causal relations between a point process and another stochastic process.

*Ann. Inst. Statist. Math.*

**34**373–387.

*Geophys. J. Int.*

**104**135–146.Ogata, Y., Imoto, M. and Katsura, K. (1991). 3-D spatial variation of b-values of magnitude-frequency distribution beneath the Kanto District, Japan.

*Geophys. J. Int.*

**104**135–146.

*J. Geophys. Res.*

**108**2318.Ogata, Y., Jones, L. and Toda, S. (2003). When and where the aftershock activity was depressed: Contrasting decay patterns of the proximate large earthquakes in southern California.

*J. Geophys. Res.*

**108**2318.

*J. Appl. Probab.*

**23A**291–310. MR803179 10.2307/3214359Ogata, Y. and Katsura, K. (1986). Point-process models with linearly parametrized intensity for application to earthquake data.

*J. Appl. Probab.*

**23A**291–310. MR803179 10.2307/3214359

*Geophys. J. Int.*

**113**727–738.Ogata, Y. and Katsura, K. (1993). Analysis of temporal and spatial heterogeneity of magnitude frequency distribution inferred from earthquake catalogues.

*Geophys. J. Int.*

**113**727–738.

*J. Roy. Statist. Soc. Ser. C*

**52**499–509. MR2012973 10.1111/1467-9876.00420Ogata, Y., Katsura, K. and Tanemura, M. (2003). Modelling heterogeneous space–time occurrences of earthquakes and its residual analysis.

*J. Roy. Statist. Soc. Ser. C*

**52**499–509. MR2012973 10.1111/1467-9876.00420

*J. Geophys. Res.*

**115**B06318.Ogata, Y. and Toda, S. (2010). Bridging great earthquake doublets through silent slip: On- and off-fault aftershocks of the 2006 Kuril Island subduction earthquake toggled by a slow slip on the outer-rise normal fault the 2007 great earthquake.

*J. Geophys. Res.*

**115**B06318.

*Geophys. J. Int.*

**121**233–254.Ogata, Y., Utsu, T. and Katsura, K. (1995). Statistical features of foreshocks in comparison with other earthquake clusters.

*Geophys. J. Int.*

**121**233–254.

*Rep. Coord. Comm. Earthq. Predict.*

**66**562–570.Ogata, Y. and Zhuang, J. (2001). Statistical examination of anomalies for the precursor to earthquakes, and the multi-element prediction formula: Hazard rate changes of strong earthquakes ($M\geqq 4$) around Beijing area based on the ultra-low frequency ground electric observation (1982–1997).

*Rep. Coord. Comm. Earthq. Predict.*

**66**562–570.

*Bull. Seismol. Soc. Amer.*

**103**3.Ogata, Y., Katsura, K., Falcone, G., Nanjo, K. and Zhuang, J. (2013). Comprehensive and topical evaluations of earthquake forecasts in terms of number, time, space, and magnitude.

*Bull. Seismol. Soc. Amer.*

**103**3.

*Scientific Reports*

**3**Article No. 2218.Omi, T., Ogata, Y., Hirata, Y. and Aihara, K. (2013). Forecasting large aftershocks within one day after the main shock.

*Scientific Reports*

**3**Article No. 2218.

*The Mechanics of the Earthquake. The California Earthquake of April*18, 1906,

*Report of the State Investigation Commission*

**2**. Carnegie Institution, Washington.Reid, H. F. (1910).

*The Mechanics of the Earthquake. The California Earthquake of April*18, 1906,

*Report of the State Investigation Commission*

**2**. Carnegie Institution, Washington.

*J. Geophys. Res.*

**99**13701–13712.Rhoades, D. A., Van Dissen, R. J. and Dowrick, D. J. (1994). On the handling of uncertainties in estimating the hazard of rupture on a fault segment.

*J. Geophys. Res.*

**99**13701–13712.

*Proc. Natl. Acad. Sci. USA*

**99**2514–2521.Rundle, J., Tiampo, K., Klein, W. and Martins, J. (2002). Self-organization in leaky threshold systems: The influence of near-mean field dynamics and its implications for earthquakes, neurobiology, and forecasting.

*Proc. Natl. Acad. Sci. USA*

**99**2514–2521.

*Seism. Res. Lett.*

**78**17–29.Schorlemmer, D., Gestenberger, M. C., Wiemer, S., Jackson, D. and Rhoades, D. (2007). Earthquake likelihood model testing.

*Seism. Res. Lett.*

**78**17–29.

*Pure Appl. Geophys.*

**167**859–876.Schorlemmer, D., Zechar, J., Werner, M., Field, E., Jackson, D., Jordan, T. and Relm Working Group (2010). First results of the regional earthquake likelihood models experiment.

*Pure Appl. Geophys.*

**167**859–876.

*Tectonophysics*

**413**63–75.Shebalin, P., Kellis-Borok, V., Gabrielov, A., Zaliapin, I. and Turcotte, D. (2006). Short-term earthquake prediction by reverse analysis of lithosphere dynamics.

*Tectonophysics*

**413**63–75.

*J. Geophys. Res.*

**110**B05S01.Steacy, S., Gomberg, J. and Cocco, M. (2005). Introduction to special section: Stress transfer, earthquake triggering, and time-dependent seismic hazard.

*J. Geophys. Res.*

**110**B05S01.

*Bull. Seismol. Soc. Amer.*

**56**185–200.Suyehiro, S. (1966). Difference between aftershocks and foreshocks in the relationship of magnitude to frequency of occurrence for the great Chilean earthquake of 1960.

*Bull. Seismol. Soc. Amer.*

**56**185–200.

*Europhys. Lett.*

**60**481–487.Tiampo, K., Rundle, J., Mcginnis, S., Gross, S. and Klein, W. (2002). Mean-field threshold systems and phase dynamics: An application to earthquake fault systems.

*Europhys. Lett.*

**60**481–487.

*Earth Planet. Sci. Lett.*

**365**17–24. MR3104696 10.1007/s10598-013-9162-5Terakawa, T., Hashimoto, C. and Matsu’ura, M (2013). Changes in seismic activity following the 2011 Tohoku-oki earthquake: Effects of pore fluid pressure.

*Earth Planet. Sci. Lett.*

**365**17–24. MR3104696 10.1007/s10598-013-9162-5

*Nature*

**419**58–61.Toda, S., Stein, R. and Sagiya, T. (2002). Evidence from the A.D. 2000 Izu Islands swarm that seismicity is governed by stressing rate.

*Nature*

**419**58–61.

*Geophys. Bull. Hokkaido Univ.*

**13**99–103.Utsu, T. (1965). A method for determining the value of $b$ in a formula $\log\mathrm{n}=a-bM$ showing the magnitude-frequency’ relation for earthquakes.

*Geophys. Bull. Hokkaido Univ.*

**13**99–103.

*J. Fac. Sci. Hokkaido Univ.*,

*Ser.*7

**3**197–266.Utsu, T. (1970). Aftershocks and earthquake statistics (II)—Further investigation of aftershocks and other earthquake sequences based on a new classification of earthquake sequences.

*J. Fac. Sci. Hokkaido Univ.*,

*Ser.*7

**3**197–266.

*Rep. Coord. Comm. Earthq. Predict.*

**21**164–166.Utsu, T. (1979). Calculation of the probability of success of an earthquake prediction (In the case of Izu-Oshima-Kinkai earthquake of 1978).

*Rep. Coord. Comm. Earthq. Predict.*

**21**164–166.

*J. Geophys. Res.*

**102**15.Wiemer, S. and Wyss, M. (1997). Mapping the frequency-magnitude distribution in asperities: An improved technique to calculate recurrence times?

*J. Geophys. Res.*

**102**15.

*Bull. Seismol. Soc. Amer.*

**100**1184–1195.Zechar, J., Gerstenberger, M. and Rhoades, D. (2010). Likelihood-based tests for evaluating space–rate-magnitude earthquake forecasts.

*Bull. Seismol. Soc. Amer.*

**100**1184–1195.

*Chinese Journal of Geophysics*(

*in Chinese with English Abstract*)

**55**1695–1709.Zhuang, J. C. and Jiang, C. S. (2012). Evaluation of the prediction performance of the Annual Cosultation Meeting on Earthquake Tendency by using the gambling score.

*Chinese Journal of Geophysics*(

*in Chinese with English Abstract*)

**55**1695–1709.

*J. Amer. Statist. Assoc.*

**97**369–380. MR1941459 10.1198/016214502760046925Zhuang, J., Ogata, Y. and Vere-Jones, D. (2002). Stochastic declustering of space–time earthquake occurrences.

*J. Amer. Statist. Assoc.*

**97**369–380. MR1941459 10.1198/016214502760046925

*Pure and Applied Geophysics*

**162**1367–1396.Zhuang, J., Vere-Jones, D., Guan, H., Ogata, Y. and Ma, L. (2005a). Preliminary analysis of observations on the ultra-low frequency electric field in a region around Beijing.

*Pure and Applied Geophysics*

**162**1367–1396.

*J. Geophys. Res.*

**110**B5.Zhuang, J., Chang, C., Ogata, Y. and Chen, Y. (2005b). A study on the background and clustering seismicity in the Taiwan region by using point process models.

*J. Geophys. Res.*

**110**B5.

*Imaging*,

*Modeling and Assimilation in Seismology*,

*Vol. II*(Y. Li, ed.) De Gruyter, Berlin.Zhuang, J., Ogata, Y., Vere-Jones, D., Ma, L. and Guan, H. (2013). Statistical modeling of earthquake occurrences based on external geophysical observations: With an illustrative application to the ultra-low frequency ground electric signals observed in the Beijing region. In

*Imaging*,

*Modeling and Assimilation in Seismology*,

*Vol. II*(Y. Li, ed.) De Gruyter, Berlin.