This paper evaluates the simulated Arctic land snow cover duration, snow water equivalent, snow cover fraction, surface albedo, and land surface temperature in the regional climate model HIRHAM5 during 2008–2010, compared with various satellite and reanalysis data and one further regional climate model (COSMO-CLM). HIRHAM5 shows a general agreement in the spatial patterns and annual course of these variables, although distinct biases for specific regions and months are obvious. The most prominent biases occur for east Siberian deciduous forest albedo, which is overestimated in the simulation for snow covered conditions in spring. This may be caused by the simplified albedo parameterization (e.g., nonconsideration of different forest types and neglecting the effect of fallen leaves and branches on snow for deciduous tree forest). The land surface temperature biases mirror the albedo biases in their spatial and temporal structures. The snow cover fraction and albedo biases can explain the simulated land surface temperature bias of ca. −3°C over the Siberian forest area in spring.

1. Introduction

The Arctic belongs to the key regions in the global climate system, evidenced by many rapid environmental changes (e.g., [1, 2]). In this context, the land surface temperature (LST) is an important variable, because it reflects the changes in the surface energy budget, the energy exchange between land and atmosphere and feedbacks with cryospheric variables like snow cover and frozen ground.

In addition to vegetation and soil moisture, the snow cover characteristics have a significant impact on the LST (e.g., [3, 4]). The isolating snow effect leads to less cooling of the ground under snow and reduced upward longwave radiative fluxes and thus to increased cooling of snow surface. But most important is the snow-albedo feedback mechanism (e.g., [5, 6]), at least when solar irradiation is present. Overall, the occurrence of a snow pack manifests in modified sensible and latent surface heat fluxes.

The processes in spring and autumn are particularly crucial here. During the transition seasons fast changes of LST occur, whereby the timing of seasonal warming/cooling is related to the onset of snow melt/snow fall. A warming climate is discussed to move the onset of snowmelt and reduce the overall snow cover duration (SCD) [7]. The observed decrease in snow cover duration over the Arctic is largely caused by earlier onset of snow melt in spring [8], which has a direct impact on air and ground temperatures. When melting occurs, the snow cover acts as a heat sink because of the large amounts of energy required for phase change.

The snow cover characteristics, namely, snow depth (or snow water equivalent, SWE) and snow cover fraction (SCF), but also SCD and other characteristics, show distinct regional heterogeneities. Regional climate models (RCMs) are suitable tools for simulating these snow characteristics and associated land surface conditions at relatively high horizontal resolution. The studies of Shkolnik et al. [9] and Klehmet et al. [10] demonstrated a reasonable agreement of RCM simulated Siberian snow characteristics compared against different reanalysis and satellite data. Further, Klehmet et al. [10] assessed the added value of their SWE simulations compared to reanalysis data.

In the present study, we evaluate the simulated LSTs of the RCM HIRHAM5 against ERA-Interim reanalysis data with respect to the following questions: is the model able to reproduce observed LSTs? Can biases in the modeled LST be explained with biases in snow cover characteristics? Does the influence of albedo on LST play a key role? To answer these questions, we analyze the observed and modeled three snow cover characteristics (SWE, SCF, and SCD). An accurate simulation of SWE is crucial in terms of hydrological processes, while SCF and SCD are important for the surface energy budget via the albedo. Thus, the simulated surface albedo and LST are simultaneously evaluated.

Different RCMs apply land surface models of quite different complexities. To present a first glimpse about the range of snow, surface albedo and temperature simulations over Siberia, some of the HIRHAM5 results have been compared with simulations from another RCM, namely, the COSMO-CLM (CCLM), which has been applied over a Siberian subdomain recently [10].

2. Observational and Modeled Data

To evaluate monthly results from HIRHAM5 simulations, the following reference datasets have been used: reanalyses-derived snow water equivalent (SWE) from Canadian Meteorological Centre (CMC), ERA-interim and satellite-derived products from GlobSnow; satellite-derived snow cover fraction (SCF) from MODIS and GlobSnow; satellite-derived surface albedo from MODIS and reanalyses data from ERA-interim; land surface temperature (LST) from ERA-interim. A description of these datasets and of the employed model is given in the following.

