A numerical model is presented, based on the local volume averaging formulation of transport phenomena in porous media, for simulating meltwater infiltration into unsaturated, frozen soil. With the defined flow and freezing boundary conditions at the snow-soil interface, using the concept of a surface local averaging volume, the time variation in profiles of temperature, liquid / ice content, infiltration / percolation rates, and rate of phase change in upper soil layers are predicted. In addition to a parametric analysis, model estimates of infiltration are compared with quantities calculated from field measurements of soil moisture changes and temperature during snow cover ablation, showing a reasonable agreement.