Modelling hydroelastic flexure of arbitrarily shaped ice shelves forced by long ocean waves

This paper presents a novel finite element method for modeling hydroelastic flexure of arbitrarily shaped Antarctic ice shelves under long ocean wave forcing, enabling the identification of resonant responses and the analysis of how shelf geometry, wave direction, and grounding proportions influence mechanical stresses and calving risks.

The Gist

Imagine the massive ice shelves floating off the coast of Antarctica not as solid, unbreakable blocks, but as giant, thin sheets of ice that act like a trampoline or a flexible diving board. When huge ocean waves crash against them, these ice sheets bend and flex. If they bend too much, they can snap, causing massive chunks of ice to break off (a process called "calving"). This is dangerous because it weakens the ice shelf's ability to hold back the massive glaciers behind it, which could eventually lead to rising sea levels.

For a long time, scientists could only model these ice shelves as simple, straight strips with uniform thickness. But real ice shelves are messy: they have weird shapes, varying thicknesses, and sit on uneven ocean floors. Modeling the bending of these complex shapes in 3D space, while also accounting for the water underneath, is like trying to solve a puzzle where every piece is a different shape and the rules keep changing. It's incredibly difficult to compute.

The New "Smart Trampoline" Model
The authors of this paper have built a new computer program that acts like a high-tech, flexible ruler. Instead of trying to force the ice shelf into a simple shape, their method uses a special kind of digital "net" (called finite elements) that can wrap around any irregular shape of ice shelf, no matter how weird.

To make the computer calculation fast enough to be useful, they used a clever trick called a "Dirichlet-to-Neumann map." Think of this like putting up a smart fence around your backyard. Instead of calculating the waves for the entire infinite ocean (which would take forever), this "smart fence" knows exactly how the waves should behave outside the fence based on what's happening right at the fence line. This allows the computer to focus its power on the ice shelf itself without getting bogged down by the rest of the ocean.

What They Discovered
Using this new tool, the researchers ran simulations to see how different factors change how much the ice shelf wiggles. Here is what they found, using simple analogies:

• The Shape Matters (The "Harbor Effect"): They tested ice shelves that were long and skinny, square, or wide and short. They found that long, skinny ice shelves (like a narrow hallway) tend to wiggle much more violently than wide ones. It's similar to how a narrow harbor can amplify waves inside it, making the water slosh higher than the waves outside. The wider the ice shelf, the more the energy spreads out, and the less it bends.
• The Angle of the Wave: If a wave hits the ice shelf straight on, it creates a specific pattern of bending. But if the wave hits at an angle (like a car hitting a curb sideways), the pattern changes completely. Some parts of the ice shelf might start shaking much harder than before, while other parts calm down. The angle of the incoming wave is a critical switch that changes which parts of the ice are in danger.
• How Much is "Stuck" to Land: Some ice shelves are mostly attached to the land (like a wide sheet), while others stick out far into the ocean like a long tongue (like the Drygalski Ice Tongue). The researchers found that the more the ice shelf sticks out into the open ocean, the less it resonates (wiggles) at the low frequencies that usually cause the most damage. However, as the "tongue" gets longer, the ice starts to shake at higher, faster frequencies.

Why This Matters
The main achievement of this paper is that they finally have a way to calculate how any shape of ice shelf will react to ocean waves, not just simple ones. They showed that the shape of the shelf, the angle of the waves, and how much of it is attached to land all dramatically change the "resonance"—the point where the ice starts to vibrate violently.

By identifying these "sweet spots" where the ice is most likely to break, this method helps scientists understand which specific ice shelves are most vulnerable to the long, rolling waves coming from the ocean. It's a step toward predicting when and where these massive ice structures might break apart.

Technical Summary

Technical Summary: Modelling Hydroelastic Flexure of Arbitrarily Shaped Ice Shelves

Problem Statement
The paper addresses the hydroelastic response of Antarctic ice shelves to long ocean waves (periods ranging from approximately 15 seconds to 4 hours), including long swell, infragravity waves, and tsunamis. These waves induce flexure in the ice shelves, generating mechanical stresses that can amplify fractures and lead to calving events. Existing mathematical models for this phenomenon are largely limited to one horizontal dimension, assuming uniform ice thickness and water depth, or rely on artificial boundary conditions that restrict the problem to finding eigenvalues (resonant frequencies) rather than solving the scattering problem for prescribed wave forcing. Furthermore, previous two-dimensional approaches have required computationally expensive "brute force" methods or complex domain mappings to handle arbitrary geometries and varying bathymetry. The challenge lies in solving a sixth-order hydroelastic partial differential equation system coupled with shallow water theory over large, irregular domains that connect to the open ocean.

