How Do Ice Shelves Calve? Peridynamic Modeling of Ice Shelf Fracture Driven by Wave Erosion, Basal Melting, and Buoyancy Flexure

Résumé technique : Modélisation peridynamique de la fracturation des plateformes glaciaires sous l'effet de l'érosion par les vagues, de la fonte basale et de la flexion de flottabilité

Énoncé du problème
Les plateformes glaciaires agissent comme des contreforts critiques pour les calottes glaciaires terrestres, ralentissant leur écoulement vers l'océan. Cependant, le détachement de grands fragments de glace, appelé vêlage, affaiblit cet effet stabilisateur, accélère l'écoulement des glaciers et contribue à l'élévation du niveau de la mer. Bien que la fréquence et l'ampleur du vêlage soient connues pour être sensibles aux forçages environnementaux, les mécanismes physiques sous-jacents, en particulier ceux pilotés par l'action des vagues océaniques, restent mal compris. Les approches numériques existantes, principalement basées sur la méthode des éléments finis (MEF), peinent à résoudre explicitement l'initiation et la propagation de fissures dans des milieux poreux non locaux comme la glace. Ces méthodes conventionnelles reposent sur des dérivées spatiales et une connectivité de maillage fixe, nécessitant souvent des critères de fracture ad hoc qui échouent à capturer la transition de l'érosion continue pilotée par les vagues vers la fracture discrète. De plus, les modèles actuels de calottes glaciaires manquent de lois de vêlage fondées sur la physique capables de représenter la fracture pilotée par les vagues et la rupture contrôlée par la flottabilité aux fronts des plateformes glaciaires.

Méthodologie
Pour pallier ces limitations, les auteurs ont développé un cadre peridynamique (PD) non local spécifiquement adapté à la simulation de la fracturation des plateformes glaciaires. Contrairement à la mécanique des milieux continus classique, la PD utilise des équations intégrales plutôt que des dérivées spatiales, accommodant naturellement les discontinuités de déplacement sans problèmes de singularité.

• Modélisation constitutive : L'étude emploie un modèle PD basé sur les liaisons, incorporant un comportement élastique-fragile avec un adoucissement dépendant du taux de déformation pour simuler la zone de processus de fracture quasi-fragile (FPZ) de la glace. Une fonction scalaire dépendante de l'histoire est introduite pour capturer la transition de l'élasticité linéaire à l'accumulation de dommages et à la rupture finale des liaisons.
• Formulation des contraintes : Pour permettre l'évaluation directe de la concentration de contraintes et de la libération d'énergie aux pointes de fissure, les auteurs ont implémenté une formulation de contrainte viriale statique de Piola-Kirchhoff première au sein du cadre PD. Cette approche, validée par rapport à des solutions analytiques et à la MEF, permet le calcul de contraintes macroscopiques à partir de forces de liaisons microscopiques.
• Scénarios de modélisation : Le cadre a été appliqué à deux scénarios idéalisés et physiquement représentatifs :
  1. Effondrement frontal : Une poutre en porte-à-faux élastique soumise à la gravité, où l'érosion induite par les vagues est simulée par le retrait progressif du support basal, conduisant à une dalle en surplomb qui cède sous son propre poids.
  2. Desserrement du pied : Une poutre élastique flottante soumise à des moments de flexion cycliques induits par l'action des vagues à l'avant, combinés à des forces de flottabilité. Ce scénario examine le mécanisme de « desserrement du pied » où les pieds de glace immergés subissent une flexion ascendante, générant des contraintes de traction qui initient et propagent des fissures.

Résultats clés
L'étude a systématiquement examiné l'influence des taux d'érosion par les vagues, de la géométrie de la plateforme, de la longueur de la fissure et de la vitesse de propagation grâce à des simulations numériques.

• Validation : Le modèle PD a été rigoureusement validé par rapport aux champs de contraintes des éléments finis, aux solutions analytiques de la théorie des poutres et à des observations récentes sur le terrain (Sartore et al., 2025). Les résultats ont démontré un excellent accord dans les champs de déplacement et les distributions de contraintes pour des géométries intactes, confirmant la capacité du modèle à reproduire les résultats de la mécanique classique tout en offrant des capacités supérieures pour la simulation de la fracture.
• Effondrement frontal : Les simulations ont révélé que, bien que la déformation globale par flexion soit relativement insensible au rayon de congé basal créé par l'érosion, la contrainte de traction longitudinale est hautement sensible à la géométrie locale. À mesure que l'encoche d'érosion s'approfondit, la concentration de contraintes augmente jusqu'à ce que la résistance à la traction de la glace soit dépassée, déclenchant la fracture.
• Dynamique du desserrement du pied : Les simulations ont clarifié la voie mécanique reliant l'érosion continue par les vagues au vêlage discret.
  • Accumulation de contraintes : La fracture n'est pas immédiate ; elle suit une période d'accumulation de contraintes pilotée par des moments de flexion cycliques.
  • Propagation de fissure : Une fois initiées, les fissures se propagent vers l'intérieur des terres, suivant des chemins de haute contrainte de traction. Le taux de propagation n'est pas constant mais varie avec le moment cyclique appliqué, avançant préférentiellement pendant les phases de contrainte de traction maximale.
  • Comportement seuil : Une longueur d'érosion critique a été identifiée. Au-delà de ce seuil, la croissance de la fissure s'accélère rapidement en raison de contraintes de flexion amplifiées. L'étude a trouvé une corrélation positive et non linéaire entre la longueur de la fissure et la longueur d'érosion par les vagues.
  • Influence du taux d'érosion : Des taux d'érosion plus faibles permettent une redistribution des contraintes sur une région plus large, conduisant à une fissuration distribuée. Des taux d'érosion plus élevés entraînent une concentration de contraintes localisée, déclenchant une initiation de fissure plus précoce et une propagation plus rapide le long de zones étroites.

Importance et affirmations
Les auteurs affirment que ce travail représente la première tentative d'employer un cadre peridynamique fondé sur la physique pour simuler les processus de vêlage de glace pilotés par l'érosion par les vagues et la flexion de flottabilité. L'importance de l'étude est triple :

1. Avancement méthodologique : Il établit la peridynamique comme un outil robuste de nouvelle génération pour la modélisation de la fracturation des plateformes glaciaires, comblant le fossé entre la mécanique de la fracture à l'échelle du laboratoire et les processus cryosphériques à grande échelle. Le cadre élimine le besoin de traitements numériques spéciaux (comme le remaillage ou des critères ad hoc) pour gérer la nucléation et la croissance de fissures sous de grandes déformations.
2. Insight mécanistique : L'étude fournit de nouvelles insights mécanistiques sur la façon dont le forçage par les vagues module la redistribution des contraintes et l'accélération de la fracture. Elle résout explicitement le mécanisme de « desserrement du pied », démontrant que la flexion induite par les vagues seule peut générer une contrainte de traction suffisante pour piloter la propagation de fissures à grande échelle et déstabiliser les plateformes glaciaires.
3. Implications pour la modélisation : En résolvant explicitement l'initiation et l'évolution du vêlage sous l'effet du forçage par les vagues océaniques, le cadre proposé offre une voie pour améliorer les paramétrisations de vêlage dans les modèles de calottes glaciaires et d'élévation du niveau de la mer. Les comportements seuils et les lois d'échelle identifiés fournissent une base quantitative pour développer des critères de vêlage fondés sur la physique qui relient l'énergie des vagues océaniques et la géométrie de la plateforme à la fréquence et à l'ampleur du vêlage, répondant ainsi à une limitation critique dans les projections actuelles de l'élévation future du niveau de la mer.