pith. sign in

arxiv: 2606.10545 · v1 · pith:IJDZWDQ4new · submitted 2026-06-09 · ⚛️ physics.class-ph

A phase-field modeling approach to sea-ice fracturing

Pith reviewed 2026-06-27 10:47 UTC · model grok-4.3

classification ⚛️ physics.class-ph
keywords phase-field modelingsea-ice fracturingfracture propagationdouble-well potentialGriffith theoryspectral methodgranular transitionbody forces
0
0 comments X

The pith

A phase-field model with double-well energy and overdamped dynamics reproduces tensile and shear fractures in sea ice, matching Griffith's linear crack-speed scaling.

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

The paper tests whether a phase-field approach can simulate how sea ice fractures under wind and current forces and transitions toward a granular state. It combines a double-well free-energy term with an overdamped displacement field and solves the equations spectrally in Fourier space while including representative body forces. Validation on an inclusion under tension and on simple shear setups shows that the model produces expected fracture patterns and recovers the linear relation between crack speed and load that Griffith's theory predicts. This matters for improving how regional and global models represent sea-ice mechanics and dynamics.

Core claim

The phase-field framework, built from a double-well free-energy formulation and an overdamped displacement response and solved by a spectral method, generates fracture patterns and displacement fields in sea ice that capture essential features of tensile and shear crack propagation, including the linear scaling of crack speed with applied load required by Griffith's theory.

What carries the argument

Phase-field framework that combines a double-well free-energy formulation with an overdamped displacement response, solved using a spectral method in Fourier space while incorporating body forces.

If this is right

  • The model remains computationally tractable for nonlinear double-well problems while including body forces typical of sea-ice forcing.
  • Fracture patterns in the tensile-inclusion benchmark and the plane-shear and Couette configurations reproduce key features of both tensile and shear propagation.
  • The framework can represent the transition from a continuous ice cover to a granular medium inside shear bands.
  • The recovered linear crack-speed scaling with load is consistent with Griffith's theory for the tested configurations.

Where Pith is reading between the lines

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

  • Embedding this local fracturing description into larger-scale sea-ice dynamics models could improve representations of deformation and ridging.
  • The same phase-field setup might be extended to include thermal or viscous effects without changing the core spectral solver.
  • Direct comparison of simulated shear-band widths against satellite observations of Arctic ice deformation would provide an independent test.

Load-bearing premise

The double-well free-energy formulation combined with an overdamped displacement response is assumed to adequately capture the physics of sea-ice fracturing and the continuous-to-granular transition under representative body forces.

What would settle it

A laboratory or field measurement showing that crack propagation speed in sea ice does not increase linearly with applied load, or that fracture patterns in simple shear deviate markedly from the simulated fields, would falsify the central claim.

read the original abstract

The thin ice that covers the polar oceans is a complex geomaterial that is constantly stressed and fractured by winds and ocean currents. In the central Arctic, this forcing produces deformations in the form of shear bands, within which individual ice plates detach, locally generating a granular medium. Capturing this transition from a continuous to a granular sea-ice cover has implications for the adequate representation of the mechanical and dynamical behavior of sea-ice in regional and large-scale models used for operational and climate prediction purposes. Our work investigates the feasibility of a phase-field approach to capture this granularization processes and focuses on fracture propagation in the material. The model combines a double-well free-energy formulation with an overdamped displacement response. The governing equations are solved using a spectral method in Fourier space. The implementation accounts for body forces, representative of the main forcings on sea-ice, and remains computationally tractable despite the highly nonlinear character of the double-well energy formulation. We first validate the framework against a benchmark problem: the opening of an inclusion embedded into an elastic matrix under tensile loading. Then, additional simple shear configurations are investigated: an inclusion solicited under plane shear and a cylindrical Couette experiment, for which the analytical solution of the displacement field is known. The resulting fracture patterns and displacement fields demonstrate that our phase-field framework captures key features of tensile and shear fracture propagation, including the linear scaling between crack speed and applied load predicted by the Griffith's theory.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 3 minor

