An earthquake prediction device comprises a predicted value calculation unit 16 that calculates, from a prediction formula below, a predicted value (MMIap) indicating a predicted intensity of a ground motion on the Modified Mercalli Intensity scale, using a maximum acceleration value (Aumax), which is a maximum absolute value among absolute values of vertical acceleration component of the ground motion, after a sensor starts detecting the ground motion caused by an earthquake. The prediction formula is: MMIap=αa log10(Aumax)+βa.