# Nonextensive Dynamics of Drifting Sea Ice

Authors: Alexandre Chmel; Victor Smirnov
Journal: ISRN Oceanography
 (2013)

---

## Abstract

Cycles of ice pack fragmentation in the Arctic Ocean are caused by the irregular drift dynamics. In February 2004, the Russian ice-research camp North Pole 32 established on a floe in the Arctic Ocean ceased its working activity and was abandoned after a catastrophic icequake. In this communication, the data collected during the last month of the field observations were used for calculating the changes in the kinetic energy of the ice floe. The energy distribution functions corresponding to periods of different drift intensity were analyzed using the Tsallis statistics, which allow one to assess a degree of deviation of an open dynamic system, such as the drifting ice, from its equilibrium state. The obtained results evidenced that the above-mentioned critical fragmentation has occurred in the period of substantially nonequilibrium dynamics of the system of ice floes. The determination of the state of the pack (in the sense of its equilibrium/nonequilibrium) could provide some useful information on forthcoming icequakes.


---

## Body

## 1. Introduction

From the viewpoint of conventional mechanics, the Arctic sea ice cover (ASIC) is the consolidated, mobile, deformable system. Prevalently shearing deformations result in regular pattern of fragmented pack in accordance with the Mohr’s mechanism of semibrittle failure of solids.At the same time, the size distribution of sea ice floes does not exhibit the random (Poissonian-like) statistics and follows the power law typical for scaling (fractal) structures [1, 2]. The conventional mechanics cannot explain this phenomenon.The scaling is a manifestation of long-range correlations between separated events in a statistical system. The correlation radius is determined by the spatial decay of the event effect. The decaying is fast (exponential) in equilibrium systems but slow (power law) in nonequilibrium ones. Therefore, the fractal structures are formed only under nonequilibrium conditions in statistical systems driven by external forcing.In recent years, a new statistical concept called the non-extensive statistical mechanics (NESM) was developed by Tsallis [3, 4] for thermodynamic description of the behavior of multiscale systems revealing fluctuations around their equilibrium state. The NESM was successively applied for assessing the degree of deviation of natural dynamic systems from their equilibrium state prior to large-scale hazards, such as earthquakes [5, 6] and floods [7]. In the ASIC, cycles of pack fragmentation and significant icequakes occur due to irregularity in the ice drift, which causes critical deformations in the ice cover. According to the NESM, the formation of fractal structures, such as fragmented ice pack, is a result of nonequilibrium dynamics of open systems. In this work, the sea ice drift dynamics during a few bursts of significant fragmentations, which could be regarded as natural hazards, were analyzed in terms of the NESM.The data on the local drift dynamics were collected in the Russian ice research camp North Pole 32 (NP 32) that was established on ice pack and drifted in 2003-2004 in the region north from Franz Josef Land from 88°N to 84°N [8]. In this work, the time series of the kinetic energy variation of an individual ice floe were analyzed using the Tsallis statistics in order to obtain a thermodynamic characteristic of the system of drifting ice prior to and during large-scale sea ice fragmentations.


## 2. Geometry and Size

A distinguishing feature of the Arctic sea-ice cover (ASIC) is a rectilinear lead pattern characterized by the constant angle ~30° between directions of intersecting leads (Figure1). This regular structure consists of parallelogram-like figures ranging from 1 to ~102 km. Marco and Thompson [9] explained this phenomenon on the base of the Mohr’s theory of failure, which suggests the following relation between the applied shear stress at failure, τ, and the normal stress, σn, on the failure plane (or line in 2D-case):(1)τ=C+σntanφ,
where C is the material cohesion strength, and φ is the angle of internal friction, φ=±(α-45°); here α is the intersection angle. The resolution of (1) relatively the intersection angle with parameters characteristic for sea ice gives one α≈30°, that is the value observed in many satellite images (Figure 1).Figure 1
AVHRR image taken from NOAA satellite on 6 December 2003 in the region around 82°N, 5°E.The parallelogram-like pattern appears only at sufficiently rapid strain rates when the semibrittle failure prevails. At relatively slow strain rates, the ductile behavior takes place, which results in formation of polygonal (diamond-like) or fully disordered structures due to cohesion flow at low yield stress [9].At the same time, there are some features in the lead pattern, which cannot be explained in the framework of classical mechanics. Rothrock and Thorndike [10] and Matsushita [11] studied the size distribution of sea ice floes and revealed a scale invariant relation between the number of floes and their size in the form(2)N(z>z′)∝z-β,
where N(z>z′) is the number of floes larger than z′; z is the dimensional characteristics, such as area, S, or length, L; β is the constant. The power law dependence means the self-similarity of the ensemble of ice floes because the power law function (2) is the unique solution of the scaling relation
(3)N(λz)=λ-βN(z),
where λ is the constant (scaling factor). In general, the power law behavior is regarded as a representation of the scale invariance of an object or process in space, time, or energy (in dependence of the measured scaling parameter).The scaling (power law) number-size distributions were found for lead patterns formed both by ice floes of regular geometry, such as parallelogram-like and diamond-like floes, and by “smoothed” (oval) floes typical for frost-free seasons.Another example of the invariance of the floe geometry was found by Chmel et al. [2] who had shown that the floe area and floe perimeter are interconnected through the following relation:(4)S∝L2/D.The spatial invariance of this kind is characterized by the parameterD called the “fractal dimension”, which, as distinct from the integer Euclidian dimension, could take a fractional value. In the “two-dimensional” case, D could take any value from 1 to 2. Highly fragmented pack with a vast of one dimensional features, such as leads and cracks is characterized by D close to 1, while the fractal dimension of consolidated, frozen ice fields should be closer to 2.Figures2 and 3 show examples of the power law dependences (2) and (4) calculated for lead-and-crack pattern in the ice cover shown in Figure 1.Figure 2
Cumulative distribution of floe sizes calculated by processing the satellite image shown in Figure1. Straight line is the best fitted power law dependence (2).Figure 3
Area versus perimeter dependences for ice-floes fragments resolved in the satellite image shown in Figure1. The straight line fits the power law dependence (4).The fragmentation of the Arctic ice pack is caused by non-uniform sea ice drift, which is accompanied with substantial stress redistributions. To understand the driving force of self-organizing in fractal structures one should consider the thermodynamic conditions of the drift dynamics.


