Coupled ensemble simulations of the Northern Hemisphere ice sheets at last two glacial maxima

Coupled climate-ice sheet models can capture important interactions between the ice sheets and the climate that can help us better understand an ice sheet's response to changes in forcings. In this respect, they are a useful tool for simulating future ice sheet and sea level changes as a result of climate change. However, such models have large uncertainties related to the choice of climate and ice sheet parameters used. The same processes that operate today, also occurred in glacial times, and previous work has shown that simulating the North American ice sheet at the Last Glacial Maximum (LGM; ~21 ka BP) provides a strong benchmark for testing coupled climate-ice sheet models and recalibrating uncertain parameters that control surface mass balance and ice flow (Gandy et al., 2023). Here, we build on this work by performing the first coupled FAMOUS-BISICLES simulations of the last two glacial maxima, including all Northern Hemisphere ice sheets interactively. The ice sheet component of this model is capable of efficiently simulating marine ice sheets, such as the Eurasian ice sheet, despite the high computational cost of higher order physics. We simulate and compare both the LGM and the Penultimate Glacial Maximum (PGM; ~140 ka BP), since both periods displayed major differences in the distribution of ice between Eurasia and North America. Uncertainty is explored by running ensembles of 120 simulations, randomly varying the uncertain parameters controlling ice sheet dynamics and climate through Latin Hypercube Sampling. We also work on improving the representation of ice streams in the model through performing internal ice temperature spin ups and sensitivity tests varying till water drainage properties. The ensemble members are evaluated against empirical data on ice sheet extent and ice stream location to find combinations of parameters that produce reasonable simulations of the North American and Eurasian ice sheets for both periods. We determine the impact of the uncertainty in these parameters on the result and whether both ice sheets show similar sensitivities to the model parameters. These simulations will provide a starting point for analysing some of the interactions between the climate and the ice sheets during glacial periods and how they may have led to different ice sheet evolutions.