Rapid neoglaciation on Ellesmere Island promoted by enhanced summer snowfall in a transient climate model simulation of the middle-late-Holocene

Arctic neoglaciation following the Holocene Thermal Maximum is an important feature of late-Holocene climate. We investigated this phenomenon using a transient 6000-year simulation with the CESM-CAM5 climate model driven by orbital forcing, greenhouse gas concentrations, and a land use reconstruction. During the first three millennia analyzed here (6–3 ka), mean Arctic snow depth increases, despite enhanced greenhouse forcing. Superimposed on this secular trend is a very abrupt increase in snow depth between 5 and 4.9 ka on Ellesmere Island and the Greenland coasts, in rough agreement with the timing of observed neoglaciation in the region. This transition is especially extreme on Ellesmere Island, where end-of-summer snow coverage jumps from nearly 0 to virtually 100% in 1 year, and snow depth increases to the model’s imposed maximum within 15 years. This climatic shift involves more than the Milankovitch-based expectation of cooler summers causing less snow melt. Coincident with the onset of the cold regime are two consecutive summers with heavy snowfall on Ellesmere Island that help to short-circuit the normal seasonal melt cycle. These heavy snow seasons are caused by synoptic-scale, cyclonic circulation anomalies over the Arctic Ocean and Canadian Archipelago, including an extremely positive phase of the Arctic Oscillation. Our study reveals that a climate model can produce sudden climatic transitions in this region prone to glacial inception and exceptional variability, due to a dynamic mechanism (more summer snowfall induced by an extreme circulation anomaly) that augments the traditional Milankovitch thermodynamic explanation of orbitally induced glacier development.