pith. sign in

arxiv: 1907.09629 · v1 · pith:QZMITZRVnew · submitted 2019-07-22 · 🌊 nlin.PS · physics.flu-dyn

Ice Spiral Patterns on the Ocean Surface

Zhi Zong , Andrei Ludu This is my paper

Pith reviewed 2026-05-24 17:15 UTC · model grok-4.3

classification 🌊 nlin.PS physics.flu-dyn
keywords ice swirlsocean surfaceNavier-Stokes equationslogarithmic spiralsBessel solutionshydrodynamic modelspiral patternsnonlinear stability
0
0 comments X

The pith

A two-dimensional compressible fluid model generates logarithmic spiral ice patterns observed on the ocean surface.

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

The paper develops a two-dimensional compressible Navier-Stokes model to explain large-scale ice swirls on the ocean. Linear solutions from this model form a basis of Bessel functions that can create various spiral patterns and allow analysis of their time evolution and stability. Restricting the nonlinear equations to quadratic terms produces solutions with logarithmic spiral geometry. These are validated through multiple mathematical approaches and match experimental observations, with stability analyzed via Arnold's convexity method showing geometric phase transitions in the Hamiltonian.

Core claim

The authors show that the linearized 2D compressible Navier-Stokes model generates a basis of Bessel solutions for spiral patterns, while restricting the nonlinear system to quadratic terms yields swirl solutions with logarithmic spiral geometry. These are analyzed via Townes solitary modes, mapping to a sine-Gordon equation, and series expansion, producing pure radial, azimuthal, and spiral modes as well as multiple-spiral combinations that match observations, with nonlinear stability confirmed by Arnold's method and Hamiltonian plots versus order parameters indicating geometric phase transitions.

What carries the argument

The 2D compressible Navier-Stokes equations restricted to quadratic nonlinear terms, which produce swirl solutions with logarithmic spiral geometry.

If this is right

  • Pure radial, azimuthal, and spiral modes arise from the fully nonlinear equations.
  • Combinations of multiple-spiral solutions can be constructed that match experimental observations.
  • Nonlinear stability of the spiral patterns follows from Arnold's convexity method.
  • Plots of the Hamiltonian versus order parameters reveal geometric phase transitions.

Where Pith is reading between the lines

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

  • The sine-Gordon mapping may link these patterns to soliton dynamics in other wave systems.
  • The approach could apply to modeling similar surface patterns in other rotating fluids.
  • Predictions of stability thresholds might guide numerical simulations of ocean surface flows under ice.

Load-bearing premise

The assumption that the two-dimensional compressible Navier-Stokes equations restricted to quadratic terms capture the dominant physics of large-scale ice swirl formation without three-dimensional effects or direct ice mechanics coupling.

What would settle it

Field measurements of ice spiral angles, formation timescales, or pattern combinations that deviate substantially from the logarithmic geometry and Bessel-based predictions of the quadratic model.

Figures

Figures reproduced from arXiv: 1907.09629 by Andrei Ludu, Zhi Zong.

