The propagation of large, storm-generated waves through sea ice has so far not been measured, limiting our understanding of how ocean waves break sea ice. Without improved knowledge of ice breakup, we are unable to understand recent changes, or predict future changes, in Arctic and Antarctic sea ice. Here we show that storm-generated ocean waves propagating through Antarctic sea ice are able to transport enough energy to break sea ice hundreds of kilometres from the ice edge. Our results, which are based on concurrent observations at multiple locations, establish that large waves break sea ice much farther from the ice edge than would be predicted by the commonly assumed exponential decay. We observed the wave height decay to be almost linear for large waves--those with a significant wave height greater than three metres--and to be exponential only for small waves. This implies a more prominent role for large ocean waves in sea-ice breakup and retreat than previously thought. We examine the wider relevance of this by comparing observed Antarctic sea-ice edge positions with changes in modelled significant wave heights for the Southern Ocean between 1997 and 2009, and find that the retreat and expansion of the sea-ice edge correlate with mean significant wave height increases and decreases, respectively. This includes capturing the spatial variability in sea-ice trends found in the Ross and Amundsen-Bellingshausen seas. Climate models fail to capture recent changes in sea ice in both polar regions. Our results suggest that the incorporation of explicit or parameterized interactions between ocean waves and sea ice may resolve this problem.

[1] We present a model for wave attenuation in the marginal ice zone (MIZ) based on a two-dimensional (one horizontal and one vertical dimension) multiple floating elastic plate solution in the frequency domain, which is solved exactly using a matched eigenfunction expansion. The only physical parameters that enter the model are length, mass, and elastic stiffness (of which, the latter two depend primarily on thickness) of the ice floes. The model neglects all nonlinear effects as well as floe collisions or ice creep and is therefore most applicable to floes which are large compared to the thickness and to wave conditions which are not extreme. The solution for a given arrangement of floes is fully coherent, and the results are therefore dependent on the exact geometry. We firstly show that this dependence can be removed by averaging over a distribution of floe lengths (we choose the Rayleigh distribution). We then show that after this averaging, the attenuation coefficient is a function of floe number and independent of floe length, provided the floe lengths are sufficiently large. The model predicts an exponential decay of energy, just as is shown experimentally. This enables us to provide explicit values for the attenuation coefficient, as a function of the average floe thickness and wave period. We compare our theoretical predictions of the wave attenuation with measured data and other scattering models. The limited data allows us to conclude that our model is applicable to large floes for short to medium wave periods (6 to 15 seconds). We also derive a floe breaking model, based on our wave attenuation model, which indicates that we are under-predicting the attenuation coefficients at long periods.Citation: Kohout, A. L., and M. H. Meylan (2008), An elastic plate model for wave attenuation and ice floe breaking in the marginal ice zone,

Key Points:• Long waves attenuate with distance traveled in proportion to the inverse of wave period squared • Significant wave height attenuates with distance traveled into the MIZ • Peak periods increase with distance traveled into the MIZ Abstract In situ measurements of ocean surface wave spectra evolution in the Antarctic marginal ice zone are described. Analysis of the measurements shows significant wave heights and peak periods do not vary appreciably in approximately the first 80 km of the ice-covered ocean. Beyond this region, significant wave heights attenuate and peak periods increase. It is shown that attenuation rates are insensitive to amplitudes for long-period waves but increase with increasing amplitude above some critical amplitude for short-period waves. Attenuation rates of the spectral components of the wavefield are calculated. It is shown that attenuation rates decrease with increasing wave period. Further, for long-period waves the decrease is shown to be proportional to the inverse of the period squared. This relationship can be used to efficiently implement wave attenuation through the marginal ice zone in ocean-scale wave models.

[1] The marginal ice zone (MIZ) is the boundary between the open ocean and ice-covered seas, where sea ice is significantly affected by the onslaught of ocean waves. Waves are responsible for the breakup of ice floes and determine the extent of the MIZ and floe size distribution. When the ice cover is highly fragmented, its behavior is qualitatively different from that of pack ice with large floes. Therefore, it is important to incorporate wave-ice interactions into sea ice-ocean models. In order to achieve this goal, two effects are considered: the role of sea ice as a dampener of wave energy and the wave-induced breakup of ice floes. These two processes act in concert to modify the incident wave spectrum and determine the main properties of the MIZ. A simple but novel parameterization for floe breaking is derived by considering alternatively ice as a flexible and rigid material and by using current estimates of ice critical flexural strain and strength. This parameterization is combined with a wave scattering model in a one-dimensional numerical framework to evaluate the floe size distribution and the extent of the MIZ. The model predicts a sharp transition between fragmented sea ice and the central pack, thus providing a natural definition for the MIZ. Reasonable values are found for the extent of the MIZ given realistic initial and boundary conditions. The numerical setting is commensurate with typical iceocean models, with the future implementation into two-dimensional sea ice models in mind.

This paper presents a wave‐in‐ice model calibration study. Data used were collected in the thin ice of the advancing autumn marginal ice zone of the western Arctic Ocean in 2015, where pancake ice was found to be prevalent. Multiple buoys were deployed in seven wave experiments; data from four of these experiments are used in the present study. Wave attenuation coefficients are calculated utilizing wave energy decay between two buoys measuring simultaneously within the ice covered region. Wavenumbers are measured in one of these experiments. Forcing parameters are obtained from simultaneous in‐situ and remote sensing observations, as well as forecast/hindcast models. Cases from three wave experiments are used to calibrate a viscoelastic model for wave attenuation/dispersion in ice cover. The calibration is done by minimizing the difference between modeled and measured complex wavenumber, using a multi‐objective genetic algorithm. The calibrated results are validated using two methods. One is to directly apply the calibrated viscoelastic parameters to one of the wave experiments not used in the calibration and then compare the attenuation from the model with measured data. The other is to use the calibrated viscoelastic model in WAVEWATCH III® over the entire western Beaufort Sea and then compare the wave spectra at two remote sites not used in the calibration. Both validations show reasonable agreement between the model and the measured data. The completed viscoelastic model is believed to be applicable to the fall marginal ice zone dominated by pancake ice.