Wind Drift, Breakdown and Pile up of the Ice Field

1. Introduction

Ice piling near the coastal and shelf structures is the characteristic phenomenon in the water areas of the freezing seas. Ice piles reach especially large volumes in the Arctic seas and impede the operation of the drilling and production platforms and the marine terminals. Therefore, for effective management of the ice conditions within these water areas, the short-term forecast of this phenomenon in according with hydrometeorological conditions is necessary.

The mechanism of the ice fields breakdown in contact with a fixed obstacle was the subject of the theoretical and model studies and had been fairly well known at present [1,2,3,4]. The drift of the ice cover in whole and the drift of the single ice fields, their pressure on the obstacle or the collision with it, accompanied by brittle destruction into small ice cakes, which form the pile owing to sinking part of ice cakes to the bottom and rising of rest part to the ice surface, are the obvious and understandable mechanisms of the phenomenon. The intensity of impact with an obstacle or the magnitude of the pressure at the point of contact depends on the drift velocity of the ice cover as a whole and the velocity of single ice field caused by the wind and the currents.

The compact ice cover drift objectives on the Arctic Basin are also quite well studied and the mathematical models have been developed to this process forecast [5,6]. However, these models are not applicable for local spatial scales of the water area adjacent to the offshore structures, and do not make it possible to predict the characteristics of the ice cover contact with the obstacles and its consequences.

There are investigations of the drift of the single ice field. The main problem is modeling the shape of the ice field, which, in real conditions, is formed under the influence of processes occurring in the drifting ice cover, including collisions, hummocking, freezing and melting. The mathematical description of the movement of the ice floes requires the parametrization of the impact of the external forces (wind and current) and the resistance to movement. For this purpose, the real irregular shape of the ice floe was simplified and represented by the flat disk for analytical description [7], and by the flat rectangular plate in the experimental studies [8].

The problem of the parametrization of the resistance to movement of the objects with the irregular shape can be solved within the theory of the ship propulsion [9]. The resistance for the ship movement is divided into separate components, each of which has its own method for the calculation. For the friction resistance evaluation, the method of the equivalent flat plate, that has same area and length as the ship hull, is applied. This approach is suitable for the single ice field with an irregular shape in plan. When assessing the friction of the wind and the current on the ice field surfaces, the equivalent flat plate approach that has the same surface area and length in the direction of movement can be applied.

This concept was previously used to assess the wind effect on the stationary ice field in contact with the fixed obstacle [10]. Comparison of the generated under the wind action pressure at the point of contact between the ice field and the obstacle with the compressive strength of the sea ice made it possible to estimate the conditions for the ice destruction at the point of contact (the dimensions of the ice field and wind speed) and, as a consequence, the occurrence of the ice pile near the obstacle and its possible volume.

In this article, same approach is applied to the ice field moving under the action of the wind and the current. The model of the dynamic destruction of the ice field is based on the notion that sea ice, as the solid material, has certain specific energy of the brittle destruction [11,12,13]. The velocity of the ice field drift determines the kinetic energy of the ice field, which is released, when the ice field suddenly stops owing to its collision with an unmovable obstacle. This kinetic energy is mostly realized in the destruction of the ice as material that is in the formation of any fragments – pieces, the volume of which corresponds to the specific energy of the ice destruction and which form the ice pile in the place of collision. It is reasonable to consider that form and dimensions of this ice pile are similar to the ice ridges, which are formed during the collision of drifting ice fields and are quite well studied [14,15].

The paper presents the model of the ice field movement under action of the wind and current that is differential equation that takes into account the effect of wind and current on the ice field, the resistance to movement from air and water masses, as well as initial conditions and a description of the variation over time in the velocity and direction of the wind and current.

This model was applied to the wind drift of an ice field with arbitrary shape over a stationary water mass. Solving the equation made it possible to estimate its kinetic energy. Based on the assumption that all the kinetic energy at collision with the obstacle and stopping the ice field is realized in the brittle destruction of the ice field, the estimates were made of the potential volume of the ice cakes pile (piling up) dependence on the wind speed and the ice field dimensions.

The results of the study were previously presented at 26^th IAHR International Symposium on Ice [16]. This publication is supplemented by the results of assessing the shape and dimensions of ice piles that can form at the place of a drifting ice field collision with a fixed obstacle.

Research results can be used to predict the formation of the ice ridges and piles near the offshore structures.

2. Problem Statement

