Numerical Modelling of Snow and Ice Thicknesses in Lake Vanajavesi, Finland
Total Page:16
File Type:pdf, Size:1020Kb
View metadata, citation and similar papers at core.ac.uk SERIES A brought to you by CORE DYNAMIC METEOROLOGY provided by Helsingin yliopiston digitaalinen arkisto AND OCEANOGRAPHY PUBLISHED BY THE INTERNATIONAL METEOROLOGICAL INSTITUTE IN STOCKHOLM Numerical modelling of snow and ice thicknesses in Lake Vanajavesi, Finland By YU YANG1,2*, MATTI LEPPA¨ RANTA2 ,BINCHENG3,1 and ZHIJUN LI1, 1State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China; 2Department of Physics, University of Helsinki, PO Box 48, FI-00014 Helsinki, Finland; 3Finnish Meteorological Institute, PO Box 503, FI-00101, Helsinki, Finland (Manuscript received 27 March 2011; in final form 7 January 2012) ABSTRACT Snow and ice thermodynamics was simulated applying a one-dimensional model for an individual ice season 2008Á2009 and for the climatological normal period 1971Á2000. Meteorological data were used as the model input. The novel model features were advanced treatment of superimposed ice and turbulent heat fluxes, coupling of snow and ice layers and snow modelled from precipitation. The simulated snow, snowÁice and ice thickness showed good agreement with observations for 2008Á2009. Modelled ice climatology was also reasonable, with 0.5 cm d1 growth in DecemberÁMarch and 2 cm d1 melting in April. Tuned heat flux from waterto ice was 0.5 W m 2. The diurnal weather cycle gave significant impact on ice thickness in spring. Ice climatology was highly sensitive to snow conditions. Surface temperature showed strong dependency on thickness of thin ice (B0.5 m), supporting the feasibility of thermal remote sensing and showing the importance of lake ice in numerical weather prediction. The lake ice season responded strongly to air temperature: a level increase by 1 or 58C decreased the mean length of the ice season by 13 or 78 d (from 152 d) and the thickness of ice by 6 or 22 cm (from 50 cm), respectively. Keywords: thermodynamic model, sensitivity tests, snow, ice, Lake Vanajavesi 1. Introduction cuts off airÁlake exchange of oxygen and reduces the production of dissolved oxygen by limiting the light pene- Freshwater lakes cover about 2% of the Earth’s land tration (Livingstone, 1993). Ice also prevents the exchange of surface. In Finland, there are ca. 190 000 lakes larger than momentum between the atmosphere and the lake water 500 m2, accounting for10% of the areaof the whole (Williams et al., 2004). Lakes affect the local weatherby country, and these lakes freeze every winter. Their ice modifying the air temperature, wind, humidity and precipi- sheets consist of congelation ice and snowÁice with snow tation in their surroundings (Ellis and Johnson, 2004; covernormallyon top. In medium-size and small lakes, the Rouse et al., 2008a, 2008b), and ice coverprotectsthe heat ice coveris usually stable, while in largelakes mechanical content of lakes by its insulation capacity and by damping displacements may take place, creating piles of ice blocks or turbulent mixing in the water body. Thus, the presence ridges. An ice cover stabilises the thermal structure of a (orabsence) of ice coverhas an impact on both regional lake. The ice bottom is at the freezing point, and there is a climate and weatherevents in the winter Áspring season. weak heat flux from the water to the ice. In spring, ice Understanding the processes and interactions of lake ice sheets gain heat mainly from solar radiation, and as the ice and atmosphere is essential for numerical weather predic- melts, all impurities contained in the ice are released into tion (NWP) (Brown and Duguay, 2010). As the spatial the waterorto the air. Lake ice plays an important role in the hydrological, resolution of NWP models becomes higher, presently biological, chemical and socio-economical regimes of boreal approaching one-kilometer scale, including lakes and lake lakes (Leppa¨ ranta, 2009). For example, a compact ice cover ice forforecastingand data assimilation in mesoscale, NWP systems have gained more attention (e.g. Eerola et al., 2010; Mironov et al., 2010; Salgado et al., 2010). *Corresponding author. Analytical, semi-empirical lake ice thickness models were email: [email protected] widely used until 1980s, and thereafter full numerical Tellus A 2012. # 2012 Y. Yang et al. This is an Open Access article distributed under the terms of the Creative Commons Attribution-Noncommercial 3.0 1 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. Citation: Tellus A 2012, 64, 17202, DOI: 10.3402/tellusa.v64i0.17202 (page numbernot forcitation purpose) 2 Y. YANG ET AL. models have also been employed (e.g. Liston and Hall, 2. Data 1995; Duguay et al., 2003). The basic principle of these models is to solve the heat conduction equation forthe 2.1. Observations during the ice season 2008Á2009 snow and ice layers in the vertical direction, as was first done in sea ice thermodynamic modelling by Maykut and Lake Vanajavesi (61.138N, 24.278E) is located in south- Untersteiner (1971). The Canadian Lake Ice Model is a western Finland (Fig. 1). It is a large, shallow and 2 snow and ice process model (Duguay et al., 2003) adapted eutrophic lake. The area is 113 km , and the mean and from a sea ice model of Flato and Brown (1996). maximum depths are 7 and 24 m, respectively. The lake has In most of the lake ice models snow thickness is a long eutrophication history with very poor water quality prescribed with a climatological growth rate or with a in the 1960s and 1970s. The ice season lasts 4Á6 months, on average from December to April, and the thickness of ice fixed snowfall scenario based on in situ data. In the boreal reaches its annual maximum value of 45Á60 cm in March. zone, however, the role of snow in the growth and melting The ice and snow thickness may show distinct spatial of lake ice is very important, and snow and ice form a variations due to snowdrift and heat flux from the lake coupled system. In analytical modelling, the influence of water. snow is usually considered by modifying the growth law In the winterof 2008 Á2009, a field programme was parameters (e.g. Leppa¨ ranta, 1993). Interactive snow and performed in Lake Vanajavesi. Hydrographical surveys snowÁice layers were included in the numerical quasi- were made regularly, an automatic ice station was set up steady ice thickness model of Leppa¨ ranta (1983). The first and optical investigations were performed (Jaatinen et al., detailed lake model study in Finland was presented by 2010; Lei et al., 2011). The average maximum snow and ice Leppa¨ ranta and Uusikivi (2002), based on a Baltic Sea thicknesses reached 15 and 41 cm, respectively, in March. model of Saloranta (2000), where snow metamorphosis and The freezing and breakup dates were 1 January and 27 snowÁice formation due to flooding were included (see also April, respectively. Shirasawa et al., 2005). In this study, a one-dimensional high-resolution thermo- dynamic snow and ice model (HIGHTSI) was applied for Lake Vanajavesi, located in southern Finland. This model contains congelation ice, snowÁice and snow layers with full heat conduction equation. Compared with earlier lake ice models, new features were superimposed ice formation, an advanced atmosphereÁice heat exchange treatment accounting forthe influence of stability of stratification to the turbulent fluxes and coupling of the snow and ice layers. Atmospheric forcing was derived from weather observations and climatology, which also drove the snow coverevolution. The simulation results were compared with measured ice and snow thickness. A case study was performed for the ice season 2008Á2009, forced by daily weather observations. Ice climatology was examined forthe 30-yearperiod1971 Á 2000; also the correlation between the observed monthly total precipitation and snow accumulation was investigated in order to understand the uncertainties of precipitation as model forcing for climatological simulation. A number of climate sensitivity simulations were carried out for the ice season. The objectives of the present work were to assess the applicability of the HIGHTSI model forlake snow and ice thermodynamics, to find out the most important factors affecting lake ice growth and melting and to evaluate the influence of climate variations on the lake ice season. Section 2 introduces data, the model is described in Fig. 1. Geographical location of Lake Vanajavesi (A); Section 3 and the results are presented and discussed surrounding observation sites Jokioinen meteorological in Section 4. Final conclusions follow in Section 5. observatory (B), Lake Kuivaja¨ rvi (C) and Lake Pa¨ a¨ ja¨ rvi (D). NUMERICAL MODELLING OF SNOW AND ICE THICKNESSES 3 Synoptic-scale weatherconditions and regionalclimate than in the winter2008 Á2009. The precipitation was over this area are represented by Jokioinen meteorological 232 mm, about the same that was observed in 2008Á2009, observatory (60.88N, 23.58E; WMO station 02863), located but during the ice growth season, the average snow some 50 km southwest of Lake Vanajavesi (Fig. 1). The accumulation was 145 mm WEQ or58% morethan in weatherforcingdata forthe lake ice model consist of wind 2008Á2009. speed (Va), air temperature (Ta), relative humidity (Rh), The average freezing and ice breakup dates of Lake cloudiness (CN) and precipitation (Prec), collected at Vanajavesi were 30 November and 30 April, respectively, in 3-hour time intervals. For the present modelling work, we 1971Á2000 (Korhonen, 2005). A large variation occurred in interpolated data to 1-hour intervals. The snow thickness these dates in individual years. The earliest freezing date on land is also available from Jokioinen. was 28 Octoberand the latest one was 28 January.Forthe The winter2008 Á2009 was mild (Fig. 2). The early winter ice breakup, the earliest and latest dates were recorded as was warm, with average air temperature of 0.68Cin 1 April and 27 May, respectively.