The paper attempts to enhance the VOF-basedfinite element method for snow deformation problem by taking into account the deposition process. The variation of snow shape due to the deposition is modeled by involving the supply from snowfall in the advection equation of the VOF function. The snowfall process is described by simplified advection and diffusion models. The developed method is verified through a comparison with a semianalytical solution for a one-dimensional deposition problem. Besides, the simulation of Nakamura's experiment for snowpack with an embedded beam is achieved. Based on these analyses the significance of snow deposition process and deformation hardening in the mechanical behavior is discussed.