Methodology
The authors develop a solution method for a hydroelastic mathematical model based on the Kirchhoff-Love plate theory for the ice shelf and linearized shallow water theory for the underlying water cavity and surrounding ocean. The key components of the methodology include:

• Governing Equations: The model couples the ice shelf flexural displacement ($\eta$) and the velocity potential in the water cavity ($\Phi$) with the velocity potential in the open ocean ($\phi$). The system assumes small amplitude motions, irrotational flow, and neglects Coriolis effects for the wave periods of interest.
• Finite Element Discretization: The computational domain is divided into the ice-shelf/cavity region ($\Omega$) and the open ocean region ($V$).
  • In $\Omega$, a specialized nonconforming hydroelastic finite element (SHEEL) is used. This is a three-node triangle element with four degrees of freedom per node: ice shelf deflection, two rotations, and the velocity potential. It employs linear Lagrange interpolation for the potential and a nonconforming Specht approximation for the plate deflection to capture high-order derivatives without requiring domain mapping.
  • In $V$, standard linear Lagrange triangles (LSWT) are used for the velocity potential.
• Boundary Conditions and Radiation:
  • The grounding line is treated as clamped (zero displacement and rotation).
  • The lateral boundaries of the ice shelf exposed to the ocean have zero normal bending moment and active shear force.
  • Crucially, to handle the connection to the open ocean without extending the computational domain to infinity, a Dirichlet-to-Neumann (DtN) radiation condition is applied on an artificial boundary ($\Gamma_R$). This allows the open ocean domain to be truncated to a relatively small size while accurately representing wave radiation.
• Variational Formulation: The problem is formulated as a monolithic system using a sesquilinear functional that incorporates the hydroelastic coupling terms. The system is solved for the scattered wave field given an incident plane wave.

Key Contributions

1. Arbitrary Geometry Handling: The method allows for the modeling of ice shelves of arbitrary shape with non-uniform thickness and variable bathymetry, overcoming the 1D limitations of previous studies.
2. Scattering Problem Solution: Unlike previous eigenproblem-focused studies, this method solves the scattering problem, quantifying the ice shelf's deflection in response to specific incident wave forcing over a continuous frequency spectrum.
3. Computational Efficiency: By combining the tailored SHEEL element with the DtN boundary condition, the method achieves high fidelity with significantly fewer degrees of freedom compared to "brute force" commercial software or methods requiring large Sommerfeld radiation boundaries.
4. Novel Parametric Studies: The efficiency of the method enables broad frequency range studies to identify resonant responses, which were previously difficult to compute for complex geometries.

Results
The method was verified against 1D models, Sommerfeld radiation conditions, and existing literature on harbour oscillations. Parametric analyses revealed:

• Aspect Ratio: Narrow, elongated ice shelves exhibit more distinct resonant peaks with higher amplitudes compared to square or wide shelves. This is linked to the "harbour paradox," where narrower openings can lead to larger internal responses.
• Incidence Angle: While the fundamental (lowest frequency) mode is largely unaffected by the angle of incidence, higher harmonics and the overall shape of the response spectrum change significantly with oblique wave angles. Some higher-order peaks are amplified for oblique incidence.
• Grounding Line Geometry (Ice Shelves vs. Ice Tongues): The proportion of the ice shelf that is grounded versus protruding into the ocean significantly alters the response. As the protruding portion increases (transitioning from a fully embedded shelf to an ice tongue), resonant peaks decrease in amplitude and shift toward higher frequencies. The unconfined portion of the shelf responds similarly to open ocean surface gravity waves, while larger responses tend to occur in the embedded region.

Significance
The paper claims that this work provides a computationally efficient, high-fidelity tool for studying the hydroelastic response of ice shelves to long ocean waves. By incorporating two-dimensional horizontal effects and arbitrary shapes, the method sets a computational basis for understanding which Antarctic ice shelves are most susceptible to wave-induced impacts. The identification of resonant responses across broad frequency ranges and varying geometries offers new insights into the mechanisms driving ice shelf flexure and potential calving, particularly for shelves that have thinned or lost sea ice barriers. The authors note that while the current model uses shallow water theory, the framework could be extended to include swell regimes (via single-mode approximation) and viscoelastic effects, though these are not implemented in the current study.