Snow and ice on Bear Lake (Alaska) sensitivity experiments with two lake ice models
ByT I D O S E M M L E R1 , 2* , B I N C H E N G3, Y U Y A N G4 , 5 a n d L A U R A R O N T U3,
1Alfred Wegener Institute for Polar and Marine Research, Am Handelshafen 12, 27570 Bremerhaven, Germany;
2Met E´ireann, Glasnevin Hill, Dublin 9, Ireland; 3Finnish Meteorological Institute, PO Box 503, FI-00101 Helsinki, Finland;4State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology,
Dalian 116024, China;5Department of Physics, University of Helsinki, FI-00014 Helsinki, Finland
(Manuscript received 8 April 2011; in final form 2 February 2012)
A B S T R A C T
Snow and ice thermodynamics of Bear Lake (Alaska) are investigated with a simple freshwater lake model (FLake) and a more complex snow and ice thermodynamic model (HIGHTSI). A number of sensitivity experiments have been carried out to investigate the influence of snow and ice parameters and of different complexity on the results. Simulation results are compared with observations from the Alaska Lake Ice and Snow Observatory Network. Adaptations of snow thermal and optical properties in FLake can largely improve accuracy of the results. Snow-to-ice transformation is important for HIGHTSI to calculate the total ice mass balance. The seasonal maximum ice depth is simulated in FLake with a bias of 0.04 m and in HIGHTSI with no bias. Correlation coefficients between ice depth measurements and simulations are high (0.74 for FLake and 0.9 for HIGHTSI). The snow depth simulation can be improved by taking into account a variable snow density. Correlation coefficients for surface temperature are 0.72 for FLake and 0.81 for HIGHTSI. Overall, HIGHTSI gives slightly more accurate surface temperature than FLake probably due to the consideration of multiple snow and ice layers and the expensive iteration calculation procedure.
Keywords: snow on lake ice, thermodynamic ice model, heat conductivity of snow, albedo, snow to ice transformation
1. Introduction
In the boreal countries, a realistic representation of freezing and melting of lakes and snow on lake ice is important for numerical weather prediction (NWP) and regional climate model applications. Due to changing climate conditions, especially at higher latitudes, the lake ice season can be strongly influenced, which in turn results in a feedback on the regional climate (Brown and Duguay, 2010). It is important to couple lake ice models to climate and NWP models to simulate the two-way feedback between lakes and regional weather/climate (MacKay et al., 2009). In such a coupled model set-up, it turns out that lakes can have a net warming effect leading to a convective precipitation increase of 2040% in late summer (Samuelsson et al., 2010). With increasing computing power, higher resolutions in limited area models become feasible. Lakes and their impact on the local weather and climate can be resolved in model grids.
Sufficient data about the lake surface state are rarely available for data assimilation in NWP or for evaluation of lake models. Detailed observations of lake ice and snow depth and temperature are sparse and mostly restricted to measurement campaigns. Therefore, a precise simulation of snow and ice cover is important but challenging.
Modelling lake ice and snow has been carried out in several studies. Heron and Woo (1994) investigated the decay of lake ice with a lake ice energy-balance model.
Vavrus et al. (1996) compare simulated freezing and melting times with climatologically observed data using a simple lake model. Duguay et al. (2003) and Jeffries et al. (2005) use the Canadian Lake Ice Model CLIMo to simulate lake ice over Alaskan ponds and compare freezing and melting times as well as ice depth to observed data. Mironov et al.
(2010) have implemented freshwater lake model (FLake) into the NWP model COSMO. They note that a quantita- tive evaluation of snow and ice in the model is challenging.
Eerola et al. (2010) have reported implementation of FLake as a parameterisation scheme into the NWP model HIgh
*Corresponding author.
email: tido.semmler@awi.de
Tellus A 2012. #2012 T. Semmler et al. This is an Open Access article distributed under the terms of the Creative Commons Attribution-Noncommercial 3.0 Unported License (http://creativecommons.org/licenses/by-nc/3.0/), permitting all non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
1
Citation: Tellus A 2012,64, 17339, DOI: 10.3402/tellusa.v64i0.17339 P UBL ISH ED B Y TH E INT ERNA TIONA L M ETEOROLOGIC AL I NSTI TUTE I N S TOC KHOL M
METEOROLOGY AND OCEANOGRAPHY
(page number not for citation purpose)
Resolution Limited Area Model (HIRLAM). They have found out that the ice cover can have an important local impact on the atmospheric surface boundary layer.
Evolution of the snow layer may add further complica- tion to the lake ice processes. For example, the change of snow depth during the melt season is strongly dependent on the surface albedo. In addition to its insulation effect, snow may contribute to the total ice mass via snow-ice and superimposed ice formation (Cheng et al., 2003). This process has largely been observed in the Northern hemi- sphere ice-covered lakes and seas (Leppa¨ranta, 1983;
Kawamura et al., 1997; Nicolaus et al., 2003). Climate change scenarios suggest increasing precipitation in the Arctic (Christensen and Christensen, 2007). Therefore, snow may become more important for ice thermodynamics.
The aim of this study is to investigate how complex a lake ice model should be in order to realistically simulate snow and ice depth and temperature, how sensitive a lake ice model is to the definition of ice and snow properties, and which further processes might be necessary to consider in the parameterisations. Investigated parameters include albedo of fresh and melting snow and ice as well as snow density and heat conductivity. Two different lake ice models are compared: a simple lake model with two water layers, one ice and one snow layer (FLake: Mironov, 2008;
Mironov et al., 2010) and a sea and lake ice model with multiple ice and snow layers (HIGHTSI: Launiainen and Cheng, 1998; Cheng et al., 2003). Both models are based on the solution of the heat conduction equations and do not consider the effects of ice dynamics.
A first step towards a lake ice model comparison is made in this study. While a lake model intercomparison project has already been initiated (Stepanenko et al., 2010), the emphasis is on lake temperature profiles rather than ice and snow processes. Our comparison is carried out for a small lake in the south-east of Alaska for which observations of snow depth, snow density, surface temperature, snowice interface temperature, conductive heat flow and snow water equivalent are measured in intervals between once every week and once every month (Jeffries and Morris, 2006; ALISON, 2011) for the first time allowing a quantitative comparison of snow and ice properties be- tween simulations and observations. Through the compar- ison of two models the most sensitive parameters can be figured out and improvements introduced into both models.
In Section 2 model descriptions of FLake and HIGHTSI are given, focusing on the snow parameterisations. Section 3 describes the observation data used for validation, section 4 the set-up of the experiments and section 5 their results.
In Section 6 conclusions are drawn.
2. Models
2.1. FLake
FLake is a thermodynamiclake model. The water module considers the heat and kinetic energy budget of two layers, that is, upper mixed-layer and basin bottom. FLake is able to predict the vertical temperature structure and mixing conditions at various depths. The lake depth is an important input parameter for FLake since it determines the heat capacity, the temperature development, the freeze- up date and the heat flux from water to ice. The concept of self-similarity (Kitaigorodskii and Miropolsky, 1970) is used to describe the temperature-depth curve for all layers in vertical. The mixing-layer depth calculation includes convective entrainment and wind-mixing effect for unstable and stable regimes, correspondingly. Both mixing regimes are treated considering the volumetric character of solar radiation heating. In order to compute momentum, sen- sible and latent heat fluxes at the lake surface, a para- meterisation scheme has been developed that accounts for specific features of the surface air layer over lakes (Mironov, 2008). In our configuration, the model considers four layers in vertical: two lake water layers, one ice and one snow layer. A detailed description of FLake is given by Mironov (2008).
We focus on improvement of snow description in FLake, which has not been thoroughly tested yet. Rather than implementing a more complex snow scheme on lake ice, special attention has been paid to parameterisation of snow properties, which have a major impact on snow and ice mass balance.
In the reference FLake, snow density (rs) is a function of snow depth:
qs¼minðqs max;qs min=ð1hsc=qwÞÞ (1) where rsmax is the maximum snow density (400 kg m3), rsminis the minimum snow density (100 kg m3) andhsis the snow depth (m), c is an empirical parameter (200 kg m4) and rw is the water density (1000 kg m3).
Snow density can vary between 100 and 400 kg m3. Snow accumulates through snowfall. With increasing snow accumulation the snow is compacting, that is, the snow density is increasing according to this parameterisation.
However, snow would actually get less dense when snow is melting, which is physically unrealistic. Furthermore, increasing snow depths from 0.01 to 1 m would increase the snow density only from 100 and 125 kg m3. Observa- tions for Bear Lake show values between 130 and 470 kg m3. Sturm and Liston (2003) noticed that snow density on the Alaska lakes is usually higher than that on the surrounding land. Average snow densities of 344 and 334 kg m3were obtained for spring 2000 and 2002 over 13
Alaskan lakes. These values are comparable to the average value of 320 kg m3 for winter snow on Arctic sea ice (Huwald et al., 2005).
The snow heat conductivity (ks) is calculated depending on the density and snow depth:
ks¼minðks max;ks minþhs1:3qs=qwÞ (2) where ksmax is the maximum snow heat conductivity (1.5 W (mK)1),ksminis the minimum snow heat conduc- tivity (0.2 W (mK)1). According to Eqs. (1) and (2), ks varies between 0.2 and 0.36 W (mK)1. This is in line with a recent investigation by Ashton (2011) who argues thatks
varies between 0.1 and 0.4 W (mK)1, with typical values from 0.2 to 0.3 W (mK)1for Northern American Arctic and Subarctic lakes. The spatial inhomogeneity, however, has an impact on the spatially averaged heat conductivity, especially when snow depth is less than 0.4 m (Sturm et al., 2002). We therefore define a variable called effective snow heat conductivity:
kseff¼ksþðkiksÞehsc (3) whereks is 0.14 W (mK)1,kiis the heat conductivity of lake ice (2.29 W (mK)1), c is an empirical constant (5 m1). Eq. (3) gives higher heat conductivity for thinner snow. For thin snow, there is likely to be bare ice in the vicinity because of wind effect. Therefore, the effective snow heat conductivity would be a mixture of heat conductivity of snow and ice within a unit area.
The surface albedo (asfc) in the reference FLake is calculated by:
asfc¼a1 1e95:6ðT0TsfcÞ=T0
þa2e95:6ðT0TsfcÞ=T0 (4) where a1refers to albedo for white ice (snow-ice) or dry snow equal to 0.6.a2is albedo for blue ice (congelation ice) or melting snow which was set to 0.1. T0 is the freezing temperature (273.15 K), Tsfc is the surface temperature.
Eq. (4) suggests the same albedo for both snow and ice.
Furthermoreasfcapproaches 0.1 whenTsfcis close toT0. The albedo values given by eq. (4) differ quite a lot from in situ measurements and may not be adequate to study snow and ice mass balance. Even observed albedo shows large variations between different research sites. For example, albedo values of 0.43 and 0.21 for white ice (snow-ice) and melting ice, respectively, were observed on the Great Lakes (Bolsenga, 1977). The albedo of melting blue ice varies from 0.2 to 0.4 (Prowse and Marsh, 1989) and from 0.4 to 0.55 (Heron and Woo, 1994) depending on ice crystal orientation. An average albedo of 0.38 was observed by Henneman and Stefan (1999) over a freshwater lake in Minnesota. The albedo for new snow can be as high as 0.81 in numerical model (Gardner and Sharp, 2010) or 0.83 from lake in situ measurements (Henneman and
Stefan, 1999) and for old dirty snow as low as 0.57 (Gardner and Sharp, 2010). For snow on sea ice, the albedo can be as high as 0.87 for dry fresh snow and 0.77 for melting snow (Perovich, 1996). There seem to be considerable uncertainties in albedo for lake snow and ice.
Adjustments of FLake model have been carried out with respect to snow density, heat conductivity and albedo.
Sensitivity experiments on depth and temperature of snow and ice are investigated (see Section 4 and Table 1 for the experiment definitions).
2.2. HIGHTSI
High resolution thermodynamicsnow and ice (HIGHTSI) model has been extensively used for sea ice studies (Launiainen and Cheng 1998, Cheng et al., 2006; Cheng et al., 2008) and has recently been applied to a lake (Yang et al., 2012, this issue). HIGHTSI focuses on snow and ice thermodynamics but does not treat water underneath.
Ice and snow depth, heat conductivity and temperature are simulated solving the heat conduction equation for multiple ice and snow layers. Apart from those basic model physics that have been used in several previous studies (e.g.
Maykut and Untersteiner, 1971; Gabison, 1987; Ebert and Curry, 1993), special attention is paid to the parameterisa- tions of the air-ice fluxes and the solar radiation penetrat- ing into the snow and ice. The atmospheric turbulent surface fluxes are parameterised taking the thermal strati- fication into account. The penetration of solar radiation into the snow and ice depends on the cloud cover, the albedo, the colour of the ice (blue or white) and optical properties of snow and ice. The global radiation penetrat- ing through the surface layer is parameterised, making the model capable of calculating subsurface melting quantita- tively (Launiainen and Cheng, 1998; Cheng et al., 2003).
Short- and long-wave radiative fluxes can either be para- meterised or prescribed based on observations or results of NWP models. The surface temperature is solved from a detailed surface heat/mass balance equation, which is defined as the upper boundary condition and also used to determine whether surface melting occurs. A heat and mass balance at the ice bottom serves as the lower boundary condition of the model.
Snow processes are considered thoroughly in HIGHTSI.
The external source of snow is precipitation. An assumed initial snow density (320 kg m3) was used to convert precipitation to snow depth. The existing snow density of 320 kg m3is modified afterwards due to ageing, according to Anderson (1976). The snow heat conductivity is para- meterised according to Sturm et al. (1997). When snow and ice are present, the Archimedes’ principle is used to calculate the ice freeboard. A slush layer is then formulated if freeboard tends to be positive. The slush layer will be
refreezing to snow-ice during cold condition. During the melting season, snow may be subject to surface and subsurface melting. The melting snow will be converted to a slush layer, which is added to the snowice interface.
This slush layer may refreeze before snow layer totally melts away. Thus, in HIGHTSI snow brings an insulation effect to reduce ice growth; snow also gives solid contribu- tion to ice mass balance. The calculation details are given in Cheng et al. (2003) and Yang et al. (2012, this issue).
HIGHTSI has incorporated 10 surface albedo parame- terisation schemes. Albedo parameterisations, with various degrees of complexity, tend to yield quite different results compared to the observations (Curry et al., 2001). The impact of albedo on large-scale sea ice simulation in the Arctic Ocean has suggested that albedo schemes with sufficient complexity can reproduce realistic basin-scale ice distributions (Liu et al., 2007). In this study we applied a more sophisticated scheme where albedo is parameterised according to temperature, snow and ice depth, solar zenith angle and atmosphericproperties (Briegleb et al., 2004).
A detailed description of HIGHTSI can be found in the study by Cheng et al. (2003).
3. Observation data over Bear Lake
In this study, we applied observations from Alaska Lake Ice and Snow Observatory Network (ALISON) (Jeffries and Morris, 2006; ALISON, 2011). The snow depth, snow density, snow and ice surface temperatures, conductive heat flux and snow water equivalent are measured weekly to monthly. Each time, the observations were carried out
every 5 m along a 100 m line. The readings were then averaged to represent seasonal time series of various parameters.
Bear Lake (60.218N, 149.368W) was selected for the sensitivity studies presented here, because it is one of the largest lakes within ALISON. Bear Lake is located 10 km north of Seward in the south-east of Alaska (Fig. 1). This lake is oriented from north to south with a maximum fetch of about 2.4 km. The mean water depth is 10 m.
There have been eight ice seasons of measurements at ALISON site so far. The extended winter season 2003/2004 can be regarded as an average winter season in terms of precipitation, temperature and snow and ice depth within the time period from 2003 to 2011, available from the ALISON (2011) dataset. It covers the longest time period among the available seasonal measurements. Additionally, the onset of ice melting is more pronounced than during the other periods. Therefore, we focus our modelling experi- ments on this winter.
Figure 2 shows the observed snow and ice depth and snow density in Bear Lake in 2003/2004. The snow depth and density varies significantly. Snow mass increases after a snowfall event. Snow depth reduces while the snow density increases due to snow redistribution and compaction. The reduction of snow depth during cold season may also be linked with snow-to-ice transformation. For example, at the beginning of January, after a snowfall event, snow depth gradually decreases while the ice depth increases.
Because the air temperature has been mostly well below the freezing point, it is unlikely that the decrease of snow depth is totally caused by melting. The decrease of snow mass is Table 1. FLake and HIGHTSI experiments setup
Flake Parameters HIGHTSI Parameters
RF S1
Reference FLake
RF setup, butrs320 kg m3 (constant)
rs(initial)100 kg m3
(unless specified differently) H1 Without snow
rs(initial)320 kg m3 for existing snow rs(t): Anderson (1976).
S2 S1 setup, butks(t): eq. (3) ks(t): eq. (2) (unless specified differently)
H2 Snow added, but without snow-to-ice
transformation. ks(t): Sturm et al. (1997) S3 S2 setup, but modify albedo to
ice: a10.5 (0.6);
a20.3 (0.1) snow: a10.87 (0.6);
a20.77 (0.1) (numbers in brackets are values used in reference version)
ki2.29 W m1K1 Tf08C
Initial lake water Temperature: 48C
H3 Snow added, with insulation and snow- to-ice transformation.
ki2.03 W m1K1 Tf 0.18C Fw0.5 W m2
S4 S3 setup, but prescribed snow density from ALISON (2011) observations
Number of layers in the snow: 1
Number of layers in the ice: 1
H4 H3 setup, but prescribed snow density from ALISON (2011) observations
Number of layers in the snow: 10
Number of layers in the ice: 20
Tf: freezing temperature;Fw: heat flux from water to ice;t: time.
connected with slush formation and further refreezing to snow-ice. The last episode of ice growth may result from refreezing of melting snow, that is, superimposed ice formation. Overall, during early winter season, the pre- sence of snow on ice does not strongly reduce the ice growth rate as expected by snow insulation effect. This also applies to other ALISON winter season observations (20042011).
4. Setup of model experiments
4.1. Forcing data
HIgh Resolution Limited Area Model forecasts were run every 6 h from 1 November 2003 to 15 June 2004 to obtain the atmospheric forcing data. Since it is planned in the future to apply both FLake and HIGHTSI as parameter- isations within operational NWP models, it is beneficial to test the stand-alone models driven by the atmosphericdata from an NWP model.
HIgh Resolution Limited Area Model is a numerical short-range weather forecasting system, used operationally within the international HIRLAM programme partners of 11 European countries. HIRLAM is based on the hydro- staticprimitive equations; the dependent variables are temperature, wind component, humidity, surface pressure,
cloud water content and turbulent kinetic energy (Unde´n et al., 2002). We applied an experimental HIRLAM version close to v.7.3, which contains improved surface parameter- isations. As a limited area model, HIRLAM requires lateral boundary conditions from global/hemispheric mod- el, normally provided by the European Centre for Medium Range Weather Forecasting (ECMWF). For the present study, the lateral boundary conditions were taken from archived ECMWF analyses.
The HIRLAM experiment domain includes the whole continental area of Alaska (Fig. 1), which will enable us to consider also other lakes of the ALISON (2011) dataset in the future. HIRLAM forecasts were run with a horizontal resolution of 0.0688 (7.5 km) on a rotated latitude longitude grid.
The 5-, 6- and 7-h forecasts of each HIRLAM simulation were saved. The instantaneous parameters (surface pressure as well as lowest model level temperature, specific humid- ity, u- and v-components) were extracted for the HIRLAM grid point nearest to Bear Lake from the 6-h forecasts. The accumulated parameters (snowfall, downwelling long-wave and direct and diffuse short-wave radiation) were calcu- lated as differences between the 5- and 7-h forecasts to get values valid at every 6 h (00, 06, 12 and 18) coordinated universal time (UTC). It should be noted that HIRLAM may predict snowfall also at near-surface temperatures above 08C; for our FLake and HIGHTSI simulations, only snowfall at temperatures below 0.18C was considered.
A comparison with observation data from a meteorolo- gical station about 15 km north of Bear Lake (Seward:
60.35398N, 149.34838W; station ID USC00508377 from the U.S. National ClimaticData Center) shows a reasonable agreement of daily temperature and precipitation values.
The HIRLAM temperature tends to be higher than the observed (bias: 1.88C) while the HIRLAM precipitation is lower than the observed (bias:1.2 mm d1). However, it has to be taken into account that the meteorological station is further away from the sea and further up in the
Nov Dec Jan Feb Mar Apr May 0.0
0.2 0.4 0.6 0.8 1.0
0 200 400 600 800 ice thickness 1000
snow thickness
Snow and ice thickness (m)
Time of year 2003-2004 snow density
Snow density (kg m -3)
Fig. 2. Snow and ice depth and snow density measured on Bear Lake during the 2003/2004 ice season.
Fig. 1. Surface elevation [m] of the HIRLAM Alaska domain.
In the south-east of Alaska the location of the Bear Lake is indicated in white.
mountains compared to the lake, and that these differences could be genuine.
Figure 3 shows the time series of HIRLAM results that are used as external forcing for the lake ice models. Starting from mid-November (13 November), HIRLAM suggests cold conditions. The cold season extends until early April (7 April) when temperature slowly rises above the freezing point. The average temperature was 5.38C during the cold season (13 November7 April) in contrast to5.88C during the warm period (8 April15 June). The HIRLAM wind speed seems to be related to temperature during these two periods. The cold temperature was accompanied by relatively large wind variations, while during the warm period the wind was weak. During the winter, the total HIRLAM snow precipitation falling when the 2-m temperature was less than 0.18C was 303 mm water
equivalent compared to 306 mm water equivalent at the Seward station.
4.2. Setup of the FLake and HIGHTSI experiments Both models were initialised on 3 November 2003 and continuously run until 15 June 2004. Due to the different nature of the models, it is necessary to initialise the water temperature in FLake while this is not the case for the ice and snow model HIGHTSI, which does not simulate the lake water. FLake simulations were initialised with a constant lake water temperature of 48C. The simulations were run in the framework of the SURFace EXternalized model (Salgado and Le Moigne, 2010). This surface model includes FLake as a lake module along with other modules for land, sea and urban areas. Since a lake fraction of 100%
is assumed for the present simulations, the results are not influenced by any module other than FLake. Since HIGH- TSI does not simulate the lake water, the sensible heat flux from water to ice cannot be calculated but is prescribed as Fw0.5 W m2. This simple approximation that has to be made in an uncoupled ice model has been further investi- gated in Yang et al. (2012, this issue). Sensitivity studies suggest that increasing the heat flux from water to ice would have the strongest impact on ice depth during the late ice season. For technical reasons a minimum snow and ice depth value has to be maintained during the whole simulation. The lake is regarded to be ice and snow free at all times at which this minimum value is not exceeded. The minimum value is exceeded, that is, ice and snow start to form, when freezing and snowfall occur. At the end of the ice season, snow and ice retain their initial values and the lake is considered to be ice free (Yang et al., 2012, this issue). The different experiments are summarised in Table 1.
5. Results
5.1. Snow and ice depth
Time series of the observed and FLake-simulated ice and snow depth are presented in Fig. 4. The reference FLake run (RF) gives the poorest ice depth compared with observations. The ice depth is significantly underestimated.
This could be caused by the fact that snow is simulated too thick because the snow density close to 100 kg m3is too low. The temporal variation of snow depth is a result of snow fall and snow melting.
When a higher snow density of 320 kg m3is introduced (S1), the accumulated snow depth decreases. The insulation effect of snow is weakened and consequently the ice depth increases but the magnitude is still well below the observa- tions. The heat conductivity in S2 is increased for thin snow 0
10 20
-20 0 20
0 50 100
0 1 2
0 400 800
Nov Dec Jan Feb Mar Apr May Jun
0 200 400
Time of year 2003-2004
Va(ms-1 ) (a)
(b)
Ta(C C)
(c)
Rh(%)
precipitation precipitation (Ta<0.1 C) (d)
Prec (mm h-1 ) (e)
Qs(Wm-2 )
-2 Ql (Wm) (f)
Fig. 3. The forcing data from HIRLAM model; (a) wind speed;
(b) air temperature; (c) relative humidity; (d) original snow precipitation from HIRLAM (blue) and snow precipitation for T2mB0.18C (red); (e) downwelling short-wave radiation; and (f) downwelling long-wave radiation.
layer. This modification leads to a much faster ice growth in early winter because snow insulation is largely reduced.
In fact, this impact is so large that when snow is getting thicker in mid-winter, the ice is still growing strongly.
Accordingly, the simulated total ice depth is in better agreement with observations until mid-February. In the second half of February, the mean air temperature is 1.68C and the downward solar radiation at the latitude of Bear Lake (60.218N) is already notably large (300 W m2). Therefore, surface melting may occur and the surface albedo will change to 0.1 according to eq.
(4). This will further enhance surface melting and even- tually lead to a snow free condition. Since the albedo of melting ice is also set to 0.1, ice melts too early in spring.
A more reasonable snow and ice surface albedo set-up in S3 improves greatly the calculated ice depth. The temporal
variability of snow is too low and the errors lie on high peaks of snowfall. The experiment S4 gives a similar ice growth than S3, with a slightly thinner ice appearing in mid-winter. The temporal snow variation is in a better agreement with the observations because observed snow density has been used.
The results of HIGHTSI sensitivity experiments are given in Fig. 5. Without snow, the seasonal maximum ice depth matches the observation but the phase considerably differs. Ice depth reaches a maximum value (0.84 m) around 9 February, while observation suggests a maximum (0.87 m) in spring (10 April). The onset of melting, defined as the start of monotonously decreasing sea ice depth, is simulated on 3 April, some 7 d earlier than observed. The simulated freeze-up date, defined as the appearance of lake ice, was 13 November and the simulated break-up date,
Nov Dec Jan Feb Mar Apr May 0.0
0.2 0.4 0.6 0.8
1.0 Hs (ob); Hi (ob);
Hs (RF) Hi (RF) Hs (S1) Hi (S1)
S no w an d ic e th ic k n es s (m )
(a)
Nov Dec Jan Feb Mar Apr May 0.0
0.2 0.4 0.6 0.8 1.0
Hs (S2) Hi (S2)
S no w an d ic e th ic k n es s (m )
Time of year 2003-2004 Time of year 2003-2004
Time of year 2003-2004 Time of year 2003-2004
(b)
Nov Dec Jan Feb Mar Apr May 0.0
0.2 0.4 0.6 0.8 1.0
Hs (S3) Hi (S3)
S no w an d ic e th ic k n es s (m )
(c)
Nov Dec Jan Feb Mar Apr May 0.0
0.2 0.4 0.6 0.8 1.0
Hs (S4) Hi (S4)
S no w an d ic e th ic k n es s (m )
(d)
Fig. 4. Observed and FLake-simulated snow and ice depth Hs and Hi. (a) Reference FLake run and sensitivity experiment S1; (bd) sensitivity experiments S2S4.
defined as the total disappearance of lake ice, was 16 May.
In reality the lake ice and snow might appear and disappear at different times in different areas of the lake, which cannot be considered in a one-dimensional lake ice model.
Unfortunately there are no observations of freeze-up and break-up dates.
The snow insulation can reduce ice growth strongly as illustrated by experiment H2. In this case, the snow depth is calculated as a function of precipitation and densification as well as surface and subsurface melting. Until 12 February, the observed snow depth shows large variations while the simulated snow depth continuously increases. The model fails to reproduce two episodes of dramatic decrease of snow depth in particular, the observed snow-free conditions on 18 and 23 February. Several factors could have led to the observed snow-free condition. The meteor- ological data suggest strong wind, which may have led to a redistribution of snow. During the second half of February,
the downwelling solar radiation is already strong. It might not only melt the snow at the surface but also internally, leading to slush at the snowice interface. This slush layer may refreeze on the colder ice layer increasing its depth.
Those effects can actually change the snow depth both from the top and from the snow-ice boundary and eventually lead to a snow-free condition.
Only small temporal snow depth variations (decrease/
increase associated respectively with surface melting/snow- fall) are simulated from mid-February onwards while the observations show a pronounced accumulation of snow at the end of February and beginning of March. The rapid melting of snow at the end of April is a result of albedo feedback mechanism, which is typically seen in high latitudes in early April when solar radiation increases significantly.
Compared to experiment H1, the ice freezing and break- up dates are about the same in experiment H2. In
Nov Dec Jan Feb Mar Apr May 0.0
0.2 0.4 0.6 0.8
1.0 Hs (ob)
Hi (ob) Hi (H1)
S no w an d ic e th ic k n es s (m )
(a)
Nov Dec Jan Feb Mar Apr May 0.0
0.2 0.4 0.6 0.8 1.0
Hs (H2) Hi (H2)
S no w an d ic e th ic k n es s (m )
Time of year 2003-2004 Time of year 2003-2004
Time of year 2003-2004 Time of year 2003-2004
(b)
Nov Dec Jan Feb Mar Apr May 0.0
0.2 0.4 0.6 0.8 1.0
Hs (H3) Hi (H3)
S no w an d ic e th ic k n es s (m )
(c)
Nov Dec Jan Feb Mar Apr May 0.0
0.2 0.4 0.6 0.8 1.0
Hs (H4) Hi (H4)
S no w an d ic e th ic k n es s (m )
(d)
Fig. 5. Observed and HIGHTSI-simulated snow and ice depth. (ad) model experiments H1H4.
experiment H3, the decrease of snow depth is caused by formation of slush layer, surface and subsurface melting as well as densification. Differently from Cheng et al. (2006), we have applied a first-order approximation to let the slush layer (formed by positive freeboard due to heavy snow load or melting snow) to be refrozen totally. In this case, the simulated snow depth is much smoother compared to the result of H2. The snow-to-ice transformation results in larger ice depth leading to a good agreement with observations, especially during the growing phase up to the first half of April when snow is still present.
In H4 the onset of ice melt is simulated only 2 d later than observed. The development of sea ice depth is very well represented compared to the other HIGHTSI simula- tions. Decreasing snow depth along with increasing ice depth is simulated at several occasions (end of December to beginning of January as well as second and third week of January) in H4. This is due to the prescribed snow density from measurements, which varies largely in contrast to the smooth change, that is, a gradual increase from 320 to 450 kg m3, in the parameterisation in experiments H2 and H3. The last observed snow growth episode in March is not comparable with the simulations (H2, H3 and H4). In reality the observed event may refer to heavy snowfall.
Since we have applied HIRLAM results as external forcing, discrepancy between the observed and simulated by HIRLAM snow precipitation may lead to errors in the snow depth simulated by HIGHTSI. The effect of wind resulting in snow redistribution is not considered in HIGHTSI, which may introduce further inaccuracy to simulated snow depth.
Simple statistics of simulated versus observed snow and ice depth are given in Table 2. The statistics confirm that the more realisticparameterisations in the FLake sensitiv- ity experiments and the inclusion of more physical pro- cesses in HIGHTSI generally improve the simulation. The FLake reference simulation shows large errors and small correlations. The ice depth correlation substantially im- proves from S2 to S3 when introducing the more realistic albedo values. The prescribed snow density in experiments S4 and H4 improves the correlation coefficients of snow
depth compared to the other experiments; however, correlation coefficients are still very low indicating that snow depth is challenging to simulate even when prescrib- ing the snow density. Except for the correlation coefficient of snow depth, the more complex HIGHTSI model shows smaller errors and higher correlations compared to the simpler FLake model.
5.2. Snow and ice temperature
A comparison between observed and simulated surface temperature (snow or ice if snow is not present) is given in Table 3. Large errors seem to appear in the FLake experiments. A number of physical reasons could lead to large errors. In early January, the too cold surface temperature of RF is linked with too small snow density.
Once the physical parameters of snow in FLake get improved with the experiments S1S4, the temperature errors are reduced but are still large. On 22 January 2004, all model runs give small errors. This is because a melting surface is simulated that matches the reality. Therefore, the surface temperature difference is confined. In early spring, the RF gives an overestimation of surface temperature.
This is associated with the unrealistically small albedo. A small albedo brings too much solar radiation down to the surface. The consequence is a too warm surface and too early onset of snow melting (Fig. 4a). With modified albedo values in S3 and S4, the surface temperature gets too cold in early spring. In spite of this, snow and ice still melt too early. Therefore, it can be concluded that an improvement of snow and ice properties in FLake without inclusion of more physical processes such as snow-to-ice transformation only improves the simulation of the surface temperature to some extent. In our case, improvements are restricted to the winter season.
HIGHTSI-simulated surface temperatures are generally closer to the observed values although the cold bias in early January still exceeds 58C. It should be noted that there are comparably little changes in the surface tem- perature errors between the four HIGHTSI experiments;
in fact bias and RMSE even increase from H1 to H4.
Table 2. Bias, root mean square error (RMSE) and correlation coefficient of simulated snow and ice depths compared to observations.
For the bias negative values indicate an underestimation by the models and positive values an overestimation
RF S1 S2 S3 S4 H1 H2 H3 H4
Bias (m) Snow 0.01 0.13 0.10 0.10 0.06 0.06 0.04 0.01
Ice 0.44 0.32 0.16 0.03 0.09 0.04 0.13 0.05 0.01
RMSE (m) Snow 0.21 0.20 0.20 0.18 0.15 0.15 0.14 0.13
Ice 0.48 0.40 0.31 0.12 0.14 0.14 0.18 0.08 0.07
Correlation coefficient Snow 0.23 0.23 0.23 0.28 0.42 0.17 0.23 0.36
Ice 0.21 0.06 0.1 0.74 0.74 0.54 0.67 0.92 0.90
Improvements from H1 to H4 are restricted to snow and ice depth.
The snowice interface temperature (Table 4) is well simulated in the FLake reference and S1 simulations.
Because of the relatively low snow heat conductivity in these two simulations the too cold temperatures at the surface do not strongly propagate to the snowice interface while this happens in S2S4. HIGHTSI simulations gen- erally show a good agreement in terms of the snowice interface temperature. The multilayer approach allows a detailed description of the heat transfer from the surface to the icewater interface.
5.3. Comparison of FLake and HIGHTSI results For both FLake and HIGHTSI, the differences between simulated and observed ice depth tend to be reduced when more physical processes are taken into consideration. In this section we focus on the comparison between best results of FLake (S4) and HIGHTSI (H4). The seasonal maximum ice depth can be simulated with a bias of 0.04 m for FLake and no bias for HIGHTSI. The calculated monthly mean snow and ice depths and ice growth rates are given in Table 5. Large differences of simulated snow depth show up in spring (March, April) Table 3. Comparison between simulated and measured surface temperatures (8C) for all model runs. The comparisons match the specific date and time of measurement (local noon). Bias, root mean square error (RMSE) and correlation coefficient are also given
Temperature difference at surface:(TcalTob)
Date Tsfc(ob) RF S1 S2 S3 S4 H1 H2 H3 H4
18 December 2003 4.3 2.74 1.26 0.01 0.18 0.73 0.31 0.31 0.22 0.52
23 December 2003 10.6 3.44 3.57 4.35 4.14 3.47 3.18 2.04 1.98 1.80
2 January 2004 7.8 18.03 13.98 7.91 8.76 8.25 4.37 6.39 6.26 6.49
7 January 2004 9.0 13.86 12.81 9.75 11.06 9.90 3.72 4.41 5.15 4.40
10 January 2004 0.7 0.09 0.08 0.66 1.20 0.64 1.75 1.62 1.39 1.60
14 January 2004 12.0 1.72 1.35 0.89 1.43 1.44 1.82 2.42 2.31 2.71
22 January 2004 0.20 0.20 0.20 0.20 0.20 0.20 0.07 0.20 0.20 0.20
29 January 2004 9.8 6.23 3.99 2.86 1.19 0.79 3.67 3.63 3.42 3.35
12 February 2004 4.7 4.47 4.51 4.61 1.09 2.34 1.96 0.45 0.26 2.45
8 March 2004 4.4 3.38 3.22 1.26 4.84 4.70 0.39 1.39 0.28 3.08
18 March 2004 5.6 5.14 5.32 4.83 4.98 5.06 2.39 1.30 2.13 0.15
Bias 1.23 0.77 0.10 2.55 2.60 0.03 0.89 0.69 1.43
RMSE 7.61 6.40 4.59 4.96 4.62 2.57 2.86 2.91 3.04
Correlation coefficient 0.55 0.62 0.69 0.68 0.72 0.82 0.81 0.82 0.81
Table 4. Comparisons between simulated and measured temperature (8C) at the snowice interface for all model runs. The comparisons match the specific date and time of measurement (local noon). Bias and root mean square error (RMSE) are also given. Note that the experiment H1 is excluded because snow is not considered
Temperature difference at snowice interface:(TcalTob)
Date Tsin RF S1 S2 S3 S4 H2 H3 H4
18 December 2003 1.29 0.56 0.42 2.17 2.33 1.37 0.59 0.61 0.26
23 December 2003 0.75 0.32 0.63 3.08 3.23 1.36 1.40 1.43 0.95
2 January 2004 2.12 1.35 0.14 3.70 3.92 8.44 1.12 1.14 0.83
7 January 2004 5.41 3.73 0.20 7.38 8.21 7.84 0.42 0.36 0.68
10 January 2004 0.22 0.17 0.09 0.57 0.89 0.22 2.22 2.26 1.75
14 January 2004 1.32 0.56 1.15 5.33 5.59 3.77 1.70 1.73 1.02
22 January 2004 0.03 0.03 0.03 0.03 0.03 0.03 1.26 1.31 1.29
29 January 2004 2.18 1.93 0.88 2.02 2.95 2.54 0.06 0.02 0.46
12 February 2004 0.15 3.42 3.31 0.03 0.04 0.03
8 March 2004 1.01 7.57 7.24 0.34 0.30 0.47
18 March 2004 0.00 8.46 7.98 0.44 0.47 0.35
Bias 1.08 0.14 3.03 4.23 4.00 0.72 0.76 0.44
RMSE 1.36 0.50 3.24 5.04 5.08 1.10 1.13 0.88
Correlation coefficient 0.75 0.86 0.90 0.73 0.72 0.83 0.83 0.86
probably due to the different treatment of downward solar radiation in FLake and HIGHTSI. FLake has only one snow layer and the solar radiation is assumed to heat the snow surface, conduct heat through the single snow layer and melt snow at the surface. HIGHTSI has 10 snow layers, the downward solar radiation is assumed to distribute exponentially down to the whole snowpack (Beer’s law). Warming and melting is accounted for through the heat balance of surface layer (first snow layer on top) and snow internal warming and melting (caused by penetrating solar radiation). In HIGHTSI it can happen that only the uppermost snow layers are affected by melting, and the snow refreezes further down in the snow pack. Another reason for the difference in spring could be higher snow heat conductivity in FLake compared to HIGHTSI, which can lead to a more efficient warming of the snow-ice pack.
The simulated monthly mean ice depths are close to each other, especially in mid-winter (JanuaryMarch). In April and May, the differences are increasing because of the different time of the onset of melting. This can be seen very clearly from ice growth rate, that is, FLake simulates the onset of melting (April) one month in advance compared with HIGHTSI (May). Nevertheless, a similar ice melting rate of the order of 0.02 m d1is seen. However, S4 and H4 provide similar ice growth rate in other months. In November, H4 is practically snow-free which means that ice grows faster than in S4. In December, the difference between S4 and H4 is narrowed. In January, S4 suggests more ice growth while in February H4 grows faster.
In March, the difference is increasing with more ice growth for S4.
The ice growth is similar between FLake (S4) and HIGHTSI (H4) for different reasons. In S4, the effective snow heat conductivity (kseff) is closely associated with the actual snow depth in FLake. Its value may be quite large (close to ki), in particular when snow is thin therefore prompting ice growth. In H4, the ice growth is largely a result of the snow-ice formation. Fig. 6 shows HIGHTSI-
simulated snow-ice and superimposed ice during the whole period. When snow is flooded (negative freeboard), the slush snow is refrozen to snow-ice. The formation of superimposed ice is more pronounced in early spring when melting snow refreezes at the snow surface and under it.
The observed time series and the time series of surface and snowice interface temperature simulated by S4 and H4 are shown in Fig. 7. FLake and HIGHTSI are able to catch temperature changes correctly, particularly for a warm surface. While the simulated surface temperature of FLake is very similar to that of HIGHTSI, the simulated snowice interface temperature is clearly colder in FLake than in HIGHTSI. In FLake, a strong heat flux through the snowpack is caused by the large effective snow heat conductivity.
6. Summary and conclusions
Two numerical models (FLake and HIGHTSI) have been used to investigate snow and ice thermodynamics. Sensi- tivity experiments have been carried out for Bear Lake, close to Seward in the south-east of Alaska. This choice was made in order to have access to a comparably detailed Table 5. Monthly mean of simulated ice depth and ice growth rate by FLake (S4) and HIGHTSI (H4)
Simulated monthly mean snow depth (m)
Simulated monthly mean ice depth (m)
Simulated monthly mean ice growth rate (cm d1)
S4 H4 S4 H4 S4 H4
November 0.06 0.02 0.12 0.21 1.0 1.3
December 0.14 0.11 0.36 0.42 0.5 0.4
January 0.16 0.19 0.58 0.58 0.6 0.4
February 0.17 0.23 0.69 0.69 0.2 0.3
March 0.09 0.24 0.76 0.81 0.4 0.2
April 0.03 0.11 0.56 0.81 2.5 0.7
May 0 0 0.03 0.31 0 2.1
-0.1 0.0 0.1 0.2 0.3
Nov Dec Jan Feb Mar Apr May
-1.5 -1.0 -0.5 0.0
-0.25 -0.20 -0.15 -0.10 -0.05
*10-3 0.00 (b)
Freeboardand snow-ice(m)
freeb oard snow-ice (a)
Melting snow(m)
T ime of year, 2003-2004 melting snow
Superimposed ice (m)
superimposed ice
Fig. 6. Snow-ice and superimposed ice simulated by HIGHTSI.
data set on lake snow and ice depth and temperature, which is available from the ALISON (2011). HIRLAM forecasts over Alaska have been created and used to drive FLake and HIGHTSI. Special attention has been paid to the physical parameterisations of the thermal and optical properties of snow as well as to the snow-to-ice transfor- mation process. We have assessed how realistically the snow accumulation and freezing and melting of the lake is treated in the models. The original parameterisations in FLake define a too small snow density, which may lead to a significant underestimation of ice depth. With the correc- tion of albedo and snow density and with the inclusion of the effective snow heat conductivity, FLake results can be substantially improved. HIGHTSI results are substantially improved when the snow-to-ice transformation is taken into account.
Several parameters and processes have been found to be important for a realisticsimulation of ice and snow:
1. the snow density, which should be variable con- sidering ageing and compacting of snow
2. the snow heat conductivity, which is closely related to the snow density and for which inhomogeneities in the snow cover due to redistribution by wind or heterogeneous snow metamorphism should be con- sidered
3. the albedo of snow and ice, which strongly influ- ences snow and ice melting in spring
4. the transformation of snow into ice through flood- ing caused by strong snowfall on thin lake ice and refreezing, or through snow melting and refreezing With the improvements the seasonal ice cover can be well simulated. The seasonal maximum ice depth is simulated with a bias of 0.04 m for Flake and with no bias for HIGHTSI. However, modelling of the snow depth is still a
challenging task, especially during the growth of snow.
This is because the increase of snow depth is strongly connected with the snowfall events, which are given as external forcing for FLake and HIGHTSI. The accuracy of snowfall prediction depends on the NWP model. On the other hand, in situ snow precipitation is also prone to errors, that is, the snow gauge tends to underestimate the local snow accumulation. Therefore, it is understandable that the models may contain large errors during the snowfall period. An additional complication may be added because of the wind drift as well as snow metamorphism affecting the snow density. Indeed, with a prescribed snow density the simulated results can be improved to some extent. In contrast, melting of snow is directly linked with thermodynamic processes handled by the snow/ice model.
Therefore, there is potential to improve the simulation of snow depth in the melting season.
Overall, HIGHTSI gives more accurate surface tempera- ture than FLake, as compared with the measurements.
HIGHTSI applies an iterative procedure to calculate the surface temperature, while FLake applies a simpler ap- proximation. Iteration can be considered an optimal method for surface temperature calculation. It may, how- ever, be too expensive to be included in the snow and ice parameterisations for the current NWP models.
In future, the advantages of both FLake and HIGH- TSI models could be combined by coupling the two models using HIGHTSI for the lake ice and snow simulation and FLake for the water underneath. This new model could be coupled to a NWP or climate model to simulate the real-time interaction between lakes and the atmosphere.
7. Acknowledgements
The third author (YY) was supported in early 2011 by a CIMO scholarship from Ministry of Education in Finland.
The HIGHTSI modelling was supported by the National Natural Science Foundation of China (grant no. 51079021, 40930848) and the Norwegian Research Council (grant no.
193592/S30). The field data are courtesy of the ALISON (http://www2.gi.alaska.edu/alison/ALISON_pubs.html).
The authors thank Dr. Ekaterina Kurzeneva for valuable discussions during the preparation of this manuscript. The suggestions by two anonymous reviewers led to significant improvements of the manuscript.
References
ALISON. 2011. Alaska Lake Ice and Snow Observatory Network.
Online at http://www.gi.alaska.edu/alison/ALISON_measures.
html.
-25 -20 -15 -10 -5 0
Nov Dec Jan Feb Mar Apr May
-15 -10 -5 0
Tsfc(C)
S4 H4 ob (b)
Tsi (C)
T ime of year, 2003-2004 (a)
Fig. 7. The time series of observed (black points) and simulated surface temperature (Tsfc), snowice interface temperature (Tsi) from S4 (blue line) and H4 (red line).
Anderson, E. A. 1976.A Point Energy and Mass Balance Model of a Snow Cover. Technical Report NWS 19, 150 pp., Natl.
Oceanic and Atmos. Admin., Washington, D. C.
Ashton, G. D. 2011. River and lake ice thickening, thinning and snow ice formation.Cold Regions Sci. Technol.68, 319.
Briegleb, B. P., Bitz, C. M., Hunke, E. C., Lipscomb, W. H., Holland, M. M. and co-authors. 2004.Scientific Description of the Sea Ice Component in the Community Climate System Model, Version Three. Technical Report NCAR/TN-463STR, National Center for AtmosphericResearch, Boulder, CO, 78 pp.
Bolsenga, S. J. 1977. Preliminary observations on the daily variation of ice albedo.J. Glaciol.18, 517521.
Brown, L. C. and Duguay, C. R. 2010. The response and role of ice cover in lake-climate interactions. Prog. Phys. Geogr. 34, 671704.
Cheng, B., Launianen, J. and Vihma, T. 2003. Modelling of Superimposed ice formation and Sub-Surface melting in the BalticSea.Geophysica39, 3150.
Cheng, B., Vihma, T., Pirazzini, R. and Granskog, M. A. 2006.
Modelling of superimposed ice formation during the spring snowmelt period in the BalticSea.Ann. Glaciol.44, 139146.
Cheng, B., Zhang, Z., Vihma, T., Johansson, M., Bian, L. and co- authors. 2008. Model experiments on snow and ice thermo- dynamics in the Arctic with CHINARE 2003 data.J. Geophys.
Res.,113, C09020.
Christensen, J. H. and Christensen, O. B. 2007. A summary of the PRUDENCE model projections of changes in European climate by the end of this century.Clim. Change81, 730.
Curry, J. A., Schramm, J., Perovich, D. and Pinto, O. 2001.
Application of SHEBA/FIRE data to evaluation of snow/ice albedo parameterizations.J. Geophys. Res.106, 1534515355.
Duguay, C. R., Flato, G. M., Jeffries, M. O., Me´nard P., Morris, K. and co-authors. 2003. Ice cover variability on shallow lakes at high latitudes: model simulations and observations.Hydrol.
Proc.17, 34653483.
Ebert, E. E. and Curry, J. A. 1993. An intermediate one- dimensional thermodynamicsea ice model for investigating ice-atmosphere interaction. J. Geophys. Res. 98(C6), 10085 10109.
Eerola, K., Rontu, L., Kourzeneva, E. and Shcherbak, E. 2010. A study on effects of lake temperature and ice cover in HIRLAM.
Boreal. EnvRes.15, 130142.
Gabison, R. 1987. A thermodynamicmodel of the formation, growth, and decay of first-year sea ice.J. Glaciol.33, 105109.
Gardner, A. S. and Sharp, M. J. 2010. A review of snow and ice albedo and the development of a new physically based broad- band albedo parameterization.J. Geophys. Res. 115, F01009, doi:10.1029/2009JF001444.
Henneman, H. E. and Stefan, H. G. 1999. Albedo models for snow and ice on a freshwater lake.Cold Regions Sci. Technol.29, 3148.
Heron, R. and Woo, M. K. 1994. Decay of a High Arctic lake-ice cover: observations and modelling.J. Glaciol.40, 283292.
Huwald, H., Tremblay, L.-B. and Blatter, H. 2005. Reconciling different observational data sets from Surface Heat Budget of the Arctic Ocean (SHEBA) for model validation purposes.
J. Geophys. Res.110, C05009.
Jeffries, M. O. and Morris, K. 2006. Instantaneous daytime conductive heat flow through snow on lake ice in Alaska.
Hydrol. Proc.20, 803815.
Jeffries, M. O., Morris, K. and Duguay, C. R. 2005. Lake ice growth and decay in central Alaska: observations and computer simulations compared.Ann. Glaciol.40, 195199.
Kawamura, T., Ohshima, K. I., Takizawa, T. and Ushio, S. 1997.
Physical, structural, and isotopic characteristics and growth processes of fast sea ice in Lutzow-Holm Bay, Antarctica.
J. Geophys. Res.102(C2), 33453355.
Kitaigorodskii, S. A. and Miropolsky, Y. Z. 1970. On the theory of the open ocean active layer. Izv. Akad. Nauk SSSR. Fizika Atmosfery I Okeana6, 178188.
Launiainen, J. and Cheng, B. 1998. Modelling of ice thermo- dynamics in natural water bodies.Cold Reg. Sci. Technol.27(3), 153178.
Leppa¨ranta, M. 1983. A growth model for black ice, snow ice and snow thickness in subantarctic basins.Nordic Hydrol.14, 5970.
Liu, J., Zhang, Z., Inoue, J. and Horton, R. M. 2007. Evaluation of snow/ice albedo parameterizations and their impacts on sea- ice simulations.Int. J. Climatol.27, 8191.
MacKay, M. D., Neale, P. J., Arp, C. D., De Senerpont Domis, L.
N., Fang, X., and co-authors. 2009. Modeling lakes and reservoirs in the climate system.Limnol. Oceanogr. 54, 2315 2329.
Maykut, G. A. and Untersteiner, N. 1971. Some results from a time-dependent thermodynamicmodel of sea ice. J. Geophys.
Res.76, 15501575.
Mironov, D. V. 2008. Parameterization of Lakes in Numerical Weather Prediction. Description of a Lake Model. Technical Report of COSMO, No. 11. Deutscher Wetterdienst, Offenbach am Main, 41 pp.
Mironov, D., Heise, E., Kourzeneva, E., Ritter, B., Schneider, N.
and co-authors. 2010. Implementation of the lake parameteriza- tion scheme FLake into the numerical weather prediction model COSMO.Boreal Env. Res.15, 218230.
Nicolaus, M., Haas, C. and Bareiss, J. 2003. Observations of superimposed ice formation at melt-onset on fast ice on Kongsfjorden, Svalbard.Phys Chem Earth28, 12411248.
Perovich, D. K. 1996. The optical properties of sea ice.CRREL Monogr.961, 25 pp.
Prowse, T. D. and Marsh, P. 1989. Thermal budget of river ice covers during break-up.Can. J. Civil Eng.16, 6271.
Salgado, R. and Le Moigne, P. 2010. Coupling of the FLake model to the Surfex externalized surface model.Boreal Env. Res.15, 231244.
Samuelsson, P., Kourzeneva, E. and Mironov, D. 2010. The impact of lakes on the European climate as simulated by a regional climate model.Boreal Env. Res.15, 113129.
Stepanenko, V. M., Goyette, S., Martynov, A., Perroud, M., Fang, X. and co-authors. 2010. First steps of a Lake Model Inter- comparison Project: LakeMIP.Boreal Env. Res.15, 191202.
Sturm, M., Holmgren, J., Ko¨nig, M. and Morris, K. 1997. The thermal conductivity of snow.J. Glaciol.43, 2641.
Sturm, M. and Liston, G. E. 2003. The snow cover on lakes of the Arctic Coastal Plain of Alaska.J. Glaciol.49(166), 370380.
Sturm, M., Perovich, D. K. and Holmgren, J. 2002. Thermal conductivity and heat transfer through the snow and ice of the Beaufort Sea. J. Geophys. Res. Oceans 107(C21), 8043, doi:10.1029/2000JC000409
Unde´n, P., Rontu, L., Ja¨rvinen, H., Lynch, P., Calvo, J. and co- authors. 2002. The HIRLAM-5 scientific documentation. 144pp.
Online at http://hirlam.org and SMHI, S-60176 Norrko¨ping, Sweden.
Vavrus, S. J., Wynne, R. H. and Foley, J. A. 1996. Measuring the sensitivity of southern Wisconsin lake ice to climate variations and lake depth using a numerical model.Limnol. Oceanogr.41, 822831.
Yang, Y., Leppa¨ranta, M., Cheng, B., Li, Z. and Rontu, L. 2012.
Numerical modelling of snow and ice thicknesses in Lake Vanajavesi, Finland.Tellus,64A, 17202.
