Drifting snow storm is an important aeolian process that reshapes alpine glaciers and polar ice shelves, and it may also affect the climate system and hydrological cycle since flying snow particles exchange considerable mass and energy with air flow. Prior studies have rarely considered the full-scale drifting snow storm in the turbulent boundary layer, thus, the transportation feature of snow flow higher in the air and its contribution are largely unknown. In this study, a large eddy simulation is combined with a subgrid scale velocity model to simulate the atmospheric turbulent boundary layer, and a Lagrangian particle tracking method is adopted to track the trajectories of snow particles. A drifting snow storm that is hundreds of meters in depth and exhibits obvious spatial structures is produced. The snow transport flux profile at high altitude, previously not observed, is quite different from that near the surface, thus, the extrapolated transport flux profile may largely underestimate the total transport flux. At the same time, the development of a drifting snow storm involves three typical stages, the rapid growth, the gentle growth and the equilibrium stages, in which the large-scale updrafts and subgrid scale fluctuating velocities basically dominate the first and second stage, respectively. This research provides an effective way to get an insight into natural drifting snow storms.