In this paper we address one of the big challenges in snow hydrology, namely the accurate simulation of the seasonal snowpack in ungauged regions. We do this by assimilating satellite observations of snow cover into a modelling framework. Importantly (and a novelty of the paper), we include a clustering approach that permits highly efficient ensemble simulations. Efficiency gains and dependency on purely global datasets, means that this method can be applied over large areas anywhere on Earth.