High‐latitude Northern Hemisphere climate during the last glacial period was characterized by a series of abrupt climate changes, known as Dansgaard‐Oeschger (DO) events, which were recorded in Greenland ice cores as shifts in the oxygen isotopic composition of the ice. These shifts in inferred Northern Hemisphere high‐latitude temperature have been linked to changes in Atlantic meridional overturning strength. The response of ocean overturning circulation to forcing is non‐linear and a hierarchy of models have suggested that it may exist in multiple steady state configurations. Here, we use a time‐dependent coarse‐resolution isopycnal model with four density classes and two basins, linked by a Southern Ocean to explore overturning states and their stability to changes in external parameters. The model exhibits hysteresis in both the steady‐state stratification and overturning strength as a function of the magnitude of North Atlantic Deep Water (NADW) formation. Hysteresis occurs as a result of two non‐linearities in the model‐‐‐the surface buoyancy distribution in the Southern Ocean and the vertical diffusivity profile in the Atlantic and Indo‐Pacific basins. We construct a metric to assess circulation configuration in the model, motivated by observations from the Last Glacial Maximum, which show a different circulation structure from the modern. We find that circulation configuration is primarily determined by NADW density. The model results are used to suggest how ocean conditions may have influenced the pattern of DO events across the last glacial cycle
