CC BY
15Reconstruction of the Ice Thickness Seasonal Evolution in the Northeastern Sea of Azov Using Different Arrays of Meteorological Data
D. D. Zavyalov
Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol, Russian Federation
e-mail: zavyalov.dd@mhi-ras.ru
Purpose. The aim of the paper was to compare the results of numerical experiments on reconstructing seasonal thermal evolution of the sea ice thickness with the data of in situ observations of the ice state in the northeastern part of the Taganrog Bay.
Methods and Results. Characteristics of the ice state in the northeastern part of the Taganrog Bay were studied using the previously developed thermodynamical model of sea ice. The data of the European Center for Medium-Range Weather Forecasts ERA-Interim, regional prognostic model SKIRON and the array of daily eight (with 3-hour intervals) observations of the basic meteorological parameters (All-Russian Research Institute of Hydrometeorological Information - World Data Center (RIHMI - WDC)) obtained at the meteorological station Taganrog, were used in the numerical experiments as the atmospheric forcing. The modeling results were compared with the in situ data for the winter seasons in 2007/2008-2010/2011. It is shown that the characteristics of the snow-ice cover resulted from application of various meteorological data as the external forcing, can be significantly different.
Conclusions. The highest similarity between the modeled ice thickness seasonal variation and the one reconstructed using the observations data was obtained at applying the RIHMI - WDC data array. In this case, both thickness and the basic stages of the snow-ice cover evolution in the Taganrog Bay were adequately reproduced in the model. As compared to the in situ data, the results of the models based on the SKIRON and ERA-Interim data were mainly overestimating and underestimating, respectively. It is related, to a great extent, to determination of the precipitation amount, the prognostic values of which in ERA-Interim are higher than those in SKIRON. However, even in the calculations taking no account of atmospheric precipitation or in those for the ice seasons when the atmospheric precipitation is very insignificant, the SKIRON based model provides the higher values of ice thickness than the values resulted from the ERA-Interim based model. Analysis of the modeling results shows that adequate reconstruction of the ice state characteristics in the Azov Sea requires preliminary setting of the thermodynamic model depending on the chosen data array used as the atmospheric forcing.
Keywords: sea ice, thermodynamics, ice thickness, atmospheric forcing, Sea of Azov, Taganrog Bay.
Acknowledgements: the investigation was carried out within the framework of the state task on theme No. 0827-2018-0003 "Fundamental studies of oceanologic processes governing state and evolution of marine environment affected by natural and anthropogenic factors, based on the observational and modeling methods".
For citation: Zavyalov, D.D., 2019. Reconstruction of the Ice Thickness Seasonal Evolution in the Northeastern Sea of Azov Using Different Arrays of Meteorological Data. Physical Oceanography, [e-journal] 26(3), pp. 247-259. doi:10.22449/1573-160X-2019-3-247-259
DOI: 10.22449/1573-160X-2019-3-247-259
© 2019, D. D. Zavyalov © 2019, Physical Oceanography
Introduction
In winter period, the Sea of Azov partially freezes, in severe winters - completely, therefore sea ice is an important component of the sea hydrological regime. The Sea of Azov ice cover is characterized by high variability
both in a single season and in a long-term period and depends on the prevailing hydrometeorological conditions. Contrasting weather and ice conditions are formed under effect of high pressure area located in the north and north-east of the Eurasian continent, and cyclonic invasions from the south-west and west [1]. Frequent eastern and northeastern winds bring continental polar and arctic air masses, cause severe frosts and intensive ice formation, and cyclonic activity, on the contrary, leads to the inflow of warm air masses and a rapid decrease in the thickness of snow-ice cover [2].
The sea ice evolution is determined by the thermal and dynamic regimes of the atmosphere and the aquatic environment [3], while the ice cover itself significantly affects the heat flux between these media. The limitations of physically based modeling of the thermal evolution of the snow-ice cover thickness are due to the availability of initial information. In order to calculate the seasonal variation of sea-ice thickness, reliable data on temperature, pressure and air humidity, wind velocity, cloudiness, the amount of precipitation and its phase are required. The quality of this data also determines the calculation quality of the snow-ice cover thickness.
A comparative analysis of calculation results of the Arctic Ocean ice and hydrological characteristics using three different data sets as atmospheric effects is presented in [4]. Comparison of NCEP/NCAR reanalysis data used as atmospheric forcing in a one-dimensional thermodynamic model of sea ice with the observations obtained at drifting stations in the Central Arctic is presented in [5].
The processes of heat and moisture transport in the Sea of Azov snow-ice system significantly differ from the processes occurring in the snow-ice cover of the northern seas. The complexity of the observations and the small number of direct measurements of the Sea of Azov snow-ice cover thickness significantly complicate the study of the winter thermodynamics of this basin. The integration of publicly available data from weather forecast models with a thermodynamic model of seasonal evolution of sea-ice thickness can be considered as a possible element of the forecast of the Sea of Azov ice conditions.
This paper presents the results of modeling the formation and melting processes of snow and ice cover in the Taganrog Bay for the winter seasons from 2007/2008 to 2010/2011. The winter of 2007/2008 in the Sea of Azov region [6] belonged to moderate winters, but its feature was contrasting weather and ice conditions, characterized by both intensive ice formation with short-term low temperatures and significant thaws. Winters of 2008/2009 and 2009/2010 were mild, with cyclonic activity predominating in the atmospheric processes. The effect of anticyclones was insignificant. The minimum air temperatures during the winter period of 2008/2009 were observed in January 2009. During the winter period of 2009/2010, two cold waves, which came from the north in mid-December 2009 and the end of January 2010, were observed. Between them, the sea near the Taganrog weather station was completely clear from the ice. In 2010/2011 ice season, two periods, differing in the general temperature background and ice cover, were clearly visible. The first - from November to December - is warm. At this time, an intensification of cyclonic activity was observed. The second - from January to March - was cold and was characterized by the invasion of the arctic cold masses into the rear of the leaving Atlantic cyclones. In the area under consideration, this winter was classified as mild one [7].
T„°c P,, kPa |KJ,m/s tp, % N, cloud points Pr, kg/m2
Dec Jan Feb Mar Dec Jan Feb Mar Dec Jan Feb Mar Dec Jan Feb Mar Dec Jan Feb Mar Dec Jan Feb Mar
- • Skiron - + Era-Interim - ♦ RIHMI
F i g. 1. Average decadal values (solid lines) of the air temperature on height 2 m Ta, the atmospheric pressure on the basin surface Pa, the wind speed on 10 m height |Va|, relative air humidity 9, general cloudiness number, as well as total precipitation for a month Pr over the northeastern part of the Taganrog Bay (47° 11'N, 38° 54'E) in December - March, 2007/2008, 2008/2009, 2009/2010 and 2010/2011. Symbols on the graphs denote minimal and maximal monthly values of the meteorological parameters
Three meteorological databases were used as information support: SKIRON prognostic atmospheric model [8] of the University of Athens (Greece) with 0.1° x 0.1° spatial resolution and 2 h time step; ERA-Interim reanalysis [9] of the European Center for Medium-Range Weather Forecasts (ECMWF) with 0.125° x 0.125° resolution and 6 h time step, as well as RIHMI-WDC* array of urgent observations over the main meteorological parameters at the Taganrog weather station (WMO index 34720) with 3 h interval.
In Fig. 1 solid lines denote ten-day mean values of air temperature at 2 m height, atmospheric pressure near the basin surface, wind velocity at 10 m height, relative humidity and general cloudiness number. The bar charts denote total precipitation above the northeastern part of the Taganrog Bay for the month from December to March 2007/2008 - 2010/2011. It can be seen from the figure that the meteorological parameters given in these arrays have a similar nature of temporal variability, and the difference in ten-day values of air temperature, pressure and wind velocity does not exceed 2 °C, 1.5 kPa and 3 m/s, respectively. The greatest differences are manifested in the forecast of precipitation amount, since precipitation is one of the most variable in time and space meteorological phenomena. Uncertainties that arise when modeling the sea ice thickness thermodynamic evolution are largely related both to the estimation of general precipitation amount getting on the sea surface and the lack of reliable information about their structure.
* RIHMI-WDS, 2018. Hydrometeorological data. [online] Available at: http://meteo.ru [Accessed: 04 September 2018] (in Russian).
Description of the model
Evolution of the Sea of Azov snow-ice cover, taking into account regional features of its formation, is described by a locally one-dimensional thermodynamic model [10, 11]. Heat distribution in the snow and ice layers is determined by thermal conductivity equations
(Pc)i,s
ÔThS(z,t)_ d(k ST{-(z,t)^ QIl(z,t)
i,s
dz
dz
dt dz
with the boundary conditions at the snow-ice cover upper
z = 0
z G
[0, h + h- ]
dz
on the surface of snow-ice interface
, dTs dTi ks—^ = ki—L, dz dz
and at the ice cover lower boundary
STis _
T = T
1 s 1 i,
ST-
-kl -T -Fb(Tf,Tw), Ti = Tf,
Sz
z = hs (t)
z = hi (t).
At the moving interfaces of different matter phases the law of energy conservation during phase transitions is fulfilled:
dh[ s d Ti s
-Pi,sZfi,s -r- = Ft (Tmi,s ) + ki,s , z = 0. dt dz
Sh ST
-Pi Lfi -r-- Fb (Tf ,Tw )- ki -r1,
St Sz
z = hi (t).
If the snow-ice cover is absent, the heating or cooling of mixed quasi-homogeneous water layer takes place:
STw
St
Pwcwhw - Fb - Ft - vPs Lfs,
z G
[o, hw ].
Here t is time; z is a vertical coordinate directed downward from the snow-ice cover upper surface (z = 0); p is a density; h is a thickness; T is a temperature; S a salinity; c is a heat capacity; k is a thermal conductivity; Lf is a heat of fusion; Ii is a solar radiation penetrating the ice; Tsfc, Tmi,s ,Tf are the temperatures on the upper surface of the snow-ice cover, the melting point of ice/snow and the water freezing temperature, respectively; v is the rate of solid precipitation; Ft is a heat flux from the atmosphere through the upper boundary of the snow-ice cover; Fb = cwpwCtb (Tw - Tf) is a heat flux from the water to the lower ice boundary,
Ctb = 10-3 m/s is a turbulent exchange coefficient. Hereinafter, i, s, w, a indices refer to the parameters of ice, snow, water and the atmosphere, respectively.
The rate of snow-ice cover thermal evolution depends on the difference between the heat flux from the atmosphere through its upper boundary and the heat
flux from the water to its lower boundary. At the upper boundary, the heat flux Ft consists of turbulent fluxes of sensible (Fs) and latent (Fl) heat, which are determined by the integral aerodynamic formulas
Fs = CpaPaStVa(Tsfc " Ta) ; F = LPaDaVa(?0(Tsfc) - 4a(Ta))
long-wave (R)
P = 4AoTsfcTa3 - AnTa4(3.765+ 0.22N3),
short-wave (F) radiation fluxes
F=Fo(1-CV)(i-a)(l-/0) ,
as well as heat fluxes (Fm) related to the cooling processes and the subsequent possible crystallization of liquid precipitation:
Fn=CwPr (Ta - 273.15)+PrLi.
Here e is the water vapor pressure; Pa is the atmospheric air pressure; Va is wind velocity; Prr is the amount of precipitation in liquid phase; St = Da = = 1.7-10"3; qa = (0.622f /Pa>10aira/(b+Ta), % = (0.622e0 / Pa ^O^^i +f is
specific air humidity at 2 m height and at the snow-ice cover upper boundary, respectively; ax = 9.5; bx = 265.5 K; e0 = 611 hPa is the pressure of saturated water vapor at 0 °C; cpa = 103 J/(kgK) is the heat capacity of air at constant pressure; f is relative humidity; L is the specific heat of sublimation; c is Stefan-Boltzmann constant; X, a is an emissivity and albedo of the underlying surface; N is the general cloudiness number; i0 is a coefficient that determining which part of the short-wave radiation penetrates deep into the snow-ice layer and is evenly distributed throughout the entire thickness;
F0 = S cos2 z0 ((cosz@ + 2.7)e -10-5 + i.085cos % + O.i)-1 is the incoming short-wave solar radiation for a cloudless sky; S is the solar constant; z0 is solar zenith angle. Assuming that the vertical profiles of ice and snow temperature can be described by a linear function, and considering the heat flux through ice and snow to be the same (according to the study of A.P. Makshtas * and [12, 13]), the thermal conductivity equations can be solved analytically. The nonlinear equation for determining the temperature of the snow-ice cover upper surface, obtained from the heat balance equation assuming the continuity of the heat flux at the snow-ice interface, is solved numerically. The parameters used in the model for physical, thermal and optical characteristics of sea water and ice are given in [14].
The basis of the model of snow cover, accumulating on the sea ice surface, is the schematization of the processes making the greatest contribution to the formation of physical and thermal snow characteristics. Due to the intermittent nature of precipitation accumulation, wind effect and changes in air temperature, snow cover is composed of separate layers that differ from each other at least in thickness, density and water content. Considering the regional features of precipitation in the winter period on the Sea of Azov coast, we take the layers of fresh and existing snow as the main structural model units.
* Makshtas, A.P., Thermal Balance of Arctic Ice in the Winter: Gidrometeoizdat, 1984. 67 p. PHYSICAL OCEANOGRAPHY VOL. 26 ISS. 3 (2019) 251
A layer of fresh snow is formed as a result of snow accumulation on the surface of the snow-ice cover during the deposition of solid precipitation. The density of freshly fallen snow ps0 is determined by COSMO algorithm [15]. This parametrization is used to calculate the characteristics of fresh snow at a short period of snow accumulation (up to 12 hours) and in the case of small amount of precipitation in the form of snow [16]. Due to the relatively small amount of solid precipitation observed over the water area of the Sea of Azov, as well as extreme instability of snow cover, the process of elastic deformation of fallen snow can be neglected. Wind compaction of snow was parametrized, based on the assumption [12] that its density rises by 20 kg/m3 at wind velocity increasing for every 1 m/s, as ps = max (ps0; 20 Va) (kg/m3). Due to snow melting or rain falling, the water, which seeps into the snow layer and leads to an increase in its density, is formed. The maximum amount of water a snow layer can contain is determined by its water retention capacity 0max *. If the liquid water amount 0w in the snow layer exceeds 0max, then its surplus moves to the bottom layer or forms water flow to the snow-ice interface, where it crystallizes if Ts < Tm. The density of each layer is calculated based on the amount of water contained in it in the liquid and solid phase. Thermal conductivity of a snow column (a set of snow layers) is determined by the Osokin's formula [17]. The snow surface albedo parameterization is taken from ECHAM5 atmosphere model [18].
Ice cover buoyancy is small and its overload occurs when the snow cover height reaches approximately 40% of the ice thickness, so the process of turning snow into ice is possible when the snow - ice section line falls below the water level. The thickness of the flooded part of the snow is calculated from the floating condition of the bodies. According to the change in the snow column height, its mass decreases and the ice mass increases.
Results of calculations
On the basis of the constructed thermodynamic model, numerical experiments for assessing the effect of external forcing choice on the reconstruction of ice thickness seasonal evolution in the north-eastern part of the Sea of Azov were carried out. The modeling results were compared with each other and with the data on the sea ice thickness taken from ice maps published by the Unified State System of Information on the Situation in the World Ocean ** (ESIMO).
In order to assess the adequacy of the forecast of the ice thickness seasonal variation, the following criteria were considered: - model error Et = hin situ - hi;
E 2
- root-mean-square deviation of hi from hin situ - o = A I—— ;
i l
- determination coefficient R2 (dispersion portion hi explained by
■ 2 _ "t
the model), R2 = 1 -,
IE,2
I(h - h )2
* Kuzmin, P.P., 1957. Fizicheskie Svoystva Snezhnogo Pokrova [Physical Properties of Snow Cover]. Leningrad: Gidrometeoizdat, 179 p. (in Russian).
** NODS, 2018. The Unified State System of Information on the World Ocean. [online] Available at: http://esimo.ru/portal/ [Accessed: 04 September 2018] (in Russian).
- Theil's inequality coefficient U (indicates the degree of time series similarity: the closer it is to zero, the closer are the compared series)
E e 2
u=_rhE_
i?^ w
- correlation coefficient K =
situ " hn situ Xhi - hi )
,E(h n situ - hn situ )\ E(h - h )2
\l i Mi
where l is a number of steps in which the measured ice thickness, corresponding to the computational time step, differed from zero.
The modeling results of ice thickness thermodynamic evolution in ice seasons from 2007/2008 to 2010/2011 are represented in Fig. 2. Since, as noted earlier, the most noticeable differences between the meteorological data of the considered arrays are manifested in the amount of precipitation, for a comparative analysis of atmospheric forcing effect on the ice thickness seasonal variation, a series of calculations were carried out both with regard to precipitation and without it [1].
Average model errors Et calculated for each of four ice seasons are shown in
Fig. 3; the season-average values of the measured hnsitu and computational hi ice
thickness, root-mean-square deviations o, determination coefficients R2, Theil's inequality coefficient U and correlation coefficients K are listed in the table. It can be seen that for the considered periods the forecast of ice thickness seasonal variation, performed without taking into account the precipitation, is overestimated. This conclusion applies to all three atmospheric forcing. The largest negative values of the model errors were obtained for the SKIRON data, and the smallest - for the RIHMI-WDC data (Fig. 3, a). However, despite the consistent ice thickness hi reassessment by the models constructed with no regard to precipitation, they adequately reconstruct (except for SKIRON, 2007/2008) the dates of the maximum ice thickness formation (see. Fig. 2). In addition, sufficiently high correlation coefficients indicate a similar nature of temporal dependences of the measured hin situ (t) and computational hi(t) ice thickness values.
It should be noted that although the ten-day values of meteorological parameters (see Fig. 1) in the given arrays are quite close to each other (except for precipitation, which were not taken into account in this series of calculations), the differences in the maximum ice thickness forecast for the same season can reach 40-70 % of the average ice thickness per season hnsitu, depending on the atmospheric forcing type.
F i g. 2. Results of modeling the ice thickness thermodynamic evolution in the 2007-2011 in the ice seasons obtained due to application of the meteorological data arrays SKIRON, ERA-Interim and RIHMI - WDC as the external forcing with the regard for precipitation (solid lines) and with no regard for precipitation (hatch lines). Grey circles correspond to the sea ice thickness data from the ESIMO ice charts, red ones - to the maximum ice thicknesses measured at the meteorological station Taganrog [1]
Ef, cm
F i g. 3. Season-average model errors e resulted from determining the ice thickness seasonal variation using meteorological data from the SKIRON, ERA-Interim and RIHMI - WDC arrays with no regard to atmospheric precipitation (a) and with regard to atmospheric precipitation (b)
In order to carry out numerical experiments taking into account the snow accumulation on the sea ice surface, it is necessary to choose a criterion by which the phase of precipitation is determined. As one of these criteria, air temperature near the sea surface can be used. The criterion for dividing precipitation into liquid and solid can be the threshold air temperature [19], below which all precipitation is classified as snow, and higher as rain. In some studies [20], the empirical dependences of rain and snow percentage on the surface temperature are given and it is assumed that there is a range of temperature values at which mixed precipitation is observed. Such dependencies are not universal and have pronounced regional in nature. The condition for dividing the total amount of prognostic precipitation into liquid and solid is to a certain extent a setup one. In this work, it was assumed that when the air temperature is below -0.5 °C, all precipitation is in the solid phase, and at air temperature above 0.3 °C, it is only in the liquid one. Within the air temperature range -0.5 °C < Ta < 0.3 °C the snow content in precipitation was determined as follows: within -0.5 ... -0.2 °C temperature interval the snow content percentage was measured linearly from 100 to 95 %, within -0.2 ... 0 °C - from 95 to 60 % and within 0 ... 0.3 °C - from 60 to 0 %.
Calculations showed that the consideration of precipitation reduces Et model error (Fig. 3, b). However, in most of the considered cases, when using SKIRON array the model remains noticeably overestimating, and when using ERA-Interim data - underestimating. From three considered arrays, ERA-Interim gives the most, and SKIRON - the least amount of predictive precipitation in both the seasonal and inter-annual cycle. It should also be noted that when modeling the thermal evolution of ice thickness for the conditions of winter 2009/2010 (which was most abundant in precipitation: their total monthly amount exceeded climatic norms by 1.5-2 times), all three models turned out to be somewhat underestimating. The most adequate values of ice thickness seasonal variation for the selected
parameterizations in the thermodynamic model were obtained using RIHMI -WDC data (Fig. 2, 3, table). In this series of calculations, the duration of the ice period, the observation time of the maximum ice thickness and its value, the error in the determination of which did not exceed 3 cm, were rather reliably reconstructed.
T a b l e
Comparison of the ice thickness model calculations with the data of in situ observations
Season h : , cm rin situ' Taking account of precipitatio n Forcing hi, cm a, cm R2 U K
SKIRON 40.4 18.1 - 0.26 0.93
00 o o Cj No ERA-Interim 29.4 7.1 0.66 0.12 0.97
23.1 RIHMI 25.5 3.1 0.93 0.06 0.98
o o SKIRON 29.8 8.6 0.50 0.15 0.89
<N Yes ERA-Interim 10.0 15.1 - 0.40 0.86
RIHMI 22.7 1.9 0.98 0.04 0.98
SKIRON 34.9 20.8 - 0.40 0.74
a\ o o JN No ERA-Interim 32.3 17.4 - 0.35 0.88
15.3 RIHMI 26.8 12.1 - 0.30 0.89
oo o o SKIRON 25.4 11.9 - 0.27 0.62
<N Yes ERA-Interim 14.9 7.1 0.04 0.22 0.26
RIHMI 17.0 4.6 0.59 0.13 0.83
SKIRON 20.4 7.9 0.37 0.19 0.95
o O JN No ERA-Interim 17.8 6.1 0.62 0.16 0.92
14.4 RIHMI 17.6 4.7 0.78 0.12 0.96
0^ o o SKIRON 11.5 7.7 0.40 0.25 0.71
<N Yes ERA-Interim 12.9 8.9 0.19 0.28 0.51
RIHMI 14.0 4.7 0.77 0.14 0.88
SKIRON 23.1 13.2 - 0.35 0.89
No ERA-Interim 20.4 9.9 - 0.29 0.92
O JN 10.9 RIHMI 15.4 5.4 0.41 0.18 0.94
O O SKIRON 17.6 7.7 - 0.24 0.87
<N Yes ERA-Interim 4.7 7.3 - 0.38 0.87
RIHMI 9.4 3.9 0.68 0.15 0.90
Conclusion
The results of modeling of the ice thickness thermodynamic evolution in 2007-2011 ice seasons, obtained using different meteorological data as external forcing, showed that at the same parameterization of physical processes in the snow-ice cover thermodynamic model, the computational snow-ice thickness values can vary significantly. The greatest similarity of the reconstructed seasonal variations in ice thickness in the Taganrog Bay northeastern part with the data on the sea ice thickness taken from ESIMO ice maps was obtained using the observations taken 8 times per day of the main meteorological parameters at the Taganrog weather station (RIHMI - WDC) in the array model. The model constructed on the basis of SKIRON data is usually overestimated, and on the ERA-Interim data - is an underestimating one. This is largely due to the determination of precipitation amount, the prognostic values of which in ERA-Interim are greater than in SKIRON. However, even in calculations without taking into account precipitation (or in ice seasons with a small amount of them), the model based on the SKIRON data gives higher ice thickness values than the one based on ERA-Interim data. The analysis of the simulation results showed that for adequate reconstruction of ice regime characteristics in the Sea of Azov, it is necessary to carry out a preliminary adjustment of the thermodynamic model depending on the selected data set used as atmospheric forcing.
REFERENCES
1. D'yakov, N.N., Timoshenko, T.Yu., Belogudov, A.A. and Gorbach, S.B., 2015. Atlas L'dov Chernogo i Azovskogo morey [Ice Atlas of the Black and Azov Seas]. Sevastopol: ECOSI-Gidrofizika, 219 p. (in Russian).
2. Dumanskaya, I.O. and Fedorenko, A.V., 2008. Analysis of the Connection of Ice Cover Parameters of the Non-Arctic Seas in the European Part of Russia with Global Atmospheric Processes. Russian Meteorology and Hydrology, [e-journal] 33(12), pp. 809-818. https://doi.org/10.3103/S106837390812008X
3. Bukatov, A.E., Zav'yalov, D.D. and Solomaha, T.A., 2016. Analiz Zavisimosti Vetrovogo Dreyfa L'da v Azovskom More ot Izmeneniy Koeffitsientov Treniya na Granitse Razdela Vozdukh-Led i Led-Voda [Analysis of the Dependence of Ice Drift in the Azov Sea on Changes of the Values of Drag Coefficients for Wind and Water]. Processes in GeoMedia, (5), pp. 28-36 (in Russian).
4. Hunke, E.C. and Holland, M.M., 2007. Global Atmospheric Forcing Data for Arctic Ice-Ocean Modeling. Journal of Geophysical Research: Oceans, [e-journal] 112(C4). C04S14. https://doi.org/10.1029/2006JC003640
5. Kulakov, M.Yu., Makshtas, A.P. and Shoutilin, S.V., 2013. Verifikatsiya Dannykh Reanaliza NCEP/NCAR po Rezul'tatam Nablyudeniy na Dreyfuyushchikh Stantsiyakh «Sevemyy Polyus» [Verification of the NCEP / NSAR Reanalysis Data by Observations at the Drifting Station "North Pole"]. Arctic and Antarctic Research, (1), pp. 88-96 (in Russian).
6. Fedorenko, A.V., 2015. Pokholodaniya na Azovskom More i Obshchaya Tsirkulyatsiya Atmosfery nad Severnym Polushariem [Cold Spells on the Sea of Azov and the General Circulation of the Atmosphere over the Northern Hemisphere]. In: E. S. Nesterov, 2015. Proceedings of Hydrometcentre of Russia. Moscow: Roshydromet. Iss. 354, pp. 138-154 (in Russian).
7. Borovskaya, R.V. and Klapan, S.N., 2011. Osobennosti Ledovykh Usloviy Kerchenskogo Proliva Zimoy 2008-2009, 2009-2010, 2010-2011 godov [Features of the Kerch Strait Ice Conditions in Winter of 2008-2009, 2009-2010, 2010-2011]. In: YugNIRO, 2011. Proceedings of Southern Scientific Research Institute of Fisheries and Oceanography. Kerch: YugNIRO. Vol. 49, pp. 123-129 (in Russian).
8. Kallos, G., Nickovic, S., Papadopoulos, A., Jovic, D., Kakaliagou, O., Misirlis, N., Boukas, L., Mimikou, N., Sakellaridis, G. [et al], 1997. The Regional Weather Forecasting System SKIRON: an Overview. In: Proceedings of the International Symposium on Regional Weather Prediction on Parallel Computer Environments. Athens, 15-17 October 1997. Athens. pp. 109-122.
9. Dee, D.P., Uppala, S.M., Simmons, A.J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M.A., Balsamo, G. [et al], 2011. The ERA-Interim Reanalysis: Configuration and Performance of the Data Assimilation System. Quarterly Journal of the Royal Meteorological Society, [e-journal] 137(656), pp. 553-597. doi:10.1002/qj.828
10. Semtner Jr., A.J., 1976. A Model for the Thermodynamic Growth of Sea Ice in Numerical Investigations of Climate. Journal of Physical Oceanography, [e-journal] 6(3), pp. 379-389. https://doi.org/10.1175/1520-0485(1976)006<0379:AMFTTG>2.0.CO;2
11. Lecomte, O., Fichefet, T., Vancoppenolle, M. and Nicolaus, M., 2011. A New Snow Thermodynamic Scheme for Large-Scale Sea-Ice Models. Annals of Glaciology, [e-journal] 52(57), pp. 337-346. doi:10.3189/172756411795931453
12. Bogorodsky, P.V., Marchenko, A.V. and Pnyushkov, A.V., 2007. Osobennosti Formirovaniya Pripaynogo L'da v Beregovoy Zone Zamerzayushchikh Morey [Features of Land Fast Ice Forming at Freezing Seas Coastal Zone]. Arctic and Antarctic Research, (3), pp. 17-27 (in Russian).
13. Klyachkin, S.V., Guzenko, R.B. and May, R.I., 2015. Chislennaya Model' Evolyutsii Ledyanogo Pokrova Arkticheskikh Morey Dlya Operativnogo Prognozirovaniya [Numerical Model of the Ice Cover Evolution in Arctic Seas for the Operational Forecasting]. Ice and Snow, [e-journal] 55(3), pp. 83-96 (in Russian). https://doi.org/10.15356/2076-6734-2015-3-83-96
14. Bukatov, A.E., Zavyalov, D.D. and Solomakha, T.A., 2017. Thermal Evolution of the Sea Ice in the Taman Bay and the Dinskoy Gulf. Physical Oceanography, [e-journal] (5), pp. 19-30. doi:10.22449/1573-160X-2017-5-19-30
15. Doms, G., Förstner, J., Heise, E., Herzog, H.-J., Mironov, D., Raschendorfer, M., Reinhardt, T., Ritter, B., Schrodin, R. [et al], 2011. A Description of the Nonhydrostatic Regional COSMO Model. Part II: Physical Parameterization. Offenbach, Germany: Deutscher Wetterdienst, 154 p. Available at: http://www.cosmo-model.org/content/model/documentation/core/cosmoPhysParamtr.pdf [Accessed: 10 May 2019].
16. Kazakova, E.V., Chumakov, M.M. and Rozinkina, I.A., 2014. Algoritm Rascheta Vysoty Svezhevypavshego Snega, Prednaznachennogo Dlya Postprotsessinga Sistem Atmosfernogo Modelirovaniya (na Primere COSMO) [Algorithm of Fresh Snow Height Calculation Based on Results of Precipitation Forecasts by Numerical Atmosphere Models]. In: E. S. Nesterov, 2015. Proceedings of Hydrometcentre of Russia. Moscow: Roshydromet. Iss. 350, pp. 164179. Available at: http://method.meteorf.ru/publ/tr/tr350/kazakova.pdf [Accessed: 10 May 2019] (in Russian).
17. Osokin, N.I., Sosnovskiy, A.V. and Chernov, R.A., 2017. Koeffitsient Teploprovodnosti Snega i Ego Izmenchivost' [Effective Thermal Conductivity of Snow and Its Variations]. Kriosfera Zemli, 21(3), pp. 60-68 (in Russian). doi:10.21782/KZ1560-7496-2017-3(60-68)
18. Roeckner, E., Bauml, G., Bonaventura, L., Brokopf, R., Esch, M., Giorgetta, M., Hagemann, S., Kirchner, I., Kornblueh, L. [et al], 2003. The Atmospheric General Circulation Model ECHAM 5. Part I: Model Description. Hamburg: Max-Planck-Institut for Meteorologie. Report № 349, 140 p. Available at: https://www.researchgate.net/publication/258437837_The_atmospheric_general_circulation_ model_ECHAM_5_PART_I_model_description [Accessed: 10 May 2019].
19. Yang, Yu., Lepparanta, M., Cheng, B. and Li, Zh., 2012. Numerical Modelling of Snow and Ice Thicknesses in Lake Vanajavesi, Finland. Tellus A: Dynamic Meteorology and Oceanography, [e-journal] 64(1). 17202. https://doi.org/10.3402/tellusa.v64i0.17202
20. Anderson, E., 2006. Snow Accumulation and Ablation Model - SNOW-17. [online] Available at: http://www.nws.noaa.gov/oh/hrl/nwsrfs/users_manual/part2/_pdf/22snow17.pdf [Accessed: 10 May 2019].
About the author:
Dmitry D. Zavyalov - Senior Research Associate, Oceanography Department, FSBSI MHI
(2 Kapitanskaya Str., Sevastopol, 299011, Russian Federation), Ph.D. (Phys.-Math.), Scopus Author
ID: 6506347014, ORCID ID: 0000-0002-7444-980X, zavyalov.dd@mhi-ras.ru
The author has read and approved the final manuscript.
The author declares that he has no conflict of interest.
