This paper presents a novel Discrete Element Method (DEM) on the GPU for efficient snow simulation. To this end, our approach employs an iterative scheme on particles that easily allows the snow density to vary vastly for simulation while still maintaining a relatively large time step. We provide computationally inexpensive ways to capture cohesion and compression in the snow that enables us to generalize the behavior of various kinds of snow (like dry, wet, etc.) by varying physical parameters within the same simulator. We achieve a speed-up of nearly eight times with one million snow particles over the existing real-time method, even while dealing with scenes containing complex object boundaries. Furthermore, our simulator not only retains stability at these large time steps but also improves upon the physical behavior of the existing method. We have also conducted a user evaluation of our approach, where a majority of the participants voted in favor of its realism value for computer games.