2.1. MODIS Satellite-Derived Data

The Moderate Resolution Imaging Spectroradiometer (MODIS; http://modis.gsfc.nasa.gov/index.php) is a passive optical and thermal sensor system aboard the Terra and Aqua satellites. The MODIS satellites provide global daily coverage. The wider range of data products derived from MODIS observations are valuable for evaluating the spatial-temporal results from model simulations.

2.1.1. MODIS Surface Albedo

The albedo data used here is the Climate Modeling Grid (CMG) Nadir Bidirectional Reflectance Distribution Function (BRDF) Adjusted Reflectance (NBAR) Product “MODIS/Terra Nadir BRDF-Adjusted Reflectance 16-Day L3 Global 0.05Deg CMG”, which is derived from white-sky albedo only. Since the white-sky albedo represents the integration over all zenith angles, the error of the high quality MODIS operational albedos at 500 m is well less than 5% at the majority of the validation sites studied thus far, and even those albedo values with low quality flags have been found to be primarily within 10% of the field data (e.g., [1113]). Note that along coastlines, the 1 km pixels that lie over shallow water are averaged into the 0.05 degree CMG pixels, which may cause biases for these areas. The albedo data are downloaded from http://modis-atmos.gsfc.nasa.gov/ALBEDO/. The data are available for every 8 days as 16-day running mean, and we interpolated the corresponding time intervals to obtain monthly means.

2.1.2. MODIS Snow Cover Fraction (SCF)

The “MODIS/Terra Snow Cover Monthly L3 Global 0.05Deg CMG, Version 5” data used in this study was obtained from the National Snow and Ice Data Center in Boulder, Colorado, USA ([14]; http://nsidc.org/data/mod10cm.html). The reason why we selected the Terra data as a primary key for snow detection is the characteristic of snow to have high visible reflectance and low reflectance in the near infrared, MODIS band 6 (1.6 μm), which is fully functional in Terra but only 30% functional in Aqua. SCF is calculated by using the regression equation of Salomonson and Appel [15], which is based on the Normalized Difference Snow Index (NDSI) values of each observing pixel. For more detailed information we refer to the MODIS snow products user guide [16]. The monthly snow product was derived from the daily 0.05 degree resolution product. It has an overall absolute accuracy of about 93%, but with remaining problems caused by snow/cloud discrimination [17]. Arsenault et al. [18] evaluated this product with in-situ measurements collected in Colorado and Washington states and identify a total SCF bias of about 7–17%. A similar study from Simic et al. [19] over the Canadian evergreen forests concludes a total error of about 10–20%.

2.2. CMC Snow Water Equivalent (SWE)

The operational analysis of snow depth data at the Canadian Meteorological Centre (CMC) is a suitable dataset for evaluating model simulations [20]. These 24 km horizontal resolution and daily datasets have been downloaded from http://nsidc.org/data/docs/daac/nsidc0447_CMC_snow_depth/. The analysis is based on station observations and a background snow depth field, which is generated by a simple snowpack model forced with European Centre for Medium-Range Weather Forecasts (ECMWF) precipitation and 2 m air temperature [20]. Thus, the analysis makes use of precipitation forecasts and an analysis of screen level temperature to estimate snowfall and snowmelt. Wherever snow depth observations are available, these are incorporated using the method of statistical interpolation, performed every six hours on a 1/3° grid. For details we refer to Brasnett [21]. In regions with relatively good station coverage the dataset is fairly reliable, while it is considered as reasonable estimate of snow depth over data-poor Arctic regions [22]. CMC derived the monthly SWE data from the snow depth using a lookup table of snow density, which is in turn based on snow course observations [23]. However, we assumed a constant snow density of 330 kg/m3 to be consistent with the approach that is used in the HIRHAM5 model (see Section 2.5).

2.3. GlobSnow Data

GlobSnow is a satellite-derived snow products produced within the GlobSnow funded by the European Space Agency (ESA) [24]. The SCF and SWE data are based on the aggregation of optical satellites for SCF and microwave satellites for SWE passing over within a calendar month, weather station observations of snow depth, and forward simulations employing a semi-empirical snow emission model. The 25 km resolution data have been downloaded from http://www.globsnow.info/. This most recent monthly L3B-product (v2.0) has major improvements over the heterogeneous landscapes and in regions with significant lake coverage [25, 26]. Concerning SWE retrieval accuracy, improvements are reported over Tundra and Northern Boreal forests [26]. The GlobSnow dataset contains satellite-retrieved information on snow extent, SCF, and SWE. The mapping of SCF is solved from a semiempirical model based on reflectance observation [27]. A root-mean-square error (RMSE) of 15% was calculated for the area of drainage basins when validating the snow cover with snow course observations [27]. Salminen et al. [28] find an average error of about 5–7% for boreal forest. The SWE data has been derived using the fuzzy mutual information (FMI) algorithm [29] in combination with a time-series melt-detection algorithm [8] that is based on passive microwave radiometer data and synoptic station data, which leads to poorer performance in mountain areas than non-mountain areas [30] due to strong orographic complexity. Wet snow also contributes to the difficulties in the retrieval algorithm [30]. Additionally, vegetation cover affects the accuracy of the data so that dense forest will for example, result in an underestimation of SWE [31]. Takala et al. [30] compared GlobSnow SWE with station data from the former Soviet Union, Finland, and Canada. The authors found that the RMSE is less than 40 mm in cases of SWE smaller than 150 mm, while the bias increases if SWE is larger than 150 mm. Even though, GlobSnow better represents the large-scale SWE pattern in comparison with MODIS, because its SWE product combines satellite-based passive microwave measurements with ground-based observations [30].

2.4. ERA-Interim Data

ERA-interim represents the ECMWF’s new generation global reanalysis [32] that is based on the Integrated Forecasting System (IFS) model. This reanalysis incorporates many important IFS improvements such as an increased model resolution, improved model physics, the use of four-dimensional variational (4D-Var) assimilation, and various other changes in the analysis methodology [33]. ERA-interim reanalysis data have a horizontal resolution of T255 (ca. 0.75°). Snow depth and SWE are first estimated by the forecast model based on Douville et al. [34] and then updated based on a Cressman analysis of station observations of snow depth and (when available) SCF by satellite data [35]. There are two reasons why the derived snow data are limited [32]: First, the surface data assimilation system uses a relatively simple data interpolation scheme and does not generate comprehensive information about data quality control and usage. Second, the incorrect introduction of geolocation during the processing of the Interactive Multisensor Snow and Ice Mapping System (IMS) product at ECMWF introduces biases. In the ECMWF forecast model, the SCF is linearly related to the SWE in meters [36]: where  m is the threshold of snow thickness in m water equivalent, above which the grid cell is fully covered with snow. Then the snow mass budget can be written as where is snow density, is the snow fall, and and are snow sublimation and melting. Snow albedo for forest-free areas is time dependent and parameterized according to the formulation of Verseghy [37] and Douville et al. [34]. The snow albedo is separately calculated for melting conditions: and nonmelting conditions:

Here, and are constants to control albedo changes and snow melting, respectively, represents the minimum snow albedo, and  s is the length of a day in seconds. If snow fall is larger than 1 kgm−2 h−1, the snow albedo is reset to its assumed maximum value of . For boreal forest-covered grid cells, the snow albedo is restricted to a maximum value of 0.2. Finally, each grid cell albedo value is a linear relation of forest and nonforest areas. For full details we refer to Viterbo and Betts [38].

2.5. HIRHAM5 Model

HIRHAM5 is an atmospheric RCM that combines the dynamical model core of HIRLAM [39] and the physical parameterizations of ECHAM5 [40]. The previous HIRHAM4 version has been applied to a wide range of Arctic climate studies (e.g., [41, 42]). Its circum-Arctic integration domain is shown in Figure 1 and has grid points at a horizontal resolution of 0.25° (about 25 km). There are 40 vertical model levels in the atmosphere, ranging from approximately 10 m above the surface up to a height of 10 hPa. A Simulation for this domain, driven by ERA-Interim data, has been carried out for the time period of 2008–2010. LST is calculated solving the energy balance equation at the surface, thus linking atmospheric, surface and soil processes. According to Roeckner et al. [40] the energy balance equation is used in the form: where is the heat capacity of the surface layer, is the latent surface heat flux, is the sensible surface heat flux, is the ground heat flux, and is the net surface radiation.

SWE is estimated by a budget equation that considers sublimation and deposition, snow fall, and snow melting. On the basis of SWE, the model derives snow cover fraction (SCF) as a function of forest fraction, interception capacity of the canopy, and slope of the terrain given by the subgrid-scale standard deviation of orography. SCF of a grid cell is given as where is the forest fraction of the grid cell, is the SCF on the canopy, and refers to the SCF of bare ground.

The surface albedo of the land-covered part of a grid cell is calculated by modifying a background albedo field according to the calculated SCF. The background albedo is derived from satellite data allocated to a high-resolution map of major ecosystem complexes [43]. Over snow-covered land, the grid mean albedo depends on a number of parameters such as the fractional forest area, the leaf area index, the bare-soil albedo, the snow albedo, and SCF of both the ground and the canopy [44]. Similar to SCF, the surface albedo is derived as where is the albedo of forest, which depends on the leaf area index, snow on canopy and canopy albedo as well as bare soil albedo, while is the ground albedo that depends on the bare soil and snow albedo as well as SCF of the ground.

2.6. CCLM Model

COSMO-CLM (CCLM, http://www.clm-community.eu/, [45]) is a nonhydrostatic atmospheric RCM that is based on the numerical weather prediction model Consortium for Small-scale Modeling (COSMO) performed in a Climate Limited-area Modelling (CLM). Recently, it has been applied for climate studies over Siberia [10]. The Siberian domain extends from the Laptev and Kara Seas to Northern Mongolia and from the West Siberian Lowland to the border of the Sea of Okhotsk. The overlapping area with the HIRHAM5 domain is indicated by the black lines in Figure 1. A hindcast simulation has been conducted over the Siberian domain with 50 km horizontal resolution for the period 2008–2010 using ERA-interim forcing data.

For the consideration of snow, a multilayer snow model within the land-surface scheme TERRA-ML is used. The multilayer snow model is introduced in a preliminary version by the Deutscher Wetterdienst. In the performed simulation, two snow layers have been chosen as more layers did not improve the results in the current setup. Each layer is described separately by temperature, water content, and porosity, according to the snow density. Further information concerning the multilayer snow model is given in Klehmet et al. [10]. SCF is calculated in CCLM using the following formula: where  m for snow [46].

The surface albedo over land is a mean of bare soil albedo, vegetation albedo and snow albedo if snow is present. In the used CCLM version, soil albedo is a function of soil-type and soil water content in the top soil layer. TERRA-ML distinguishes between seven soil types based on the Digital Soil Map of the World (http://www.fao.org/home/en/). Vegetation is not explicitly included. Rather an external dataset of land cover is used and specific plant characteristics are derived from the dominant land cover and land use type. These characteristics are fractional area covered by plants (plant cover), leaf area index, and roughness length. Deciduous or evergreen forest is considered separately to better account for the influence on the snow albedo of forests. The snow albedo of forests is given by where is the snow albedo on bare soil, and are the forest fractions of evergreen or deciduous forest, and and the snow albedo of both forest types.

In addition, a time-dependent snow albedo is included by using an ageing condition. A detailed description of CCLM is presented by [46].

3. Analysis Methods

All the gridded observational data have been interpolated onto the HIRHAM5 grid. This allows us to compare the datasets to one another and to evaluate the HIRHAM5 simulation by reference to these various observations. Further, monthly area-averages over a Siberian domain (see Section 2.6) have been calculated for all the datasets.

According to the simulation period, the paper discusses monthly and seasonal (spring-MAM, autumn-SON) SWE, snow cover duration (SCD), SCF, albedo, and LST averages over 2008–2010. On the one hand, we present the spatial maps for spring and autumn and on the other hand, we discuss the annual cycle of the variables averaged over the Siberian domain.

As an additional measure for snow evaluation we calculated the SCD from daily CMC data and the HIRHAM5 model. Discussing SCD on a seasonal scale allows a simplified approach on analyzing earlier onset of snow melt or later onset of snow fall, as these occur in the transition seasons. Spring SCD counts the number of days for each grid cell during the period from March 1st to May 31st for which the cell was snow-covered. Analogous, autumn SCD considers the period from September 1st to November 30th. In accordance with other studies (e.g., [4749]), a threshold of daily 2 cm SWE was applied to decide if the grid cell is snow-covered or not. Thus, spring SCD characterizes the snow free start, while autumn SCD characterizes the timing of snow accumulation onset after summer and before permanent snow cover occurs.

Snow cover and depth, albedo, and LST are highly dependent on land cover or vegetation type. Therefore, we also separately discuss them with respect to forest and nonforest grid cells. Based on a global 1 km-resolution dataset from Hagemann [43], a threshold of 50% forest fraction has been applied according to Figure 1.

4. Results and Discussion

4.1. Snow Cover and Albedo
4.1.1. Spatial Patterns in Circum-Arctic Domain

Snow Cover Duration (SCD). The spatial patterns of spring and autumn SCD are shown in Figure 2. The spring patterns show a strong zonally symmetric structure, but orographically influenced, as SCD is largely affected by surface temperature during the snow melt period. The shortest SCD is found over lower latitudes (e.g., over Northern Europe and west Russia); it ranges there from 10 days to 60 days. The longest SCD of more than 70 days occurs over the West Siberian Plains, the Taymir Peninsula, and the Canadian archipelago and in mountain areas. The timing of snow cover onset in autumn is temperature and precipitation driven. Therefore, the autumn SCD spatial patterns are more regionalized than those in spring. Largest SCD of more than 60 days occur here over the mountainous regions which experience very early snow fall.

Generally, the simulated spatial SCD patterns agree with those from CMC, but biases in absolute numbers occur. In spring, the model underestimates SCD of up to 25 days. Thus, the model shows the largest underestimation of up to 40% in the sub-Arctic areas. Smallest differences (of up to −4 days) are over the regions with the largest SCD (Taymir Peninsula, Lena river basin area, and Canadian archipelago). The largest underestimation of SCD in sub-Arctic areas in spring is caused by an underestimated snow amount, which might be related to an insufficient cyclone transported moisture and the localization of precipitation amount [10]. The general underestimation of SCD is further associated with biased melting onset and melting rate, which could be caused by the simplified snow parameterizations (such as the one-layer snow model or nonconsideration of snow aging) in HIRHAM5. In autumn, the model under- and overestimates (depending on the region) the SCD compared to CMC. The differences are of small-scale and manifested over mountainous regions, where the model tends to show higher SCD compared to CMC. The reasons could be, on the one hand, the higher model resolution and associated higher precipitation and snow depth and, on the other hand, the lower station density in CMC.

Snow Water Equivalent (SWE). The simulated spatial patterns of SWE are comparable with those from CMC but show some disagreement with GlobSnow and ERA-interim (Figures 3 and 4). Largest CMC-derived amount of SWE occur over the West Siberia plains and the Rocky Mountains, with SWE of more than 18 cm in spring and more than 4 cm in autumn. The model tends to overestimate the SWE in those regions by ca. 20–30% compared to CMC. The model biases compared to GlobSnow and ERA-interim are larger than CMC. GlobSnow shows the same large differences relative to CMC. Particularly, GlobSnow does not show a large amount of SWE over Alaska, the Canadian archipelago and the mountain areas. The differences between GlobSnow and CMC can reach more than 8 cm. This may be due to the sparse measurements and the simplicity of the retrieval algorithm in GlobSnow (the algorithm is designed particularly for boreal forest; [50]).

Snow Cover Fraction (SCF). Spring and autumn SCF spatial patterns are shown in Figures 5 and 6. For autumn, GlobSnow has only been averaged over September and October due to lack of optical satellite acquisitions which is caused by reduced incoming and reflected radiation at high latitudes.

In spring, clear zonally symmetric and orography-dependent structures appear, similarly to those seen in SCD. The largest snow cover in HIRHAM5 is found over the west Siberian plains and the Taymir Peninsula, the Kolyma River basin, and the Canadian archipelago. SCF reaches there up to more than 0.7. Other regions (e.g., Lena river basin, Northern Europe, and Alaska) show smaller SCF of 0.3–0.6. MODIS and GlobSnow SCF spring patterns are in good agreement to each other. Both show high values of SCF of 0.7–0.9 over the whole Arctic. Compared with MODIS, HIRHAM5 calculates generally similar spatial patterns but much smaller amplitudes of SCF within the whole domain in both seasons. The regionally dependent bias ranges between −0.4 and −0.1. The modeled spatial SCF patterns follow closely those of SWE because SCF is derived from SWE (see Section 2.5). Compared with GlobSnow, the HIRHAM5 biases are quite similar with the comparison against MODIS in both seasons. However in autumn, the HIRHAM5 bias compared to GlobSnow shows smaller scale biases of higher amplitude. This is caused by the lack of GlobSnow data due to reduced incoming and reflected radiation at high latitudes in this season.

Surface Albedo. The simulated surface albedo shows biases in spring and autumn, compared to both ERA-interim and MODIS (Figure 7). In spring, both “HIRHAM5 minus ERA-interim” and “HIRHAM5 minus MODIS” show a main feature: HIRHAM5 has a positive albedo bias compared to both datasets (of up to 0.4 compared with ERA-interim and of up to 0.3 compared with MODIS) over the areas where the forest cover fraction is relatively large (see Figure 1) and forest is mainly comprised of deciduous trees, that is, over the West Russian Arctic, Lena River basin, and parts of river basins in Alaska and Canada. In autumn, a similar spatial bias pattern is obvious in comparison with ERA-interim, while the bias in comparison with MODIS is largely latitude dependent. The results indicate that HIRHAM5 overestimates the surface albedo in these forest areas for both seasons. It is worth noting that this does not happen in the Scandinavia area where forest cover fraction is even higher, because the forest type is different: deciduous tree versus evergreen tree. Also there is no snow cover in Scandinavia during spring months despite being at high altitudes. Fallen leaves and branches on snow coverage, which are not considered in the albedo parameterization, would cause a much darker albedo despite snow coverage. This could be responsible for the positive bias over the Siberian region. This result is in agreement with Roesch and Roeckner [51], even though HIRHAM5 shows shorter spring SCD and underestimated SCF. Further, in spring, there are relatively large differences between the two biases “HIRHAM5 minus ERA-interim” and “HIRHAM5 minus MODIS” in specific areas like the Canadian archipelago and the Kola Peninsula. There HIRHAM5 shows large negative albedo biases (of up to −0.3) with respect to MODIS, but only small negative or even positive biases compared to ERA-interim. In autumn, large positive albedo biases (compared to MODIS) are located over the West Siberia Plain, the Taymir Peninsula, and the Canadian archipelago, which may come from the longer autumn SCD (Figure 2). Biases with respect to ERA-interim derived albedo are a result of the different parameterizations compared to HIRHAM5. As for MODIS-derived albedo, cloud cover affection is one important limitation. On the one hand, satellite observation is not continuous due to cloud cover; on the other hand, when there is snow cover, optical satellite products may contain undetected and unmasked cloud cover for snow and ice land surface [52].

4.1.2. Annual Cycle over the Siberian Subdomain

As shown, snow cover and albedo are influenced by the vegetation cover. Therefore, their annual cycle has been analyzed separately for forest and nonforest regions. To enable a comparison with CCLM simulations (see Section 2.6), the evaluation is restricted here to the Siberian subdomain (see Figure 1). Due to the difference of forest fraction between HIRHAM5 and CCLM, we consider the forest cover fraction between 0.5–1.0 in both model as forest area and forest cover fraction between 0–0.5 in both model as nonforest area.

Snow Water Equivalent (SWE). Figure 8 shows significant differences in the annual cycle of SWE, particularly for forest areas. Over forest, the observational datasets ERA-interim, CMC, and GlobSnow show differences from November to May. While CMC and GlobSnow largely agree, ERA-interim shows systematically higher SWE by ca. 5 cm. For nonforest areas, the observational data are closer to each other, but GlobSnow shows smaller SWE compared to CMC and ERA-interim data from February to April.

For both regions, HIRHAM5 tends to underestimate SWE, but particularly in forest areas. The underestimation of SWE in the forest region is also presented by CCLM, but mainly during the winter months. Except for May, both models are better in reproducing SWE over nonforest areas than over forest areas. This may be caused by the complexity of forest vegetation.

Snow Cover Fraction (SCF). Figure 9 shows that for both forest and nonforest areas, MODIS and GlobSnow-derived SCF agree well in spring (March–May), but differ significantly in autumn (October-November). A systematic underestimation of HIRHAM5-simulated SCF especially in spring has already been discussed (see Section 4.1.1; Figures 5 and 6), although the model calculation of SCF considers various surface properties (e.g., orography, surface type, and vegetation cover). Over both areas, forest and nonforest areas, HIRHAM5 strongly underestimates SCF by 0.2–0.5 in April and May compared to both MODIS and GlobSnow. In autumn, the models agreement with the observations is better, particularly with GlobSnow. Compared to MODIS, the model again shows an underestimated SCF. In October and November, HIRHAM5 agrees with GlobSnow, while MODIS shows higher SCF especially for nonforest areas. Takala et al. [30] argue that GlobSnow is of better quality than MODIS in the Russian area. Compared with HIRHAM, CCLM shows a better simulation of the observed annual cycle of SCF. The results agree well with MODIS and GlobSnow in spring, though a small overestimation in autumn compared with GlobSnow exists. The better CCLM performance could be connected with the separate consideration of deciduous and evergreen forest, which better accounts for their influence on the snow and albedo over forests.

Surface Albedo. Though HIRHAM5 underestimates SWE and SCF, the surface albedo is largely overestimated for the Siberian forest region in spring by more than 0.2 (about 60%–70%), compared to MODIS and ERA-interim (Figure 10). Over forest regions, the snow-masking effect of fallen leaves and branches that darkens the snow albedo seems to be an important factor to affect the surface albedo simulation for snow covered areas. Note that lack of solar radiation lessens the importance of the snow-masking effect on the albedo during the heavy snow season from October to March [52]. However, the snow-masking effect causes a darker observed albedo. This may play an important role during the transition months of May and September when solar radiation is available. Over nonforest region, the HIRHAM5 albedo shows better agreement with MODIS, especially for the period from April to October. The albedo overestimation in spring is much smaller than over forest. HIRHAM5 reproduces the albedo decrease from April to June with the same rate as MODIS over nonforest. Both for forest and for nonforest areas, the CCLM albedo shows better agreement with MODIS than HIRHAM5, associated with the sophisticated, snow-dependent albedo scheme in CCLM (Section 2.6).

During the snow-free months from June to September, the HIRHAM5 and CCLM albedo data agree with the observational data of MODIS and ERA-interim both for forest and for nonforest areas. During the snow free period, the model albedo depends only on temperature and forest fraction, and the results show that the applied parameterizations are adequate.

It is worth noting that the ERA-interim albedo over the Siberian forest region shows practically no annual cycle but is more or less fixed to 0.2. This is due to the snow albedo scheme in the ECMWF forecast model, which limits the deep-snow albedo for boreal forest to 0.2 [38].

4.2. Land Surface Temperature (LST)

Figure 11 shows the comparison of HIRHAM5 LST and ERA-Interim LST. The spring and autumn LST spatial bias patterns (Figure 11) resemble the albedo bias patterns (Figure 7). In spring, HIRHAM5 is up to 4 K colder over the Siberian forest region. The autumn LST bias ranges between −3 K and 3 K and is obviously not related to the vegetation cover.

By comparing the simulated LST annual cycle with ERA-interim over the Siberian region (Figure 12), we recognize that both models simulate the LST annual cycle well and the deviations from March to October are smaller than in the other months. In winter, the models tend to be warmer than ERA-interim, particularly CCLM. The LST agreement in spring and autumn is good, except for HIRHAM5 over forest areas in March and April. Here, the better SCF and albedo simulation in CCLM affect the better CCLM LST results. Though the HIRHAM5 forest albedo is higher for most of the months, the temperature bias is only limited to March and April, because the snow albedo affection to LST is obvious when the solar radiation becomes stronger (e.g., March and April), while the solar radiation (and albedo effect) in winter is quite small. And this effect will decrease with progressive snow melting.

At the beginning, we aimed to compare the simulated LST with the MODIS-derived LST. However, the following limitations have to be considered such that we decided not to do so: LST MODIS observational data are only for clear-sky cases. Thus, it is difficult to compare satellite data with modeled results for whole time period. Even if the comparison would be focused on clear-sky LST, due to different temporal and spatial samplings between satellite and model, such a comparison is not firm. Also because our HIRHAM5 simulation is relatively short, if one compares the two datasets, the sampling of clear-sky versus cloudy-sky data needs to be considered carefully in order to keep the temporal uniformity of the comparison. Therefore, we cannot directly compare LST-derived by MODIS and modeled LST. With the state of observational data, comparison between model and satellite data has to be limited to LST-driving parameters like SCF and SWE. Therefore, the modeled LST is only compared with ERA-interim data due to limitations of satellite-derived LST.

5. Summary and Conclusion

In this study, HIRHAM5 modeled Arctic land snow cover; albedo and LST are evaluated in spring and autumn using various satellite and reanalysis data. In addition, the HIRHAM data is compared with CCLM as second RCM.

HIRHAM5 can generally capture the main characteristics of the spatial patterns and the annual evolution of SWE, SCF, albedo, and LST, although significant biases are detected.

It has been shown that the simulation of snow cover characteristics (SWE, SCF, and albedo) over forest areas is difficult for the model, particularly in spring during the snow melt. The albedo and temperature bias in spring for Siberian deciduous tree forest is obvious in HIRHAM5. This could be caused by neglecting the snow-masking effect of fallen leaves and branches for deciduous tree forest in the SCF parameterization, since SCF is a key variable to control the albedo simulation in spring. The effect of improved SCF and albedo simulations on the improved LST is shown by the CCLM results. Less snow cover leads to lower albedo, which causes more solar radiation absorption and this further accelerates snow melting. This process is a positive feedback and small biases can be amplified and then reflected in the LST bias. The implication of this result is that the simulated spring and autumn LST can be mainly improved by improving the regional snow cover and thus the surface albedo.

Our finding of the unsatisfactory nature of the Siberian forest albedo parameterization is in agreement with Roesch and Roeckner [51], who pointed out this weakness for the applied ECHAM5 parameterization. Thus we suggest for the albedo parameterization a separate consideration of forest types (e.g., such as in CCLM) and the implementation of the effect of fallen leaves and branches on snow. The simulated nonforest SWE, SCF, and albedo are in reasonable agreement with the observations. The applied temperature and forest fraction dependent albedo parameterization works well.

There are also limitations in the use of the selected datasets. Differences between the different data sets have been presented. For instance, SWE derived from CMC and GlobSnow and autumn SCF derived from MODIS and GlobSnow differ considerably. For such cases it is difficult to evaluate the model simulations and the comparisons need to take all available observational data sets into account. In addition, optical satellite retrievals are difficult for dense cloud cover and for cloud cover over snow and ice surfaces, which will cause contamination of data [52].

The evaluation of our study is restricted to only 3 years: a time span where all data products overlap. However, this short period is sufficient when looking at mean values and the evaluation results do not differ when more years are included (e.g., 10 years is tested). Nevertheless, for further climate statistics the data of long-term simulations should be included in future work.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This work has been funded by the Helmholtz Climate Initiative REKLIM (Regional Climate Change), a joint research project of the Helmholtz Association of German research centers (HGF). The authors thank Ines Hebestadt’s support in the HIRHAM5 simulation.