## 3. Nonextensive Dynamics

The thermodynamic state determines the principal trends in the behavior of every statistical system. The most significant perturbations occur when the system deviates from its equilibrium state. Under nonequilibrium conditions, individual “events” in the system are not independent but affect each other at distances much longer than the size of an excited locality. A substantial effect of local event on adjacent sites disturbs the entropic additivity of the process, which in the simplest case of two independent equilibrium subsystemsA and B can be expressed as:
(5)S(A+B)=S(A)+S(B),
where S is the entropy.The thermodynamic description of processes in nonequilibrium statistical systems was developed by Tsallis [3, 4]. The main idea forwarded in his approach was to introduce into additive expression (5) a cross-member that would take into account some interactions between subsystems through a “parameter of nonextensivity” q:(6)S(A+B)=S(A)+S(B)+S(A)S(B)(1-q)k.
Here k is the Boltzmann constant.Passing over details of the Tsallis’ consideration (see [3, 4] for details), we note that the classical statistical mechanics responds the limit q→1 ((6) transforms into (5)); for q<1 the formalism imposes a limited number of events in equilibrium system; finally, the value q>1 indicates the nonequilibrium state of the dynamic system—this is the case of non-extensive dynamics, which is incompatible with the entropic additivity.The value of the parameterq serves as a measure of deviation of the statistical system from its equilibrium state. The stronger deviation, the higher q. Accordingly, the higher q, the larger correlation radius of interactions, and the cooperative effects are more pronounced.Abe and Suzuki [12, 13] used the non-extensive paradigm in their analysis of the distance and waiting time distributions between successive earthquakes and derived a universal expression for the probability of events in nonequilibrium system:
(7)P(x>x′)=[1-(1-q)x′x0]1/(1-q),
where P(x>x′) is the probability of being the value x higher than a threshold x′, x0 is the mean value of x in the series.Equation (7) was used by Abe and Suzuki for fitting the parameters q and x0 to known distances [12] and waiting time [13] distributions between earthquakes. The expression for the cumulative energy distribution N(ε>ε′) versus ε′ can written in the following form:
(8)N(ε>ε′)=P(ε>ε′)×N=N[1-(1-q)ε′ε0]1/(1-q),
where N(ε>ε′) is the number of events satisfying the inequality ε>ε′, and N is the total number of events. Some opportunities given by the applying of (8) in studies of sea ice drift dynamics were demonstrated recently by Chmel et al. [14].In this work we used (8) for analyzing the time series of changes in drift speed of an individual ice floe. In this context, an “event” should be regarded as a change in the kinetic energy of the floe; ε0 is the mean value of ε in a selected time interval.


## 4. Energy Exchange

