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.

A large collaborative program has studied the coupled air‐ice‐ocean‐wave processes occurring in the Arctic during the autumn ice advance. The program included a field campaign in the western Arctic during the autumn of 2015, with in situ data collection and both aerial and satellite remote sensing. Many of the analyses have focused on using and improving forecast models. Summarizing and synthesizing the results from a series of separate papers, the overall view is of an Arctic shifting to a more seasonal system. The dramatic increase in open water extent and duration in the autumn means that large surface waves and significant surface heat fluxes are now common. When refreezing finally does occur, it is a highly variable process in space and time. Wind and wave events drive episodic advances and retreats of the ice edge, with associated variations in sea ice formation types (e.g., pancakes, nilas). This variability becomes imprinted on the winter ice cover, which in turn affects the melt season the following year.

Continuum‐based models that describe the propagation of ocean waves in ice‐infested seas are considered, where the surface ocean layer (including ice floes, brash ice, etc.) is modeled by a homogeneous viscoelastic material which causes waves to attenuate as they travel through the medium. Three ice layer models are compared, namely a viscoelastic fluid layer model currently being trialed in the spectral wave model WAVEWATCH III® and two simpler viscoelastic thin beam models. All three models are two dimensional. A comparative analysis shows that one of the beam models provides similar predictions for wave attenuation and wavelength to the viscoelastic fluid model. The three models are calibrated using wave attenuation data recently collected in the Antarctic marginal ice zone as an example. Although agreement with the data is obtained with all three models, several important issues related to the viscoelastic fluid model are identified that raise questions about its suitability to characterize wave attenuation in ice‐covered seas. Viscoelastic beam models appear to provide a more robust parameterization of the phenomenon being modeled, but still remain questionable as a valid characterization of wave‐ice interactions generally.

A theoretical model is used to study wave energy attenuation and directional spreading of ocean wave spectra in the marginal ice zone (MIZ). The MIZ is constructed as an array of tens of thousands of compliant circular ice floes, with randomly selected positions and radii determined by an empirical floe size distribution. Linear potential flow and thin elastic plate theories model the coupled water-ice system. A new method is proposed to solve the time-harmonic multiple scattering problem under a multi-directional incident wave forcing with random phases. It provides a natural framework for tracking the evolution of the directional properties of a wave field through the MIZ. The attenuation and directional spreading are extracted from ensembles of the wave field with respect to realizations of the MIZ and incident forcing randomly generated from prescribed distributions. The averaging procedure is shown to converge rapidly so that only a small number of simulations need to be performed. Far field approximations are investigated, allowing efficiency improvements with negligible loss of accuracy. A case study is conducted for a particular MIZ configuration. Observed exponential attenuation of wave energy through the MIZ is reproduced by the model, while the directional spread is found to grow linearly with distance. Directional spreading is shown to weaken when the wavelength becomes larger than the maximum floe size.

Quantifying the rate of wave attenuation in sea ice is key to understanding trends in the Antarctic marginal ice zone extent. However, a paucity of observations of waves in sea ice limits progress on this front. We deployed 14 waves-in-ice observation systems (WIIOS) on Antarctic sea ice during the Polynyas, Ice Production, and seasonal Evolution in the Ross Sea expedition (PIPERS) in 2017. The WIIOS provide in situ measurement of surface wave characteristics. Two experiments were conducted, one while the ship was inbound and one outbound. The sea ice throughout the experiments generally consisted of pancake and young ice <0.5 m thick. The WIIOS survived a minimum of 4 d and a maximum of 6 weeks. Several large-wave events were captured, with the largest recorded significant wave height over 9 m. We find that the total wave energy measured by the WIIOS generally decays exponentially in the ice and the rate of decay depends on ice concentration.