A phase-field modeling approach to sea-ice fracturing
Pith reviewed 2026-06-27 10:47 UTC · model grok-4.3
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.
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
- 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.
Referee Report
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)
- [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.
- [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)
- [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.
- [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.
- [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
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
-
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
-
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
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
axioms (1)
- domain assumption Double-well free-energy formulation with overdamped displacement accurately represents sea-ice fracture physics.
Reference graph
Works this paper leans on
-
[1]
2025 , note =
NASA Worldview , author =. 2025 , note =
2025
-
[2]
2025 , note =
SASIP Project , author =. 2025 , note =
2025
-
[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
2001
-
[4]
Thorndike and R
A.S. Thorndike and R. Colony. Arctic Ocean Buoy Program. DTIC. 1980
1980
-
[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
2023
-
[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
2013
-
[7]
Fracture Mechanics: Fundamentals and Applications
Anderson, T.L. Fracture Mechanics: Fundamentals and Applications. 2005
2005
-
[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
1996
-
[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
2016
-
[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]
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
2009
-
[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
2016
-
[13]
Surveys in Geophysics
Scaling of fracture and faulting of ice on earth. Surveys in Geophysics. 2002
2002
-
[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
2019
-
[15]
A Dynamic Thermodynamic Sea Ice Model
Hibler, W.D. A Dynamic Thermodynamic Sea Ice Model. Journal of Physical Oceanography. 1979
1979
-
[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
1997
-
[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
2011
-
[18]
Fracture
Irwin, G.R. Fracture. In: Handbuch der Physik, Band VI, Elastizität und Plastizität. 1958
1958
-
[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
2015
-
[20]
Global Fiducials Library Image Gallery
USGS. Global Fiducials Library Image Gallery. 2014
2014
-
[21]
International Arctic Buoy Program (IABP)
NCEI NOAA. International Arctic Buoy Program (IABP)
-
[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
2022
-
[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
2023
-
[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
2025
-
[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
1957
-
[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
2018
-
[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
1989
-
[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
1998
-
[29]
and Levine, H
Collins, J. and Levine, H. Diffuse interface model of diffusion-limited crystal growth. Phys. Rev. B. 1985
1985
-
[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
2019
-
[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
2020
-
[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]
Elbanna and Jinhyun Choo , title =
Fan Fei and Md Shumon Mia and Ahmed E. Elbanna and Jinhyun Choo , title =. 2025 , howpublished =
2025
-
[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
2018
-
[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
2010
-
[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
2020
-
[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
2009
-
[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
2011
-
[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
2022
-
[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]
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
2016
-
[42]
Daniel Feltham. Sea ice rheology. Annual Review of Fluid Mechanics. 2008. doi:https://doi.org/10.1146/annurev.fluid.40.111406.102151
-
[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]
Spatschek and E
R. Spatschek and E. Brener and A. Karma. Phase field modeling of crack propagation. Philosophical magazine. 2010
2010
-
[45]
Phase field theory in materials physics, the holograph equation
Pierre Galenko. Phase field theory in materials physics, the holograph equation. 2024
2024
-
[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
1921
-
[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
1972
-
[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]
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]
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
2023
-
[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
2004
-
[52]
Creep and fracture of ice
Erland Schulson & Paul Duval. Creep and fracture of ice. 2009
2009
-
[53]
Freund, L. B. Dynamic fracture mechanics. 1998
1998
-
[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
2014
-
[55]
Landau and E
L.D. Landau and E. M. Lifshitz. Theory of Elasticity: volume 7. Vol. 7. 2012
2012
-
[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
2023
-
[57]
Journal of Open Source Software , year =
Vidar Skogvoll and Jonas Ronning , title =. Journal of Open Source Software , year =
-
[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]
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
1990
-
[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
2000
-
[61]
and Francfort, G
Bourdin, B. and Francfort, G. and Marigo, J.-J. The variational approach to fracture. Journal of Elasticity. 2008
2008
-
[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
2007
-
[63]
Serreze, Roger G
Mark C. Serreze, Roger G. Barry , title =. Global and Planetary Change , year =
-
[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]
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]
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]
and Glowacki, O
Herman, A. and Glowacki, O. , TITLE =. The Cryosphere , VOLUME =. 2012 , NUMBER =
2012
-
[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 =
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.