Abstract
A model of ice floe break-up under ocean wave forcing in the marginal ice zone (MIZ) is proposed to investigate how floe size distribution (FSD) evolves under repeated wave break-up events. A three-dimensional linear model of ocean wave scattering by a finite array of compliant circular ice floes is coupled to a flexural failure model, which breaks a floe into two floes provided the two-dimensional stress field satisfies a break-up criterion. A closed-feedback loop algorithm is devised, which (i) solves the wave-scattering problem for a given FSD under time-harmonic plane wave forcing, (ii) computes the stress field in all the floes, (iii) fractures the floes satisfying the break-up criterion, and (iv) generates an updated FSD, initializing the geometry for the next iteration of the loop. The FSD after 50 break-up events is unimodal and near normal, or bimodal, suggesting waves alone do not govern the power law observed in some field studies. Multiple scattering is found to enhance break-up for long waves and thin ice, but to reduce break-up for short waves and thick ice. A break-up front marches forward in the latter regime, as wave-induced fracture weakens the ice cover, allowing waves to travel deeper into the MIZ.