Pith. sign in

REVIEW 2 minor 48 references

Analysis and numerical simulations of a landfast ice model

T0 review · 0 major / 2 minor · reviewed 2026-05-08 · grok-4.3

Challenge this review Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A model for landfast ice has local strong solutions, global solutions near constant equilibria without external forces, and time-periodic solutions, with numerics showing stationary states of vanishing velocity.

desk verdict The paper gives new local and global strong well-posedness plus time-periodic results for the landfast ice extension of the Hibler model, backed by simple numerics on stationary states. read the letter →

arxiv 2604.23596 v1 submitted 2026-04-26 math.AP

Felix Brandt , Carolin Mehlmann This is my paper · ORCID
classification math.AP
keywords landfasticeseadynamicswell-posednessviscous-plasticmodelpartialdifferentialequationsstationaryequilibriatime-periodicsolutionsnumericalsimulation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper studies equations that describe sea ice fixed to coastlines and therefore nearly motionless. It proves that these equations admit unique local-in-time strong solutions from suitable starting data. When external forces are removed and the initial state is close to a uniform equilibrium, the solutions extend to all future times. The same equations also possess solutions that repeat after a fixed period. Numerical runs illustrate the ice velocity dropping to exactly zero and remaining there, matching the analytical guarantees of stability near equilibria.

What carries the argument

The coupled system of nonlinear partial differential equations that evolves ice thickness, concentration, and velocity under a viscous-plastic constitutive law modified to allow vanishing motion.

What would settle it

An explicit initial condition near a constant equilibrium, run without external forces, that develops large velocities or ceases to exist after finite time would contradict the global well-posedness and stability statements.

Watch

Extended reading notes

Core claim

The system of nonlinear partial differential equations for landfast ice thickness, concentration, and velocity possesses strong local solutions for general initial data, strong global solutions when external forces vanish and data lies near constant equilibria, and time-periodic solutions; numerical experiments confirm the formation of stationary equilibria in which ice velocity is identically zero.

Load-bearing premise

The ice's resistance to deformation and its interaction with the coast are captured exactly by the chosen viscous-plastic relations and forcing terms.

Editorial extensions

If this is right

  • Strong solutions exist at least locally in time for initial data in appropriate Sobolev spaces.
  • Global-in-time strong solutions exist when external forces are absent and initial data is sufficiently close to a constant equilibrium.
  • Time-periodic solutions exist for the forced system.
  • Numerical solutions starting near equilibrium converge to stationary states with exactly zero velocity.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The stability near equilibria could be used to construct numerical schemes that exactly preserve zero-velocity states once reached.
  • The same analytical approach may extend to other grounded-ice models whose rheology permits zero motion.
  • Periodic solutions suggest that coastal ice could exhibit repeating cycles under seasonally varying winds or currents.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

0 major / 2 minor

Summary. The manuscript analyzes a landfast ice model extending the Hibler viscous-plastic sea-ice model with landfast modifications from Lemieux et al. It establishes local strong well-posedness, global strong well-posedness without external forces for initial data near constant equilibria (Theorem 4.1), and existence of time-periodic solutions (Theorem 5.2). Numerical simulations demonstrate the formation of stationary equilibrium states with vanishing ice velocity, consistent with the theoretical results.

Significance. If the well-posedness theorems hold, the paper makes a valuable contribution by providing rigorous mathematical analysis for a model relevant to Arctic climate simulations. The combination of existence results for local, global, and periodic solutions, along with numerical illustrations of key physical behaviors like stationary states, offers new insights into the stability and long-term dynamics of landfast ice. This bridges applied mathematics and geophysical modeling effectively.

minor comments (2)
  1. [Numerical simulations] The numerical illustrations show qualitative agreement with the theorems but lack quantitative error metrics or convergence rates, which would help confirm the accuracy of the simulations.
  2. [Introduction] A more detailed comparison with existing mathematical analyses of the classical Hibler model would contextualize the novelty of the landfast extensions.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of our manuscript, which accurately summarizes the local and global strong well-posedness results, the existence of time-periodic solutions, and the supporting numerical simulations for the landfast ice model. We appreciate the recognition of the work's relevance to Arctic climate modeling and the bridging of mathematical analysis with geophysical applications.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity detected in analytical claims

