Land models require evaluation in order to understand results and guide future development. Examining functional relationships between model variables can provide insight into the ability of models to capture fundamental processes and aid in minimizing uncertainties or deficiencies in model forcing. This study quantifies the proficiency of land models to appropriately transfer heat from the soil through a snowpack to the atmosphere during the cooling season (Northern Hemisphere: October-March). Using the basic physics of heat diffusion, we investigate the relationship between seasonal amplitudes of soil versus air temperatures due to insulation from seasonal snow. Observations demonstrate the anticipated exponential relationship of attenuated soil temperature amplitude with increasing snow depth and indicate that the marginal influence of snow insulation diminishes beyond an "effective snow depth" of about 50 cm. A snow and heat transfer metric (SHTM) is developed to quantify model skill compared to observations. Land models within the CMIP5 experiment vary widely in SHTM scores, and deficiencies can often be traced to model structural weaknesses. The SHTM value for individual models is stable over 150 years of climate, 1850-2005, indicating that the metric is insensitive to climate forcing and can be used to evaluate each model's representation of the insulation process.

The Project for the Intercomparison of Land-Surface Parameterization Schemes (PILPS) aims to improve understanding and modeling of land surface processes. PILPS phase 2(d) uses a set of meteorological and hydrological data spanning 18 yr (1966-83) from a grassland catchment at the Valdai water-balance research site in Russia. A suite of stand-alone simulations is performed by 21 land surface schemes (LSSs) to explore the LSSs' sensitivity to downward longwave radiative forcing, timescales of simulated hydrologic variability, and biases resulting from single-year simulations that use recursive spinup. These simulations are the first in PILPS to investigate the performance of LSSs at a site with a well-defined seasonal snow cover and frozen soil. Considerable model scatter for the control simulations exists. However, nearly all the LSS scatter in simulated root-zone soil moisture is contained within the spatial variability observed inside the catchment. In addition, all models show a considerable sensitivity to longwave forcing for the simulation of the snowpack, which during the spring melt affects runoff, meltwater infiltration, and subsequent evapotranspiration. A greater sensitivity of the ablation, compared to the accumulation, of the winter snowpaek to the choice of snow parameterization is found. Sensitivity simulations starting at prescribed conditions with no spinup demonstrate that the treatment of frozen soil (moisture) processes can affect the long-term variability of the models. The single-year recursive runs show large biases, compared to the corresponding year of the control run, that can persist through the entire year and underscore the importance of performing multiyear simulations.