A new national daily snowfall record was measured in Finland on
8 January 2016 when it snowed 73 cm (31 mm as liquid water) in less than a
day in Merikarvia on the western coast of Finland. The area of the most
intense snowfall was very small, which is common in convective precipitation.
In this work we used hourly weather radar images to identify the sea-effect
snowfall case and to qualitatively estimate the performance of HARMONIE, a
non-hydrostatic convection-permitting weather prediction model, in simulating
the spatial and temporal evolution of the snowbands. The model simulation,
including data assimilation, was run at 2.5 km horizontal resolution and
65 levels in vertical. HARMONIE was found to capture the overall sea-effect
snowfall situation quite well, as both the timing and the location of the
most intense snowstorm were properly simulated. Based on our preliminary
analysis, the snowband case was triggered by atmospheric instability above
the mostly ice-free sea and a low-level convergence zone almost perpendicular
to the coastline. The simulated convective available potential energy (CAPE)
reached a value of 87 J kg
The year 2016 started with imposing a new national snowfall record in
Finland, located in the Baltic Sea region in northern Europe. The snow depth
increased by 73 cm in Merikarvia (61.85
The Merikarvia snowfall case was caused by a lower tropospheric mesoscale phenomenon, lake-effect (here sea-effect) convection. If a cold air mass flows over a warm ice-free water surface, the water area acts as a source of heat and moisture. This produces an unstable boundary layer over the water body. As a result, the turbulent heat and moisture fluxes from the water surface are large and generate shallow convection that induces small and intensive convective precipitation which can drift to the coast as snowbands (Markowski and Richardson, 2010; Mazon et al., 2015).
Surface heating of cold-air advection over relatively warm water is the most
important factor affecting the occurrence of snowbands (Mazon et al., 2015,
and references therein). As a rule of thumb, the temperature difference
between the 850 hPa temperature and the water surface should be at least
13
Temperature (in
To form and reach the coast, intense sea-effect snowbands require a
sufficient fetch over open water and a suitable wind direction and speed.
According to Markowski and Richardson (2010), the former is usually over
75 km and the latter is 7–15 m s
Predicting and understanding of sea-effect snowfall is important as these events can cause serious problems for the fluency and safety of coastal transportation (road/train/aircraft) and electricity transmission lines. Several factors affect the formation of sea-effect snowfall. The aims of this case study are to identify the main triggering factors that induced these snowbands and to qualitatively verify weather prediction simulations with the aid of hourly radar images. In this paper only preliminary results are presented as the study is still ongoing. In the future, more snowband cases will be selected and evaluated.
Hourly radar reflectivity images from the weather radar network of the Finnish Meteorological Institute (FMI, Saltikoff et al., 2010; Gregow et al., 2017) were used to preliminary identify the sea-effect snow case on 8 January 2016. Moreover, radar-based hourly accumulated precipitation images were qualitatively compared to simulations of precipitation with a weather prediction model (HARMONIE; see below). Observational information about the sea ice cover extent and sea surface temperatures over the Gulf of Bothnia were provided by ice charts of FMI (Fig. 1). In order to explore the synoptic situation during the Merikarvia event, air pressure and temperature fields from the HARMONIE simulations and from the operational ECMWF medium range forecasts were used together with analyses made by FMI meteorologists (fronts, temperature, precipitating areas, clouds, weather symbols, and wind barbs).
HARMONIE (ALADIN-HIRLAM system, version 38h1.2) is a non-hydrostatic convection-permitting limited area Numerical Weather Prediction (NWP) model (Bénard et al., 2010; Brousseau et al., 2011). It has been developed in collaboration by ALADIN and HIRLAM NWP model consortia including multiple European weather centres. Physical parametrizations of the HARMONIE model include radiation, surface, shallow convection, turbulence and microphysics schemes (Seity et al., 2011). In HARMONIE, the radiation parametrizations of ECMWF are used (Iacono et al., 2008; Mlawer et al., 1997; Morcrette et al., 2008). The surface scheme SURFEX simulates fluxes of heat, moisture and momentum between the surface and the atmosphere (Masson et al., 2013). Shallow convection parametrization uses scheme by Pergaud et al. (2009). The turbulence parametrization scheme, which is based on the TKE equation, was developed by Cuxart et al. (2000). For microphysics, three-class ice parametrization scheme ICE3 (Pinty and Jabouille, 1998) is applied. The data assimilation system used in HARMONIE for upper air observations is 3D-Var. In surface data assimilation the Optimal Interpolation (OI) method is used.
Ice chart on 8 January 2016. Water temperature isotherms are shown with black solid lines along with the sea surface mean temperatures (1971–2000) encircled. Purple/red shades indicate different ice formations (new ice, grey ice, very close ice). Ice thickness is shown in squares (cm). Merikarvia is marked with a green dot, Pori with a blue dot, Rauma with a black dot, and Helsinki with an orange dot.
In this work the AROME configuration of the HARMONIE system was used. We run the model for the Merikarvia case following the operational practise at FMI, i.e. 2.5 km horizontal resolution and 65 model levels in vertical, with the model top at 10 hPa, and a domain that covers Finland, Scandinavia and the Baltic countries. The lateral boundary conditions were obtained from the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS) global NWP model. The data assimilation system used in HARMONIE for upper air observations is 3D-Var. In surface data assimilation the Optimal Interpolation (OI) method is used. The assimilated observations included surface synoptic observations both on land and sea (SYNOP and SHIP), aircraft reports (AMDAR, AIREP, ACARS), buoy observations (BUOY), radiosondes (TEMP) and wind profiler (PILOT) observations.
It would be possible to run the HARMONIE system also as a climate model. The HARMONIE Climate system (HCLIM, Lindstedt et al., 2015) is an active development branch in the ALADIN-HIRLAM community.
During the turn of the year from 2015 to 2016, a high-pressure area started to approach the Fennoscandia from the east, which enabled a flow of very cold and dry air masses. High sea-level pressure values (1040–1050 hPa) were observed over Finland a week before the Merikarvia event but those had decreased to 1010–1015 hPa over the Bothnian Sea, a part of the Gulf of Bothnia in front of Merikarvia, until 8 January. A ridge of high 850 hPa geopotential heights extended from southern Finland to central Sweden, while there were centres of low geopotential heights in the north-eastern and southwestern corners of the domain (Fig. 1). At the level of 500 hPa westerly winds prevailed.
One-hour accumulated precipitation (mm h
Vertical cross section from 1000 to 700 hPa of the equivalent
potential temperature (K) along 62
A thermodynamic diagram taken over the sea at 62
Vertical cross section from 1000 to 500 hPa of cloud water and
ice (g m
The autumn and early winter had been exceptionally warm in 2015. The mean
autumn temperatures in western coast of Finland were 2–2.5
On 8 January the snowband event had several hours to form. According to the radar images (FMI, 2016), snowfall started in Merikarvia at about 09:30 UTC and finally ended after about 18 h when the westerly flow turned to southerly. Small intense snowbands were detected by radars still on 10 January over the Gulf of Bothnia. On 10 January, a low pressure system started to approach from the southwest, causing southeasterly flow with warm advection, which led to less favorable conditions for convection to form.
HARMONIE could capture the snowband event quite well when compared to the radar observations (Fig. 3). According to the simulation, separate cells of convective precipitation were formed over the Gulf of Bothnia during the night between 7 and 8 January. The convective cells evolved to one narrow snowband which hit the Finnish coastline perpendicularly on 8 January. The hourly accumulated precipitation intensity is comparable to radar images and the timing and shape of the snowband was correctly simulated (Fig. 3). The simulated total 24 h accumulation (30 mm) corresponded very well the observed (31 mm) value at Merikarvia. The location of the snowband was simulated well although slightly north to the observed snowband.
The simulated 2 m temperatures were low: the minimum over the western coast
of Finland was around
Over the sea near Merikarvia, convective available potential energy (CAPE;
calculated as pseudoadiabatic with ice) and the lifted condensation level
(LCL, in hPa) increased in time from early morning values of 10 J kg
Resulting from colliding winds from the south and north, there was a zone of low-level convergence extending from the sea to Merikarvia (13:00 and 21:00 UTC in Figs. 3 and 6). The convergence zone was almost perpendicular to the Finnish coast. The air mass south of it was distinctly drier (relative humidity (RH) of 60–70 % at 2 m height and 20–70 % at 850 hPa) than the air mass north of the zone (RH of 80–100 % at 2 m and 80–100 % at 850 hPa). At the height of 2 m the strongest RH gradient was found in the same area as the snowband and at the 850 hPa level the RH gradient was tilted towards south. The snowband was simulated to form along the convergence zone and it drifted towards the coastline as the southerlies and southwesterlies were dominating during the first half of the day. Towards the evening the southerlies turned to southeasterlies, the northerlies intensified (Figs. 3–6), and the area of maximum wind velocities moved from the Finnish coastline towards the Swedish coast. Directional wind shear with height varied along the day over the Gulf of Bothnia. During the most intense snowfall, the wind direction was almost the same at 850 and 1000 hPa i.e. from the southwest (Fig. 5).
In Merikarvia the maximum fetch was
The offshore snowfall over the Bothnian Bay during the morning (05:00 UTC in Fig. 3) was found to be related to scattered zones of low-level convergence. Besides, in addition to the very intense snowfall on the western coast, some weaker sea-effect snowfall areas were also observed on the southwestern and southern coastline of Finland. HARMONIE simulated some low-level convergence there but no snowfall.
Weather radars observed several snowbands over the Gulf of Finland near the southern coast for almost two weeks before the Merikarvia event. However, at that time the wind direction had not been favourable for the snowbands to hit the coast. With a more south-easterly flow, the new snow record could have taken place in the Helsinki metropolitan area on the southern coast of Finland and severe weather might have lasted for several days there. Snowfalls had also been detected over the Gulf of Bothnia near the western coast of Finland already three days before the Merikarvia event, with almost the same intensity, but those snowfalls remained over the sea.
Fortunately, the impacts of the record snowfall in Merikarvia remained
small. The intense snowfall area was very limited (
It is intriguing that the two most intense daily snowfalls have happened almost at the same place; the distance between Rauma and Merikarvia is only 100 km (Fig. 2). So, what are the factors favouring intense sea-effect snowfalls in this area? The Gulf of Bothnia is 80–240 km wide and has an average depth of 60 m. In addition, the topography around the gulf is flat. Usually the westerly and southwesterly airflow is too warm for snowbands to form in this area but when the air mass is sufficiently cold and the gulf is still ice-free, the fetch over the open Gulf of Bothnia can be much longer than on the narrower Gulf of Finland.
The roles of colliding land breezes from the coasts of Finland and Sweden or Finland and Estonia; orographic lifting; and a concave shape of a coastline in promoting offshore convergence and formation of snowbands have recently studied by Savijärvi (2012, 2015) and Mazon et al. (2015). Based on our preliminary analysis, the Merikarvia snowfall case was triggered by atmospheric instability above the mostly ice-free sea and a low-level convergence zone. The zone was located almost perpendicular to the Finnish coastline and to the east of a concave part of the Swedish coast. Whether synoptic, mesoscale or smaller-scale factors, or all them, contributed to the formation of the convergence zone remains as a topic beyond the scope of this paper.
Finally it may be noted that the simulated CAPE of 87 J kg
The weather conditions were favourable for snowband to form over the whole Bothnian Sea (a part of the Gulf of Bothnia) as the sea was relatively warm and ice free, and arctic air mass prevailed over most of Fennoscandia. Due to colliding winds over the Bothnian Sea a very localised extreme of snowfall was produced and drifted towards Merikarvia as the maximum winds blew from southwest. The narrow snowband perpendicular to coast line caused record snowfalls to the small area of Merikarvia in the western coast of Finland. The non-hydrostatic convection-permitting weather prediction model HARMONIE captured the overall sea-effect snowfall situation quite well. Both the timing and the location of the most intense snowstorm was correctly simulated. In addition the accumulated precipitation was simulated correctly. Nevertheless, the simulated weaker snowfall did not spread as broadly along the coastline of Finland as was observed.
The preliminary results discussed here suggest that HARMONIE is able to simulate sea-effect snowfall sufficiently accurately and thus can be used to examine prominent weather situations to form snowbands. According to our plan, more sea-effect snowfall cases will be selected and evaluated.
The modelled data (in grib format) and the hourly radar reflectivity images used in this study can be provided upon request to the corresponding author. The source code of the HARMONIE model is not publicly available.
The authors declare that they have no conflict of interest.
This article is part of the special issue “16th EMS Annual Meeting & 11th European Conference on Applied Climatology (ECAC)”. It is a result of the 16th EMS Annual Meeting & 11th European Conference on Applied Climatology (ECAC), Trieste, Italy, 12–16 September 2016.
This work was supported financially by VYR (State Nuclear Waste Management
Fund) through SAFIR2018 (The Finnish Research Programme on Nuclear Power
Plant Safety 2015–2018,