Ice Floe Breaking in Contemporary Third Generation Operational Wave Models

The dynamical zone observed at the edge of the consolidated ice area where are found the wave-fractured floes (i.e. marginal ice zone or MIZ) has become an important topic in ocean modeling. As both operational and climate ocean models now seek to reproduce the complex atmosphere-ice-ocean system with realistic coupling processes, many theoretical and numerical studies have focused on understanding and modeling this zone. Few attempts have been made to embed wave-ice interactions specific to the MIZ within a two-dimensional model, giving the possibility to calculate both the attenuation of surface waves by sea ice and the concomitant breaking of the sea ice-cover into smaller floes. One of the first challenges consists in improving the parameterization of wave-ice dynamics in contemporary third generation operational wave models. A simple waves-in-ice model (WIM) similar to the one proposed by Williams et al. (2013a,b) was implemented in WAVEWATCH III. This WIM considers ice floes as floating elastic plates and predicts the dimensionless attenuation coefficient by the use of a lookup-table-based, wave scattering scheme. As in Dumont et al. (2011), the different frequencies are treated individually and floe breaking occurs for a particular frequency when the expected wave amplitude exceeds the allowed strain amplitude, which considers ice floes properties and wavelength in ice field. The model is here further refined and tested in idealized two-dimensional cases, giving preliminary results of the performance and sensitivity of the parameterization to initial wave and ice conditions. The effects of the wave-ice coupling over the incident wave spectrum are analyzed as well as the resulting floe size distribution. The model gives prognostic values of the lateral extent of the marginal ice zone with maximum ice floe diameter that progressively increases with distance from the ice edge.