Summary. The manuscript presents a phase-field model for sea-ice fracturing that combines a double-well free-energy formulation with an overdamped displacement response, solved via a Fourier spectral method that incorporates body forces. It validates the framework on three 2D benchmark problems—an inclusion under tensile loading, an inclusion under plane shear, and a cylindrical Couette flow—claiming that the resulting fracture patterns and displacement fields reproduce key features of tensile and shear propagation, including the linear scaling of crack speed with applied load predicted by Griffith theory. The work is positioned as a feasibility study for capturing the continuous-to-granular transition in sea ice.

Significance. If the validations hold, the approach offers a computationally tractable route to represent fracturing and the emergence of granular behavior in sea-ice models, which is relevant for improving mechanical representations in regional and climate-scale simulations. Credit is due for the explicit inclusion of body forces representative of wind and current forcing and for demonstrating that the spectral solver remains tractable despite the nonlinearity of the double-well term. The current scope remains limited to simplified 2D elastic benchmarks without direct sea-ice observations or inertial dynamics.

major comments (2)
  1. [Validation section (Couette and shear benchmarks)] Validation section (Couette and shear benchmarks): the manuscript states that the displacement field matches the known analytical solution, but does not report any quantitative error measure (e.g., L2 norm or pointwise maximum deviation) between the simulated and analytical fields; without this, the accuracy of the elastic response underlying the fracture patterns cannot be assessed.
  2. [Results on Griffith scaling] Results on Griffith scaling: the abstract and results claim reproduction of the linear scaling between crack speed and applied load, yet the text provides neither the procedure used to extract crack speed from the phase-field evolution nor the set of load values and corresponding speeds; this quantitative support is load-bearing for the central claim that the model captures Griffith features.
minor comments (3)
  1. [Methods] The definition and normalization of the phase-field variable (order parameter) and the specific form of the double-well potential are introduced without an explicit equation reference in the methods; adding the governing equations with numbered labels would improve clarity.
  2. [Figures] Figure captions for the fracture patterns do not indicate the color scale for the phase-field variable or the time at which the snapshot is taken; this hinders direct comparison with the claimed propagation features.
  3. [Introduction] The introduction cites general sea-ice mechanics literature but omits key phase-field fracture references (e.g., Bourdin et al. 2000 or subsequent spectral implementations); adding these would contextualize the numerical approach.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment and constructive feedback on our manuscript. We address the major comments point by point below. Both points raised are valid and we will incorporate the requested information in a revised version.

read point-by-point responses
  1. Referee: Validation section (Couette and shear benchmarks): the manuscript states that the displacement field matches the known analytical solution, but does not report any quantitative error measure (e.g., L2 norm or pointwise maximum deviation) between the simulated and analytical fields; without this, the accuracy of the elastic response underlying the fracture patterns cannot be assessed.

    Authors: We agree that quantitative error measures would strengthen the validation. In the revised manuscript, we will compute and report the L2 norm and maximum pointwise deviation between the simulated displacement fields and the analytical solutions for both the plane shear and Couette flow benchmarks. revision: yes

  2. Referee: Results on Griffith scaling: the abstract and results claim reproduction of the linear scaling between crack speed and applied load, yet the text provides neither the procedure used to extract crack speed from the phase-field evolution nor the set of load values and corresponding speeds; this quantitative support is load-bearing for the central claim that the model captures Griffith features.

    Authors: We acknowledge that the procedure for extracting crack speed and the specific load-speed data pairs are not detailed in the current text. We will add a description of the crack speed extraction method (based on tracking the position of the phase-field interface over time) and include a table or figure showing the applied loads and corresponding measured crack speeds to support the linear scaling claim. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper presents a phase-field model for sea-ice fracturing and validates it on independent benchmark problems (tensile inclusion, plane shear, Couette flow) whose analytical displacement solutions are known from linear elasticity and are external to the model. The reproduction of Griffith-theory scaling is a consistency check against a pre-existing analytical prediction rather than a fit or self-referential step. No load-bearing derivation reduces to a fitted parameter renamed as prediction, a self-citation chain, or an ansatz smuggled via prior work by the same authors. The central claim therefore rests on external benchmarks and does not collapse to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based on abstract only; the model rests on standard phase-field assumptions and numerical methods without explicit free parameters or invented entities listed.

