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.

A new satellite passive microwave (PMW) melt onset retrieval algorithm based on temporal variations in the differences of the brightness temperature between 19 and 37 GHz is shown to be as effective as radar (e.g., QuikScat) measurements. The PMW technique shows improved melt estimates that are more closely linked to observed snow‐off dates than previous studies. An integrated pan‐Arctic (north of 50°N) melt onset date (MOD) dataset is produced by combining estimates on land and sea ice for the entire satellite PMW record. During the 1979–2011 period, significant trends of 2~3 days (decade)⁻¹ to earlier MOD are mainly concentrated over the Eurasian land sector of the Arctic, consistent with changes in spring snow cover extent observed with visible satellite data. The variability and change in melt onset are largely driven by spring surface air temperature, with insignificant influence from low‐frequency modes of atmospheric circulation.

The Canadian Earth System Model version 5.0 (CanESM5.0), the most recent major version of the global climate model developed at the Canadian Centre for Climate Modelling and Analysis (CCCma) at Environment and Climate Change Canada (ECCC), has been used extensively in climate research and for providing future climate projections in the context of climate services. Previous studies have shown that CanESM5.0 performs well compared to other models and have revealed several model biases. To address these biases, the CCCma has recently initiated the “Analysis for Development” (A4D) activity, a coordinated analysis activity in support of CanESM development. Here we describe the goals and organization of this effort and introduce two variants (“p1” and “p2”) of a new CanESM version, CanESM5.1, which features important improvements as a result of the A4D activity. These improvements include the elimination of spurious stratospheric temperature spikes and an improved simulation of tropospheric dust. Other climate aspects of the p1 variant of CanESM5.1 are similar to those of CanESM5.0, while the p2 variant of CanESM5.1 features reduced equilibrium climate sensitivity and improved El Niño–Southern Oscillation (ENSO) variability as a result of intentional tuning of the atmospheric component. The A4D activity has also led to the improved understanding of other notable CanESM5.0 and CanESM5.1 biases, including the overestimation of North Atlantic sea ice, a cold bias over sea ice, biases in the stratospheric circulation and a cold bias over the Himalayas. It provides a potential framework for the broader climate community to contribute to CanESM development, which will facilitate further model improvements and ultimately lead to improved climate change information.

The Canadian Land Surface Scheme Including Biogeochemical Cycles (CLASSIC) is an open-source community model designed to address research questions that explore the role of the land surface in the global climate system. Here, we evaluate how well CLASSIC reproduces the energy, water, and carbon cycle when forced with quasi-observed meteorological data. Model skill scores summarize how well model output agrees with observation-based reference data across multiple statistical metrics. A lack of agreement may be due to deficiencies in the model, its forcing data, and/or reference data. To address uncertainties in the forcing, we evaluate an ensemble of CLASSIC runs that is based on three meteorological data sets. To account for observational uncertainty, we compute benchmark skill scores that quantify the level of agreement among independent reference data sets. The benchmark scores demonstrate what score values a model may realistically achieve given the uncertainties in the observations. Our results show that uncertainties associated with the forcing and observations are considerably large. For instance, for 10 out of 19 variables assessed in this study, the sign of the bias changes depending on what forcing and reference data are used. Benchmark scores are much lower than expected, implying large observational uncertainties. Model and benchmark score values are mostly similar, indicating that CLASSIC performs well when considering observational uncertainty. Future model development should address (i) a positive albedo bias and resulting shortwave radiation bias in parts of the Northern Hemisphere (NH) extratropics and Tibetan Plateau, (ii) an out-of-phase seasonal gross primary productivity cycle in the humid tropics of South America and Africa, (iii) a lacking spatial correlation of annual mean net ecosystem exchange with site-level measurements, (iv) an underestimation of fractional area burned and corresponding emissions in the boreal forests, (v) a negative soil organic carbon bias in high latitudes, and (vi) a time lag in seasonal leaf area index maxima in parts of the NH extratropics. Our results will serve as a baseline for guiding and monitoring future CLASSIC development.

This study investigates the role of leaf/plant area index (LAI/PAI) specification on the large spread of winter albedo simulated by climate models. To examine the sensitivity of winter albedo to LAI, we perform a sensitivity analysis using two methods commonly used to compute albedo in snow‐covered forests as well as diagnostic calculations within version 4.2 of the Canadian Atmospheric Model for which PAI is systematically varied. The results show that the simulated albedo is very sensitive to negative PAI biases, especially for smaller PAI values. The LAI and surface albedo of boreal forests in the presence of snow simulated by the Coupled Model Intercomparison Project Phase 5 (CMIP5) models are evaluated using satellite observations. The evaluation of CMIP5 models suggest that inaccurate tree cover fraction due to improper plant functional type specification or erroneous LAI parameterization in some models explains, in part, an observed positive bias in winter albedo over boreal forest regions of the Northern Hemisphere. This contributes to a large intermodel spread in simulated surface albedo in the presence of snow over these regions and is largely responsible for uncertainties in simulated snow‐albedo feedback strength. Errors are largest (+20–40%) in models with large underestimation of LAI but are typically within ±15% when simulated LAI is within the observed range. This study underscores the importance of accurate representation of vegetation distribution and parameters in realistic simulation of surface albedo.