pith. sign in

arxiv: 1907.03490 · v1 · pith:4R25FNUAnew · submitted 2019-07-08 · ⚛️ physics.flu-dyn · physics.ao-ph· physics.data-an

Single-spectrum prediction of kurtosis of water waves in a non-conservative model

D. Eeltink , A. Armaroli , Y.M. Ducimeti\`ere , J. Kasparian , M. Brunetti This is my paper

Pith reviewed 2026-05-25 01:07 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn physics.ao-phphysics.data-an
keywords water waveskurtosisspectral bandwidthrogue wavesnonlinear Schrödinger equationBenjamin-Feir indexforced-damped modelwave statistics
0
0 comments X

The pith

A quadratic relation between kurtosis and spectral bandwidth holds robustly for water waves under wind forcing and damping.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper models one-directional water waves with random initial conditions using a forced-damped higher-order nonlinear Schrödinger equation. It tracks how statistics change during and after a sudden wind episode, with wave action increasing then decreasing. The kurtosis of wave heights, a key indicator for rogue wave probability, shows a quadratic dependence on the Benjamin-Feir index that varies with forcing details. In contrast, a consistent quadratic relation emerges between kurtosis and the spectral bandwidth across different forcing and damping strengths. This suggests that measuring the spectrum alone could indicate the chance of extreme waves in a sea state.

Core claim

In simulations of the forced-damped higher-order NLS equation, the kurtosis of the wave height distribution maintains a simple quadratic dependence on the spectral bandwidth that does not depend on the specific forcing and damping coefficients, unlike its relation to the Benjamin-Feir index. This relation permits estimation of rogue wave likelihood from a single measured spectrum. In addition, the evolution of kurtosis after the wind episode can be predicted by combining the bandwidth measurement with the damping coefficient.

What carries the argument

The quadratic relation between kurtosis and bandwidth extracted from the forced-damped higher-order NLS simulations.

If this is right

  • Single-spectrum assessment of rogue wave likelihood becomes possible.
  • Post-wind evolution of kurtosis can be predicted by combining bandwidth with the damping coefficient.
  • The relation provides a practical tool for sea state analysis independent of exact wind details.
  • Kurtosis increases with bandwidth during wind forcing and decreases afterward.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The relation might extend to other wave models or real ocean data if the bandwidth-kurtosis link proves model-independent.
  • Measurements of wave spectra in the field could directly inform rogue wave forecasts without needing full time series.
  • Combining this with existing BFI methods could improve hybrid predictions for extreme wave risk.

Load-bearing premise

The quadratic kurtosis-bandwidth relation stays the same regardless of the exact forcing and damping values used in the higher-order NLS model.

What would settle it

Running simulations with substantially different forcing and damping coefficients and checking whether the fitted quadratic coefficients remain unchanged would test the claim.

Figures

Figures reproduced from arXiv: 1907.03490 by A. Armaroli, D. Eeltink, J. Kasparian, M. Brunetti, Y.M. Ducimeti\`ere.

Figure 1
Figure 1. Figure 1: FIG. 1. Temporal evolution. Solid lines: [PITH_FULL_IMAGE:figures/full_fig_p012_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Surface elevation for [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Spectrum for [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Exceedance probability distribution (EPD) of the envelope [PITH_FULL_IMAGE:figures/full_fig_p016_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Kurtosis as a function of BFI after the wind episode. Circles: [PITH_FULL_IMAGE:figures/full_fig_p017_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Kurtosis as a function of bandwidth [PITH_FULL_IMAGE:figures/full_fig_p018_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. (a) Coefficient [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Time evolution of a) kurtosis b) bandwidth c) spectral mean envelope [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Spectrum at [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10 [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Power spectrum for (a) the surface elevation [PITH_FULL_IMAGE:figures/full_fig_p026_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Kurtosis as a function of BFI, bandwidth, steepness, spectral mean for (a) full spectrum, [PITH_FULL_IMAGE:figures/full_fig_p027_12.png] view at source ↗
read the original abstract

We study statistical properties after a sudden episode of wind for water waves propagating in one direction. A wave with random initial conditions is propagated using a forced-damped higher order Nonlinear Schr\"odinger equation (NLS). During the wind episode, the wave action increases, the spectrum broadens, the spectral mean shifts up and the Benjamin-Feir index (BFI) and the kurtosis increase. Conversely, after the wind episode, the opposite occurs for each quantity. The kurtosis of the wave height distribution is considered the main parameter that can indicate whether rogue waves are likely to occur in a sea state, and the BFI is often mentioned as a means to predict the kurtosis. However, we find that while there is indeed a quadratic relation between these two, this relationship is dependent on the details of the forcing and damping. Instead, a simple and robust quadratic relation does exist between the kurtosis and the bandwidth. This could allow for a single-spectrum assessment of the likelihood of rogue waves in a given sea state. In addition, as the kurtosis depends strongly on the damping and forcing coefficients, by combining the bandwidth measurement with the damping coefficient, the evolution of the kurtosis after the wind episode can be predicted.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript studies statistical properties of unidirectional water waves after a sudden wind episode, propagated via a forced-damped higher-order Nonlinear Schrödinger equation. During the wind episode wave action and spectral bandwidth increase while the Benjamin-Feir index (BFI) and kurtosis rise; the opposite occurs after the episode ends. The authors report a quadratic kurtosis-BFI relation that depends on forcing and damping details, but claim a simple, robust quadratic relation exists between kurtosis and bandwidth that is independent of those details and therefore permits single-spectrum assessment of rogue-wave likelihood. They further suggest that combining the measured bandwidth with the damping coefficient allows prediction of post-episode kurtosis evolution.

Significance. If the claimed robustness of the kurtosis-bandwidth quadratic holds across forcing/damping variations, the result would supply a practical, spectrum-only diagnostic for rogue-wave probability that is less sensitive to model-specific dissipation than the BFI. The work usefully distinguishes the two predictors in a non-conservative setting and could inform operational wave forecasting.

major comments (2)
  1. [Abstract and Results] Abstract and Results section: the headline claim that a single quadratic k = a + b·BW + c·BW² remains independent of forcing and damping coefficients is load-bearing yet only asserted, not demonstrated. No quantitative ranges for forcing amplitudes, damping rates, or wind-episode durations are stated, nor is any table or figure shown that tracks stability of the fitted coefficients a, b, c when those parameters are varied by O(1) factors.
  2. [Results] Results section: the empirical quadratic fits are extracted from simulation output; the manuscript must report the number of independent realizations, the statistical uncertainty on each kurtosis and bandwidth value, and the goodness-of-fit metrics (e.g., R² or residual standard error) to allow assessment of whether the reported relation is statistically distinguishable from noise.
minor comments (2)
  1. [Methods] Clarify the precise operational definition of spectral bandwidth (e.g., second-moment width, half-width at half-maximum, or other) used in the quadratic regression.
  2. [Figures] Figure captions should state the specific forcing and damping coefficients employed for each plotted data set so readers can judge the explored parameter domain.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and for recognizing the potential utility of a spectrum-only rogue-wave diagnostic in non-conservative settings. We address each major comment below. Where the comments correctly identify missing quantitative details, we will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [Abstract and Results] Abstract and Results section: the headline claim that a single quadratic k = a + b·BW + c·BW² remains independent of forcing and damping coefficients is load-bearing yet only asserted, not demonstrated. No quantitative ranges for forcing amplitudes, damping rates, or wind-episode durations are stated, nor is any table or figure shown that tracks stability of the fitted coefficients a, b, c when those parameters are varied by O(1) factors.

    Authors: We agree that the independence claim requires explicit support. Although the simulations underlying the manuscript did employ multiple values of the forcing amplitude, damping rate, and episode duration, the manuscript did not tabulate those ranges or display the resulting coefficient stability. In the revised manuscript we will add a table (or supplementary figure) in the Results section that lists the specific parameter ranges explored and the corresponding fitted quadratic coefficients, thereby demonstrating that a, b, and c remain stable under O(1) changes. revision: yes

  2. Referee: [Results] Results section: the empirical quadratic fits are extracted from simulation output; the manuscript must report the number of independent realizations, the statistical uncertainty on each kurtosis and bandwidth value, and the goodness-of-fit metrics (e.g., R² or residual standard error) to allow assessment of whether the reported relation is statistically distinguishable from noise.

    Authors: We accept that these statistical diagnostics should have been included. The revised manuscript will state the number of independent realizations used for each reported data point, attach standard-error estimates to the kurtosis and bandwidth values, and report the R² (or equivalent residual metric) for each quadratic fit. revision: yes

Circularity Check

0 steps flagged

No significant circularity; empirical relation extracted from simulations

full rationale

The paper's central claim is an empirical quadratic relation between kurtosis and bandwidth identified from numerical integrations of the forced-damped higher-order NLS equation with random initial conditions. This relation is observed directly in the simulation outputs during and after wind episodes and is not defined in terms of itself, nor is any fitted parameter renamed as a first-principles prediction. The asserted independence from specific forcing/damping coefficients is checked by parameter variation inside the same model, which constitutes an explicit computational test rather than a self-referential reduction. No self-citation load-bearing steps, uniqueness theorems, or smuggled ansatzes appear in the derivation chain. The result is therefore self-contained as a model-specific empirical finding.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on numerical propagation of the forced-damped higher-order NLS equation and the empirical observation that the kurtosis-bandwidth quadratic is insensitive to the specific values of the forcing and damping coefficients.

free parameters (1)
  • forcing and damping coefficients
    Varied across simulations to test and demonstrate independence of the kurtosis-bandwidth relation from their specific values.
axioms (1)
  • domain assumption The forced-damped higher-order nonlinear Schrödinger equation provides a sufficiently accurate description of unidirectional water-wave statistics under wind and damping.
    All reported statistics and relations are obtained exclusively from numerical solutions of this equation.

pith-pipeline@v0.9.0 · 5777 in / 1304 out tokens · 31249 ms · 2026-05-25T01:07:29.775770+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages · 8 internal anchors

  1. [1]

    It contributes to improve the numerical stability of our solver, as it prevents unstable growth for this mode and all modes below

    This is attributed to the Taylor expansion of the growth rate. It contributes to improve the numerical stability of our solver, as it prevents unstable growth for this mode and all modes below. C. Comparison with Plant growth rate Plant and Wright [19] derived the following linear growth rate for the evolution of the intensity: βP = δξu2 ∗ cpK2k− 4νk2 (II...

    work page

  2. [2]

    The Rayleigh probability distribution corresponds to the envelope height distribution of a Gaussian, linear, process

    In general, the probability of finding a wave that exceeds x times Hs, the exceedance probability P (|a| > Hs), is a measure of how dangerous the sea state is. The Rayleigh probability distribution corresponds to the envelope height distribution of a Gaussian, linear, process. The Tayfun ([26]) and Fedele-Tayfun ([27]) distributions, take into account seco...

    work page

  3. [3]

    The quadratic dependence occurs in our parametric range, even for norm values sampled just after the wind episode

    to derive the quadratic dependence of the kurtosis on the bandwidth. The quadratic dependence occurs in our parametric range, even for norm values sampled just after the wind episode. Appendix C shows that for the non-truncated spectrum, for the spectral mean too, a quadratic relationship with the kurtosis exists. However, when the spectrum is truncated t...

    work page

  4. [4]

    S. Y. Annenkov and V. I. Shrira, Physica D: Nonlinear Phenomena 152-153, 665 (2001)

    work page 2001

  5. [5]

    On the kurtosis of deep-water gravity waves

    F. Fedele, Journal of Fluid Mechanics 782, 25 (2015), arXiv:1412.8231

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  6. [6]

    Fedele, J

    F. Fedele, J. Brennan, S. Ponce de Le´ on, J. Dudley, and F. Dias, Scientific Reports 6, 27715 (2016)

    work page 2016

  7. [7]

    Fedele, C

    F. Fedele, C. Lugni, and A. Chawla, Scientific Reports 7, 1 (2017)

    work page 2017

  8. [8]

    Onorato, L

    M. Onorato, L. Cavaleri, S. Fouques, O. Gramstad, P. A. Janssen, J. Monbaliu, A. R. Osborne, C. Pakozdi, M. Serio, C. T. Stansberg, A. Toffoli, and K. Trulsen, Journal of Fluid Mechanics 627, 235 (2009)

    work page 2009

  9. [9]

    Waseda, T

    T. Waseda, T. Kinoshita, and H. Tamura, Journal of Physical Oceanography 39, 621 (2009)

    work page 2009

  10. [10]

    Toffoli, O

    A. Toffoli, O. Gramstad, K. Trulsen, J. Monbaliu, E. Bitner-Gregersen, and M. Onorato, Journal of Fluid Mechanics 664, 313 (2010)

    work page 2010

  11. [11]

    N. Mori, M. Onorato, and P. A. E. M. Janssen, Journal of Physical Oceanography 41, 1484 (2011)

    work page 2011

  12. [12]

    Onorato, S

    M. Onorato, S. Residori, U. Bortolozzo, A. Montina, and F. T. Arecchi, Physics Reports Rogue waves and their generating mechanisms in different physical contexts, 528, 47 (2013)

    work page 2013

  13. [13]

    Spectral up- and downshifting of Akhmediev breathers under wind forcing

    D. Eeltink, A. Lemoine, H. Branger, O. Kimmoun, A. Chabchoub, M. Brunetti, J. Kasparian, C. Kharif, J. D. Carter, A. Chabchoub, M. Brunetti, and J. Kasparian, Physics of Fluids 29, 107103 (2017), arXiv:1709.09381

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  14. [14]

    Wave amplification in the framework of forced nonlinear Schrodinger equation: the rogue wave context

    A. Slunyaev, A. Sergeeva, and E. Pelinovsky, Physica D: Nonlinear Phenomena 303, 18 (2015), arXiv:1407.2443 [physics.flu-dyn]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  15. [15]

    P. A. E. M. Janssen, Journal of Physical Oceanography 33, 863 (2003)

    work page 2003

  16. [16]

    K. B. Dysthe, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 369, 105 (1979)

    work page 1979

  17. [17]

    Zakharov, Journal of Applied Mechanics and Technical Physics 9, 190 (1968)

    V. Zakharov, Journal of Applied Mechanics and Technical Physics 9, 190 (1968)

    work page 1968

  18. [18]

    Origin of heavy tail statistics in equations of the Nonlinear Schr\"odinger type: an exact result

    M. Onorato, D. Proment, G. El, S. Randoux, and P. Suret, Physics Letters A 380, 3173 (2016), arXiv:1601.04317

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  19. [19]

    Shemer and A

    L. Shemer and A. Sergeeva, Journal of Geophysical Research: Oceans 114, 1 (2009)

    work page 2009

  20. [20]

    Hara and C

    T. Hara and C. C. Mei, Journal of Fluid Mechanics 230, 429 (1991)

    work page 1991

  21. [21]

    Fabrikant, Wave Motion 2, 355 (1980)

    A. Fabrikant, Wave Motion 2, 355 (1980). 29

    work page 1980

  22. [22]

    W. J. Plant and J. W. Wright, Journal of Fluid Mechanics 82, 767 (1977)

    work page 1977

  23. [23]

    Song and M

    J.-B. Song and M. L. Banner, Journal of Physical Oceanography 32, 2541 (2002)

    work page 2002

  24. [24]

    S. Y. Annenkov and V. I. Shrira, Geophysical Research Letters 36, 1 (2009)

    work page 2009

  25. [25]

    Balac and F

    S. Balac and F. Mah´ e, Computer Physics Communications 184, 1211 (2013)

    work page 2013

  26. [26]

    Toffoli, A

    A. Toffoli, A. Babanin, M. Onorato, and T. Waseda, Geophysical Research Letters 37, L05603 (2010)

    work page 2010

  27. [27]

    Fedele, Journal of Fluid Mechanics 748, 692 (2014)

    F. Fedele, Journal of Fluid Mechanics 748, 692 (2014)

    work page 2014

  28. [28]

    V. E. Zakharov, S. I. Badulin, P. A. Hwang, and G. Caulliez, Journal of Fluid Mechanics 780, 503 (2015), arXiv:1411.7235

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  29. [29]

    M. A. Tayfun, Journal of Geophysical Research 85, 1548 (1980)

    work page 1980

  30. [30]

    M. A. Tayfun and F. Fedele, Ocean Engineering 34, 1631 (2007)

    work page 2007

  31. [31]

    Mori and P

    N. Mori and P. A. E. M. Janssen, Journal of Physical Oceanography 36, 1471 (2006)

    work page 2006

  32. [32]

    F. E. Snodgrass, G. W. Groves, K. F. Hasselmann, G. R. Miller, W. H. Munk, and W. H. Powers, Philosophical Transactions of the Royal Society of London A: Mathematical, Physical and Engineering Sciences 259, 431 (1966)

    work page 1966

  33. [33]

    Jiang, J

    H. Jiang, J. E. Stopa, H. Wang, R. Husson, A. Mouche, B. Chapron, and G. Chen, Journal of Geophysical Research: Oceans 121, 1446 (2016)

    work page 2016

  34. [34]

    Recurrence in the high-order nonlinear Schr\"odinger equation: a low dimensional analysis

    A. Armaroli, M. Brunetti, and J. Kasparian, Physical Review E 96, 1 (2017), arXiv:1703.09482

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  35. [35]

    Nonlinear stage of Benjamin-Feir instability in forced/damped deep water waves

    A. Armaroli, D. Eeltink, M. Brunetti, and J. Kasparian, Physics of Fluids 30, 17102 (2018), arXiv:1709.07850

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  36. [36]

    Onorato, A

    M. Onorato, A. R. Osborne, M. Serio, and L. Cavaleri, Physics of Fluids 17, 1 (2005)

    work page 2005

  37. [37]

    El Koussaifi, A

    R. El Koussaifi, A. Tikan, A. Toffoli, S. Randoux, P. Suret, and M. Onorato, Physical Review E 97, 012208 (2018)

    work page 2018

  38. [38]

    F. Dias, A. I. Dyachenko, and V. E. Zakharov, Physics Letters, Section A: General, Atomic and Solid State Physics 372, 1297 (2008), arXiv:0704.3352

    work page internal anchor Pith review Pith/arXiv arXiv 2008

  39. [39]

    Trulsen, I

    K. Trulsen, I. Kliakhandler, K. B. Dysthe, and M. G. Velarde, Physics of Fluids 12, 2432 (2000). 30

    work page 2000