The problem under consideration is quite complicated. Therefore, in order to assess the consequences of the collision of an ice field with an obstacle, this process had to be schematized. Three processes are considered and three tasks (subproblems) are solved, namely, the drift and the destruction of the ice field, and then the formation of a pile from the ice fragments. These tasks are solved by different methods, but they are connected in such way: the results of the first one solving are the basis for solving the second task. To solve the third task, which is the purpose of the study, the results of solving the second task are used.

The first task is to describe the movement of the ice field under the influence of the external effects: the wind and the currents that act on the ice field, setting it in motion through the friction forces and the velocity pressure on the end surfaces of the ice field. To do this, it is required to construct the equations of the ice field motion under the action of the wind and the current in some sufficiently large water area so without the necessity to consider the influence of the boundaries and the shallow water. The solution of this equation makes it possible to find the speed and the trajectory of the ice field drift under the action of the time-varying wind and the currents and to estimate its kinetic energy.

The solution of the second problem - the assessment of the volume of the ice field fragments as result of the breakdown under the collision with the fixed obstacle - does not involve the study of the process of the brittle fracture itself, which requires the special analysis of the processes occurring in the crystal structure of the sea ice. For the practical application of under consideration problem, it suffices to estimate "from above", which consists in comparing the kinetic energy of the ice field, which depends on its dimensions and the drift velocity, and the specific energy of the brittle destruction of the ice. The main assumption is that all the kinetic energy is converted and spent on the destruction of the ice field. This makes it possible to estimate the potential volume of the ice fragments depending on the wind speed and the ice field dimensions and its thickness.

There are sufficient grounds to assume that the pattern of formation of the pile from ice fragments hear an obstacle after collision with the drifting ice field are similar to the patterns of formation of the ice ridges (hummocks) during the collision of the ice fields. This assumption allows us to consider that the configuration (geometric form) of the ice ridge and the ice pile are similar and on the basis of generalized data on the structure of the ice ridges to construct a method to estimating the height of the pile (like a hummock sail) above the ice cover.

3. Equation of the Ice Field Drift

The movement of the ice field with the specific dimensions and mass under the effect of the wind and the current is considered. It is assumed that the ice field does not contact other ice fields and drifts far enough from the shores of the water area to have no contact with the coastline, and the change in the depth of the sea does not affect the current speed. The wind is assumed to be uniform over the area of the ice field and its characteristic is the velocity vector, which does not change in its direction and modulus. The same assumption is made for the water flow under the ice field. The friction of the water mass and the atmospheric air on the surface of the ice field occurs in the turbulent regime, and the coefficient of friction is determined by the regularities that were determined for a flat turbulent boundary layer [9,17].

Under the stated assumptions for the coordinate system, where OX axis is directed along the parallel, and OY axis along the meridian, the equation of the motion of the ice field in the vector form is following (Figure 1).

$$\left(M+m_w\right)\hspace{0.17em}\left(\frac{d\overrightarrow{V}}{dt}+f\hspace{0.17em}\overrightarrow{V}\right)={\gamma }_a{\zeta }_{as}\left|\overrightarrow{W}-\overrightarrow{V}\right|\hspace{0.17em}\left(\overrightarrow{W}-\overrightarrow{V}\right)\hspace{0.17em}{S}_{aw}+{\gamma }_w{\zeta }_{ws}\left|\overrightarrow{U}-\overrightarrow{V}\right|\hspace{0.17em}\left(\overrightarrow{U}-\overrightarrow{V}\right)\hspace{0.17em}{S}_{uw}\hspace{0.17em}.$$

(1)

In this equation, M is the mass of the ice field, m_w is the added mass of sea water, S_aw is the characteristic surface area of the ice field for its part rising above the water surface, S_uw is the characteristic surface area of the ice field under the water surface, V is the speed of the ice field motion, W – the wind velocity, U – the current velocity, t – the time, f – the Coriolis parameter, γ_a – the atmospheric air density, γ_w – the sea water density, ζ_a =f_a(W-V) – the air resistance coefficient on the ice field surface, depending on the velocity of its drift and the wind speed, ζ_w = f_w(U-V) is the coefficient of the water mass resistance to the movement of the ice field.

The initial conditions for the Equation (1) are as follows

$$\begin{array}{l}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}t=0\hspace{1em}x=0,\hspace{0.33em}y=0,\hspace{0.33em}\overrightarrow{V}=0;\hspace{1em}\\ \overrightarrow{W}=f_{W1}\left(0\right)\times cos\left(\overrightarrow{W},OX\right)=f_{W2}\left(0\right);\hspace{1em}\overrightarrow{U}=f_{U1}\left(0\right)\times cos\left(\overrightarrow{U},OX\right)=f_{U2}\left(0\right).\end{array}$$

