This paper presents theoretical foundations, numerical implementation and examples of application of the two-dimensional Discrete-Element bonded-particle Sea Ice model – DESIgn. In the model, sea ice is represented as an assemblage of objects of two types: disk-shaped "grains" and semi-elastic bonds connecting them. Grains move on the sea surface under the influence of forces from the atmosphere and the ocean, as well as interactions with surrounding grains through direct contact (Hertzian contact mechanics) and/or through bonds. The model has an experimental option of taking into account quasi-three-dimensional effects related to the space- and time-varying curvature of the sea surface, thus enabling simulation of ice breaking due to stresses resulting from bending moments associated with surface waves. Examples of the model's application to simple sea ice deformation and breaking problems are presented, with an analysis of the influence of the basic model parameters ("microscopic" properties of grains and bonds) on the large-scale response of the modeled material. The model is written as a toolbox suitable for usage with the open-source numerical library LIGGGHTS. The code, together with full technical documentation and example input files, is freely available with this paper and on the Internet.

This paper presents the analysis of floe-size distribution (FSD) data obtained in laboratory experiments of ice breaking by waves. The experiments, performed at the Large Ice Model Basin (LIMB) of the Hamburg Ship Model Basin (Hamburgische Schiffbau-Versuchsanstalt, HSVA), consisted of a number of tests in which an initially continuous, uniform ice sheet was broken by regular waves with prescribed characteristics. The floes' characteristics (surface area; minor and major axis, and orientation of equivalent ellipse) were obtained from digital images of the ice sheets after five tests. The analysis shows that although the floe sizes cover a wide range of values (up to 5 orders of magnitude in the case of floe surface area), their probability density functions (PDFs) do not have heavy tails, but exhibit a clear cut-off at large floe sizes. Moreover, the PDFs have a maximum that can be attributed to wave-induced flexural strain, producing preferred floe sizes. It is demonstrated that the observed FSD data can be described by theoretical PDFs expressed as a weighted sum of two components, a tapered power law and a Gaussian, reflecting multiple fracture mechanisms contributing to the FSD as it evolves in time. The results are discussed in the context of theoretical and numerical research on fragmentation of sea ice and other brittle materials.

Among many mechanisms potentially contributing to wave energy attenuation in sea ice are wave‐induced ice floe collisions. At present, little is known about collision patterns and their phase‐averaged effects under different combinations of sea ice properties (ice thickness, floe size, etc.) and wave forcing (wavelength and steepness). The existing parameterizations of collision‐related effects are therefore based on several simplifying, unverified assumptions. In this work, wave‐induced motion and collisions of ice floes are analyzed numerically with a model based on momentum equations for an arbitrary number of floes, with source terms computed by integrating local forcing (wave‐induced dynamic pressure, surface drag, etc.) over the surface area/volume of each floe. It is shown that this simple model, with prescribed wave forcing (i.e., no wave‐ice interactions), is capable of reproducing observed surge amplitudes up to floe sizes comparable with wavelength. A full Hertzian contact model is used instead of a simple hard‐disk algorithm, which makes the model suitable for simulating both rapid collisions and prolonged contact between floes. The model equations are used to formulate heuristic collision criteria based on relative floe size, ice concentration, and wave steepness. The model is then run for different combinations of those three parameters, together with different restitution and drag coefficients, in order to analyze possible motion/collision patterns within the multidimensional parameter space, and phase‐averaged effects of collisions: kinetic and contact stress, granular temperature, and work done by forces acting on the ice.

Sea-ice floe-size distribution (FSD) in ice-pack covered seas influences many aspects of ocean-atmosphere interactions. However, data concerning FSD in the polar oceans are still sparse and processes shaping the observed FSD properties are poorly understood. Typically, power-law FSDs are assumed although no feasible explanation has been provided neither for this one nor for other properties of the observed distributions. Consequently, no model exists capable of predicting FSD parameters in any particular situation. Here I show that the observed FSDs can be well represented by a truncated Pareto distribution P(x)=x$-${}1$-${}$\alpha${}\,exp[(1$-${}$\alpha${})/x], which is an emergent property of a certain group of multiplicative stochastic systems, described by the generalized Lotka-Volterra (GLV) equation. Building upon this recognition, a possibility of developing a simple agent-based GLV-type sea-ice model is considered. Contrary to simple power-law FSDs, GLV gives consistent estimates of the total floe perimeter, as well as floe-area distribution in agreement with observations.

In this paper, a coupled sea ice–wave model is developed and used to analyze wave-induced stress and breaking in sea ice for a range of wave and ice conditions. The sea ice module is a discrete-element bonded-particle model, in which ice is represented as cuboid grains floating on the water surface that can be connected to their neighbors by elastic joints. The joints may break if instantaneous stresses acting on them exceed their strength. The wave module is based on an open-source version of the Non-Hydrostatic WAVE model (NHWAVE). The two modules are coupled with proper boundary conditions for pressure and velocity, exchanged at every wave model time step. In the present version, the model operates in two dimensions (one vertical and one horizontal) and is suitable for simulating compact ice in which heave and pitch motion dominates over surge. In a series of simulations with varying sea ice properties and incoming wavelength it is shown that wave-induced stress reaches maximum values at a certain distance from the ice edge. The value of maximum stress depends on both ice properties and characteristics of incoming waves, but, crucially for ice breaking, the location at which the maximum occurs does not change with the incoming wavelength. Consequently, both regular and random (Jonswap spectrum) waves break the ice into floes with almost identical sizes. The width of the zone of broken ice depends on ice strength and wave attenuation rates in the ice.