1. An earthquake prediction method wherein: telluric current induction field vectors and telluric currents are estimated based on vector differences between observed magnetic field vectors of magnetic fields that are observed within an observation area and a geomagnetic vector; and change over time of the telluric currents and a state of telluric currents within said observation area are compared and collated with past patterns of change over time of the telluric currents and a state of the telluric currents to estimate a seismofocal zone, time of occurrence, and seismic intensity of a seismic event.