Abstract: A data driven, near photospheric magnetohydrodynamic model predicts spikes in the horizontal current density, and associated resistive heating rate per unit volume Q. The spikes appear as increases by orders of magnitude above background values in neutral line regions (NLRs) of active regions (ARs). The largest spikes typically occur a few hours to a few days prior to M or X flares. The spikes correspond to large vertical derivatives of the horizontal magnetic field. The model takes as input …