full rationale

The paper establishes local strong well-posedness, global strong well-posedness near equilibria without forcing, and existence of time-periodic solutions for the landfast ice extension of the Hibler viscous-plastic system. These are proved using standard fixed-point and continuation arguments for quasilinear parabolic-hyperbolic PDEs once the landfast modifications (vanishing velocity near coast/grounding) are incorporated into the constitutive law. The modeling framework is taken from the external citation to Lemieux et al. [27], but the existence theorems themselves are independent mathematical results with no reduction to fitted parameters, self-definitional loops, or load-bearing self-citations. Numerical illustrations of equilibria are presented as consistent with the theorems rather than as substitutes for them. The derivation chain is therefore self-contained against external benchmarks.

Assumptions & free parameters 0 free parameters · 1 assumptions · 0 invented entities

The paper relies on the standard viscous-plastic rheology and landfast ice parameterizations from the cited references; no new free parameters, ad-hoc axioms, or invented physical entities are introduced in the abstract.

assumptions (1)
  • domain assumption Standard assumptions on initial data regularity and forcing terms required for strong solutions of the viscous-plastic sea-ice PDE system
    Invoked to obtain local and global well-posedness

how reviews work

0 comments
Cite this review

Pith. "Pith review of Analysis and numerical simulations of a landfast ice model." pith.science (2026). https://pith.science/paper/2604.23596

@misc{pith2026260423596,
  author       = {Pith},
  title        = {Pith review of: Analysis and numerical simulations of a landfast ice model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2604.23596}},
  note         = {Machine review of arXiv:2604.23596}
}
read the original abstract

In this manuscript, we consider a common modeling framework for Arctic landfast ice based on the work of Lemieux et al. [27], which is designed for use in large-scale climate models. This approach extends the classical viscous-plastic sea-ice model introduced by Hibler [18], which remains the most used model for simulating large-scale sea-ice dynamics in climate science. In particular, landfast ice refers to sea-ice that is attached to the coastline or grounded and therefore exhibits nearly vanishing motion. We present a rigorous analytical and numerical study of this landfast ice model. The main analytical contributions are the local strong well-posedness, the global strong well-posedness in the absence of external forces and for initial data close to constant equilibrium solutions, and the existence of time-periodic solutions. Complementing the analysis, we perform numerical simulations that illustrate key qualitative differences between landfast ice and classical viscous-plastic sea-ice models. In particular, the simulations reveal the formation of stationary equilibrium states characterized by vanishing ice velocity. These observations are consistent with the global-in-time existence result close to equilibria established in Theorem 4.1 as well as the time-periodic result in Theorem 5.2. The combined analytical and numerical results provide new insight into the structure, stability, and long-term behavior of landfast ice dynamics.

Figures

Figures reproduced from arXiv: 2604.23596 by the authors.

