Recognition: unknown
A mathematical study of an elastic-viscous-plastic sea-ice model with the Kelvin-Voigt rheology
Pith reviewed 2026-05-07 13:17 UTC · model grok-4.3
The pith
Introducing Kelvin-Voigt viscosity into the momentum balance of an elastic-viscous-plastic sea-ice model yields local well-posedness with advection and global well-posedness without it through a new stress tensor bound.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We formulate an elastic-viscous-plastic sea-ice model with Kelvin-Voigt regularisation in the momentum balance and prove the local well-posedness for the model with the advection term and the global well-posedness in the absence of the advection term. A crucial component of the proof is a new L^∞-estimate for the stress tensor which relies on the damping structure. We handle viscosity coefficients without a cutoff from above and prove existence for much less regular initial data.
What carries the argument
The damping structure provided by the Kelvin-Voigt term in the momentum balance, enabling the L^∞ estimate for the stress tensor.
If this is right
- Solutions exist locally in time even when the full advection term is retained in the momentum equation.
- Global-in-time solutions exist when the advection term is omitted from the momentum balance.
- Well-posedness holds for viscosity coefficients that lack any artificial upper bound.
- Existence results apply to initial data of lower regularity than those treated in the authors' prior Voigt-EVP work.
Where Pith is reading between the lines
- The damping-derived stress bound may transfer to other non-Newtonian fluid systems that combine elastic and viscous effects.
- Numerical climate codes could adopt this regularization to bypass cutoff restrictions that currently limit fidelity.
- Long-time stability or attractor properties of the solutions remain open for further analysis.
Load-bearing premise
The Kelvin-Voigt damping term in the momentum balance produces a uniform L^∞ bound on the stress tensor under the given model assumptions.
What would settle it
A concrete initial datum for which the stress tensor norm becomes unbounded in finite time despite the damping term, or for which no solution exists even locally.
read the original abstract
Motivated by the elastic-viscous-plastic (EVP) sea-ice model [E. C. Hunke and J. K. Dukowicz, J. Phys. Oceanogr., 27, 9 (1997), 1849--1867], which is used in large-scale numerical climate simulations, we proposed in [D. W. Boutros, X. Liu, M. Thomas and E. S. Titi, arXiv:2505.03080 (2025)] the use of the inviscid Voigt regularisation for the constitutive (stress-tensor) relation and proved the global well-posedness of the resulting model. The EVP model treats sea ice as a non-Newtonian fluid. In turn, elastic-viscous-plastic solids often involve a Kelvin-Voigt viscosity in terms of the strain rate. Therefore, in the present work we formulate an elastic-viscous-plastic sea-ice model with a Kelvin-Voigt regularisation in terms of the strain rate. In other words, we introduce the Voigt regularisation in the momentum balance rather than in the constitutive relation (for the stress tensor). We then prove the local well-posedness for the Kelvin-Voigt EVP model with the advection term, in the momentum balance, and the global well-posedness in the absence of the advection term (following a very standard approximation in the latter case). A crucial component of the proof of these results, is a new $L^\infty$-estimate for the stress tensor which relies on the damping structure. Note that, both with and without the advection term, we are able to handle the case of viscosity coefficients without a cutoff from above, which remains a major open problem for the closely related Hibler sea-ice model. We are also able to prove the existence of solutions for much less regular initial data compared to our previous paper on the Voigt-EVP model.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper formulates an elastic-viscous-plastic sea-ice model with Kelvin-Voigt regularization placed in the momentum balance (rather than the constitutive relation) and proves local well-posedness of the system including the advection term together with global well-posedness when advection is omitted. Both results rest on a new L^∞ estimate for the stress tensor that exploits the damping structure of the Kelvin-Voigt term; the analysis accommodates viscosity coefficients without an upper cutoff and initial data of lower regularity than in the authors' prior Voigt-EVP work.
Significance. If the claimed L^∞ bound on the stress tensor is valid, the results constitute a concrete advance for the mathematical theory of sea-ice models used in climate simulations. They remove the artificial upper cutoff on viscosity coefficients that remains open for the related Hibler model and lower the regularity threshold on initial data, while the global existence result without advection follows standard approximation techniques.
major comments (1)
- [Proof of the a priori L^∞ estimate for the stress tensor (in the section containing the energy estimates and maximum-prc] The new L^∞ estimate for the stress tensor (central to both Theorem statements on local and global well-posedness) must be shown to control the nonlinear plastic contribution and any transport effects while using only the given constitutive law and the Kelvin-Voigt damping term. The derivation should explicitly confirm that no implicit upper bound on the viscosity coefficients or extra regularity on the data is required, as this is the load-bearing step for the main claims.
minor comments (2)
- [Abstract] The abstract and introduction would benefit from an explicit statement of the precise function spaces in which local and global solutions are obtained.
- [Introduction] A short comparison table or paragraph contrasting the regularity assumptions here with those in the authors' previous Voigt-EVP paper (arXiv:2505.03080) would clarify the improvement.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable comments on our paper. We are pleased that the significance of the results is recognized. Below we provide a point-by-point response to the major comment.
read point-by-point responses
-
Referee: [Proof of the a priori L^∞ estimate for the stress tensor (in the section containing the energy estimates and maximum-prc] The new L^∞ estimate for the stress tensor (central to both Theorem statements on local and global well-posedness) must be shown to control the nonlinear plastic contribution and any transport effects while using only the given constitutive law and the Kelvin-Voigt damping term. The derivation should explicitly confirm that no implicit upper bound on the viscosity coefficients or extra regularity on the data is required, as this is the load-bearing step for the main claims.
Authors: We appreciate the referee's emphasis on the importance of this estimate. The L^∞ bound is derived in the energy estimates section by testing the constitutive relation with the stress tensor itself, leveraging the positive damping term from the Kelvin-Voigt regularization in the momentum equation to control the growth. This approach directly handles the nonlinear plastic term through the structure of the model and accounts for transport effects via integration by parts or commutator estimates that do not require additional regularity. The derivation relies solely on the constitutive law and the damping, without any upper bound on the viscosity coefficients (which is a key feature distinguishing our result from the Hibler model) and with the initial data regularity as stated in the theorems. To address the referee's request for explicit confirmation, we will revise the manuscript by adding a remark immediately following the estimate that summarizes these properties and verifies the absence of implicit assumptions. revision: yes
Circularity Check
No circularity: direct PDE existence proof via a priori estimates
full rationale
The paper establishes local well-posedness (with advection) and global well-posedness (without advection) for the Kelvin-Voigt EVP sea-ice model through standard mathematical techniques: derivation of a new L^∞ bound on the stress tensor from the damping term in the momentum equation, followed by approximation arguments and energy estimates. This chain relies on the constitutive law and the Kelvin-Voigt regularization without any reduction to fitted parameters, self-definitional loops, or load-bearing self-citations. The cited prior work (arXiv:2505.03080) supplies only motivational context for the model choice and is not invoked to justify the current estimates or uniqueness. The derivation is self-contained against external PDE benchmarks and does not rename known results or smuggle ansatzes.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The sea-ice dynamics follow the standard elastic-viscous-plastic constitutive framework with the Kelvin-Voigt term inserted into the momentum equation.
- standard math Standard assumptions on the spatial domain, boundary conditions, and Sobolev-type function spaces hold for the well-posedness analysis.
Reference graph
Works this paper leans on
-
[1]
Adams, R. A. Sobolev Spaces . Academic Press, 1975
1975
-
[2]
Andreu, F., Ballester, C., Caselles, V., and Maz \'o n, J. M. The Dirichlet problem for the total variation flow. Journal of Functional Analysis 180 , 2 (2001), 347--403
2001
-
[3]
Andreu, F., Ballester, C., Caselles, V., and Maz \'o n, J. M. Minimizing total variation flow. Differential Integral Equations 14 , 3 (2001), 321--–360
2001
-
[4]
The classical Kelvin--Voigt problem for incompressible fluids with unknown non-constant density: Existence, uniqueness and regularity
Antontsev, S., De Oliveira, H., and Khompysh, K. The classical Kelvin--Voigt problem for incompressible fluids with unknown non-constant density: Existence, uniqueness and regularity. Nonlinearity 34 , 5 (2021), 3083
2021
-
[5]
The total variation flow in R ^N
Bellettini, G., Caselles, V., and Novaga, M. The total variation flow in R ^N . Journal of Differential Equations 184 , 2 (2002), 475--525
2002
-
[6]
Rigorous analysis of the interaction problem of sea ice with a rigid body
Binz, T., Brandt, F., and Hieber, M. Rigorous analysis of the interaction problem of sea ice with a rigid body. Mathematische Annalen 389 , 1 (2024), 591--625
2024
-
[7]
Interaction of geophysical flows with sea ice dynamics
Binz, T., Brandt, F., and Hieber, M. Interaction of geophysical flows with sea ice dynamics. Nonlinear Differential Equations and Applications NoDEA 33 , 2 (2026), 35
2026
-
[8]
The elastic--viscous--plastic method revisited
Bouillon, S., Fichefet, T., Legat, V., and Madec, G. The elastic--viscous--plastic method revisited. Ocean Modelling 71\/ (2013), 2--12
2013
-
[9]
Bouillon, S., Maqueda, M. A. M., Legat, V., and Fichefet, T. An elastic--viscous--plastic sea ice model formulated on Arakawa B and C grids. Ocean Modelling 27 , 3-4 (2009), 174--184
2009
-
[10]
W., Liu, X., Thomas, M., and Titi, E
Boutros, D. W., Liu, X., Thomas, M., and Titi, E. S. Global well-posedness of the elastic-viscous-plastic sea-ice model with the inviscid Voigt-regularisation . arXiv:2505.03080\/ (2025)
-
[11]
Well-posedness of Hibler’s parabolic-hyperbolic sea ice model
Brandt, F. Well-posedness of Hibler’s parabolic-hyperbolic sea ice model. Journal of Evolution Equations 25 , 3 (2025), 1--31
2025
-
[12]
Rigorous analysis and dynamics of Hibler’s sea ice model
Brandt, F., Disser, K., Haller-Dintelmann, R., and Hieber, M. Rigorous analysis and dynamics of Hibler’s sea ice model. Journal of Nonlinear Science 32 , 4 (2022), 50
2022
-
[13]
Time periodic solutions to Hibler’s sea ice model
Brandt, F., and Hieber, M. Time periodic solutions to Hibler’s sea ice model. Nonlinearity 36 , 6 (2023), 3109
2023
-
[14]
Well posedness of the regularized-Hibler model of sea-ice dynamics
Chatta, S., and Khouider, B. Well posedness of the regularized-Hibler model of sea-ice dynamics. Authorea preprint\/ (2025). https://doi.org/10.22541/essoar.175157581.19836982/v1
-
[15]
Linear well posedness of regularized equations of sea-ice dynamics
Chatta, S., Khouider, B., and Kesri, M. Linear well posedness of regularized equations of sea-ice dynamics. Journal of Mathematical Physics 64 , 5 (2023)
2023
-
[16]
About lifespan of regular solutions of equations related to viscoelastic fluids
Chemin, J.-Y., and Masmoudi, N. About lifespan of regular solutions of equations related to viscoelastic fluids. SIAM Journal on Mathematical Analysis 33 , 1 (2001), 84--112
2001
-
[17]
Chen, N., Deng, Q., and Stechmann, S. N. Superfloe parameterization with physics constraints for uncertainty quantification of sea ice floes. SIAM/ASA Journal on Uncertainty Quantification 10 , 4 (2022), 1384--1409
2022
-
[18]
Note on global regularity for 2D Oldroyd-B fluids with diffusive stress
Constantin, P., and Kliegl, M. Note on global regularity for 2D Oldroyd-B fluids with diffusive stress. Archive for Rational Mechanics and Analysis 206\/ (2012), 725–--740
2012
-
[19]
Magnetic relaxation of a Voigt--MHD system
Constantin, P., and Pasqualotto, F. Magnetic relaxation of a Voigt--MHD system. Communications in Mathematical Physics 402 , 2 (2023), 1931--1952
2023
-
[20]
High Reynolds number and high Weissenberg number Oldroyd-B model with dissipation
Constantin, P., Wu, J., Zhao, J., and Zhu, Y. High Reynolds number and high Weissenberg number Oldroyd-B model with dissipation. Journal of Evolution Equations 21\/ (2021), 2787--2806
2021
-
[21]
A., and Thorndike, A
Coon, M., Maykut, G., Pritchard, R., Rothrock, D. A., and Thorndike, A. S. Modeling the pack ice as an elastic-plastic material. AIDJEX Bulletin , 24 (1974), 1--105
1974
-
[22]
A Maxwell elasto-brittle rheology for sea ice modelling
Dansereau, V., Weiss, J., Saramito, P., and Lattes, P. A Maxwell elasto-brittle rheology for sea ice modelling. The Cryosphere 10 , 3 (2016), 1339--1359
2016
-
[23]
N., and Chen, N
Deng, Q., Stechmann, S. N., and Chen, N. Particle-continuum multiscale modeling of sea ice floes. Multiscale Modeling & Simulation 22 , 1 (2024), 230--255
2024
-
[24]
On the singular limit in Hibler's sea ice model
Denk, R., Gmeineder, F., and Hieber, M. On the singular limit in Hibler's sea ice model. arXiv:2511.09327\/ (2025)
-
[25]
Global existence and uniqueness for Hibler's visco-plastic sea-ice model
Dingel, S., and Disser, K. Global existence and uniqueness for Hibler's visco-plastic sea-ice model. arXiv:2508.16537\/ (2025)
-
[26]
M., and Rousset, F
Elgindi, T. M., and Rousset, F. Global regularity for some Oldroyd-B type models. Communications on Pure and Applied Mathematics 68 , 11 (2015), 2005--2021
2015
-
[27]
Very singular diffusion equations: second and fourth order problems
Giga, M.-H., and Giga, Y. Very singular diffusion equations: second and fourth order problems. Japan Journal of Industrial and Applied Mathematics 27 , 3 (2010), 323--345
2010
-
[28]
Very singular diffusion equations
Giga, M.-H., Giga, Y., and Kobayashi, R. Very singular diffusion equations. In Taniguchi Conference on Mathematics Nara'98\/ (2001), vol. 31, Mathematical Society of Japan, pp. 93--126
2001
-
[29]
Existence results for the flow of viscoelastic fluids with a differential constitutive law
Guillop \'e , C., and Saut, J.-C. Existence results for the flow of viscoelastic fluids with a differential constitutive law. Nonlinear Analysis: Theory, Methods & Applications 15 , 9 (1989), 849--869
1989
-
[30]
D., Feltham, D
Heorton, H. D., Feltham, D. L., and Tsamados, M. Stress and deformation characteristics of sea ice in a high-resolution, anisotropic sea ice model. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences 376 , 2129 (2018), 20170349
2018
-
[31]
Hibler, W. D. A dynamic thermodynamic sea ice model. Journal of Physical Oceanography 9 , 4 (1979), 815--846
1979
-
[32]
Hunke, E. C. Viscous--plastic sea ice dynamics with the EVP model: Linearization issues. Journal of Computational Physics 170 , 1 (2001), 18--38
2001
-
[33]
C., and Dukowicz, J
Hunke, E. C., and Dukowicz, J. K. An elastic--viscous--plastic model for sea ice dynamics. Journal of Physical Oceanography 27 , 9 (1997), 1849--1867
1997
-
[34]
C., and Dukowicz, J
Hunke, E. C., and Dukowicz, J. K. The elastic--viscous--plastic sea ice dynamics model in general orthogonal curvilinear coordinates on a sphere—incorporation of metric terms. Monthly Weather Review 130 , 7 (2002), 1848--1865
2002
-
[35]
2D Voigt Boussinesq equations
Ignatova, M. 2D Voigt Boussinesq equations. Journal of Mathematical Fluid Mechanics 26 , 1 (2024), 15
2024
-
[36]
F., Hibler III, W
Ip, C. F., Hibler III, W. D., and Flato, G. M. On the effect of rheology on seasonal sea-ice simulations. Annals of Glaciology 15\/ (1991), 17--25
1991
-
[37]
V., Danilov, S., Sidorenko, D., Hutter, N., Losch, M., Goessling, H., Rakowsky, N., Scholz, P., Sein, D., Wang, Q., et al
Koldunov, N. V., Danilov, S., Sidorenko, D., Hutter, N., Losch, M., Goessling, H., Rakowsky, N., Scholz, P., Sein, D., Wang, Q., et al. Fast EVP solutions in a high-resolution sea ice model. Journal of Advances in Modeling Earth Systems 11 , 5 (2019), 1269--1284
2019
-
[38]
Kreyscher, M., Harder, M., Lemke, P., and Flato, G. M. Results of the sea ice model intercomparison project: Evaluation of sea ice rheology schemes for use in climate simulations. Journal of Geophysical Research: Oceans 105 , C5 (2000), 11299--11320
2000
-
[39]
Kupferman, R., Mangoubi, C., and Titi, E. S. A Beale-Kato-Madja breakdown criterion for an Oldroyd-B fluid in the creeping flow regime. Communications in Mathematical Sciences 6 , 1 (2008), 235--256
2008
- [40]
-
[41]
Larios, A., and Titi, E. S. On the higher-order global regularity of the inviscid Voigt-regularization of three-dimensional hydrodynamic models. Discrete and Continuous Dynamical Systems - Series B 14 , 2 (2010), 603--627
2010
-
[42]
Larios, A., and Titi, E. S. Higher-order global regularity of an inviscid Voigt-regularization of the three-dimensional inviscid resistive magnetohydrodynamic equations. Journal of Mathematical Fluid Mechanics 16 , 1 (2014), 59--76
2014
-
[43]
A., Tremblay, B., Holland, D
Lemieux, J.-F., Knoll, D. A., Tremblay, B., Holland, D. M., and Losch, M. A comparison of the Jacobian-free Newton--Krylov method and the EVP model for solving the sea ice momentum equation with a viscous-plastic formulation: A serial algorithm study. Journal of Computational Physics 231 , 17 (2012), 5926--5944
2012
-
[44]
A First Course in Sobolev Spaces
Leoni, G. A First Course in Sobolev Spaces . American Mathematical Society, 2017
2017
-
[45]
The Drift of Sea Ice
Lepp \"a ranta, M. The Drift of Sea Ice . Springer, 2011
2011
-
[46]
On the initial-boundary value problem of the incompressible viscoelastic fluid system
Lin, F., and Zhang, P. On the initial-boundary value problem of the incompressible viscoelastic fluid system. Communications on Pure and Applied Mathematics 61 , 4 (2008), 539--558
2008
-
[47]
Global solutions for some Oldroyd models of non-Newtonian flows
Lions, P., and Masmoudi, N. Global solutions for some Oldroyd models of non-Newtonian flows. Chinese Annals of Mathematics 21 , 02 (2000), 131--146
2000
-
[48]
Liu, X., Thomas, M., and Titi, E. S. Well-posedness of Hibler’s dynamical sea-ice model. Journal of Nonlinear Science 32 , 4 (2022), 49
2022
- [49]
-
[50]
On solving the momentum equations of dynamic sea ice models with implicit solvers and the elastic--viscous--plastic technique
Losch, M., and Danilov, S. On solving the momentum equations of dynamic sea ice models with implicit solvers and the elastic--viscous--plastic technique. Ocean Modelling 41\/ (2012), 42--52
2012
-
[51]
A parallel Jacobian-free Newton--Krylov solver for a coupled sea ice-ocean model
Losch, M., Fuchs, A., Lemieux, J.-F., and Vanselow, A. A parallel Jacobian-free Newton--Krylov solver for a coupled sea ice-ocean model. Journal of Computational Physics 257\/ (2014), 901--911
2014
-
[52]
Sea-ice dynamics on triangular grids
Mehlmann, C., and Korn, P. Sea-ice dynamics on triangular grids. Journal of Computational Physics 428\/ (2021), 110086
2021
-
[53]
The uniqueness and global solvability of boundary-value problems for the equations of motion for aqueous solutions of polymers
Oskolkov, A. The uniqueness and global solvability of boundary-value problems for the equations of motion for aqueous solutions of polymers. Journal of Mathematical Sciences 8 , 4 (1977), 427--455
1977
-
[54]
Nonlocal problems for the equations of Kelvin-Voight fluids and their -approximations
Oskolkov, A. Nonlocal problems for the equations of Kelvin-Voight fluids and their -approximations. Journal of Mathematical Sciences 87\/ (1997), 3393--3408
1997
-
[55]
On the dynamics of grounded shallow ice sheets: Modeling and analysis
Piersanti, P., and Temam, R. On the dynamics of grounded shallow ice sheets: Modeling and analysis. Advances in Nonlinear Analysis 12 , 1 (2023), 20220280
2023
-
[56]
Improving the Jacobian free Newton--Krylov method for the viscous--plastic sea ice momentum equation
Seinen, C., and Khouider, B. Improving the Jacobian free Newton--Krylov method for the viscous--plastic sea ice momentum equation. Physica D: Nonlinear Phenomena 376\/ (2018), 78--93
2018
-
[57]
Showalter, R. E. Well-posed problems for a partial differential equation of order 2m+1 . SIAM Journal on Mathematical Analysis 1 , 2 (1970), 214--231
1970
-
[58]
Showalter, R. E. Nonlinear degenerate evolution equations and partial differential equations of mixed type. SIAM Journal on Mathematical Analysis 6 , 1 (1975), 25--42
1975
-
[59]
E., and Ting, T
Showalter, R. E., and Ting, T. W. Pseudoparabolic partial differential equations. SIAM Journal on Mathematical Analysis 1 , 1 (1970), 1--26
1970
-
[60]
A viscous continuum theory of sea ice motion based on stochastic floe dynamics
Toppaladoddi, S. A viscous continuum theory of sea ice motion based on stochastic floe dynamics. Journal of Fluid Mechanics 1014\/ (2025), A6
2025
-
[61]
L., and Wilchinsky, A
Tsamados, M., Feltham, D. L., and Wilchinsky, A. V. Impact of a new anisotropic rheology on simulations of Arctic sea ice. Journal of Geophysical Research: Oceans 118 , 1 (2013), 91--107
2013
-
[62]
V., and Feltham, D
Wilchinsky, A. V., and Feltham, D. L. Modelling the rheology of sea ice as a collection of diamond-shaped floes. Journal of Non-Newtonian Fluid Mechanics 138 , 1 (2006), 22--32
2006
-
[63]
On an efficient numerical method for modeling sea ice dynamics
Zhang, J., and Hibler III, W. On an efficient numerical method for modeling sea ice dynamics. Journal of Geophysical Research: Oceans 102 , C4 (1997), 8691--8702
1997
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.