(2)

The equation of the movement should be supplemented with the dependences of the velocity and the direction of the wind and the current on the time

$$\left|\overrightarrow{W}\right|=f_{W1}(t),\hspace{1em}\cos \left(\overrightarrow{W},OX\right)=f_{W2}\left(t\right);\hspace{1em}\hspace{0.33em}\hspace{1em}\left|\overrightarrow{U}\right|=f_{U1}(t),\hspace{1em}\cos \left(\overrightarrow{U},OX\right)=f_{U2}\left(t\right)\hspace{0.17em}.$$

(3)

The components included in the right side of the Equation (1) characterize the effects of the wind and the current on the ice field. If the wind velocity exceeds the drift velocity of the ice field, that is (W - V)_{x, y} > 0, then the wind "moves" the ice field, otherwise: (W - V)_x,y < 0, the air medium "slows down" the drift of the ice field. The action of the current is similar: at (U – V)_x,y > 0, the current sets the ice field in motion, and at (U – V)_x,y, the water mass prevents its drift. The magnitude of the resistance forces depends on the velocity of the ice field relative to the wind velocity and relative to the water surface. The wave resistance is assumed to be negligible due to the relatively low ice drift velocity.

The ice field has the irregular shape, which was formed during the drift as result of the interaction with other ice fields. It is reasonable to represent the ice field as the rectangular parallelepiped with following dimensions: L_if - the length, В_if - the width and h_ice - the ice thickness (Figure 1). In order not to consider the orientation of the ice field relative to the trajectory of its motion, as the characteristic linear dimension for calculating the friction and the Reynolds number Re, it is reasonable to take the following value, which is close to the average size for the “non-elongated” shape of the ice field in plan and is determined by the following formula (S_if – the ice field area):

$$L_f=\sqrt{L_{if}\hspace{0.17em}{B}_{if}}=\sqrt{S_{if}}\hspace{0.17em}.$$

(4)

4. Steady Wind Drift of the Ice Field

It is considered the drift of the ice field along the stationary water mass under the action of the wind with the constant velocity and direction. The curvature of the drift trajectory due to the Coriolis effect is not took into account. In this case, the equation of the motion (1) takes the following form (the first term is the effect of the wind causing drift, the second term is the resistance of the water mass to the movement of the ice field.)

$${\gamma }_a{\zeta }_{as}{\left(w-v\right)}^2\hspace{0.17em}{S}_{aw}-{\gamma }_w{\zeta }_{ws}\hspace{0.17em}{v}^2\hspace{0.17em}{S}_{uw}=0\hspace{0.17em}.$$

(5)

The coefficients of the air resistance ζ_as and the water one ζ_ws included in the equation consist of two components each: ζ^a_vh – the head (frontal) resistance of the part of the ice field rising above the water level (the head resistance for the submerged under water part - ζ^w_vh) and the frictional resistance of the air ζ^a_fr and the water ζ^w_fr on the flat surfaces of the ice field. The scheme of action of these forces is shown in Figure 2.

In this case, it is possible to present the Reynolds number and the coefficient of the friction resistance on the interface atmospheric air – ice with following form

$${\mathrm{Re}}_{ai}=\frac{w_{sf}}{{\nu }_a}{L}_f,\hspace{1em}\hspace{1em}\hspace{1em}{\zeta }_{fr}^a=\frac{0.455}{{\left(\lg {\mathrm{Re}}_{ai}\right)}^{2.58}}\hspace{0.17em}.$$

(6)

Here, w_sf – the wind velocity on the upper surface of the ice field, which is determined in respect to the surface rise over the water surface h_ia, ν_a – the kinematic viscosity of the atmospheric air. Value h_ia is connected with the ice field thickness h_ice and the ice density γ_ice and the sea water density γ_w by following form

$$h_{ia}=\frac{{\gamma }_w-{\gamma }_{ice}}{{\gamma }_w}{h}_{ice}.$$

(6)

Applying known formula for the wind velocity within the boundary layer over the sea surface [18] the following formula was deduced

$$w_{sf}=w(h_{ia})=W_{10}\left(1+\frac{\sqrt{{с}_w}}{{\chi }_0}\ln \frac{h_{ia}}{z_{10}}\right)\hspace{0.17em}.$$

(8)

