Effect of orographic gravity wave drag on Northern Hemisphere climate in transient simulations of the last deglaciation

Long transient simulations of the last deglaciation are increasingly being performed to identify the drivers of multiple rapid Earth system changes that occurred in this period. Such simulations frequently prescribe temporal variations in ice sheet properties, which can play an important role in controlling atmospheric and surface climate. To preserve a model’s standard performance in simulating climate, it is common to apply time dependent orographic variations, including parameterised sub-grid scale orographic variances, as anomalies from the pre-industrial state. This study investigates the causes of two abrupt climate change events in the Northern Hemisphere extratropics occurring between 16 and 14 thousand years ago in transient simulations of the last deglaciation from the Hadley Centre coupled general circulation model (HadCM3). One event is characterized by regional Northern Hemisphere changes comprising a centennial scale cooling of ~ 10 °C across Fennoscandia followed by rapid warming in less than 50 years as well as synchronous shifts in the Northern Annular Mode. The second event has comparable but temporally reversed characteristics. Sensitivity experiments reveal the climate anomalies are exclusively caused by artificially large values of orographic gravity wave drag, resulting from the combined use of the orographic anomaly method along with a unique inclusion of transient orography that linearly interpolates between timesteps in the ice sheet reconstruction. Palaeoclimate modelling groups should therefore carefully check the effects of their sub-grid scale orographic terms in transient palaeoclimate simulations with prescribed topographic evolution.