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).

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.

Leaf area index (LAI) and its seasonal dynamics are key determinants of vegetation productivity in nature and as represented in terrestrial biosphere models seeking to understand land surface atmosphere flux dynamics and its response to climate change. Non-structural carbohydrates (NSCs) and their seasonal variability are known to play a crucial role in seasonal variation in leaf phenology and growth and functioning of plants. The carbon stored in NSC pools provides a buffer during times when supply and demand of carbon are asynchronous. An example of this role is illustrated when NSCs from previous years are used to initiate leaf onset at the arrival of favourable weather conditions. In this study, we incorporate NSC pools and associated parameterizations of new processes in the modelling framework of the Canadian Land Surface Scheme-Canadian Terrestrial Ecosystem Model (CLASS–CTEM) with an aim to improve the seasonality of simulated LAI. The performance of these new parameterizations is evaluated by comparing simulated LAI and atmosphere–land CO₂ fluxes to their observation-based estimates, at three sites characterized by broadleaf cold deciduous trees selected from the FLUXNET database. Results show an improvement in leaf onset and offset times with about 2 weeks shift towards earlier times during the year in better agreement with observations. These improvements in simulated LAI help to improve the simulated seasonal cycle of gross primary productivity (GPP) and as a result simulated net ecosystem productivity (NEP) as well.

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.

Global land cover information is required to initialize land surface and Earth system models. In recent years, new land cover (LC) datasets at finer spatial resolutions have become available while those currently implemented in most models are outdated. This study assesses the applicability of the Climate Change Initiative (CCI) LC product for use in the Canadian Land Surface Scheme (CLASS) through comparison with finer resolution datasets over Canada, assisted with reference sample data and a vegetation continuous field tree cover fraction dataset. The results show that in comparison with the finer resolution maps over Canada, the 300 m CCI product provides much improved LC distribution over that from the 1 km GLC2000 dataset currently used to provide initial surface conditions in CLASS. However, the CCI dataset appears to overestimate needleleaf forest cover especially in the taiga-tundra transition zone of northwestern Canada. This may have partly resulted from limited availability of clear sky MEdium Resolution Imaging Spectrometer (MERIS) images used to generate the CCI classification maps due to the long snow cover season in Canada. In addition, changes based on the CCI time series are not always consistent with those from the MODIS or a Landsat-based forest cover change dataset, especially prior to 2003 when only coarse spatial resolution satellite data were available for change detection in the CCI product. It will be helpful for application in global simulations to determine whether these results also apply to other regions with similar landscapes, such as Eurasia. Nevertheless, the detailed LC classes and finer spatial resolution in the CCI dataset provide an improved reference map for use in land surface models in Canada. The results also suggest that uncertainties in the current cross-walking tables are a major source of the often large differences in the plant functional types (PFT) maps, and should be an area of focus in future work.