In this case W₁₀ (= w) – the wind velocity on the height z₁₀ = 10м above the sea level, χ₀ = 0.4 – the Karman constant and c_w – the coefficient of the sea surface friction. This coefficient depends on the degree of the wave roughness of the sea surface. For the sea surface covered with drifting ice, this coefficient can be estimated by its minimum value c_w = 1.11⋅10^-3.

The friction resistance of the ice field on the water mass depends on the drift velocity v and is determined similarly to the formulas for air resistance (6), that is

$${\mathrm{Re}}_{wi}=\frac{v}{{\nu }_w}{L}_f,\hspace{1em}\hspace{1em}\hspace{0.33em}\hspace{1em}{\zeta }_{fr}^w=\frac{0.455}{{\left(lgR{e}_{wi}\right)}^{2.58}}\hspace{0.17em}.$$

(9)

Head resistance coefficients can be assumed to be same for the over surface and under surface parts of the ice field. The front and back faces should be considered as the extended barriers, and the following value ζ^a_vh = ζ^w_vh = 2.1 can be accepted [19].

Substituting the obtained results into Equation (4), we get the following transcendental equation, which allows estimating the drift velocity of single ice field

$${\gamma }_a\left[{\zeta }_{vh}^a\hspace{0.17em}{h}_{ia}+{\zeta }_{fr}^a\right.\hspace{0.17em}\left.L_{if}\right]{\left(w-v\right)}^2\hspace{0.17em}-{\gamma }_w\left[{\zeta }_{vh}^w\hspace{0.17em}\left(h_{ice}-h_{ia}\right)+{\zeta }_{fr}^w\hspace{0.17em}{L}_{if}\right]\hspace{0.17em}\hspace{0.17em}{v}^2=0\hspace{0.17em}.$$

(10)

Figure 3 presents the dependence of the ice field (dimension L_f = 200 m) drift velocity on the wind velocity for various thickness of an ice. There is also the dependence of the compact ice cover drift velocity on the wind velocity (the low: 2% of the wind speed) that was established by the observations in the Arctic Ocean [5,6], and the velocity of the wind drift current at the latitude ϕ = 70°N calculated by following formula [19]

$$v_{wd}=\frac{0.0127\hspace{0.17em}{W}_{10}}{\sqrt{\sin \hspace{0.17em}\phi }}.$$

(11)

The presented materials show that the velocity of the wind drift of single ice field exceeds significantly the wind surface current velocity and the velocity of the compact ice cover wind drift in the same water area. These materials also show that the drift velocity of the ice field decreases with increasing an ice thickness. It is possible to explain by following way: the contribution of the head resistance of the submerged part of the ice field into the total resistance to the wind drift increases as the ice thickness growths.

5. Evaluation of the Kinetic Energy of Drifting Ice Field

The value of the kinetic energy of a body is determined by its mass and velocity. For the drifting ice field with thickness h_ice and characteristic dimension L_f, the kinetic energy is the function of the wind speed w and is given by following formula

$$E_{kif}(w)=\frac{1}{2}{\gamma }_{ice}\hspace{0.17em}{h}_{ice}\hspace{0.17em}{L}_f^2\hspace{0.17em}{\left[v\left(w\right)\right]}^2.$$

(12)

Expression (12) shows that the kinetic energy of single ice field increases with the growth of its thickness and area and increases in proportion to the square of the drift velocity v(w). At the same time, the drift velocity of the ice field v depends itself on the dimensions and the thickness of the ice field, therefore, to assess the specific nature of this dependence, the special numerical simulation was performed with variation the included in (12) parameters by solving Equation (10).

Figure 4 shows the dependence of the kinetic energy of the drifting ice field on the wind speed for various ice thicknesses. The presented results of modeling give possibility to approximate the dependence of kinetic energy of the ice field on the wind speed by the power dependence

$$E_{kif}(w)\sim w^{1.97} .$$

(13)

Figure 5 shows the dependence of the kinetic energy of the drifting ice field on its thickness for various wind velocities. These data give possibility to approximate the dependence of the kinetic energy on the thickness of ice field h_ice by following power formula

$$E_{kif}(h_{ice})\sim h_{ice}^{0.748}.$$

(14)

The results of computations of the dependence of the kinetic energy of the drifting ice field on its characteristic dimension for various wind velocities (presented on a logarithmic scales) show the existence of a power-law dependence between the kinetic energy and dimension. This dependence has the following form

$$E_{kif}(L_f)\sim L_f^{2.18}.$$

(15)