The energy exchange between a drifting ice floe and its environment can be found from changes in the speed of the floe caused by impact or shear interactions with neighboring floes. The kinetic energy change,ε, is directly proportional to the speed change squared of the ice floe, ε∝ΔV2 (where ΔV2=|V′2-V′′2|, V′ and V′′ are the ice floe speeds measured in two subsequent time intervals). Therefore, the distribution function, N(ε>ε′), is equivalent to the distribution of the values ΔV2.The time series of changes in the kinetic energy of the ice floe were recorded in the Russian ice research camp NP 32.In February 2004, a series of icequakes caused the intensive ice cover fragmentation in the region of the NP 32 activity. The GPS monitoring of the camp position showed a series of local pack displacements with the most intensive shocks on 18 February. In addition, multiple fragmentations of the ice floe, on which the NP 32 drifted, were detected as beginning from the day of the main lead propagation and during subsequent days of February. As a result, the scientific observations were ceased on 26 February; on 2 March 2004 the research station NP 32 was abandoned in connection with multiple breakage of the ice floe. The time series of the floe kinetic energy changes were now used for assessing the disequilibrium of the pack state prior to, during, and after the occurred icequake.The valuesV were found from the changes in the geopositioning data determined by a GPS transducer. The changes in the ice floe speed were measured in 10 minutes intervals. The measurement procedure and the accuracy of the determination of changes in the value of the ice floe speed using GPS transducer are given in [15]. A time series of kinetic energy changes (proportional to the measured ΔV2 values) in the time period of drift from 1 February to 26 February 2004 is depicted in Figure 4(a). The whole period of observations was divided in four “time windows” in accordance the varied intensity of floe displacements in different periods of time.Figure 4
(a) Amplitudes of kinetic energy changes found from the ice floe speed data in the time interval 1 February to 26 February 2004. A few gaps are due to instrumental failures caused by local breakings in the ice floe. (b) and (c) Histograms of the mean value of kinetic energy changeε0  and the parameter of nonextensivity  q, respectively.Figure5(a) shows the energy distributions N(ε>ε′) versus ε′ in each time window. The measured distributions were approximated by the function (8) using the procedure of best fitting applied to the parameters ε0 and q. The calculated analytical curves are shown in Figure 5(a) by lines. The found best-fitted values of the mean value of the energy change and the parameter of nonextensivity are shown as histograms in Figures 4(b) and 4(c), respectively. One can see that the parameter ε0 (Figure 4(b)) reproduces adequately the intensity of ice floe mobility (Figure 4(a)) in different time intervals. In a similar way, the q-value is minimal in the window A that is prior to the beginning of the sea ice fragmentation. Its closeness to unity evidences the almost equilibrium state of the drifting sea ice. In the subsequent interval B characterized by the moderate ice floe mobility, the  q-value increased from 1.11 to 1.19 thus indicating an increased deviation of the ice pack from the equilibrium state. The highest value of the parameter q (1.31) was achieved in window C when the largest changes in the floe kinetic energy were recorded. At the final stage of observations (window D), the value q returned to its initial value.Figure 5
Experimental (signs) and calculated (curves) energy distributionsN(ε>ε′) versus ε′ in time windows A to D (see Figure 4) in the periods of time covering the icequakes in February 2004. The lines fit (8).The variation of theq-value reflects trends in the behavior of ice pack influenced by the external forcing. A deviation of the parameter q from unity means the nonequilibrium state of drifting sea ice. Under this condition, the fragmentation of the pack would tend to form fractal structures.


## 5. Discussion

Classical statistical mechanics are constructed on the implication of the independence of individual events involved in the mechanical process; each event produces an additive contribution. The additivity (extensivity) of this kind is the basic property of closed systems.The mechanical behavior of open, nonequilibrium thermodynamic systems, such as the ASIC, is affected by the energy exchange between its components. It is the open character of the ASIC that causes nonlinear, synergic effects, such as the fractal geometry of lead-and-crack pattern. The flow of the external energy through the system increases the correlation radius of the event cross impact. The effect of local perturbations decays as the power law with covering much more significant neighborhood than that in the case of exponential decay inherent in the system under equilibrium condition. The principle of independence of individual events becomes disturbed, and the system exhibits a trend to self-organization.Intuitively, one could expect that the scaling would be more pronounced under more nonequilibrium conditions. However, during a long time after the discovery of fractals by Mandelbrot no quantitative parameters were suggested for assessing the deviation of a dynamic system from its equilibrium state. Meanwhile, this problem grows out of abstract discussions because of its significance for forecasting the behavior of open natural systems. Apparently, the occurrence of natural hazards depends on the actual state of the given system. It appears that it is for this reason first attempts to estimate in numerical parameters the role of nonequilibrium state were made in seismology [5, 6, 16].In a similar way, the parameter of nonextensivity could serve as a characteristics of the state of drifting sea ice, which would provide information on potential changes in the pack motion and connectedness.The example given in this communication demonstrates the drift dynamics of an individual ice floe during a substantial ice pack fragmentation, which had catastrophic subsequences for the research camp established on the pack. The terminal stage of the floe fracturing was preceded by the steady rise of the degree of nonequilibrium. A similar scenario one could expect for any object of marine engineering, which undergoes irregular impacts from the drifting ice. A response of the sea ice fragmentation on fluctuations in the thermodynamic state of the ice pack (equilibrium/nonequilibrium) contain a certain potential for forecasting icequakes in the Arctic region.


## 6. Conclusion

Multilevel fractal structures, such as the crack-and-lead pattern of fragmented sea ice, are formed in open dynamic systems driven by the external forcing. The statistical analysis of the drift dynamics in the Arctic Ocean in the period of time covering a large-scale ice pack fragmentation allowed us to assess the fluctuations in the thermodynamic state of the mobile sea ice. It was revealed that the most significant sea ice fragmentation characterized by the strongest fluctuations in the ice floe kinetic energy occurred in highly nonequilibrium state of the pack. The deviation of the drifting pack from its equilibrium state could be regarded as a precursor of forthcoming icequakes.



---

*Source: 540407-2012-12-02.xml*