Jul 03, 2023
The dynamics of unstable waves in sea ice
Scientific Reports volume 13, Article number: 13654 (2023) Cite this article 527 Accesses 7 Altmetric Metrics details Wave and sea ice properties in the Arctic and Southern Oceans are linked by
Scientific Reports volume 13, Article number: 13654 (2023) Cite this article
527 Accesses
7 Altmetric
Metrics details
Wave and sea ice properties in the Arctic and Southern Oceans are linked by feedback mechanisms, therefore the understanding of wave propagation in these regions is essential to model this key component of the Earth climate system. The most striking effect of sea ice is the attenuation of waves at a rate proportional to their frequency. The nonlinear Schrödinger equation (NLS), a fundamental model for ocean waves, describes the full growth-decay cycles of unstable modes, also known as modulational instability (MI). Here, a dissipative NLS (d-NLS) with characteristic sea ice attenuation is used to model the evolution of unstable waves. The MI in sea ice is preserved, however, in its phase-shifted form. The frequency-dependent dissipation breaks the symmetry between the dominant left and right sideband. We anticipate that this work may motivate analogous studies and experiments in wave systems subject to frequency-dependent energy attenuation.
Arctic and Antarctic sea ice play a prominent role in the Earth system by regulating heat and momentum exchanges over large spatial scales1,2,3,4. The sea ice properties are intimately linked to ocean wave properties via feedback mechanisms in the marginal ice zone (MIZ)5,6,7 which around Antarctica, fed by intense Southern Ocean waves all year round8, extends for hundreds of kilometers9,10,11. Rapid evolution of the polar regions driven by climate change12,14,14 have revived and energised research activities in understanding waves properties and feedback in the MIZ7,15, including in the emerging Arctic MIZ16.
Example of Southern Ocean waves (wave height \(\approx 5\) m and peak period \(\approx 12\) s) propagating in a MIZ comprised of small ice floes (1–10 m) as seen from the icebreaker S.A. Agulhas II (beam 21.7 m, for visual reference) on the 24 July 2022 at 59\(^\circ\)S and 1\(^\circ\)E, and schematic of exponential dissipation for a monochromatic wave of unit amplitude propagating from left to right. In the schematic the green line denotes the surface elevation and the red line the wave envelope which undergoes exponential attenuation with distance.
In the MIZ exterior, where the sea ice cover is a mixture of small floes (much shorter than the wavelength) and interstitial frazil ice17,18, as shown in Fig. 1, viscous-like losses have been identified as the main wave attenuation mechanism19,21,21. In the MIZ interior, where floes are larger and comparable to the wavelength, wave attenuation by scattering dominates20. At leading order in wave steepness, i.e. the wave nonlinearity parameter, each wave component attenuates exponentially with distance, see schematic in Fig. 1, and at a frequency-dependent attenuation rate19,22,23. That is, shorter waves are attenuated faster than their longer counterparts. For a comprehensive review of waves in sea ice we refer the reader to Meylan et al.19 and Squire20, and references therein.
Narrowband ocean wave dynamics can be accurately described by the nonlinear Schrödinger equation (NLS). One intriguing dynamical phenomenon, which is responsible for the formation of large-amplitude and coherent waves and has attracted the scientific interest since the late 60s, is the modulation instability (MI)24. In fact , and in contrast to the linear stability analysis of Stoke waves, the complete growth and decay cycles can be described within the NLS framework25. More recently, several studies have been devoted to investigate the effect of wave dissipation on the phase-shifted recurrent MI focusing cycles26,28,29,30,31,31. The latter studies highlight the phase-shifted recurrence in the long-term evolution of nonlinear and unstable waves, when constant, weak, and linear dissipation effects are at play. That said, the NLS can be also adapted to accommodate the influence of sea ice attenuation on the waves by including viscous-like losses as a dissipative term that matches the decay rate of the linear amplitude, as shown by32,33. Within this context, it as been shown that it is important to account for the ice-induced frequency dependency in the attenuation of ocean waves in the MIZ34.
We will focus our attention on the evolution dynamics in infinite depth regime since most of the sea ice processes are relevant in deep-water regimes. In this case the weakly nonlinear spatial evolution of damped hydrodynamic waves can be described by the d-NLS [e.g26,34]:
where \(\psi\) is the complex wave amplitude, which is evolving along the space co-ordinate x, while t denotes the time in the frame of reference moving with at the group speed (\(t=t'-x/c_{g}\) where \(c_g\) is the group speed and \(t'\) the time in the fixed frame of reference), i is the imaginary unit, g the gravitational acceleration, \(k_0\) the wavenumber of the carrier wave, \(\omega _0=\sqrt{gk_0}\) the angular frequency, and \(\mathfrak {D}(\omega )\) the frequency-dependent linear attenuation rate. Eq. (1) is therefore written in the frame of reference moving with at the group speed. We remark that the damping term is introduced using an heuristic approach, and is equivalent to the one introduced in34 with the incorporation of a whole range of dissipative parameters while accounting for the group velocity as reference celerity.
Following19, the viscous-like wave dissipation due to sea ice follows a frequency power law with a particular feature that higher frequency components undergo stronger attenuation:
where we use \(\omega =\omega _0+\delta \omega\) in the formulation of the dissipation coefficient to explicitly highlight the difference in frequency between the carrier wave angular frequency \(\omega _0\) and other spectral components. The power law exponent n depends on the physical mechanisms at play20 (e.g. viscous losses16, basal friction35, scattering21,36, overwash37,38, floe slamming39), and usually falls in the range 2–419. The same applies for the real coefficients \(\alpha _n\), which also depends on the same dynamical and complex interactions19 and, therefore, on effective sea ice properties (e.g. thickness, density, effective viscosity)40. The power law dissipation with exponent \(n=3\), as mentioned in Eq. (2), was found to agree well with field observations and the mathematical modelling19, and will be therefore considered herein. The exponent n is derived from the classical linearised water wave problem when an extra pressure term is added to the free surface dynamic boundary condition to model sea ice. When the pressure term is assumed proportional to the vertical velocity the resulting wave dispersion relation has imaginary part (equivalent to attenuation rate) proportional to \(\omega ^3\)19. Note that the model is equivalent to the one discussed in34, but in this case the dissipative coefficients \(\alpha\) incorporates all the physical parameters and the evolution of the wave envelope is written in the frame of reference moving at group speed. The model can be easily modified to include other forms of frequency-dependent dissipation.
We investigate the classical MI problem of a monochromatic wave train with frequency \(\omega _0\) and initial amplitude \(a_{i}\) subjected to initially small symmetric sideband perturbations \(a_{i}^l\) and \(a_{i}^r\):
Note that no phase shift between the carrier and sidebands is imposed. The initial amplitude of the left and right sideband, \(a_{i}^l\) and \(a_{i}^r\), respectively, is set as 1% of the carrier wave amplitude, i.e. \(a_{i}^l=a_{i}^r=0.01a_{i}\), and the frequency difference as \(\delta \omega =0.1\omega _0\). Note that the chosen \(\delta \omega\) gives the maximum growth rate of the sidebands in ice free waters24,41.
The wave steepness, defined using wavenumber of the carrier wave and the sum of all amplitudes, [e.g.42]:
is initially set to 0.1, also defining \(T_p=2\pi /\omega _0=12\) s. The corresponding wavelength and wavenumber are \(L=225\) m and \(k_0=2.8\times 10^{-2}\) \(\hbox {m}^{-1}\) respectively. Wave period and steepness are representative of intense storm waves at the edge of the Antarctic MIZ18, wave conditions few tens of kilometers away from the sea ice edge and within the MIZ are illustrated in Fig. 1. Besides the conservative case used as a reference, 4 dissipation levels defined by \(\alpha _3\) value which span a large range are analysed (summarise in Table 1). The dissipative cases are arbitrarily defined with respect to each other as low, medium, high and very high.
This work will focus on the physical effects on the unstable wave dynamic when considering frequency-dependent dissipation in the wave modelling to third-order in wave nonlinearity. We particularly show that the shifted MI recurrence is retained, but with a noticeable decrease in the recurrent focusing period with the increase of the sea ice dissipation value. Moreover, the asymmetric damping in wave energy components leads to an instrinsic behaviour of the dominant sidebands in the respective phase space. We anticipate that this study will motivate numerical and experimental studies in several nonlinear wave systems governed by a frequency dependent forced/damped NLS, e.g. optical cavities43, nonlinear optics44, exciton-polariton Bose–Einstein condensates45, plasma physics46, and metamaterials47.
The spatial evolution of the unstable dimensionless, normalized, envelope \(|\Psi |=|\psi |/a_{i}\) in absence of dissipation, i.e. when considering the conservative case, is shown in Fig. 2a. Indeed, this is a well-anticipated and know intrinsic dynamic, which involves recurrent focusing cycles with same wave amplification factor along the dimensionless space co-ordinate \(X=x/L\). Moreover, when considering adjusting the wave packet motion with respect to the group speed, all periodic wave amplification maxima occur at the same dimensionless time \(T=t/T_p\). This corresponds to the pulsation dynamics of a B–Type doubly-periodic breather48 for which the recurrence distance is \(\approx\)160 wavelengths for the chosen wave parameters. Indeed, this is in good agreement with the theoretical value predicted using:
Spatial evolution of the unstable wave envelope in the time domain for increasing dissipation (top to bottom) from none (a) to very high (e). The corresponding dissipation values \(\mathfrak {D}\) for the carrier wave are reported in Table 1. In panels (a), i.e. the conservative B-Type breather, and (b), i.e. the dissipative A-Type breather, the recurrence distances (coloured bars at the bottom) and the recurrence patterns (dashed white lines) are highlighted.
When accounting for a gradual increase of dissipation, as defined by the parameter \(\mathfrak {D}\), phase-shifted focusing recurrence, which is a characteristic A–Type periodic feature26,49, emerges. The evolution of the respective wave packets are shown in Fig 2b–e. This occurs already for very small values in wave dissipation. Moreover, the distance between the cycles diminishes compared to the conservative case, with the shortening being more evident for growing degrees of dissipation for low and medium dissipation. The first pseudo-recurrence cycle is 150 wavelengths for low dissipation and 130 wavelengths for medium dissipation. The picture is more complex in the two most dissipative cases in which the length of the first pseudo-focusing cycle increases marginally, while immediately followed by the next slight wave focusing. We attribute this behaviour to the rapid decrease of wave amplitude with the increase of wave dissipation, and thus the lower wave nonlinearity results in a slower growth-decay cycle for the MI. This can be also traced in Eq. (5) in which the wave steepness appears in the denominator. As such, the rapid energy loss contributes to the increase of the recurrence period.
By using as a boundary condition the numerical results at the end of each pseudo-recurrence cycle we can obtain an updated value for the recurrence length from Eq. (5). Nevertheless, we note that the predictive ability of Eq. (5) deteriorates at each cycle and for increasing degree of dissipation in low and medium dissipation. The formula usually underpredicts the recurrence length compared to the numerical simulations, e.g. the recurrence length of the third cycle is underpredicted by \(\approx 10\%\) in low dissipation and by \(\approx 30\%\) in medium dissipation. For the case of high dissipation, the prediction of the recurrence length in fact improves. Compared to the numerical simulations, the pseudo-recurrence length predicted using Eq. (5) is shorter by \(\approx 10\%\) at the third cycle and longer by \(\approx 20\%\) at the subsequent cycle. This underlines the fact that the wave prediction become unreliable in the case of very high dissipation.
The spatial evolution of the mean wave amplitude \(\langle \Psi \rangle\), computed with respect to the time variable, for the different dissipation levels is summarised in Fig. 3. Note that the vertical axes is in logarithmic scale to highlight deviation from the benchmark exponential decay. In the conservative case, energy is naturally conserved along the propagation in the in space coordinate, i.e. \(\langle \Psi (X)\rangle =1\). With the increase of the level of the exponential dissipation rate, i.e. the negative slope is linear in Fig. 3, the mean wave amplitude at the end of the computational domain is summarised in Table 1. These values are within 1% of the energy level of the carrier wave component subjected to linear attenuation, i.e. \(\exp {[-\mathfrak {D}(\omega _0)x]}\). Therefore, the presence of modulation instability cycles does not alter the overall energy that is carried by waves in the dissipative sea ice domain. The negligible difference contrasts with simulations for random waves in which the nonlinear cases remained significantly more energetic in the sea ice cover compared to corresponding linear cases and the energy decayed less than exponentially [cf.34]. However, it is worth noting that for a continuous spectrum the carrier wave component progressively shifted to longer wave period50 which are also less dissipative. On the other hand, the three wave system has a more symmetric spectrum and the period of the carrier wave component remains constant during the long-term evolution in space. It is also noteworthy that for the lowest dissipation level, at the end of the spatial domain in our numerical simulations (approximately three recurrence cycles), the mean wave amplitude is reduced by less than 1%. Nonetheless, the spatial pattern of the recurrence cycles still becomes phase-shifted, see Fig. 2b.
Spatial evolution of the mean wave amplitude in logarithmic scale for the dissipative cases, from low (yellow) to very high (blue), as listed in Table 1.
In each of the simulations performed, we tracked the maximum amplification in space, irrespective of time, i.e. \(\max |\Psi |\) shown as thick red line in Fig. 4, and compared to the conservative case, depicted in thin red line in Fig. 4. As shown in the surface plots in Fig. 2, the shortening of the recurrence cycle is evident even for low dissipation and as a result the maxima occur at location different from the ones for the conservative case but amplification is almost unaltered. Strikingly, in the medium dissipation case, see Fig. 4b, the spatial frequency of local maxima is almost doubled compared to the conservative case, i.e. the locations in which rogue waves occur, i.e. \(\max |\Psi |\ge 2\), is increased for medium level of dissipation. On the first cycle of recurrence the maximum amplification is only slightly diminished and at the end of the numerical domain (after 5 cycles of phase-shifted recurrence) the maximum amplification is \(\approx 2\). For greater levels of dissipation, that is, high and very high and as in Fig. 4c–d, the shortening of recurrence cycle is even more pronounced. In the high dissipative case, only the first two cycles of phase-shifted recurrences have amplification larger than one, since the first cycle has an amplification \(\approx 1.5\). For very high dissipation, the damping is predominant and the amplitude amplification never exceeds one.
Spatial evolution of the wave amplitudes for low to very high dissipation (top to bottom) for the carrier peak energy (black), left (blue) and right sideband (green), and the maximum amplitude (red). Thin lines denote the conservative case. The dash-dotted line denotes linear decay of the carrier wave amplitude component.
The evolution of the carrier wave (\(A=a/a_{i}\)) and first left and right order sideband \(\bigg (A^l=a^{l}/a_{i}\), \(A^r=a^{r}/a_{i}\bigg )\) is also shown in Fig. 4. All the cases show the energy exchange from the carrier to the sidebands during the growing phase of the MI cycle and the opposite in the decaying phase. The frequency dependent dissipation would imply that the left (right) sideband undergoes 30% lower (greater) dissipation than the carrier. That said, the difference is small but noticeable only in the high dissipation case difference, in all the other cases the difference is negligible, i.e. the blue and green lines overlap. For low dissipation, only an accelerated cycle in the dynamics of the carrier and the first order sidebands is observed but their amplitude remain almost unaltered compared to the conservative case (thin lines in Fig. 4). For medium dissipation the sidebands never return to the initial amplitude level but at the end of each decay cycle are more energetic, i.e. each local minimum has a higher value of the previous one. For high dissipation rate a similar behaviour is noted, but the higher dissipation rate means that the sidebands are eventually damped. In the high dissipative case, the energy exchange between the carrier and the sidebands is almost suppressed after the first cycle, and the carrier band decays exponentially (the thick black line overlaps the dash dotted line in Fig. 4d).
To study the dynamics of the first order left sideband we construct its phase space diagram as51:
where \(\Delta \phi ^l\) denotes the phase difference between the carrier wave mode and the dominant left sideband, i.e. \(\Delta \phi ^{l} = \phi -\phi ^{l}\). The phase space diagram of the right sideband is obtained in similar manner.
Figure 5 reports the conservative case (in black) in which the subsequent cycles of the MI repeat identically in time, i.e. they repeat the same track, and are confined in the right side of the phase space diagram (right and left sideband behave exactly in the same way). For the selected initial condition, the trajectory in the phase space diagram is confined within the separatrix, and recurrence cycle are in-phase. This is consistent with the expected B–Type recurrence. For low to medium dissipation (Fig.5a–d), trajectories switch to an outer trajectory (shaped like an eight) and result in phase-shifted recurrence cycles, i.e. A–Type recurrence. Subsequent cycles for the left (right) sideband move clockwise (anti-clockwise) and at a higher degree for low dissipation. For frequency-independent attenuation no rotation of the main axes in the phase space diagram is observed and the left and right sideband behaves in the same manner, see dotted line in Fig. 5a–d. For high to very-high dissipation the switch to an eight-like shape of the trajectory that spans the left and right hand side of the phase space is also observed, but subsequent loops rapidly degenerate into spiralling cycles due to the substantial attenuation, see Fig. 5e–h.
Phase space diagram of the left (a, c, e, g) and right sideband (b, d, f, h) for low to very high dissipation (top to bottom). The black lines depict the conservative case and the dotted line the frequency independent dissipation. Note that the axis’ limits change in each panel.
The different behaviour of the frequency-dependent and frequency-independent and constant attenuation is clearly shown by examining the phase difference between the sidebands and the carrier, see Fig. 6. For the conservative case, i.e. the case of zero dissipation, the phase difference is the same for the dominant left and right sideband and confined between \(-\pi /2\) and \(\pi /2\), see continuous black line in Fig. 6. For the frequency-dependent dissipation left and right sideband span angles between \(-\pi\) to \(\pi\) and behave differently, particularly, for low dissipation levels (Fig. 6a), with the right sideband lagging behind the left one in the rotation. It is worth noting that for frequency-independent attenuation, as denoted by dotted line in Fig. 6, the right and left sideband behave the same and have an intermediate behaviour between the sidebands in the frequency dependent case. Differences between the right and left sidebands in the dissipative case tend to disappear for the higher attenuation levels, see Fig. 6c–d.
Phase-difference for low to very high dissipation (top to bottom) for the left (blue) and right sideband (green). The black line shows the conservative case and the dotted line the frequency independent dissipation.
The classical dynamics of the modulational instability is studied in the presence of a frequency-dependent attenuation, this is in contrast with previous works in which forcing and damping were considered frequency independent, e.g.26,29,31,51. Nevertheless, similarly to previous works, a small dissipation is capable of altering the recurrence cycle from a B–Type (non-shifted) to A–Type (shifted) breather-type evolution. The most important effect of the differential attenuation rate modifies the dynamics of the left/right sidebands and there is a lag between the two. However, compared to frequency-independent attenuation, differences in energy and recurrence pattern appear unaffected. The large initial nonlinearity of the system triggers pseudo-recurrence cycles with large amplification of the initial sea state, particularly, for low and medium dissipation. Interestingly, for medium dissipation rates the shortening of the cycle means that there are more instances in the domain with amplifications greater than two, compared to the conservative case. For comparison, random sea states, see Ref.34, reverted towards Gaussian statistics in presence of dissipation.
Despite the complex dynamics in the nonlinear evolution, the total energy decays in agreement with linear theory, this contrasts the random waves (cf.34) in which for large dissipation make a considerable difference between the linear and nonlinear case. This difference is attributed to the shift of the peak frequency in a continuous wave spectrum, whereas for the classical MI case the carrier wave frequency remains unshifted.
Ultimately, the results show that an energetic swell system propagating in the Arctic and Antarctic sea ice can generate exceptionally large waves where none were expected, especially, when dissipation is mild to medium, and could cause unexpected hazardous conditions for vessels operating in the marginal ice zone. The present research highlights the need to further explore the dynamics of high-order NLS-type systems for broadband processes and in the presence of a frequency dependent damping, which could be applied to other physical systems, e.g. nonlinear optics, Bose–Einstein condensates, and metamaterials, in addition to hydrodynamics waves.
The Benjamin–Feir Instability index52,53 is
The instability condition for the framework evolving in time is \(0<K<2\sqrt{2}k^2a\), so the one in space becomes \(0<\Omega <\sqrt{2}ka\omega\), since \({\Omega }/{\omega }={K}/{(2k)}\). That is, for the same initial long-wave perturbation and assuming a typical exponential energy attenuation under the sea ice the following holds:
For the chosen wave properties the sea ice Benjamin–Feir index at the ice edge \(BFI^{SI}\) is 1.41, i.e. \(BFI^{SI}(x=0)=1.41\), for all considered \(\mathfrak {D}\) values. This is a sufficient condition to trigger growth-decay cycles, typical to the MI54,55.
The theoretical recurrence period for the growth and decay cycle of the unstable wave train is given by56. In the framework of a wave propagating we apply the canonical transformation using the group velocity, to obtain the recurrence length (Eq. 5)). Compared to56, the ratio between the carrier and the sidebands is modified to account for possible asymmetry in their amplitude. For the chosen boundary conditions the recurrence cycle corresponds to 155 wavelengths using the carrier wavelength L as normalising factor.
Various level of dissipation on unstable wave dynamics are analysed in this work: from the conservative case (used as benchmark) to very high attenuation values. Recalling that \(n=3\) is used as the power law exponent for the dissipation, the corresponding scaling parameter \(\alpha _3\) are reported in Table 1. The linear attenuation rate for the dominant wave frequency is also given, i.e. \(\mathfrak {D}(\omega _0)=\alpha _3\omega _0^3\). The dissipation lengthscale can be expressed by the ratio \(k/\mathfrak {D}(\omega _0)\) and ranges from infinity for the conservative case, i.e. \(\mathfrak {D}=0\), to approximately 1400 wavelengths for the very high dissipative case.
It is worth noting that dissipation scales \(\omega ^3\) in our model, therefore the carrier wave component and the sidebands are subjected to different attenuation rates. In particular, by using \(\delta \omega = 0.1\omega _0\), we obtain that the left sideband attenuates at a rate \(0.73\mathfrak {D}(\omega _0)\) and the right sideband at a rate \(1.33\mathfrak {D}(\omega _0)\). Therefore, for sidebands with the maximum growth rate and cubic attenuation with respect with frequency, the difference in attenuation rate between the left and right sideband is almost twice stronger, when considering \(\mathfrak {D}(\omega _0+\delta \omega )/\mathfrak {D}(\omega _0-\delta \omega )=1.83\).
The d-NLS (1) is solved numerically, by advancing \(\psi (x,t)\) in space using the fourth order Runge-Kutta method. A time-periodic domain which encompasses a full cycle of modulation, including 10 wave periods for the chosen \(\delta \omega\), is used. As such, this makes an efficient and accurate computation of time derivatives in the Fourier domain possible. The frequency-dependent dissipative term can be computed in Fourier space in a straightforward manner using the Fourier and inverse Fourier transform (\(\mathcal {F}\) and \(\mathcal {F}^{-1}\) respectively) as the following:
The time domain is discretized using \(2^6\) elements and the resulting timestep is \(\delta t=\)1.875 s. Thus, \(T/\delta t=6.4\) and the unstable wave envelope is propagated over 100 km in space using \(\delta x = 1\) m (in dimensionless form the computational domain corresponds to \(X=x/L=445\) wavelengths). Temporal and spatial resolutions guarantee numerical stability.
All data generated or analysed during this study are included in this published article.
Vihma, T. Effects of arctic sea ice decline on weather and climate: A review. Surv. Geophys. 35(5), 1175–1214 (2014).
Article ADS Google Scholar
Massom, R. A. & Stammerjohn, S. E. Antarctic sea ice change and variability - physical and ecological implications. Polar Sci. 4(2), 149–186. https://doi.org/10.1016/j.polar.2010.05.001 (2010).
Article ADS Google Scholar
Vichi, M. et al. Effects of an explosive polar cyclone crossing the Antarctic marginal ice zone. Geophys. Res. Lett. 46(11), 5948–5958. https://doi.org/10.1029/2019GL082457 (2019).
Article ADS Google Scholar
Hepworth, E., Messori, G. & Vichi, M. Association between extreme atmospheric anomalies over antarctic sea ice, southern ocean polar cyclones and atmospheric rivers. J. Geophys. Res. Atmos. 127(7), 036121. https://doi.org/10.1029/2021JD036121 (2021).
Article Google Scholar
Wadhams, P. In The Seasonal Ice Zone (ed. Untersteiner, N.) 825–991 (Springer, 1986).
Eayrs, C. et al. Understanding the seasonal cycle of antarctic sea ice extent in the context of longer-term variability. Rev. Geophys. 57(3), 1037–1064. https://doi.org/10.1029/2018RG000631 (2019).
Article ADS Google Scholar
Squire, V. A. Marginal ice zone dynamics. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 380(2235), 20210266. https://doi.org/10.1098/rsta.2021.0266 (2022).
Article Google Scholar
Derkani, M. H. et al. Wind, waves, and surface currents in the southern ocean: Observations from the antarctic circumnavigation expedition. Earth Syst. Sci. Data 13(3), 1189–1209. https://doi.org/10.5194/essd-13-1189-2021 (2021).
Article ADS Google Scholar
Stopa, J. E., Sutherland, P. & Ardhuin, F. Strong and highly variable push of ocean waves on southern ocean sea ice. Proc. Natl. Acad. Sci. 115(23), 5861–5865. https://doi.org/10.1073/pnas.1802011115 (2018).
Article ADS CAS PubMed PubMed Central Google Scholar
Brouwer, J. et al. Altimetric observation of wave attenuation through the Antarctic marginal ice zone using ICESat-2. Cryosph. Discuss. 2021, 1–40. https://doi.org/10.5194/tc-2021-367 (2021).
Article Google Scholar
Womack, A., Vichi, M., Alberello, A. & Toffoli, A. Atmospheric drivers of a winter-to-spring lagrangian sea-ice drift in the eastern antarctic marginal ice zone. J. Glaciol. 68, 999–1013. https://doi.org/10.1017/jog.2022.14 (2022).
Article Google Scholar
Stammerjohn, S., Massom, R., Rind, D. & Martinson, D. Regions of rapid sea ice change: An inter-hemispheric seasonal comparison. Geophys. Res. Lett.https://doi.org/10.1029/2012GL050874 (2012).
Article Google Scholar
Turner, J. & Comiso, J. Solve Antarctica’s sea-ice puzzle. Nature 547(7663), 275–277. https://doi.org/10.1038/547275a (2017).
Article ADS CAS PubMed Google Scholar
Eayrs, C., Li, X., Raphael, M. N. & Holland, D. M. Rapid decline in Antarctic sea ice in recent years hints at future change. Nat. Geosci. 14(7), 460–464. https://doi.org/10.1038/s41561-021-00768-3 (2021).
Article ADS CAS Google Scholar
Bennetts, L. G., Bitz, C. M., Feltham, D. L., Kohout, A. L. & Meylan, M. H. Theory, modelling and observations of marginal ice zone dynamics: Multidisciplinary perspectives and outlooks. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 380(2235), 20210265. https://doi.org/10.1098/rsta.2021.0265 (2022).
Article Google Scholar
Thomson, J. et al. Overview of the arctic sea state and boundary layer physics program. J. Geophys. Res. Oceans 123(12), 8674–8687. https://doi.org/10.1002/2018JC013766 (2018).
Article ADS Google Scholar
Alberello, A. et al. Brief communication: Pancake ice floe size distribution during the winter expansion of the antarctic marginal ice zone. Cryosphere 13(1), 41–48. https://doi.org/10.5194/tc-13-41-2019 (2019).
Article ADS Google Scholar
Alberello, A. et al. Three-dimensional imaging of waves and floes in the marginal ice zone during a cyclone. Nat. Commun. 13, 4590 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Meylan, M. H. et al. Dispersion relations, power laws, and energy loss for waves in the marginal ice zone. J. Geophys. Res. Oceans 123(5), 3322–3335. https://doi.org/10.1002/2018JC013776 (2018).
Article ADS Google Scholar
Squire, V. A. Ocean wave interactions with sea ice: A reappraisal. Annu. Rev. Fluid Mech. 52(1), 37–60. https://doi.org/10.1146/annurev-fluid-010719-060301 (2020).
Article ADS MathSciNet MATH Google Scholar
Toffoli, A., Pitt, J. P. A., Alberello, A. & Bennetts, L. G. Modelling attenuation of irregular wave fields by artificial ice floes in the laboratory. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 380(2235), 20210255. https://doi.org/10.1098/rsta.2021.0255 (2022).
Article ADS CAS Google Scholar
Kohout, A. L. et al. Observations of exponential wave attenuation in Antarctic sea ice during the PIPERS campaign. Ann. Glaciol. 61(82), 196–209. https://doi.org/10.1017/aog.2020.36 (2020).
Article ADS Google Scholar
Waseda, T. et al. Observation of anomalous spectral downshifting of waves in the okhotsk sea marginal ice zone. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 380(2235), 20210256. https://doi.org/10.1098/rsta.2021.0256 (2022).
Article Google Scholar
Benjamin, T. B. & Feir, J. E. The disintegration of wave trains on deep water part 1. theory. J. Fluid Mech. 27(3), 417–430. https://doi.org/10.1017/S002211206700045X (1967).
Article ADS MATH Google Scholar
Yuen, H. C. & Lake, B. M. Nonlinear dynamics of deep-water gravity waves. Adv. Appl. Mech. 22, 67–229 (1982).
Article MathSciNet MATH Google Scholar
Kimmoun, O. et al. Modulation instability and phase-shifted fermi-pasta-ulam recurrence. Sci. Rep. 6, 28516 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Mussot, A. et al. Fibre multi-wave mixing combs reveal the broken symmetry of fermi-pasta-ulam recurrence. Nat. Photon. 12(5), 303–308. https://doi.org/10.1038/s41566-018-0136-1 (2018).
Article ADS CAS Google Scholar
Pierangeli, D. et al. Observation of fermi-pasta-ulam-tsingou recurrence and its exact dynamics. Phys. Rev. X 8, 041017. https://doi.org/10.1103/PhysRevX.8.041017 (2018).
Article Google Scholar
Eeltink, D. et al. Separatrix crossing and symmetry breaking in NLSE-like systems due to forcing and damping. Nonlinear Dyn. 102(4), 2385–2398 (2020).
Article CAS PubMed PubMed Central Google Scholar
Coppini, F., Grinevich, P. & Santini, P. Effect of a small loss or gain in the periodic nonlinear Schrödinger anomalous wave dynamics. Phys. Rev. E 101(3), 032204 (2020).
Article ADS MathSciNet CAS PubMed Google Scholar
Zaug, C. R. & Carter, J. D. Dissipative models of swell propagation across the pacific. Stud. Appl. Math. 147(4), 1519–1537. https://doi.org/10.1111/sapm.12423 (2021).
Article MathSciNet MATH Google Scholar
Wu, G., Liu, Y. & Yue, D. K. P. A note on stabilizing the Benjamin–Feir instability. J. Fluid Mech. 556, 45–54. https://doi.org/10.1017/S0022112005008293 (2006).
Article ADS MathSciNet MATH Google Scholar
Dias, F., Dyachenko, A. I. & Zakharov, V. E. Theory of weakly damped free-surface flows: A new formulation based on potential flow solutions. Phys. Lett. A 372(8), 1297–1302. https://doi.org/10.1016/j.physleta.2007.09.027 (2008).
Article ADS CAS MATH Google Scholar
Alberello, A. & Părău, E. A dissipative nonlinear Schrödinger model for wave propagation in the marginal ice zone. Phys. Fluidshttps://doi.org/10.1063/5.0089866 (2022).
Article Google Scholar
Boutin, G. et al. Floe size effect on wave-ice interactions: Possible effects, implementation in wave model, and evaluation. J. Geophys. Res. Oceans 123(7), 4779–4805. https://doi.org/10.1029/2017JC013622 (2018).
Article ADS Google Scholar
Squire, V. A., Dugan, J. P., Wadhams, P., Rottier, P. J. & Liu, A. K. Of ocean waves and sea ice. Annu. Rev. Fluid Mech. 27(1), 115–168. https://doi.org/10.1146/annurev.fl.27.010195.000555 (1995).
Article ADS MathSciNet Google Scholar
Bennetts, L. G. et al. An idealised experimental model of ocean surface wave transmission by an ice floe. Ocean Modell. 96, 85–92. https://doi.org/10.1016/j.ocemod.2015.03.001 (2015).
Article ADS Google Scholar
Toffoli, A. et al. Sea ice floes dissipate the energy of steep ocean waves. Geophys. Res. Lett. 42(20), 8547–8554. https://doi.org/10.1002/2015GL065937 (2015).
Article ADS Google Scholar
Nelli, F., Bennetts, L., Skene, D. & Toffoli, A. Water wave transmission and energy dissipation by a floating plate in the presence of overwash. J. Fluid Mech. 889, 19. https://doi.org/10.1017/jfm.2020.75 (2020).
Article ADS MathSciNet Google Scholar
Rabault, J., Sutherland, G., Jensen, A., Christensen, K. H. & Marchenko, A. Experiments on wave propagation in grease ice: Combined wave gauges and particle image velocimetry measurements. J. Fluid Mech. 864, 876–898. https://doi.org/10.1017/jfm.2019.16 (2019).
Article ADS CAS Google Scholar
Tulin, M. P. & Waseda, T. Laboratory observations of wave group evolution, including breaking effects. J. Fluid Mech. 378, 197–232. https://doi.org/10.1017/S0022112098003255 (1999).
Article ADS Google Scholar
Alberello, A. & Iafrati, A. The velocity field underneath a breaking rogue wave: Laboratory experiments versus numerical simulations. Fluids 4(2), 68. https://doi.org/10.3390/fluids4020068 (2019).
Article ADS CAS Google Scholar
Herr, T. et al. Temporal solitons in optical microresonators. Nat. Photon. 8(2), 145–152 (2014).
Article ADS CAS Google Scholar
Naveau, C. et al. Experimental characterization of recurrences and separatrix crossing in modulational instability. Opt. Lett. 44(22), 5426–5429. https://doi.org/10.1364/OL.44.005426 (2019).
Article ADS PubMed Google Scholar
Deng, H., Haug, H. & Yamamoto, Y. Exciton-polariton Bose–Einstein condensation. Rev. Mod. Phys. 82, 1489–1537. https://doi.org/10.1103/RevModPhys.82.1489 (2010).
Article ADS CAS Google Scholar
Shukla, P. K. & Eliasson, B. Nonlinear aspects of quantum plasma physics. Phys. Usp. 53(1), 51. https://doi.org/10.3367/UFNe.0180.201001b.0055 (2010).
Article ADS CAS Google Scholar
Demiquel, A., Achilleos, V., Theocharis, G. & Tournat, V. Modulation instability in nonlinear flexible mechanical metamaterials. Phys. Rev. E 107, 054212. https://doi.org/10.1103/PhysRevE.107.054212 (2023).
Article ADS MathSciNet CAS PubMed Google Scholar
Conforti, M., Mussot, A., Kudlinski, A., Trillo, S. & Akhmediev, N. Doubly periodic solutions of the focusing nonlinear Schrödinger equation: Recurrence, period doubling, and amplification outside the conventional modulation-instability band. Phys. Rev. A 101(2), 023843 (2020).
Article ADS MathSciNet CAS Google Scholar
Akhmediev, N. N. et al. Nonlinear Pulses and Beams (Springer, 1997).
MATH Google Scholar
Hasselmann, K. F., Barnett, T. P., Bouws, E., Carlson, H. C., Cartwright, D. E., Enke, K., Ewing, J. A., Gienapp, H., Hasselmann, D. E., Kruseman, P., Meerburg, A., Müller, P. M., Olbers, D. J., Richter, K., Sell, W. & Walden, H. Measurements of wind-wave growth and swell decay during the joint north sea wave project (jonswap). (1973)
Trillo, S. & Wabnitz, S. Dynamics of the nonlinear modulational instability in optical fibers. Opt. Lett. 16(13), 986–988. https://doi.org/10.1364/OL.16.000986 (1991).
Article ADS CAS PubMed Google Scholar
Janssen, P. A. Nonlinear four-wave interactions and freak waves. J. Phys. Oceanogr. 33(4), 863–884 (2003).
2.0.CO;2" data-track-action="article reference" href="https://doi.org/10.1175%2F1520-0485%282003%2933%3C863%3ANFIAFW%3E2.0.CO%3B2" aria-label="Article reference 52" data-doi="10.1175/1520-0485(2003)332.0.CO;2">Article ADS MathSciNet Google Scholar
Serio, M., Onorato, M., Osborne, A. R. & Janssen, P. A. On the computation of the Benjamin–Feir index. Il Nuovo Cimento C 28(6), 893–903 (2005).
Google Scholar
Tulin, M. P. & Waseda, T. Laboratory observations of wave group evolution, including breaking effects. J. Fluid Mech. 378, 197–232 (1999).
Article ADS Google Scholar
Onorato, M. et al. Extreme waves, modulational instability and second order theory: Wave flume experiments on irregular waves. Eur. J. Mech. B. Fluids 25(5), 586–601. https://doi.org/10.1016/j.euromechflu.2006.01.002 (2006) (Rogue waves).
Article ADS MathSciNet MATH Google Scholar
Stiassnie, M. & Kroszynski, U. I. Long-time evolution of an unstable water-wave train. J. Fluid Mech. 116, 207–225. https://doi.org/10.1017/S0022112082000433 (1982).
Article ADS MathSciNet MATH Google Scholar
Download references
The research presented in this paper was carried out on the High Performance Computing Cluster supported by the Research and Specialist Computing Support service at the University of East Anglia. AA acknowledges Dr. Davide Proment and Dr. Alberto Villois for interesting discussions. AC acknowledges support from Kyoto University’s Hakubi Center for Advanced Research.
School of Mathematics, University of East Anglia, Norwich, NR4 7TJ, UK
Alberto Alberello & Emilian Părău
School of Civil Engineering, The University of Sydney, Sydney, NSW, 2006, Australia
Amin Chabchoub
Hakubi Center for Advanced Research, Disaster Prevention Research Institute, Kyoto University, Kyoto, 606-8501, Japan
Amin Chabchoub
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
A.A. designed the research. A.A., E.P. and A.C. conducted the analysis and contributed to the writing of the article.
Correspondence to Alberto Alberello.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Alberello, A., Părău, E. & Chabchoub, A. The dynamics of unstable waves in sea ice. Sci Rep 13, 13654 (2023). https://doi.org/10.1038/s41598-023-40696-3
Download citation
Received: 05 July 2023
Accepted: 16 August 2023
Published: 22 August 2023
DOI: https://doi.org/10.1038/s41598-023-40696-3
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.