Thus, the multifactorial dependence of the kinetic energy of the drifting ice field on its size and the wind velocity appears, which is contained in equation (10) implicitly form. Formulas (13)–(15) obtained as result of the computer modeling make it possible to estimate qualitatively the change in the kinetic energy of the drifting ice field with its dimensions and the wind speed variation.

6. Possible Volume of the Ice Fragments after Collision the Ice Field with an Obstacle

The sea ice is the brittle material; therefore, it fails as result of the mechanical impact, and small ice pieces (fragments, ice cakes) are formed. The ratio of the impact energy and the volume of the formed fragments is characterized by the impact energy of the destruction. According to the available measurement results [11,12,13], the value of the specific energy of the sea ice destruction lies within the fairly wide range: e_F = 0.34 ÷ 4.5 kJ/m³. These results were obtained when testing the ice samples that had the salinity from 0 to 6.0‰ and the temperature from -1°С to -60°С.

The obtained solution for the kinetic energy value of the drifting ice field (12) allows to estimate the volume of the ice destruction (volume of the ice fragments) after the drifting ice field collision with the fixed obstacle. It is assumed the obstacle is absolutely solid and its length along the waterline exceeds the dimensions of the ice field. The estimation “from above” obtained in this case has the following form

$$Q_{if}(w)=\frac{{\gamma }_{ice}}{2\hspace{0.17em}{e}_F}\hspace{0.17em}{h}_{ice}\hspace{0.17em}{L}_f^2\hspace{0.17em}{\left[v(w)\right]}^2.$$

(16)

Using this formula, it is possible to evaluate the volume of the ice fragments that can appear after the drifting ice field collision for some average value of the specific energy of the ice destruction, for example, e_F = 1.5 kJ/m³. Figure 6 shows the results of calculations: the dependence of the volume of the ice fragments on the wind speed for the ice field with the characteristic dimension L_f = 200 m for several ice thicknesses, which are typical for the first-year ice. These materials demonstrate, when the drifting ice field breaks down in the collision with the fixed obstacle, the large (hundreds of the cubic meters) volumes of the ice fragments can be formed, which will form the pile with significant volume that increases with the ice thickness increasing.

The volume of the ice fragments can be comparable to the total volume of the ice field. Figure 7 shows, for the wind speeds above 20 m/s, the relative volume of the ice fragments may exceed 10% of the pre-collision ice field volume.

7. Estimation of the Possible Dimensions of the Ice Fragments Piles

Accumulations of the ice fragments near the hulls of the stationary drilling and mining platforms make difficult their operation, as the ice piles impede the mooring for the workboats and tankers. Above stated results of the performed investigations of the wind drift and breakdown of the ice field give possibilities for evaluations of the dimensions of the ice piles near board of a stationary platform. For this purpose it is assumed the collision of drifting ice field with a platform is the main mechanism for the formation of the ice accumulation near the boards of the platforms and other offshore structures [4,20,21].

There are sufficient reasons to believe that the ice fragments formed as result of the collision of a drifting ice field with an obstacle form a pile or ridge just as it happens when drifting ice fields collide with each other and hummocks and ridges of hummocks are formed.

Therefore, it is possible to apply the results of studying the shape of the hummocks and the ridges of ice fragments to estimate the sizes of the pile of ice fragments in front of the obstacle [14,15,22,23]. The main characteristics that is reasonable apply for analysis are following [15]:

• the ratio of the sail height H_s to the keel depth H_k is equal to Hk / H_s = β₁ = 4.4,
• the ratio of the keel width W_k to the sail height H_s is equal to W_k / H_s = β₃ = 15.1,
• the slope angle of the sail α_s = 20.7° and slope angle of the keel α_k = 26.6°.

Figure 8 shows the schema of the pile of ice fragments in front of the obstacle wall, which is half of the ridge the average cross section of which is presented in [15]. The adopted schema of the ice pile makes it possible to represent the cross-sectional area of the pile up as the sum of the sail in the form of a triangle and the area of the keel in the form of a trapezoid.

The parameters of the ice pile presented above [15] are sufficient to determine unambiguously the cross-sectional area A_ip of the ice section pile in front of wall using a single parameter: the sail height H_s. Accordingly, the height of the sail of pile Hs can be expressed in term of the cross section of pile area A_ip. The corresponding formula has the following form

$$H_s=\sqrt{\frac{2\hspace{0.17em}{A}_{ip}}{ctg({\alpha }_s)+{\beta }_1\left[{\beta }_3-{\beta }_1\hspace{0.17em}ctg\left({\alpha }_k\right)\right]}}.$$

