Global-scale modelling of mountain glacier evolution since the last interglacial

REFERENCE: Barndon, S., Lima, A.C., Chandler, D.M., Wiersma, A.T., Leger, T., Prats, R.P., Rentier, E.S. & Flantua, S.G.A. Global-scale modelling of mountain glacier evolution since the last interglacial. DESCRIPTION: This study is the result of the Instructed Glacier Model (IGM) transient paleoglacier simulations in 8 mountain ranges through the Late Pleistocene ( ∼130 ka BP) up to the present-day (2000 CE) at 500 m resolution. Mountain regions were subdivided into 14 simulation domains for modelling purposes. The simulation domains consists of four simulation domains covering the Northern Andes (Sierra Nevada de Santa Marta, Cordillera Oriental and Cordillera de Merida, Cordillera Central and Cordillera Occidental (Colombia), and Cordillera Real and Cordillera Occidental (Ecuador)); three simulation domains in the Ethiopian Highlands (Semien Mountains, Bale Mountains, and Arsi Mountains); two simulation domains in the Eastern Rift Mountains (Mount Elgon, and Mount Kenya and Mount Kilimanjaro). The Rwenzori Mountains, Sierra Nevada (US), Greater Caucasus, Pyrenees, and European Alps were kept as individual simulation domains. In this study, we assess and discuss IGM capabilities and performance in paleoglacier distribution, area and volume in different climatic-topographic conditions using both local processing and super-computer interface. Additionally, we propose a novel method for model evaluation that considers equifinality in model parameters-outputs, facilitating the visualization of robust agreement of glaciated areas independent of parameter calibration from glaciated areas sensitive and highly sensitive to parameter variations. Modelled outputs of paleoglacier area, ice thickness and surface velocity are going to be released upon the final publication of this manuscript, as well as the model evaluation outputs. DOCUMENTS: manuscript.pdf contains the preprint manuscript. Supplementary Table 1: S1 - Regional scaling factors.xlsx contains the list of regional proxies and modelled reconstructions for temperature used to calculate the regional scaling factors.