axioms (1)
  • domain assumption Double-well free-energy formulation with overdamped displacement accurately represents sea-ice fracture physics.
    Invoked as the core of the governing equations in the abstract.

pith-pipeline@v0.9.1-grok · 5819 in / 1167 out tokens · 26188 ms · 2026-06-27T10:47:38.367585+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

68 extracted references · 9 canonical work pages

  1. [1]

    2025 , note =

    NASA Worldview , author =. 2025 , note =

  2. [2]

    2025 , note =

    SASIP Project , author =. 2025 , note =

  3. [3]

    Phase-Field Model of Mode III Dynamic Fracture

    Alain Karma and Daniel Kessler and Herbert Levine. Phase-Field Model of Mode III Dynamic Fracture. Physical Review Letters. 2001

  4. [4]

    Thorndike and R

    A.S. Thorndike and R. Colony. Arctic Ocean Buoy Program. DTIC. 1980

  5. [5]

    Hughes and Chad M

    Yu-Sheng Lo and Thomas J.R. Hughes and Chad M. Landis. Phase-field fracture modeling for large structures. Journal of the Mechanics and Physics of solids. 2023

  6. [6]

    Drift, Deformation and Fracture of Sea Ice – A perspective across scales

    Jérôme Weiss. Drift, Deformation and Fracture of Sea Ice – A perspective across scales. 2013

  7. [7]

    Fracture Mechanics: Fundamentals and Applications

    Anderson, T.L. Fracture Mechanics: Fundamentals and Applications. 2005

  8. [8]

    Deformation of the Arctic Ocean Sea Ice Cover between November 1996 and April 1997: A Qualitative Survey

    Kwok, Ronald. Deformation of the Arctic Ocean Sea Ice Cover between November 1996 and April 1997: A Qualitative Survey. 2001

  9. [9]

    Hegyi, B. M. and Taylor, P. C. The unprecedented 2016–2017 Arctic sea ice growth season: The crucial role of atmospheric rivers and longwave fluxes. Geophysical Research Letters,. 2018

  10. [10]

    Journal of Geophysical Research: Oceans , doi =

    Michael Steele and Jinlun Zhang and Drew Rothrock and Harry Stern , title =. Journal of Geophysical Research: Oceans , doi =

  11. [11]

    Arctic sea ice velocity field: General circulation and turbulent‐like fluctuations

    Pierre Rampal. Arctic sea ice velocity field: General circulation and turbulent‐like fluctuations. Journal of Geophysical Research. 2009

  12. [12]

    Journal of Geophysical Research: Solid Earth

    Cohesion versus friction in controlling the long‐term strength of a self‐healing experimental fault. Journal of Geophysical Research: Solid Earth. 2016

  13. [13]

    Surveys in Geophysics

    Scaling of fracture and faulting of ice on earth. Surveys in Geophysics. 2002

  14. [14]

    Micro-Seismic Monitoring of a Shear Fault Within a Floating Ice Plate

    Cédric Lachaud and David Marsan and Maurine Montagnat and Jérôme Weiss and Ludovic Moreau and Florent Gimbert. Micro-Seismic Monitoring of a Shear Fault Within a Floating Ice Plate. Journal of Geophysical Research: Solid Earth. 2019

  15. [15]

    A Dynamic Thermodynamic Sea Ice Model

    Hibler, W.D. A Dynamic Thermodynamic Sea Ice Model. Journal of Physical Oceanography. 1979

  16. [16]

    and Dukowicz, J.K

    Hunke, E.C. and Dukowicz, J.K. An Elastic–Viscous–Plastic Model for Sea Ice Dynamics. Journal of physical Oceanography. 1997

  17. [17]

    and Bouillon, S

    Girard, L. and Bouillon, S. and Weiss, J. and Amitrano, D. and Fichefet, T. and Legat, V. A new modeling framework for sea-ice mechanics based on elasto-brittle rheology. Annals of Glaciology. 2011

  18. [18]

    Fracture

    Irwin, G.R. Fracture. In: Handbuch der Physik, Band VI, Elastizität und Plastizität. 1958

  19. [19]

    High-Order C onservative R emapping with a posteriori MOOD stabilization on polygonal meshes

    Blanchard, G and Loubere, R. High-Order C onservative R emapping with a posteriori MOOD stabilization on polygonal meshes. 2015

  20. [20]

    Global Fiducials Library Image Gallery

    USGS. Global Fiducials Library Image Gallery. 2014

  21. [21]

    International Arctic Buoy Program (IABP)

    NCEI NOAA. International Arctic Buoy Program (IABP)

  22. [22]

    A New Brittle Rheology and Numerical Framework for Large-Scale Sea-Ice Models

    Einar Olason and Guillaume Boutin and Anton Korosov and Pierre Rampal and Timothy Williams and Madlen Kimmritz and Véronique Dansereau and Abdoulaye Samaké. A New Brittle Rheology and Numerical Framework for Large-Scale Sea-Ice Models. Journal of Advances in Modeling Earth Systems. 2022

  23. [23]

    The variational approach to fracture: an example of fruitful collaboration between mechanics and applied mathematics

    Jean-Jacques Marigo. The variational approach to fracture: an example of fruitful collaboration between mechanics and applied mathematics. Comptes Rendus Mécanique. 2023

  24. [24]

    and Haapala, J

    Uusinoka, M. and Haapala, J. and Åström, J. and Lensu, M. and Polojärvi, A. Scale invariance in kilometer-scale sea ice deformation. The Cryosphere. 2025

  25. [25]

    Analysis of Stresses and Strains Near the End of a Crack Traversing a Plate

    Irwin, G.R. Analysis of Stresses and Strains Near the End of a Crack Traversing a Plate. Journal of applied mechanics. 1957

  26. [26]

    and Kay, E

    Goosse, H. and Kay, E. K and Armour, K. C. and Bodas-Salcedo, A. and Chepfer, H. and Docquier, D. and Jonko, A. and Kushner, P. and Lecomte, O. and Massonnet, F. and Park, H.-S. and Pithan, F. and Svensson, G. and Vancoppenolle, M. Quantifying climate feedbacks in polar regions. Nature communications. 2018

  27. [27]

    and Jayant Shah

    Mumford, David. and Jayant Shah. Optimal Approximations by Piecewise Smooth Functions and Associated Variational Problems. Communications on Pure and Applied Mathematics,. 1989

  28. [28]

    and Rappel, W.-J

    Karma, A. and Rappel, W.-J. Quantitative phase-field modeling of dendritic growth in two and three dimensions. Physical Review. E, Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics. 1998

  29. [29]

    and Levine, H

    Collins, J. and Levine, H. Diffuse interface model of diffusion-limited crystal growth. Phys. Rev. B. 1985

  30. [30]

    A phase-field method for modeling cracks with frictional contact

    Fan Fei and Jinhyun Choo. A phase-field method for modeling cracks with frictional contact. International Journal for Numerical methods in engineering. 2019

  31. [31]

    A phase-field model of frictional shear fracture in geologic materials

    Fan Fei and Jinhyun Choo. A phase-field model of frictional shear fracture in geologic materials. Computer methods in applied mechanics and engineering. 2020

  32. [32]

    Double-phase-field formulation for mixed-mode fracture in rocks

    Fan Fei and Jinhyun Choo. Double-phase-field formulation for mixed-mode fracture in rocks. Computer Methods in Applied Mechanics and Engineering , volume =. 2021 , issn =. doi:https://doi.org/10.1016/j.cma.2020.113655 , url =

  33. [33]

    Elbanna and Jinhyun Choo , title =

    Fan Fei and Md Shumon Mia and Ahmed E. Elbanna and Jinhyun Choo , title =. 2025 , howpublished =

  34. [34]

    High-accuracy phase-field models for brittle fracture based on a new family of degradation functions

    Juan Michael Sargado. High-accuracy phase-field models for brittle fracture based on a new family of degradation functions. Journal of the Mechanics and Physics of Solids. 2018

  35. [35]

    A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits

    Christian Miehe and Martina Hofacker and Fabian Welschinger. A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits. Computer methods in applied mechanics and engineering. 2010

  36. [36]

    van Dijk and J.J

    N.P. van Dijk and J.J. Espadas-Escalante and P. Isaksson. Strain energy density decompositions in phase-field fracture theories for orthotropy and anisotropy. International Journal of Solids and Structures. 2020

  37. [37]

    Regularized formulation of the variational brittle fracture with unilateral contact: Numerical experiments

    Hanen Amor and Jean-Jacques Marigo and Corrado Maurini. Regularized formulation of the variational brittle fracture with unilateral contact: Numerical experiments. Journal of Mechanics and Physics of Solids. 2009

  38. [38]

    IPCC climate models do not capture Arctic sea ice drift acceleration : consequences in terms of projected sea ice thinning and decline

    Pierre Rampal and Jérôme Weiss and Clotilde Dubois and Jean Michel Campin. IPCC climate models do not capture Arctic sea ice drift acceleration : consequences in terms of projected sea ice thinning and decline. Journal of geophysical research. 2011

  39. [39]

    Dinh and D

    H. Dinh and D. Giannakis and J. Slawinska and G. Stadler. Phase-field models of floe fracture in sea ice. EGUsphere. 2022

  40. [40]

    Laws of crack motion and phase-field models of fracture

    Vincent Hakim and Alain Karma. Laws of crack motion and phase-field models of fracture. Journal of the Mechanics and Physics of Solids. 2009. doi:https://doi.org/10.1016/j.jmps.2008.10.012

  41. [41]

    A Maxwell elasto-brittle rheology for sea ice modelling

    Véronique Dansereau and Jérôme Weiss and Pierre Saramito and Philippe Lattes. A Maxwell elasto-brittle rheology for sea ice modelling. The Cryosphere. 2016

  42. [42]

    Sea ice rheology

    Daniel Feltham. Sea ice rheology. Annual Review of Fluid Mechanics. 2008. doi:https://doi.org/10.1146/annurev.fluid.40.111406.102151

  43. [43]

    Phase field modeling of brittle compressive-shear fractures in rock-like materials: A new driving force and a hybrid formulation

    Shuwei Zhou and Xiaoying Zhuang and Timon Rabczuk. Phase field modeling of brittle compressive-shear fractures in rock-like materials: A new driving force and a hybrid formulation. Computer methods in applied mechanics and engineering

  44. [44]

    Spatschek and E

    R. Spatschek and E. Brener and A. Karma. Phase field modeling of crack propagation. Philosophical magazine. 2010

  45. [45]

    Phase field theory in materials physics, the holograph equation

    Pierre Galenko. Phase field theory in materials physics, the holograph equation. 2024

  46. [46]

    The phenomena of rupture and flow in solids

    Alan Arnold Griffith. The phenomena of rupture and flow in solids. Philosophical Transactions of the Royal Society of London. 1921

  47. [47]

    Ground state structures in ordered binary alloys with second neighbor interactions

    S.M Allen and J.W Cahn. Ground state structures in ordered binary alloys with second neighbor interactions. Acta Metallurgica. 1972

  48. [48]

    Physical Review Letters , volume =

    Dynamic Instabilities of Fracture under Biaxial Strain Using a Phase Field Model , author =. Physical Review Letters , volume =. 2004 , month =. doi:10.1103/PhysRevLett.93.105504 , url =

  49. [49]

    Physical Reviex Letters , volume =

    Unsteady Crack Motion and Branching in a Phase-Field Model of Brittle Fracture , author =. Physical Reviex Letters , volume =. 2004 , month =. doi:10.1103/PhysRevLett.92.245510 , url =

  50. [50]

    A phase field model for fractures in ice shelves

    Sondershaus, Rabea and Humbert, Angelika and Müller, Ralf. A phase field model for fractures in ice shelves. Proceedings in Applied Mathematics and Mechanics. 2023

  51. [51]

    Compressive shear faults within arctic sea ice: Fracture on scales large and small

    Erland Schulson. Compressive shear faults within arctic sea ice: Fracture on scales large and small. Journal of Geophysical Research. 2004

  52. [52]

    Creep and fracture of ice

    Erland Schulson & Paul Duval. Creep and fracture of ice. 2009

  53. [53]

    Freund, L. B. Dynamic fracture mechanics. 1998

  54. [54]

    Swell and sea in the emerging Arctic Ocean

    Jim Thomson and Erick Rogers. Swell and sea in the emerging Arctic Ocean. Geophysical Research Letters. 2014

  55. [55]

    Landau and E

    L.D. Landau and E. M. Lifshitz. Theory of Elasticity: volume 7. Vol. 7. 2012

  56. [56]

    A climate data record of year-round global sea-ice drift from the EUMETSAT Ocean and Sea Ice Satellite Application Facility (OSI SAF)

    Thomas Lavergne and Emily Down. A climate data record of year-round global sea-ice drift from the EUMETSAT Ocean and Sea Ice Satellite Application Facility (OSI SAF). Earth System Science Data. 2023

  57. [57]

    Journal of Open Source Software , year =

    Vidar Skogvoll and Jonas Ronning , title =. Journal of Open Source Software , year =

  58. [58]

    Revisiting brittle fracture as an energy minimization problem

    Francfort, G.A. and Marigo, J.-J , title = "Revisiting brittle fracture as an energy minimization problem", journal = "Journal of the Mechanics and Physics of Solids", year =

  59. [59]

    Tortorelli

    Luigi Ambrosio and Vincenzon M. Tortorelli. Approximation of functional depending on jumps by elliptic functional via t‐convergence. Communications on Pure and Applied Mathematics. 1990

  60. [60]

    and Francfort, G

    Bourdin, B. and Francfort, G. and Marigo, J.-J. Numerical experiments in revisited brittle fracture. Journal of the Mechanics and Physics of Solids. 2000

  61. [61]

    and Francfort, G

    Bourdin, B. and Francfort, G. and Marigo, J.-J. The variational approach to fracture. Journal of Elasticity. 2008

  62. [62]

    Numerical implementation of the variational formulation of quasi-static brittle fracture

    Bourdin, Blaise. Numerical implementation of the variational formulation of quasi-static brittle fracture. 2007

  63. [63]

    Serreze, Roger G

    Mark C. Serreze, Roger G. Barry , title =. Global and Planetary Change , year =

  64. [64]

    npj Computational Materials , volume =

    A Unified Field Theory of Topological Defects and Non-Linear Local Excitations , author =. npj Computational Materials , volume =. doi:10.1038/s41524-023-01077-6 , urldate =

  65. [65]

    and Chen, Ying and Qu, Meng and Cheng, Xuhua , year=

    Liu, Minghao and Lei, Ruibo and Lin, Long and Hutchings, Jennifer K. and Chen, Ying and Qu, Meng and Cheng, Xuhua , year=. Seasonal and spatial variations in sea ice kinematics and their response to storms in the Arctic Transpolar Drift region in 2021–2022 , volume=. doi:10.1017/jog.2025.10099 , journal=

  66. [66]

    Journal of Geophysical Research: Oceans , volume =

    Oikkonen, Annu and Haapala, Jari and Lensu, Mikko and Karvonen, Juha and Itkin, Polona , title =. Journal of Geophysical Research: Oceans , volume =. doi:https://doi.org/10.1002/2016JC012387 , url =. https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1002/2016JC012387 , abstract =

  67. [67]

    and Glowacki, O

    Herman, A. and Glowacki, O. , TITLE =. The Cryosphere , VOLUME =. 2012 , NUMBER =

  68. [68]

    Journal of Geophysical Research: Oceans , volume =

    Weiss, Jérôme , title =. Journal of Geophysical Research: Oceans , volume =. doi:https://doi.org/10.1002/2017JC013236 , url =. https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1002/2017JC013236 , abstract =