(17)

In order to estimate the elevation of the pile above the ice cover surface H_sel, the thickness of the ice field h_ice and the density of the sea ice γ_ice should be taken into account.

$$H_{s el}=\sqrt{\frac{2\hspace{0.17em}{A}_{ip}}{ctg({\alpha }_s)+{\beta }_1\left[{\beta }_3-{\beta }_1\hspace{0.17em}ctg\left({\alpha }_k\right)\right]}}-\left(1-{\gamma }_{ice}\right)\hspace{0.17em}{h}_{ice}.$$

(18)

The value of the cross-sectional area of the ice pile A_ip can be represented as the volume of ice fragments in the following to versions:

• the volume of the ice fragments per unit length of the pileup, if the width of the ice field colliding with the obstacle is less than the length of the platform board L_pl, or
• the volume of the ice fragments per unit length of the platform board, if the width of the drifting ice field is greater than the length of the platform board.

$$A_{i p}=\frac{Q}{L_f},\hspace{1em}if\hspace{0.33em}{L}_f\le L_{pl};\hspace{1em}\hspace{1em}\hspace{1em}{A}_{i p}=\frac{Q}{L_{pl}},\hspace{1em}if\hspace{0.33em}{L}_f>L_{pl}.$$

(19)

Based on this assumption, the estimations of the height of the ice accumulation near the board of the Prirazlomnaya platform were performed depending on the dimensions of the ice fields and the wind velocity using Formula (18). The computation results are presented in Figure 9.

The estimates of the elevation of the ice pile near the platform board showed the developed model of the ice pile formation made it possible to obtain dimensions close to those observed in real conditions. In particular, by the look of Photo [24] the height of the ice piles observed near the Platform Prirazlomnaya is comparable to the height of the foredeck of the work vessel and to the freeboard of the platform caisson. That is the height of the ice piles in Photo [24] was about 3 m. Results of modeling in Fig. 9 demonstrate that such ice piles could be forming as result of single collision with the platform of the ice fields with dimensions of more than L_f > 500 m in plan and thickness of about h_i = 1 m, that drifted at a wind speed more than w = 15 m/s. These hydrometeorological and ice conditions are observed in the Pechora Sea water area in reality.

8. Discussion

As result of the study, the mathematical model of the isolated ice field drift under the influence of the wind and the current was derived, which takes into account:

• the ice field mass, its dimensions and the attached water mass,
• the action of the air and the water masses, which, if their velocities exceed the drift speed, set the ice field in motion, and otherwise brake owing to resistance to the ice field motion,
• the main components of the impact of the air and the water masses on the ice field that are parameterized in the form of friction on the upper and lower surfaces of the ice field and in the form of the head pressure (resistance) on the end surfaces of the ice field rising above the water and head resistance (pressure) on the submerged part of its forward surface.

The ice field drift equation was applied to the simplest case of the wind drift of the ice field over the unmovable sea water mass, and modeling of this process was performed, which made it possible to understand the dependence of the ice field drift velocity on the wind speed, its horizontal dimensions and thickness. These results made it possible to estimate the kinetic energy of the ice field and the probable volume of the ice fragments, when it collides with the fixed obstacle.

It was assumed that the patterns of the hummocks and ice ridges, on the one hand, and the piles of ice fragments near the stationary platform boards, on the other, as well as their shape and structure, are similar. This made it possible to develop the method and estimate the elevation of the ice piles near the stationary platform board in dependence on the dimension of the ice field and the wind speed.

The height of the ice piles near the Platform Prirazlomnaya estimated by comparing its dimensions with the height of foredeck of work vessel and platform’s caisson freeboard in the photo turns out to be close to the results of calculations based on the developed ice pile model.

9. Conclusions

The performed study is based on the application of methodology developed in the Naval Architect Theory for estimation the water resistance to the ship hull movement, which have not previously been applied in the studies of the dynamic of the sea ice cover.

The results of estimating the elevations of the piles of ice fragments (cakes) in relation to the Platform Prirazlomnaya are in good agreement with the observed data.

These outcomes provide a basis for future studies, in particular, the ice field drift dynamic under varying wind and tidal currents.

The results of the study (developed models) are recommended to be applied to solve the problems of the ice conditions management in the water areas, where the platforms and terminals are located in the Arctic seas.

Developed approach can be applied for investigations the consequences of the ice field and the fixed obstacle collisions under action of another variants of the wind and the currents impacts on the ice field.