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.

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.

National audienceValidation of a linear numerical model of wave interactions with floating compliant discs is sought using data obtained from the wave basin experiments reported in Part 1 (Montiel et al. J. Fluid Mech., vol. 723, 2013, pp. 604–628). Comparisons are made for both single-disc tests and the two-disc tests in which wave interactions between discs are observed. The deflection of the disc or discs is separated into the natural modes of vibration in vacuo. The decomposition allows the rigid-body motions and flexural motions to be analysed separately. Rigid-body motions are accurately replicated by the numerical model but, although passable agreement is found, the amplitudes of flexural modes are consistently overestimated. Extensions of the numerical model are used to discount the experimental configuration as a source of the discrepancies. An enhanced viscoelastic model for the discs is also proposed, which results in improved model/data agreement for the flexural motions but cannot account for all of the disagreement

This paper investigates the attenuation and directional spreading of large amplitude waves traveling through pancake ice. Directional spectral density is analyzed from in situ wave buoy data collected during a 3‐day storm event in October 2015 in the Beaufort Sea. Two proxy metrics for wave amplitude obtained from energy density spectra, namely, spectral amplitude and significant wave height, are used to track the waves as they propagate along transects through the array of buoys in the predominantly pancake ice field. Two types of wave buoys are used in the analysis and compared, exhibiting significant differences in the wave energy density and directionality estimates. Although exponential decay is observed predominantly, one of the two buoy types indicates a potential positive correlation between wave energy density and the occurrence of linear wave decay, as opposed to exponential decay, in accord with recent observations in the Antarctic marginal ice zone. Factors affecting the validity of this observation are discussed. An empirical power law with exponent 2.2 is also found to hold between the exponential attenuation coefficient and wave frequency. The directional content of the wave spectrum appears to decrease consistently along the wave transects, confirming that wave energy is being dissipated by the pancake ice as opposed to being scattered by ice cakes.