Figure 1
Figure 1. Comparison of sea-ice concentration, thickness and velocity after 2 days of simulation using the standard viscous–plastic model (I, II) and its landfast ice extension (III, IV), both under dominant rightward wind forcing. While in the standard case (I, II), the ice drifts with the wind, the landfast ice extension (III, IV) leads to the formation of stationary sea-ice along the coast and an opening of sea-ice (polyny… view at source ↗
Figure 2
Figure 2. Sea-ice evolution without external forcing. The kinetic energy is presented in Panel I. Panel II and Panel III show the sea-ice concentration after 6 and 12 simulated days respectively. In the absence of external forces the sea-ice velocity approaches zero over time, which results in a stable sea-ice cover over time. (compare Panels II and IV in view at source ↗
Figure 3
Figure 3. Sea-ice evolution under constant rightward wind forcing over 27 simulated days. Panels I–IV show the velocity field scaled by the sea-ice concentration. In Panel I, the scaling factor is 1, as the entire domain is initially ice-covered. Over time, the system approaches a stationary state, indicated by near-zero velocities in ice-covered regions (gray plateau in Panels III and IV), while non-zero velocities occur onl… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

48 extracted references · 48 canonical work pages

  1. [1]

    Amann,Linear and Quasilinear Parabolic Problems

    H. Amann,Linear and Quasilinear Parabolic Problems. Monographs in Mathematics, vol. 89, Birkh¨ auser, 1995

    work page 1995

  2. [2]

    Arendt, S

    W. Arendt, S. Bu, The operator-valued Marcinkiewicz multiplier theorem and maximal regularity.Math. Z.240(2002), 311–343

    work page 2002

  3. [3]

    T. Binz, F. Brandt, M. Hieber, Interaction of geophysical flows with sea ice dynamics.NoDEA Nonlinear Differential Equations Appl.33(2026), Paper No. 35

    work page 2026

  4. [4]

    Blockley, M

    E. Blockley, M. Vancoppenolle, E. Hunke, C. Bitz, D. Feltham, J.-F. Lemieux, M. Losch, E. Maisonnave, D. Notz, P. Rampal, S. Tietsche, B. Tremblay, A. Turner, F. Massonnet, E. ´Olason, A. Roberts, Y. Aksenov, T. Fichefet, G. Garric, D. Iovino, G. Madec, C. Rousset, D. Salas y Melia, D. Schroeder, The future of sea ice modeling: where do we go from here?.B...

    work page 2020

  5. [5]

    Bothe, J

    D. Bothe, J. Pr¨ uss,Lp-theory for a class of non-Newtonian fluids.SIAM J. Math. Anal.39(2007), 379–421

    work page 2007

  6. [6]

    Boutros, X

    D. Boutros, X. Liu, M. Thomas, E.S. Titi, Global well-posedness of the elastic-viscous-plastic sea- ice model with the inviscid Voigt-regularisation.Math. Models Methods Appl. Sci.To appear, 2026. https://doi.org/10.1142/S0218202526500156

    work page doi:10.1142/s0218202526500156 2026

  7. [7]

    Geophysical Flow Models: An Approach by Quasilinear Evolution Equations , url=

    F.C.H.L. Brandt,Geophysical Flow Models: An Approach by Quasilinear Evolution Equations. PhD thesis, Technische Universit¨ at Darmstadt, 2024. https://doi.org/10.26083/tuprints-00027378

    work page doi:10.26083/tuprints-00027378 2024

  8. [8]

    Brandt, Well-posedness of Hibler’s parabolic-hyperbolic sea ice model.J

    F. Brandt, Well-posedness of Hibler’s parabolic-hyperbolic sea ice model.J. Evol. Equ.25(2025), Paper No. 82

    work page 2025

Show all 48 references
  1. [9]

    Brandt, K

    F. Brandt, K. Disser, R. Haller-Dintelmann, M. Hieber, Rigorous analysis and dynamics of Hibler’s sea ice model. J. Nonlinear Sci.32(2022), Paper No. 50

    2022

  2. [10]

    Brandt, M

    F. Brandt, M. Hieber, Time periodic solutions to Hibler’s sea ice model.Nonlinearity36(2023), 3109–3124

    2023

  3. [11]

    Chatta, B

    S. Chatta, B. Khouider, Well posedness of the regularized-Hibler model of sea-ice dynamics. Preprint, 2025

    2025

  4. [12]

    Chatta, B

    S. Chatta, B. Khouider, M. Kesri, Linear well-posedness of regularized equations of sea-ice dynamics.J. Math. Phys.64 (2023), Paper 051504

    2023

  5. [13]

    R. Denk, G. Dore, M. Hieber, J. Pr¨ uss, A. Venni, New thoughts on old results of R.T. Seeley.Math. Ann.328(2004), 545–583

    2004

  6. [14]

    R. Denk, F. Gmeineder, M. Hieber, On the singular limit in Hibler’s sea ice model. arXiv:2511.09327. ANALYSIS AND NUMERICAL SIMULATIONS OF A LANDFAST ICE MODEL 31

  7. [15]

    R. Denk, M. Hieber, J. Pr¨ uss,R-Boundedness, Fourier Multipliers and Problems of Elliptic and Parabolic Type. Mem. Amer. Math. Soc., vol. 166, 2003, no. 788

    2003

  8. [16]

    Dingel, K

    S. Dingel, K. Disser, Global existence and uniqueness for Hibler’s visco-plastic sea-ice model. arXiv:2508.16537

  9. [17]

    Amann,One-Parameter Semigroups for Linear Evolution Equations

    H. Amann,One-Parameter Semigroups for Linear Evolution Equations. Graduate Texts in Mathematics, vol. 194, Springer-Verlag, New York, 2000

    2000

  10. [18]

    Hibler, A dynamic thermodynamic sea ice model.J

    W.D. Hibler, A dynamic thermodynamic sea ice model.J. Phys. Oceanogr.9(1979), 815–846

    1979

  11. [19]

    Hieber, C

    M. Hieber, C. Stinner, Strong time periodic solutions to Keller-Segel systems: an approach by the quasilinear Arendt-Bu theorem.J. Differential Equations269(2020), 1636–1655

    2020

  12. [20]

    Hunke, J.K

    E.C. Hunke, J.K. Dukowicz, An elastic-viscous-plastic model for sea ice dynamics.J. Phys. Oceanogr.27(1997), 1849–1867

    1997

  13. [21]

    Kimmrich, S

    M. Kimmrich, S. Danilov, M. Lorsch, On the convergence of the modified elastic-viscous-plastic method for solving the sea ice momentum equation.J. Comput. Phys.296(2015), 90–100

    2015

  14. [22]

    Kreyscher, M

    M. Kreyscher, M. Harder, P. Lemke, G. Flato, M. Gregory, Results of the sea ice model intercomparison project: evaluation of sea ice rheology schemes for use in climate simulations.J. Geophys. Res.105(2000), 11299–11320

    2000

  15. [23]

    K¨ onig Beatty, D.M

    C. K¨ onig Beatty, D.M. Holland, Modeling landfast sea ice by adding tensile strength.J. Phys. Oceanogr.40(2010), 185–198

    2010

  16. [24]

    Kunstmann, L

    P.C. Kunstmann, L. Weis, Perturbation theorems for maximal L p-regularity.Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4)30(2001), 415–435

    2001

  17. [25]

    Kunstmann, L

    P.C. Kunstmann, L. Weis, MaximalL p-regularity for Parabolic Equations, Fourier Multiplier Theorems andH ∞- functional Calculus. In:Functional Analytic Methods for Evolution Equations. M. Iannelli, R. Nagel and S. Piazzera (eds.), Springer, 2004, 65–311

    2004

  18. [26]

    LeCrone, J

    J. LeCrone, J. Pr¨ uss, M. Wilke, On quasilinear parabolic evolution equations in weighted Lp-spaces II.J. Evol. Equ.14 (2014), 509–533

    2014

  19. [27]

    Lemieux, F

    J.-F. Lemieux, F. Dupont, P. Blain, F. Roy, G.C. Smith, G.M. Flato, Improving the simulation of landfast ice by combining tensile strength and a parameterization for grounded ridges.J. Geophys. Res. Oceans121(2016), 7354– 7368

    2016

  20. [28]

    Lemieux, B

    J.-F. Lemieux, B. Tremblay, Numerical convergence of viscous-plastic sea ice models.J. Geophys. Res.114(2009), C05009

    2009

  21. [29]

    X. Liu, M. Thomas, E.S. Titi, Well-posedness of Hibler’s dynamical sea-ice model.J. Nonlinear Sci.32(2022), Paper No. 49

    2022

  22. [30]

    X. Liu, M. Thomas, E.S. Titi, Plastic limit of a viscoplastic Burgers equation – a toy model for sea-ice dynamics. arXiv:2601.06489

  23. [31]

    Mahoney, Landfast sea ice in a changing Arctic

    A.R. Mahoney, Landfast sea ice in a changing Arctic. In:Arctic Report Card 2018. E. Osborne, J.A. Richter-Menge, M.O. Jeffries (Eds.), NOAA, 2014

    2018

  24. [32]

    Mehlmann, S

    C. Mehlmann, S. Kahl, A hybrid particle-continuum method for simulating landfast sea ice via subgrid iceberg inter- action.J. Comput. Phys.554, Paper No. 114752

  25. [33]

    Mehlmann, P

    C. Mehlmann, P. Korn, Sea-ice on triangular grids.J. Comput. Phys.428(2021), Paper No. 110086

    2021

  26. [34]

    Mehlmann, T

    C. Mehlmann, T. Richter, A modified global Newton solver for viscous-plastic sea ice models.Ocean Model.116(2017), 96–117

    2017

  27. [35]

    Mehlmann, T

    C. Mehlmann, T. Richter, A finite element multigrid framework to solve the sea ice momentum equation.J. Com- put. Phys.348(2017), 847–861

    2017

  28. [36]

    Piersanti, R

    P. Piersanti, R. Temam, On the dynamics of shallow ice sheets: modeling and analysis.Adv. Nonlinear Anal.12 (2023), Paper No. 20220280

    2023

  29. [37]

    Pr¨ uss, G

    J. Pr¨ uss, G. Simonett, Maximal regularity for evolution equations in weighted L p-spaces.Arch. Math.82(2004), 415–431

    2004

  30. [38]

    Pr¨ uss, G

    J. Pr¨ uss, G. Simonett,Moving Interfaces and Quasilinear Parabolic Evolution Equations. Monographs in Mathematics, vol. 105, Birkh¨ auser, 2016

    2016

  31. [39]

    Pr¨ uss, G

    J. Pr¨ uss, G. Simonett, M. Wilke, Critical spaces for quasilinear parabolic evolution equations and applications.J. Dif- ferential Equations264(2018), 2028–2074

    2018

  32. [40]

    Pr¨ uss, G

    J. Pr¨ uss, G. Simonett, R. Zacher, On convergence of solutions to equilibria for quasilinear parabolic problems.J. Dif- ferential Equations246(2009), 3902–3931

    2009

  33. [41]

    Seinen, B

    C. Seinen, B. Khouider, Improving the Jacobian free Newton-Krylov method for the viscous-plastic sea ice momentum equation.Physica D376–377(2018), 78–93

    2018

  34. [42]

    Y. Shih, C. Mehlmann, M. Losch, G. Stadler, Robust and efficient primal-dual Newton-Krylov solvers for viscous-plastic sea-ice models.J. Comput. Phys.474(2023), Paper No. 111802

    2023

  35. [43]

    Stroeve, D

    J. Stroeve, D. Notz, Changing state of Arctic sea ice across all seasons.Environ. Res. Lett.13(2018), Paper No. 103001

    2018

  36. [44]

    Strong, I.G

    C. Strong, I.G. Rigor, Arctic marginal ice zone trending wider in summer and narrower in winter.Geophys. Res. Lett.40 (2013), 4864–4868

    2013

  37. [45]

    Triebel,Interpolation Theory, Function Spaces, Differential Operators

    H. Triebel,Interpolation Theory, Function Spaces, Differential Operators. North-Holland, 1978

    1978

  38. [46]

    Weis, Operator-valued Fourier multiplier theorems and maximal L p-regularity.Math

    L. Weis, Operator-valued Fourier multiplier theorems and maximal L p-regularity.Math. Ann.319(2001), 735–758

    2001

  39. [47]

    Yaremchuk, G

    M. Yaremchuk, G. Panteleev, On the Jacobian approximation in sea ice models with viscous-plastic rheology.Ocean Model.177(2022), Paper 102078. 32 FELIX BRANDT AND CAROLIN MEHLMANN

    2022

  40. [48]

    Zhang, W.D

    J. Zhang, W.D. Hibler, On an efficient numerical method for modeling sea ice dynamics.J. Geophys. Res.102(1991), 8691–8702. Department of Mathematics, University of California at Berkeley, Berkeley, 94720, CA, USA. Email address:fbrandt@berkeley.edu Institute of Analysis and N...

    1991

Pith tools

Reviewed May 8, 2026 · model on record in the stance chip above.