The time evolution of drifting snow under a steady wind is estimated using a new numerical model of drifting snow. In the model, Lagrangian stochastic theory is used to incorporate the effect of turbulence on the motion of drifting-snow particles. This method enables us to discuss both the saltation and the suspension process. Aerodynamic entrainment, grain/bed collision (splash process), wind modification and particle size distribution are also taken into account. The calculations show that the time needed by the total mass flux to reach a steady state appears to be 3–5 s. Vertical profiles of horizontal mass flux near the surface show a similar tendency. In contrast, it takes >50 s for the wind speed and the whole mass-flux profile to reach a steady state. This longer time depends on the time-scale of the turbulent diffusion, which is responsible for the mass flux extending to an order of a few meters in height. Applying Taylor’s hypothesis, the estimated length scale at which drifting snow reaches equilibrium is around 400 m. This result is comparable with previously reported field observations.