Abstract. Snowpacks modulate water storage over extended land regions and at the same time play a central role in the surface albedo feedback, impacting the climate system energy balance. Despite the complexity of snow processes and their importance for both land hydrology and global climate, several state-of-the-art land surface models and Earth System Models still employ relatively simple descriptions of the snowpack dynamics. In this study we present a newly-developed snow scheme tailored to the Geophysical Fluid Dynamics Laboratory (GFDL) land model version 4.1. This new snowpack model, named GLASS (Global LAnd–Snow Scheme), includes a refined and dynamical vertical-layering snow structure that allows us to track the temporal evolution of snow grain properties in each snow layer, while at the same time limiting the model computational expense, as is necessary for a model suited to global-scale climate simulations. In GLASS, the evolution of snow grain size and shape is explicitly resolved, with implications for predicted bulk snow properties, as they directly impact snow depth, snow thermal conductivity, and optical properties. Here we describe the physical processes in GLASS and their implementation, as well as the interactions with other surface processes and the land–atmosphere coupling in the GFDL Earth System Model. The performance of GLASS is tested over 10 experimental sites, where in situ observations allow for a comprehensive model evaluation. We find that when compared to the current GFDL snow model, GLASS improves predictions of seasonal snow water equivalent, primarily as a consequence of improved snow albedo. The simulated soil temperature under the snowpack also improves by about 1.5 K on average across the sites, while a negative bias of about 1 K in snow surface temperature is observed.