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

Resumen Técnico: Modelado Peridinámico de la Fractura de Plataformas de Hielo Impulsada por Erosión Ondulatoria, Fusión Basal y Flexión por Flotabilidad

Planteamiento del Problema
Las plataformas de hielo actúan como contrafuertes críticos para las capas de hielo terrestres, frenando su flujo hacia el océano. Sin embargo, la desprendimiento de hielo (calving), es decir, la separación de grandes fragmentos de hielo, debilita este efecto estabilizador, acelerando el flujo de los glaciares y contribuyendo al aumento del nivel del mar. Si bien se sabe que la frecuencia y magnitud del desprendimiento de hielo son sensibles a las fuerzas ambientales, los mecanismos subyacentes basados en la física, particularmente aquellos impulsados por la acción de las olas oceánicas, permanecen poco comprendidos. Los enfoques numéricos existentes, basados principalmente en el Método de Elementos Finitos (MEF), luchan por resolver explícitamente la iniciación y propagación de grietas en medios porosos no locales como el hielo. Estos métodos convencionales dependen de derivadas espaciales y conectividad de malla fija, requiriendo a menudo criterios de fractura ad hoc que no logran capturar la transición desde la erosión continua impulsada por olas hasta la fractura discreta. Además, los modelos actuales de capas de hielo carecen de leyes de desprendimiento de hielo basadas en la física capaces de representar la fractura impulsada por olas y la falla controlada por flotabilidad en los frentes de las plataformas de hielo.

Metodología
Para abordar estas limitaciones, los autores desarrollaron un marco de Peridinámica (PD) no local específicamente adaptado para simular la fractura de plataformas de hielo. A diferencia de la mecánica de medios continuos clásica, la PD utiliza ecuaciones integrales en lugar de derivadas espaciales, acomodando naturalmente las discontinuidades de desplazamiento sin problemas de singularidad.

• Modelado Constitutivo: El estudio emplea un modelo PD basado en enlaces que incorpora un comportamiento elástico-frágil con ablandamiento dependiente de la tasa de deformación para simular la zona de proceso de fractura cuasi-frágil (FPZ) del hielo. Se introduce una función escalar dependiente de la historia para capturar la transición desde la elasticidad lineal hasta la acumulación de daño y la eventual ruptura de enlaces.
• Formulación de Esfuerzos: Para permitir la evaluación directa de la concentración de esfuerzos y la liberación de energía en las puntas de las grietas, los autores implementaron una formulación de esfuerzos viriales estáticos de primer Piola-Kirchhoff dentro del marco de la PD. Este enfoque, validado frente a soluciones analíticas y MEF, permite el cálculo del esfuerzo macroscópico a partir de fuerzas de enlace microscópicas.
• Escenarios del Modelo: El marco se aplicó a dos escenarios idealizados y físicamente representativos:
  1. Colapso Frontal: Una viga en voladizo elástica sometida a la gravedad, donde la erosión inducida por olas se simula mediante la eliminación progresiva del soporte basal, lo que conduce a una losa en voladizo que falla bajo su propio peso.
  2. Aflojamiento del Pie: Una viga elástica flotante sometida a momentos de flexión cíclicos inducidos por la acción de las olas en el frente, combinados con fuerzas de flotabilidad. Este escenario investiga el mecanismo de "aflojamiento del pie" donde los pies de hielo sumergidos experimentan una flexión ascendente, generando esfuerzos tensiles que inician y propagan fracturas.

Resultados Clave
El estudio investigó sistemáticamente la influencia de las tasas de erosión por olas, la geometría de la plataforma, la longitud de la grieta y la velocidad de propagación mediante simulaciones numéricas.

• Validación: El modelo PD fue validado rigurosamente frente a campos de esfuerzos de elementos finitos, soluciones analíticas de teoría de vigas y observaciones de campo recientes (Sartore et al., 2025). Los resultados demostraron un excelente acuerdo en los campos de desplazamiento y distribuciones de esfuerzo para geometrías intactas, confirmando la capacidad del modelo para replicar resultados de mecánica clásica mientras ofrece capacidades superiores para la simulación de fracturas.
• Colapso Frontal: Las simulaciones revelaron que, si bien la deformación global por flexión es relativamente insensible al radio de filete basal creado por la erosión, el esfuerzo tensil longitudinal es altamente sensible a la geometría local. A medida que la muesca de erosión se profundiza, la concentración de esfuerzos aumenta hasta que se supera la resistencia tensil del hielo, desencadenando la fractura.
• Dinámica de Aflojamiento del Pie: Las simulaciones aclararon la vía mecánica que vincula la erosión continua por olas con el desprendimiento de hielo discreto.
  • Acumulación de Esfuerzos: La fractura no es inmediata; sigue un período de acumulación de esfuerzos impulsado por momentos de flexión cíclicos.
  • Propagación de Grietas: Una vez iniciadas, las grietas se propagan tierra adentro, siguiendo caminos de alto esfuerzo tensil. La tasa de propagación no es constante, sino que varía con el momento cíclico aplicado, avanzando preferentemente durante las fases de esfuerzo tensil máximo.
  • Comportamiento Umbral: Se identificó una longitud de erosión crítica. Más allá de este umbral, el crecimiento de la grieta se acelera rápidamente debido a los esfuerzos de flexión amplificados. El estudio encontró una correlación positiva y no lineal entre la longitud de la grieta y la longitud de erosión por olas.
  • Influencia de la Tasa de Erosión: Las tasas de erosión más bajas permiten una redistribución de esfuerzos sobre una región más amplia, lo que conduce a una fisuración distribuida. Las tasas de erosión más altas resultan en una concentración de esfuerzos localizada, desencadenando una iniciación de grieta más temprana y una propagación más rápida a lo largo de zonas estrechas.

Significado y Afirmaciones
Los autores afirman que este trabajo representa el primer intento de emplear un marco de peridinámica basado en la física para simular procesos de desprendimiento de hielo impulsados por erosión ondulatoria y flexión por flotabilidad. El significado del estudio es triple:

1. Avance Metodológico: Establece la peridinámica como una herramienta robusta de próxima generación para el modelado de fracturas en plataformas de hielo, cerrando la brecha entre la mecánica de fracturas a escala de laboratorio y los procesos criosféricos a gran escala. El marco elimina la necesidad de tratamientos numéricos especiales (como remallado o criterios ad hoc) para manejar la nucleación y crecimiento de grietas bajo grandes deformaciones.
2. Perspectiva Mecanística: El estudio proporciona nuevas perspectivas mecanísticas sobre cómo la forzamiento ondulatorio modula la redistribución de esfuerzos y la aceleración de fracturas. Resuelve explícitamente el mecanismo de "aflojamiento del pie", demostrando que la flexión inducida por olas por sí sola puede generar suficiente esfuerzo tensil para impulsar la propagación de grietas a gran escala y desestabilizar las plataformas de hielo.
3. Implicaciones de Modelado: Al resolver explícitamente la iniciación y evolución del desprendimiento de hielo bajo forzamiento de olas oceánicas, el marco propuesto ofrece una vía para mejorar las parametrizaciones de desprendimiento de hielo en modelos de capas de hielo y aumento del nivel del mar. Los comportamientos umbral identificados y las leyes de escala proporcionan una base cuantitativa para desarrollar criterios de desprendimiento de hielo basados en la física que vinculen la energía de las olas oceánicas y la geometría de la plataforma con la frecuencia y magnitud del desprendimiento de hielo, abordando una limitación crítica en las proyecciones actuales del aumento futuro del nivel del mar.