We developed a new modeling framework within an Earth system model to study mountain glacier evolution under different climate scenarios, applied here to the European Alps. Substantial Alpine glacier mass loss is projected under current climate conditions, with near-total loss under further warming. This is the first use of a 3D, higher-order ice-flow model for regional glacier simulations, enabling assessment of coupled land ice–Earth system processes.