This paper describes ESM-SnowMIP, an international coordinated modelling effort to evaluate current snow schemes, including snow schemes that are included in Earth system models, in a wide variety of settings against local and global observations. The project aims to identify crucial processes and characteristics that need to be improved in snow models in the context of local- and global-scale modelling. A further objective of ESM-SnowMIP is to better quantify snow-related feedbacks in the Earth system. Although it is not part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6), ESM-SnowMIP is tightly linked to the CMIP6-endorsed Land Surface, Snow and Soil Moisture Model Intercomparison (LS3MIP).

Twenty-seven models participated in the Earth System Model–Snow Model Intercomparison Project (ESM-SnowMIP), the most data-rich MIP dedicated to snow modeling. Our findings do not support the hypothesis advanced by previous snow MIPs: evaluating models against more variables and providing evaluation datasets extended temporally and spatially does not facilitate identification of key new processes requiring improvement to model snow mass and energy budgets, even at point scales. In fact, the same modeling issues identified by previous snow MIPs arose: albedo is a major source of uncertainty, surface exchange parameterizations are problematic, and individual model performance is inconsistent. This lack of progress is attributed partly to the large number of human errors that led to anomalous model behavior and to numerous resubmissions. It is unclear how widespread such errors are in our field and others; dedicated time and resources will be needed to tackle this issue to prevent highly sophisticated models and their research outputs from being vulnerable because of avoidable human mistakes. The design of and the data available to successive snow MIPs were also questioned. Evaluation of models against bulk snow properties was found to be sufficient for some but inappropriate for more complex snow models whose skills at simulating internal snow properties remained untested. Discussions between the authors of this paper on the purpose of MIPs revealed varied, and sometimes contradictory, motivations behind their participation. These findings started a collaborative effort to adapt future snow MIPs to respond to the diverse needs of the community.

This paper describes in situ meteorological forcing and evaluation data, and bias-corrected reanalysis forcing data, for cold regions' modelling at 10 sites. The long-term datasets (one maritime, one arctic, three boreal, and five mid-latitude alpine) are the reference sites chosen for evaluating models participating in the Earth System Model-Snow Model Intercomparison Project. Periods covered by the in situ data vary between 7 and 20 years of hourly meteorological data, with evaluation data (snow depth, snow water equivalent, albedo, soil temperature, and surface temperature) available at varying temporal intervals. Thirty-year (1980–2010) time series have been extracted from a global gridded surface meteorology dataset (Global Soil Wetness Project Phase 3) for the grid cells containing the reference sites, interpolated to 1 h time steps and bias-corrected. Although the correction was applied to all sites, it was most important for mountain sites hundreds of metres higher than the grid elevations and for which uncorrected air temperatures were too high and snowfall amounts too low. The discussion considers the importance of data sharing to the identification of errors and how the publication of these datasets contributes to good practice, consistency, and reproducibility in geosciences. The Supplement provides information on instrumentation, an estimate of the percentages of missing values, and gap-filling methods at each site. It is hoped that these datasets will be used as benchmarks for future model development and that their ease of use and availability will help model developers quantify model uncertainties and reduce model errors. The data are published in the repository PANGAEA and are available at https://doi.pangaea.de/10.1594/PANGAEA.897575.

Datasets derived from measurements at Sodankylä, Finland, for driving and evaluating snow models are presented. This is the first time that such complete datasets have been made available for a site in the Arctic. The continuous October 2007–September 2014 driving data comprise all of the meteorological variables required as inputs for physically based snow models at hourly intervals: incoming solar and longwave radiation, snowfall and rainfall rates, air temperature, humidity, wind speed and atmospheric pressure. Two versions of the driving data are provided: one using radiation and wind speed measurements made above the height of the trees around the clearing where the evaluation data were measured and one with adjustments for the influence of the trees on conditions close to the ground. The available evaluation data include automatic and manual measurements of bulk snow depth and snow water equivalent, and profiles of snow temperature, snow density and soil temperature. A physically based snow model is driven and evaluated with the datasets to illustrate their utility. Shading by trees is found to extend the duration of both modelled and observed snow cover on the ground by several days a year.

Recent studies show that shrubs are colonizing higher latitudes and altitudes in the Arctic. Shrubs affect the wind transport, accumulation and melt of snow, but there have been few sensitivity studies of how shrub expansion might affect snowmelt rates and timing. Here, a three-source energy balance model (3SOM), which calculates vertical and horizontal energy fluxes – thus allowing within-cell advection – between the atmosphere, snow, snow-free ground and vegetation, is introduced. The three-source structure was specifically adopted to investigate shrub–tundra processes associated with patchy snow cover that single- or two-source models fail to address. The ability of the model to simulate the snow regime of an upland tundra valley is evaluated; a blowing snow transport and sublimation model is used to simulate premelt snow distributions and 3SOM is used to simulate melt. Some success at simulating turbulent fluxes in point simulations and broad spatial pattern in distributed runs is shown even if the lack of advection between cells causes melt rates to be underestimated. The models are then used to investigate the sensitivity of the snow regime in the valley to varying shrub cover and topography. Results show that, for domain average shrub fractional cover ≤0.4, topography dominates the pre- and early melt energy budget but has little influence for higher shrub cover. The increase in domain average sensible heat fluxes and net radiation with increasing shrub cover is more marked without topography where shrubs introduce wind-induced spatial variability of snow and snow-free patches. As snowmelt evolves, differences in the energy budget between simulations with and without topography remain relatively constant and are independent of shrub cover. These results suggest that, to avoid overestimating the effect of shrub expansion on the energy budget of the Arctic, future large-scale investigations should consider wind redistribution of snow, shrub bending and emergence, and sub-grid topography as they affect the variability of snow cover.

The 30-year simulations of seasonal snow cover in 22 physically based models driven with bias-corrected meteorological reanalyses are examined at four sites with long records of snow observations. Annual snow cover durations differ widely between models, but interannual variations are strongly correlated because of the common driving data. No significant trends are observed in starting dates for seasonal snow cover, but there are significant trends towards snow cover ending earlier at two of the sites in observations and most of the models. A simplified model with just two parameters controlling solar radiation and sensible heat contributions to snowmelt spans the ranges of snow cover durations and trends. This model predicts that sites where snow persists beyond annual peaks in solar radiation and air temperature will experience rapid decreases in snow cover duration with warming as snow begins to melt earlier and at times of year with more energy available for melting.