The climate of the Last Glacial showed much higher variability than the present day climate. Some of this larger variability is due to the interaction of the larger ice sheets with the other components of the climate system. To gain further understanding of the glacial climate variability, we perform coupled ice sheet climate model experiments with mPISM and ECHAM5/MPIOM/LPJ under Last Glacial Maximum (LGM) boundary conditions. We focus on one of the most dominant features of glacial climate variability: the Heinrich events.