Comment on “Evaluation of surface albedo and snow cover in AR4 coupled models” by A. Roesch

[1] We believe that one of the main conclusions reached by Roesch [2006, paragraph 1] (hereinafter referred to as R2006) from his evaluation of the performance of the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Review (AR4) climate model snow cover simulations, i.e., that ‘‘most IPCC AR4 climate models predict excessive snow mass in spring and suffer from a delayed spring snow melt . ’’, is incorrect for North America and probably for Eurasia as well. This erroneous conclusion is a result of deficiencies in the U.S. Air Force (USAF) global snow depth climatology of Foster and Davy [1988] (hereinafter referred to as USAF) which were previously documented by Brown et al. [2003] (hereinafter referred to as B2003) and is exacerbated by the application of an inappropriate method for estimating snow-water equivalent (SWE) from snow depth. [2] R2006 (paragraph 13) stated that the USAF data set was ‘‘generally considered to be one of the most reliable and accurate snow depth climatologies available. ’’ This statement is supported neither by the conclusions of B2003 nor by the USAF data set documentation. The USAF data set is a hand-drawn analysis based on manual interpretation of a wide range of mainly hard copy sources (e.g., atlases and climate summaries) covering various time periods. Foster and Davy [1988] clearly indicate in the documentation that they had varying levels of confidence in their data sources with low confidence in most of the highlatitude data. In mountainous areas, Foster and Davy [1988] used mean monthly precipitation totals from World Meteorological Organization climatic atlases (1970 and 1981) and assumed a mean snow density of 300 kg m 3 to arrive at a snow depth. The manual analysis was digitized on a 47 km Northern Hemisphere polar stereographic (PS) grid, and it is not clear where R2006 obtained the 1 1 equal-area grid version of the data. In the comparisons included in this comment the USAF data were used in their original gridded format. [3] B2003 carried out a detailed daily snow depth analysis on a 1/3 latitude-longitude grid over North America (NA) for the second phase of the Atmospheric Model Intercomparison Project (AMIP-2) period (1979–1996) that incorporated 8000 snow depth observations per day from U.S. cooperative stations and Canadian climate stations. The snow depth analysis is based on the scheme developed by Brasnett [1999] and employed operationally at the Canadian Meteorological Centre. The first-guess field used a simple snow accumulation, aging, and melt model driven by 6-hourly values of air temperature and precipitation from the European Centre for Medium-Range Weather Forecasting Reanalysis, December 1978 to February 1994 (ERA-15) with extensions from the Tropical Ocean–Global Atmosphere operational data archive. The gridded snow depth and estimated SWE results were found to agree well with available independent in situ and satellite data over midlatitudinal regions of the continent, and the snow depth climatology exhibited several improvements over USAF, namely, a more realistic representation of snow depths over the western cordillera and improved snow cover extent in the fall and spring periods when compared to the NOAA weekly snow cover product [Robinson et al., 1993]. The problems with the USAF product during the spring melt period are highlighted in Figure 1, which demonstrates that it ablates the snowpack too rapidly in the spring, especially over the western cordillera, resulting in significant underestimations of SWE as well as snow-covered area. For the comparisons included in this paper the B2003 data set was reinterpolated to the same 47 km PS grid as the USAF data set. [4] The B2003 data set was used to evaluate snow cover simulations of the atmospheric general circulation models participating in the AMIP-2 project by Frei et al. [2005] (several of these models were included in the R2006 evaluation), and the conclusion, contrary to R2006, was that most of the models simulated the seasonal timing and relative spatial patterns of continental-scale SWE fairly well, with a tendency to overestimate the rate of ablation during spring (compare Figure 1 of R2006 to Figure 1a of Frei et al. [2005]). This tendency can also be observed in AMIP-2–simulated snow cover extent for NA [Frei et al., 2003, Figure 3a]. It will now be shown that the different conclusions are attributable to (1) underestimated snow depths over NA during the spring period in the USAF data set and (2) the SWE-density relationship used by R2006 which systematically underestimates density in the spring period. JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112, D22102, doi:10.1029/2006JD008339, 2007 Click Here for Full Article