Heinrich events, related to large‐scale surges of the Laurentide ice sheet, represent one of the most dramatic types of abrupt climate change occurring during the last glacial. Here, using a coupled atmosphere‐ocean‐biosphere‐ice sheet model, we simulate quasi‐periodic large‐scale surges from the Laurentide ice sheet. The average time between simulated events is about 7,000 yrs, while the surging phase of each event lasts only several hundred years, with a total ice volume discharge corresponding to 5–10 m of sea level rise. In our model the simulated ice surges represent internal oscillations of the ice sheet. At the same time, our results suggest the possibility of a synchronization between instabilities of different ice sheets, as indicated in paleoclimate records.