Subgrid snow cover is one of the key parameters in global land models since snow cover has large impacts on the surface energy and moisture budgets, and hence the surface temperature. In this study, the Subgrid Snow Distribution (SSNOWD) snow cover parameterization was incorporated into the Minimal Advanced Treatments of Surface Interaction and Runoff (MATSIRO) land surface model. SSNOWD assumes that the subgrid snow water equivalent (SWE) distribution follows a lognormal distribution function, and its parameters are physically derived from geoclimatic information. Two 29-yr global offline simulations, with and without SSNOWD, were performed while forced with the Japanese 25-yr Reanalysis (JRA-25) dataset combined with an observed precipitation dataset. The simulated spatial patterns of mean monthly snow cover fraction were compared with satellite-based Moderate Resolution Imaging Spectroradiometer (MODIS) observations. The snow cover fraction was improved by the inclusion of SSNOWD, particularly for the accumulation season and/or regions with relatively small amounts of snowfall; snow cover fraction was typically underestimated in the simulation without SSNOWD. In the Northern Hemisphere, the daily snow-covered area was validated using Interactive Multisensor Snow and Ice Mapping System (IMS) snow analysis datasets. In the simulation with SSNOWD, snow-covered area largely agreed with the IMS snow analysis and the seasonal cycle in the Northern Hemisphere was improved. This was because SSNOWD formulates the snow cover fraction differently for the accumulation season and ablation season, and represents the hysteresis of the snow cover fraction between different seasons. The effects of including SSNOWD on hydrological properties and snow mass were also examined.