Despite the growing spatiotemporal density of geophysical observations at subduction zones, predicting the timing and size of future earthquakes remains a challenge. Here we simulate multiple seismic cycles in a laboratory‐scale subduction zone. The model creates both partial and full margin ruptures, simulating magnitude M w 6.2–8.3 earthquakes with a coefficient of variation in recurrence intervals of 0.5, similar to real subduction zones. We show that the common procedure of estimating the next earthquake size from slip‐deficit is unreliable. On the contrary, machine learning predicts well the timing and size of laboratory earthquakes by reconstructing and properly interpreting the spatiotemporally complex loading history of the system. These results promise substantial progress in real earthquake forecasting, as they suggest that the complex motion recorded by geodesists at subduction zones might be diagnostic of earthquake imminence.