Figure 1
Figure 1. Figure 1: Large sea-ice swirl observed in the arctic ocean from airplane [9, 10]. [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Cross section in the θ = 0 vertical plane for the functions Φ1 (r, 0, 0) as solutions of the initial condition problems Φ0 = J1(3r/L) in solid red, r −0.3J1(3r/L) dotted curve, r −1/2J1(3r/L) dashed curve, r −1J1(3r/L) solid red, and r −2J1(3r/L) dotted-dashed curve. All curves represent Archimedean spiral with the same pitch but various decay laws with r [PITH_FULL_IMAGE:figures/full_fig_p016_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The same cross section and the same initial conditions as in Fig. 2, [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Solution Eq. (61) for m = 1 and two values for a0 = h/g = v/u for the nonlinear mass conservation equation describing coexisting logarithmic spiral patterns of ice in water. Corresponding spirals’ equations G(ln r) = mθ are plotted in the inset. This solution generates a distribution of logarithmic spirals in the plane with centers along a straight line, [PITH_FULL_IMAGE:figures/full_fig_p025_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Ice swirls in the Labrador current. Courtesy of NASA and [11]. [PITH_FULL_IMAGE:figures/full_fig_p026_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Example of time evolution of a spiral whose phase [PITH_FULL_IMAGE:figures/full_fig_p028_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Spiral shapes as locus of points in plane with Φ = 1 (ice crests) [PITH_FULL_IMAGE:figures/full_fig_p032_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Density plot of the IFF function Φ(r, θ) as a logarithmic spiral solution obtained from Eq. (80). ×ExpZ r 0 Γ(1 + m)  Jm+1 rωr √ a c  − Jm−1  rωr √ a c  √ aωr 2c  a1Γ(1 − m)J−m  rωr √ a c  + Γ(1 + m)Jm  rωr √ a c dr+a0  , (80) a0,1 being constants and with the phase obtained by implementing this ampli￾tude f in Eq. (76). One example of such solution is presented in [PITH_FULL_IMAGE:figures/f… view at source ↗
Figure 9
Figure 9. Figure 9: Fully nonlinear system stability diagrams for large and slow spinning [PITH_FULL_IMAGE:figures/full_fig_p039_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Same stability diagram as in Fig. 9 for medium size spirals, diameter [PITH_FULL_IMAGE:figures/full_fig_p039_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Dependence of the total energy of the spiral system in relative units [PITH_FULL_IMAGE:figures/full_fig_p040_11.png] view at source ↗
read the original abstract

We investigate a new two-dimensional compressible Navier-Stokes hydrodynamic model design to explain and study large scale ice swirls formation at the surface of the ocean. The linearized model generates a basis of Bessel solutions from where various types of spiral patterns can be generated and their evolution and stability in time analyzed. By restricting the nonlinear system of equations to its quadratic terms we obtain swirl solutions emphasizing logarithmic spiral geometry. The resulting solutions are analyzed and validated using three mathematical approaches: one predicting the formation of patterns as Townes solitary modes, another approach mapping the nonlinear system into a sine-Gordon equation, and a third approach uses a series expansion. Pure radial, azimuthal and spiral modes are obtained from the fully nonlinear equations. Combinations of multiple-spiral solutions are also obtained, matching the experimental observations. The nonlinear stability of the spiral patterns is analyzed by Arnold's convexity method, and the Hamiltonian of the solutions is plotted versus some order parameters showing the existence of geometric phase transitions.

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

3 major / 1 minor

Summary. The paper develops a 2D compressible Navier-Stokes model for large-scale ice swirls on the ocean surface. Linearized equations produce a Bessel-function basis for spiral patterns whose time evolution and stability are analyzed; truncation to quadratic nonlinear terms yields swirl solutions with logarithmic spiral geometry. These are examined via three routes (Townes solitary modes, mapping to the sine-Gordon equation, and series expansion), pure radial/azimuthal/spiral modes are extracted from the full nonlinear system, multiple-spiral combinations are stated to match observations, and nonlinear stability is assessed by Arnold convexity with Hamiltonian plots versus order parameters indicating geometric phase transitions.

Significance. If the quadratic truncation and resulting spirals were shown to reproduce measured ice-swirl geometry and scales, the work would supply a compact hydrodynamic framework linking Bessel bases, Townes modes, and sine-Gordon structure to ocean-surface patterns. At present the internal mathematical consistency is of interest to nonlinear hydrodynamics, but the absence of any quantitative comparison to field data keeps the physical significance low.

major comments (3)
  1. [Abstract] Abstract: the claim that 'combinations of multiple-spiral solutions are also obtained, matching the experimental observations' is unsupported; the manuscript presents no observational datasets, measured parameters (pitch, radius, wavelength), error bars, or side-by-side metrics, so the central assertion of explanatory power rests on unshown steps.
  2. [Model reduction] Model design and quadratic reduction: the premise that restricting the 2D compressible Navier-Stokes equations to quadratic nonlinear terms 'preserves the essential spiral geometry' is asserted without explicit demonstration (e.g., comparison of full versus truncated solutions or bounds on higher-order contributions), yet this truncation is load-bearing for the logarithmic-spiral claim.
  3. [Validation approaches] Validation section: all three mathematical routes (Townes modes, sine-Gordon mapping, series expansion) together with the Arnold stability analysis remain internal to the truncated model; none supplies an external test against real ice-swirl geometry or scales, leaving the physical applicability unverified.
minor comments (1)
  1. [Stability analysis] The notation for the order parameters used in the Hamiltonian plots versus geometric phase transitions should be defined explicitly in the text or a table.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments on the manuscript. We address each major comment below, indicating where revisions will be made.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that 'combinations of multiple-spiral solutions are also obtained, matching the experimental observations' is unsupported; the manuscript presents no observational datasets, measured parameters (pitch, radius, wavelength), error bars, or side-by-side metrics, so the central assertion of explanatory power rests on unshown steps.

    Authors: We agree that the abstract phrasing overstates the connection to observations. The manuscript is a theoretical study, and the reference to matching observations is based on qualitative geometric similarity rather than quantitative comparison. We will revise the abstract to remove the claim of matching experimental observations and instead note that the solutions produce spiral geometries of the type seen in ice swirls. revision: yes

  2. Referee: [Model reduction] Model design and quadratic reduction: the premise that restricting the 2D compressible Navier-Stokes equations to quadratic nonlinear terms 'preserves the essential spiral geometry' is asserted without explicit demonstration (e.g., comparison of full versus truncated solutions or bounds on higher-order contributions), yet this truncation is load-bearing for the logarithmic-spiral claim.

    Authors: The quadratic truncation is introduced to isolate the leading nonlinear terms responsible for the logarithmic spiral solutions while permitting analytical treatment. We acknowledge that the manuscript does not supply an explicit comparison or error bound relative to the full system. We will add a short justification, including a scaling argument showing that higher-order terms remain perturbative for the amplitudes and wavenumbers of interest in the spiral regime. revision: yes

  3. Referee: [Validation approaches] Validation section: all three mathematical routes (Townes modes, sine-Gordon mapping, series expansion) together with the Arnold stability analysis remain internal to the truncated model; none supplies an external test against real ice-swirl geometry or scales, leaving the physical applicability unverified.

    Authors: The three routes and the Arnold analysis constitute internal consistency checks within the model. The work is framed as a hydrodynamic framework rather than a data-validated prediction. We will add a limitations paragraph that explicitly states the absence of quantitative field comparisons and outlines observable signatures (e.g., predicted pitch-angle dependence on rotation rate) that could be tested against future observations. revision: partial

Circularity Check

0 steps flagged

Derivation uses independent mathematical reductions with no tautological steps.

full rationale

The abstract and provided context describe a derivation starting from 2D compressible Navier-Stokes equations, linearizing to obtain a Bessel function basis, truncating to quadratic nonlinear terms to produce logarithmic spiral solutions, then applying Townes solitary modes, a sine-Gordon mapping, and series expansions for validation, plus Arnold stability analysis. None of these steps are shown to reduce by construction to fitted parameters, self-citations, or renamed inputs; they are standard analytic techniques applied to the model. The qualitative statement that multi-spiral combinations 'match the experimental observations' is an unsupported claim rather than a load-bearing derivation step that equates output to input. No self-citation chains or ansatzes smuggled via prior work are referenced. The paper is therefore self-contained in its mathematical chain against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Only the abstract is available, so the ledger is populated from the stated model assumptions; the central claim rests on the compressible 2D NS equations plus the validity of the quadratic truncation and the three named mathematical mappings.

axioms (2)
  • domain assumption The ocean surface ice layer can be modeled by a two-dimensional compressible Navier-Stokes system
    Stated in the model design sentence of the abstract
  • ad hoc to paper Restricting the nonlinear terms to quadratic order preserves the essential spiral geometry
    Invoked when obtaining swirl solutions from the quadratic system

pith-pipeline@v0.9.0 · 5685 in / 1479 out tokens · 26600 ms · 2026-05-24T17:15:03.058195+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

43 extracted references · 43 canonical work pages

  1. [2]

    Guyenne, and Pˇ arˇ au, Numerical Simulation of Solitary-Wave Scattering and Damping in Fragmented Sea Ice,The 27th Int

    P. Guyenne, and Pˇ arˇ au, Numerical Simulation of Solitary-Wave Scattering and Damping in Fragmented Sea Ice,The 27th Int. Ocean Polar Eng. Conf. (2017)

    work page 2017

  2. [3]

    E. I. Pˇ arˇ au, Solitary interfacial hydroelastic waves, Phil. Trans. A 376 (2017) 0099

    work page 2017

  3. [4]

    Li, Ocean surface waves in an ice-free Arctic Ocean, Ocean Dynamics 66 (2016)

    J.-G. Li, Ocean surface waves in an ice-free Arctic Ocean, Ocean Dynamics 66 (2016). 8989-1004

    work page 2016

  4. [5]

    Huang, K

    L. Huang, K. Ren, M. Li, Z. Tukovi´ c, P. Cardiff, and G. Thomas, Fluid- structure interaction of a large ice sheet in waves, Ocean Engineering 182 (2019) 102-111

    work page 2019

  5. [6]

    Ningbo Zhang, Xing Zheng, Qingwei Ma, Study on wave-induced kinematic responses and flexures of ice floe by Smoothed Particle Hydrodynamics, Computers & Fluids (2019)

    work page 2019

  6. [7]

    L. G. Bennetts, T. D. Williams, Water wave transmission by an array of floating discs, Proc. R. Soc. Lond. Math. Phys. Eng. Sci. 471 (2015) 20140698

    work page 2015

  7. [8]

    Toffoli, L

    A. Toffoli, L. G. Bennetts, M. H. Meylan, C. Cavaliere, A. Alberello, J. El- snab, et al, Sea ice floes dissipate the energy of steep ocean waves, Geophys. Res. Lett. 42 (2015) 85478554

    work page 2015

  8. [9]

    The Asahi Shimbun Weekly ( http : //gettyimages.com))

    work page

  9. [10]

    5, 2016 (coasthttps://mashable.com/2016/07/04/hurricane- sea-ice-swirl-pictures/#1e1J5FeTEiq6)

    Andrew Freedman, Mini hurricane of sea ice spotted of the Canadian, Mashable, Jul. 5, 2016 (coasthttps://mashable.com/2016/07/04/hurricane- sea-ice-swirl-pictures/#1e1J5FeTEiq6)

    work page 2016

  10. [11]

    Petrie and A

    B. Petrie and A. Isenor, The near-surface circulation and exchange in the Newfoundland Grand Banks region, Atmospheric-Ocean, Atmosphere- Ocean 23, 3 (1987) 209-227

    work page 1987

  11. [12]

    B. C. Justin, J. M. Amundson, R. Cassotto, C.-C. Kuod, and M. Dennind, Quantifying flow and stress in ice m´ elange, the worlds largest granular material, Proc. Natl. Acad. Sci. 115, 20 (2018) 5105-5110. 42

    work page 2018

  12. [13]

    Lin, and F

    C. Lin, and F. H. Shu, On the spiral structure of disk galaxies, Fluid Mechanics 2 (1987) 561-570

    work page 1987

  13. [14]

    V. A. Squire, Of ocean waves and sea-ice revisited, Cold Reg. Sci. Technol. 49 (2007) 110133

    work page 2007

  14. [15]

    Fily, and D

    M. Fily, and D. A. Rothrock, Sea ice tracking by nested correlations, IEEE transactions on geoscience and remote sensing 5 (1987) 570-580

    work page 1987

  15. [16]

    Abramowitz and I

    M. Abramowitz and I. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (Cambridge University Press, Cambridge 1964)

    work page 1964

  16. [17]

    Peng and K.-Q

    J. Peng and K.-Q. Zhu, J. Fluid Mech. , 512 (2004) 21-45

    work page 2004

  17. [18]

    M. G. McPhee, Ice-ocean momentum transfer for the AIDJEX ice model, Aidjex Bull 29 (1975) 93-111

    work page 1975

  18. [19]

    H.-J. Dai, J. C. McWilliams, and J.-H. Liang, Wave-driven mesoscale cur- rents in a marginal ice zone, Ocean Modelling 134 (2019) 1-17

    work page 2019

  19. [20]

    Guyenne, and E

    P. Guyenne, and E. I. Pˇ arˇ au, Computations of fully nonlinear hydroelastic solitary waves on deep water, J. Fluid Mechanics 713 (2012) 307-329

    work page 2012

  20. [21]

    M. W. Mueller and W. D. Arnett, Propagating Star formation and Irregular Structure in Spiral Galaxies, Astrophys. J. 210 (1976) 670-678

    work page 1976

  21. [22]

    A. E. Karperaki, FEM hydroelastic models with application to the nonlin- ear response of large floating bodies in shallow wave conditions, Procedia Computer Science 66 (2015) 122-31

    work page 2015

  22. [23]

    Landau and E

    L. Landau and E. M. Lifshitz, Fluid Mechanics (Pergamon Press, Oxford 1987)

    work page 1987

  23. [24]

    G. L. Lamb, Jr., Elements of Soliton Theory (John Willey & Sons, New York 1980)

    work page 1980

  24. [25]

    M. S. Longuet-Higgins, Statistical properties of wave groups in a random sea state, Phil. Trans. Roy. Soc. London A 312, 1521 (1984) 219-250

    work page 1984

  25. [26]

    M. J. Ablowitz, Nonlinear Dispersive Waves. Asymptotic Analysis and Soli- tons (Cambridge University Press, Cambridge 2011)

    work page 2011

  26. [27]

    D. J. Benny, and G. J. Roskes, Wave instabilities, Stud. Appl. Math. 48 (1969) 377-385

    work page 1969

  27. [28]

    V. D. Djordjevic, and L. G. Redekopp, One two-dimensional packets of capillary-gravity waves, J. Fluid Mech. 79 (1977) 703-714

    work page 1977

  28. [29]

    Sulem, and P

    C. Sulem, and P. Sulem, The Nonlinear Schr¨ odinger equation: Self- Focusing and Wave Collapse (Springer-Verlag, Berlin 1999). 43

    work page 1999

  29. [30]

    Rosenau, On non-analytic solitary waves formed by a nonlinear dis- persion, Phys

    P. Rosenau, On non-analytic solitary waves formed by a nonlinear dis- persion, Phys. Lett. A 230, 56, (1997) 305318; On a class of nonlinear dispersive-dissipative interactions, Physica D 230, 56 (1998) 535546; Com- pact and noncompact dispersive structures, Phys. Lett. A 275, 3 (2000) 193203

    work page 1997

  30. [31]

    A. Ludu, G. Stoitcheva and J. P. Draaye, Similarity Analysis of Nonlin- ear Equations and Bases of Finite Wavelength Solitons, International, J. Modern Phys. E 9 (2000) 263-278

    work page 2000

  31. [32]

    A. D. Polyanin and V. F. Zaitsev, Handbook of Ordinary Differential Equa- tions (CRC Press, Boca Raton 2018). G. Fibich, and G. Papanicolau, Self- focusing in the perturbed and unperturbed nonlinear Schr¨ odinger equation in critical dimension, SIAM J. Appl. Math. 60 (1999) 183-240

    work page 2018

  32. [33]

    Weinstein, Nonlinear Schr¨ odinger equations and sharp interpolation estimates, Comm

    M. Weinstein, Nonlinear Schr¨ odinger equations and sharp interpolation estimates, Comm. Math. Phys. 87 (1983) 567-576

    work page 1983

  33. [34]

    Papanicolau, C

    G. Papanicolau, C. Sulem, P. L. Sulem, and X. P. Wang, The focusing singularity of the Davey-Stewartson equations for gravity-capillary surface waves, Physica D 72 (1994) 61-86

    work page 1994

  34. [35]

    M. J. Ablowitz, and Z. Musslimani, Spectral renormalization method for computing self-localized solutions to nonlinear systems, Opt. Lett. 30 (2005) 2140-2142

    work page 2005

  35. [36]

    M. J. Ablowitz, and H. Segur, Solitons and The Inverse Scattering Trans- form (SIAM, Philadelphia 1981)

    work page 1981

  36. [37]

    M. J. Ablowitz, G. Biondini, S. Blair, Nonlinear Schr¨ odinger equations with mean terms in non-resonant multi-dimensional quadratic materials, Phys. Rev. E 63 (2001) 605-620

    work page 2001

  37. [38]

    L. C. Crasovan, J. P. Torres, D. Mihalache, and L. Tomer, Arresting wave collapse by self-rectification, Phys. Rev. Lett. 9 (2003) 063904

    work page 2003

  38. [39]

    Ludu, Boundaries of a Complex World, Sect

    A. Ludu, Boundaries of a Complex World, Sect. 4.4 (Springer-Verlag, Hei- delberg 2016)

    work page 2016

  39. [40]

    Sukhinin, A

    A. Sukhinin, A. B. Aceves, J.-C. Diels, L. Arissian, Collapse events of two- color optical beams, Phys. Rev. A 95 (2017) 031801

    work page 2017

  40. [41]

    V. E. Zakharov and E. A. Kuznetsov, Physica, 18D (1986) 455-463

    work page 1986

  41. [42]

    D. D. Holm, J. E. Marsden, T. Ratiu, and A. Weinstein, Nonlinear stability conditions and a priori estimates for barotropic hydrodynamics, Phys. Lett. A 98, 1-2 (1983) 15-21

    work page 1983

  42. [43]

    Zlotnik, On the energy dissipative spatial discretization of the barotropic quasi-gasdynamic and compressible NavierStokes equations in polar coor- dinates, Russ

    A. Zlotnik, On the energy dissipative spatial discretization of the barotropic quasi-gasdynamic and compressible NavierStokes equations in polar coor- dinates, Russ. J. Num. Analysis Math. Modelling 33, 3 (2018) 199-210. 44

    work page 2018

  43. [44]

    P. G. Drazin and W. H. Reid, Hydrodynamic stability (Cambridge Univ. Press, London, 1981). 45

    work page 1981