Global- and continental-scale snow simulation has been used to produce snow estimates and attribute the change of snow into hydrological variables. However, uncertainties in global- and continental-scale snow simulation due to model structure, model parameters, and meteorological forcing has not been well documented. In the present study, we used MATSIRO land surface model and conducted ensemble simulations for snow, i.e., perturbed parameter ensemble simulation and multi precipitation ensemble simulation. The results show that uncertainties in parameters are larger in the melting season than accumulation season and uncertainties in precipitation are large in accumulation seasons. The uncertainties of parameters in melting season are comparable to that of precipitation.