A snow model forced by temperature and precipitation is used to simulate the spatial distribution of snow water equivalent (SWE) over a 600 000 km² portion of the province ofQuebec, Canada. We propose to improve model simulations by assimilating SWE data from sporadic manual snow surveys with a particle filter. A temporally and spatially correlated perturbation of the meteorological forcing is used to generate the set of particles. The magnitude of the perturbations is fixed objectively. First, the particle filter and direct insertion were both applied on 88 sites for which measured SWE consisted of more or less five values per year over a period of 17 years. The temporal correlation of perturbations enables us to improve the accuracy and the ensemble dispersion of the particle filter, while the spatial correlation leads to a spatial coherence in the particle weights. The spatial estimates of SWE obtained with the particle filter are compared with those obtained through optimal interpolation of the snow survey data, which is the current operational practice in Quebec. Cross-validation results as well as validation against an independent dataset show that the proposed particle filter enables us to improve the spatial distribution of the snow water equivalent compared with optimal interpolation.