To provide estimates of past glacier mass changes, an adequate initial state is required. However, information about past glacier states at regional or global scales is largely incomplete. Our study presents a new way to initialize the Open Global Glacier Model from past climate information and present-day geometries. We show that even with perfectly known but incomplete boundary conditions, the problem of model initialization leads to nonunique solutions, and we propose an ensemble approach.