pith. machine review for the scientific record. sign in

arxiv: 2604.26295 · v1 · submitted 2026-04-29 · 🧮 math.AP · physics.ao-ph· physics.geo-ph

Recognition: unknown

A mathematical study of an elastic-viscous-plastic sea-ice model with the Kelvin-Voigt rheology

Authors on Pith no claims yet

Pith reviewed 2026-05-07 13:17 UTC · model grok-4.3

classification 🧮 math.AP physics.ao-phphysics.geo-ph
keywords elastic-viscous-plastic sea-ice modelKelvin-Voigt rheologywell-posednessstress tensor estimateadvection termviscosity regularizationnon-Newtonian fluid
0
0 comments X

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.

The authors study a version of the elastic-viscous-plastic sea-ice model that incorporates Kelvin-Voigt rheology directly in the momentum equation. They establish local-in-time existence and uniqueness of solutions when the advection term is present and global existence when it is omitted. The key technical step is deriving a new supremum bound on the stress tensor that exploits the damping provided by the viscosity term. This approach succeeds without imposing an artificial upper limit on the viscosity coefficients and works for initial conditions that are less smooth than those required in related prior results. Such mathematical guarantees are relevant because these models underpin large-scale simulations of sea ice in climate studies.

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

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

  • 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.

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

1 major / 2 minor

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)
  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)
  1. [Abstract] The abstract and introduction would benefit from an explicit statement of the precise function spaces in which local and global solutions are obtained.
  2. [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

1 responses · 0 unresolved

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
  1. 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

0 steps flagged

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

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the standard EVP sea-ice equations plus the added Kelvin-Voigt damping term in momentum balance, together with classical PDE theory for existence; no free parameters or new physical entities are introduced.

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.
    The model is constructed by modifying the established EVP equations as described in the abstract.
  • standard math Standard assumptions on the spatial domain, boundary conditions, and Sobolev-type function spaces hold for the well-posedness analysis.
    Implicit in any local/global existence proof for evolutionary PDE systems.

pith-pipeline@v0.9.0 · 5679 in / 1504 out tokens · 115690 ms · 2026-05-07T13:17:08.354816+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

63 extracted references · 6 canonical work pages

  1. [1]

    Adams, R. A. Sobolev Spaces . Academic Press, 1975

  2. [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

  3. [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

  4. [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

  5. [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

  6. [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

  7. [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

  8. [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

  9. [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

  10. [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. [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

  12. [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

  13. [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

  14. [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. [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)

  16. [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

  17. [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

  18. [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

  19. [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

  20. [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

  21. [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

  22. [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

  23. [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

  24. [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. [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. [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

  27. [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

  28. [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

  29. [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

  30. [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

  31. [31]

    Hibler, W. D. A dynamic thermodynamic sea ice model. Journal of Physical Oceanography 9 , 4 (1979), 815--846

  32. [32]

    Hunke, E. C. Viscous--plastic sea ice dynamics with the EVP model: Linearization issues. Journal of Computational Physics 170 , 1 (2001), 18--38

  33. [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

  34. [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

  35. [35]

    2D Voigt Boussinesq equations

    Ignatova, M. 2D Voigt Boussinesq equations. Journal of Mathematical Fluid Mechanics 26 , 1 (2024), 15

  36. [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

  37. [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

  38. [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

  39. [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

  40. [40]

    Larios, A., Lunasin, E., and Titi, E. S. Global well-posedness for the 2D Boussinesq system without heat diffusion and with either anisotropic viscosity or inviscid Voigt- regularization. arXiv:1010.5024\/ (2010)

  41. [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

  42. [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

  43. [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

  44. [44]

    A First Course in Sobolev Spaces

    Leoni, G. A First Course in Sobolev Spaces . American Mathematical Society, 2017

  45. [45]

    The Drift of Sea Ice

    Lepp \"a ranta, M. The Drift of Sea Ice . Springer, 2011

  46. [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

  47. [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

  48. [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

  49. [49]

    Liu, X., Thomas, M., and Titi, E. S. Plastic limit of a viscoplastic Burgers equation -- a toy model for sea-ice dynamics. arXiv:2601.06489\/ (2026)

  50. [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

  51. [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

  52. [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

  53. [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

  54. [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

  55. [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

  56. [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

  57. [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

  58. [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

  59. [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

  60. [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

  61. [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

  62. [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

  63. [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