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

Resumo Técnico: Modelagem Peridinâmica da Fratura de Plataformas de Gelo Impulsionada por Erosão Ondulatória, Derretimento Basal e Flexão por Flutuação

Declaração do Problema
As plataformas de gelo atuam como contrafortes críticos para as camadas de gelo continentais, retardando seu fluxo para o oceano. No entanto, a calota de gelo — o desprendimento de grandes fragmentos de gelo — enfraquece esse efeito estabilizador, acelerando o fluxo dos glaciares e contribuindo para a elevação do nível do mar. Embora se saiba que a frequência e a magnitude da calota são sensíveis ao forçamento ambiental, os mecanismos subjacentes baseados na física, particularmente aqueles impulsionados pela ação das ondas oceânicas, permanecem pouco compreendidos. As abordagens numéricas existentes, baseadas principalmente em Métodos de Elementos Finitos (MEF), lutam para resolver explicitamente a iniciação e propagação de trincas em meios porosos não locais como o gelo. Esses métodos convencionais dependem de derivadas espaciais e conectividade de malha fixa, frequentemente exigindo critérios de fratura ad hoc que falham em capturar a transição da erosão contínua impulsionada por ondas para a fratura discreta. Além disso, os modelos atuais de camadas de gelo carecem de leis de calota baseadas na física capazes de representar a fratura impulsionada por ondas e a falha controlada pela flutuação nas fronteiras das plataformas de gelo.

Metodologia
Para abordar essas limitações, os autores desenvolveram uma estrutura Peridinâmica (PD) não local especificamente adaptada para simular a fratura de plataformas de gelo. Diferentemente da mecânica do contínuo clássica, a PD utiliza equações integrais em vez de derivadas espaciais, acomodando naturalmente descontinuidades de deslocamento sem problemas de singularidade.

• Modelagem Constitutiva: O estudo emprega um modelo PD baseado em ligações que incorpora comportamento elástico-frágil com amolecimento dependente da taxa de deformação para simular a zona de processo de fratura quasi-frágil (ZPF) do gelo. Uma função escalar dependente da história é introduzida para capturar a transição da elasticidade linear para o acúmulo de danos e eventual ruptura das ligações.
• Formulação de Tensão: Para permitir a avaliação direta da concentração de tensão e liberação de energia nas pontas de trincas, os autores implementaram uma formulação de tensão virial estática de Piola-Kirchhoff de primeira ordem dentro da estrutura PD. Essa abordagem, validada contra soluções analíticas e MEF, permite o cálculo da tensão macroscópica a partir de forças microscópicas de ligação.
• Cenários do Modelo: A estrutura foi aplicada a dois cenários idealizados e fisicamente representativos:
  1. Colapso Frontal: Uma viga em balanço elástica submetida à gravidade, onde a erosão induzida por ondas é simulada pela remoção progressiva do suporte basal, levando a uma laje em balanço que falha sob seu próprio peso.
  2. Afrouxamento do Pé: Uma viga elástica flutuante submetida a momentos de flexão cíclicos induzidos pela ação das ondas na frente, combinados com forças de flutuação. Este cenário investiga o mecanismo de "afrouxamento do pé", onde pés de gelo submersos experimentam flexão ascendente, gerando tensões de tração que iniciam e propagam fraturas.

Principais Resultados
O estudo investigou sistematicamente a influência das taxas de erosão ondulatória, geometria da plataforma, comprimento da trinca e velocidade de propagação por meio de simulações numéricas.

• Validação: O modelo PD foi rigorosamente validado contra campos de tensão de elementos finitos, soluções analíticas da teoria de vigas e observações de campo recentes (Sartore et al., 2025). Os resultados demonstraram excelente concordância nos campos de deslocamento e distribuições de tensão para geometrias intactas, confirmando a capacidade do modelo de replicar resultados da mecânica clássica enquanto oferece capacidades superiores para simulação de fratura.
• Colapso Frontal: As simulações revelaram que, embora a deformação global de flexão seja relativamente insensível ao raio do filete basal criado pela erosão, a tensão de tração longitudinal é altamente sensível à geometria local. À medida que o entalhe de erosão se aprofunda, a concentração de tensão aumenta até que a resistência à tração do gelo seja excedida, desencadeando a fratura.
• Dinâmica de Afrouxamento do Pé: As simulações esclareceram o caminho mecânico que liga a erosão ondulatória contínua à calota discreta.
  • Acúmulo de Tensão: A fratura não é imediata; segue um período de acúmulo de tensão impulsionado por momentos de flexão cíclicos.
  • Propagação de Trincas: Uma vez iniciadas, as trincas propagam-se para o interior, seguindo caminhos de alta tensão de tração. A taxa de propagação não é constante, mas varia com o momento cíclico aplicado, avançando preferencialmente durante fases de tensão de tração máxima.
  • Comportamento de Limiar: Um comprimento crítico de erosão foi identificado. Além desse limiar, o crescimento da trinca acelera rapidamente devido às tensões de flexão amplificadas. O estudo encontrou uma correlação positiva e não linear entre o comprimento da trinca e o comprimento da erosão ondulatória.
  • Influência da Taxa de Erosão: Taxas de erosão mais baixas permitem a redistribuição de tensão sobre uma região mais ampla, levando a trincas distribuídas. Taxas de erosão mais altas resultam em concentração de tensão localizada, desencadeando iniciação de trincas mais precoce e propagação mais rápida ao longo de zonas estreitas.

Significado e Alegações
Os autores afirmam que este trabalho representa a primeira tentativa de empregar uma estrutura peridinâmica baseada na física para simular processos de calota de gelo impulsionados por erosão ondulatória e flexão por flutuação. O significado do estudo é triplo:

1. Avanço Metodológico: Estabelece a peridinâmica como uma ferramenta robusta e de próxima geração para modelagem de fratura de plataformas de gelo, preenchendo a lacuna entre a mecânica da fratura em escala de laboratório e processos criosféricos em grande escala. A estrutura elimina a necessidade de tratamentos numéricos especiais (como remalhamento ou critérios ad hoc) para lidar com a nucleação e crescimento de trincas sob grandes deformações.
2. Insight Mecanístico: O estudo fornece novos insights mecanísticos sobre como o forçamento ondulatório modula a redistribuição de tensão e a aceleração da fratura. Resolve explicitamente o mecanismo de "afrouxamento do pé", demonstrando que a flexão induzida por ondas sozinha pode gerar tensão de tração suficiente para impulsionar a propagação de trincas em grande escala e desestabilizar plataformas de gelo.
3. Implicações de Modelagem: Ao resolver explicitamente a iniciação e evolução da calota sob forçamento de ondas oceânicas, a estrutura proposta oferece um caminho para melhorar as parametrizações de calota em modelos de camadas de gelo e elevação do nível do mar. Os comportamentos de limiar e leis de escala identificados fornecem uma base quantitativa para o desenvolvimento de critérios de calota baseados na física que ligam a energia das ondas oceânicas e a geometria da plataforma à frequência e magnitude da calota, abordando uma limitação crítica nas projeções atuais da elevação futura do nível do mar.