Attenuation of long waves through regions of irregular floating ice and bathymetry

This paper presents a revised, energy-conserving theoretical model for the attenuation of long waves through regions of irregular floating ice and random bathymetry, which corrects previous over-predictions by utilizing ensemble averaging of transfer matrix eigenvalues and successfully reproduces key features of field data, including frequency-dependent attenuation rates and high-frequency roll-over effects.

The Gist

Here is an explanation of the paper "Attenuation of long waves through regions of irregular floating ice and bathymetry" using simple language and creative analogies.

The Big Picture: Waves, Ice, and the "Ghost" Decay

Imagine you are standing on a beach watching a long, rolling ocean wave travel toward the shore. Now, imagine that instead of a smooth sandy bottom, the ocean floor is bumpy and uneven (like a rocky road). Or, imagine the surface of the water is covered in a chaotic mess of broken ice floes, like a giant, floating puzzle with pieces of different thicknesses.

When a wave hits these irregularities, it doesn't just keep going smoothly. It scatters. Some of it bounces back, some gets trapped, and the wave loses energy as it travels. This loss of energy is called attenuation.

For decades, scientists have tried to write math equations to predict exactly how much energy a wave loses in these messy environments. However, the old math had a big problem: it was too pessimistic. It predicted that waves would die out much faster than they actually do in real life.

This paper by Lloyd Dafydd and Richard Porter fixes that math. They discovered that the old way of calculating the average energy loss was accidentally counting "ghost" energy loss—energy that disappears only because of how we do the math, not because of physics.

---

The Problem: The "Ghost" in the Machine

To understand the fix, let's use an analogy.

Imagine you have a choir of 100 singers. Each singer is slightly out of tune with the others.

• The Old Math: The scientists tried to predict the volume of the choir by taking the average of every possible combination of singers. Because the singers are out of tune, their voices cancel each other out in the average. The math concluded, "Wow, the choir is almost silent!"
• The Reality: In any single performance (a single wave), the singers are all singing at once. They might be slightly out of tune, but they don't cancel each other out completely. The sound is loud, just a bit "fuzzy."

The old mathematical models were averaging the "silence" caused by the voices canceling each other out (phase cancellation). This created a "fictitious decay"—a fake prediction that the wave was dying out faster than it really was.

The Paper's Fix: The authors realized they needed to separate the "real" energy loss (caused by the wave hitting rocks or ice) from the "fake" loss (caused by the math averaging out the noise). They developed a new method that removes this "ghost" decay, ensuring the model only counts energy that is actually being scattered or reflected.

The New Model: A Bumpy Road and a Heavy Blanket

The authors built a new, simpler model to test their theory. They looked at two scenarios:

1. The Bumpy Road (Random Bathymetry): Imagine the ocean floor is a road with random potholes and bumps. As a wave rolls over it, it gets jiggled and slowed down.
2. The Heavy Blanket (Floating Ice): Imagine the ocean surface is covered by a blanket of broken ice. Some parts of the blanket are thick, some are thin. As the wave tries to push through, the heavy ice resists the motion, slowing the wave down.

The "Long Wave" Assumption:
The authors focused on "long waves" (like tsunamis or very long swells). Think of these waves as a giant, slow-moving snake. Because the snake is so long, it doesn't care about tiny pebbles on the road; it only cares about the big hills and valleys. This allows them to use a simpler set of equations (Ordinary Differential Equations) rather than incredibly complex 3D simulations.

The Results: The "Rollover" Effect

When they ran their new, corrected math and compared it to computer simulations, they found some fascinating things:

1. The Peak: The wave loses the most energy when the size of the bumps (or ice chunks) matches a specific "sweet spot" relative to the wave's length. This is similar to Bragg Resonance. Imagine pushing a child on a swing; if you push at the exact right rhythm, the swing goes high. Here, if the ice bumps match the wave rhythm, the wave gets "jammed" and loses energy rapidly.
2. The Rollover: This is the most exciting part. In many real-world measurements of waves in the Arctic, scientists noticed something weird: as the waves get higher frequency (shorter, faster ripples), the energy loss stops increasing and actually starts to decrease. This is called the "rollover effect."
   • Previous theories couldn't explain this. They thought waves would just keep losing more and more energy as they got faster.
   • This new model naturally predicts this rollover. It shows that once the waves get too fast for the ice bumps to "grab" onto them, the scattering becomes less effective, and the waves actually travel further.

Why Does This Matter?

This isn't just about math puzzles. It has real-world consequences for our changing planet:

• Climate Change: As the Arctic warms, sea ice is breaking up into smaller, thinner chunks. This changes how ocean waves travel through the ice.
• Prediction: If we want to predict how far a storm surge will travel into the Arctic, or how much ice will be broken up by a storm, we need accurate models.
• The Takeaway: By fixing the "ghost decay" in the math, this paper gives us a more accurate tool to understand how waves interact with a chaotic, broken ice world. It suggests that the "rollover" effect seen in nature is likely a real physical phenomenon caused by the randomness of the ice, not just a measurement error.

Summary in One Sentence

The authors fixed a flaw in old wave math that was overestimating how fast waves die out in rough water and broken ice, revealing that waves actually survive longer than we thought and explaining a mysterious "rollover" effect where fast waves stop losing energy as quickly.

Technical Summary

Here is a detailed technical summary of the paper "Attenuation of long waves through regions of irregular floating ice and bathymetry" by Lloyd Dafydd and Richard Porter.

1. Problem Statement

The paper addresses the attenuation (spatial decay) of surface waves propagating through inhomogeneous, disordered media, specifically focusing on two scenarios:

1. Random Bathymetry: Waves traveling over a seabed with randomly varying depth.
2. Random Floating Ice: Waves traveling under a continuous cover of fragmented floating ice with randomly varying thickness.

The Core Issue: Existing theoretical models for wave attenuation in such random media (based on multiple scattering and ensemble averaging) have historically over-predicted the rate of decay. Previous studies (e.g., Devillard et al., Mei et al.) suggested that attenuation scales with the square of the frequency ($\omega^2$) across all frequencies or decays exponentially in short-wave regimes. However, numerical simulations and field data often show different behaviors, including a "rollover effect" where attenuation peaks and then decreases at higher frequencies. The authors argue that previous theoretical failures stem from how ensemble averaging is performed, specifically the inclusion of "fictitious decay" caused by phase cancellations that do not occur in individual wave realizations.

2. Methodology

A. Mathematical Model

The authors utilize a linearized long-wavelength (shallow water) model.

• Derivation: The model is derived from the Euler equations for inviscid, incompressible flow. It involves a depth-averaging process and an expansion to second order in the small parameter $\epsilon = H/L$ (ratio of depth to wavelength).
• Governing Equation: The scattering process is reduced to solving a single Ordinary Differential Equation (ODE) for a "pseudo-potential" $\Omega(x)$:
  $$(\hat{h}(x)\Omega'(x))' + K\Omega(x) = 0$$
  where $K = \omega^2/g$.
  • For random bathymetry (no ice), $\hat{h}(x)$ depends on the fluid depth $h(x)$.
  • For random ice, the model is extended to include the inertia of the floating ice. The depth $h(x)$ is replaced by an effective depth term involving the ice submergence $d(x)$. This second-order extension is crucial as it captures vertical acceleration effects neglected in first-order shallow water models.

B. Analytical Approach: Multiple Scales with Corrected Averaging

The authors employ a multiple-scales analysis assuming small randomness amplitude ($\sigma \ll 1$).

• The Correction: In standard multiple-scales analysis, ensemble averaging of the first-order solution often introduces terms that represent "fictitious decay." These terms arise from phase cancellations between waves when averaging over different realizations of the random medium, rather than from actual energy dissipation or scattering.
• The Fix: The authors explicitly identify and remove these fictitious decay terms from the ensemble averaging process. They enforce a condition where the leading-order solution is deterministic, while the statistical properties are captured in the slow-varying amplitudes of left- and right-propagating waves.
• Energy Conservation: The revised approach is proven to conserve energy, ensuring that the calculated attenuation represents true physical scattering rather than mathematical artifacts.

C. Numerical Simulations

To validate the theory, the authors perform numerical simulations:

• Random Generation: Random functions for bed elevation or ice thickness are generated using a weighted moving average method to satisfy specific Gaussian correlation properties.
• Transfer Matrix Method: Instead of calculating reflection/transmission coefficients directly (which can be affected by boundary reflections), they use a transfer matrix to propagate the wave solution across the random region.
• Decay Measurement: The attenuation rate is accurately measured by computing the eigenvalues of the transfer matrix. If the eigenvalues are real, the decay rate is derived from their logarithm. This method isolates the decay due to the random medium itself, avoiding interference from boundary junctions.
• Ensemble Averaging: Results are averaged over hundreds of random realizations ($N=500$) to compare with the theoretical predictions.

3. Key Contributions

1. Correction of Ensemble Averaging: The primary theoretical contribution is the identification and removal of "fictitious decay" terms in the ensemble averaging process. This resolves the long-standing discrepancy between theoretical predictions (which over-predicted decay) and numerical simulations.
2. Energy Conservation: The revised model explicitly demonstrates energy conservation, a property often violated or obscured in previous multiple-scales approaches to random media.
3. Extension to Floating Ice: The authors successfully extend the shallow-water ODE model to include a continuous cover of fragmented floating ice with variable thickness. This includes the necessary second-order terms to account for ice inertia.
4. Prediction of the "Rollover Effect": The model predicts a peak in attenuation at a specific wavenumber-correlation length relationship ($k_0\Lambda \approx 1$), followed by an exponential decay at higher frequencies. This "rollover" matches features observed in field data, which previous models failed to capture without introducing non-physical damping terms.

4. Key Results

• Attenuation Scaling:
  • Low Frequencies ($k_0\Lambda \ll 1$): Attenuation scales with the square of the wavenumber ($k_i \propto k_0^2$), which corresponds to $\omega^2$ in shallow water.
  • Peak Attenuation: The maximum attenuation occurs near Bragg resonance ($k_0\Lambda \approx 1$), where the wavelength is roughly twice the correlation length of the randomness.
  • High Frequencies ($k_0\Lambda > 1$): Attenuation decays exponentially as $k_0$ increases. Crucially, the authors show this decay is a feature of the long-wavelength model and not a finite-depth effect, contradicting earlier interpretations by Devillard et al. and Mei et al.
• Comparison with Simulations: The corrected theoretical predictions align almost perfectly with numerical simulations of finite random sections, whereas the uncorrected (previous) theories significantly over-predicted the decay.
• Ice vs. Bathymetry: The mathematical structure for random ice thickness is identical to random bathymetry, differing only in the definition of the scaling coefficients ($C_1, C_2, C_3$) and the dispersion relation.
• Rollover Effect: The model naturally produces the high-frequency rollover (attenuation peak followed by decrease) observed in Marginal Ice Zone (MIZ) field data, suggesting that random variations in ice thickness are a plausible physical mechanism for this phenomenon, without needing ad-hoc damping terms.

5. Significance and Limitations

Significance:

• The paper provides a rigorous mathematical explanation for why previous theories over-predicted wave decay in random media.
• It offers a unified framework for analyzing wave scattering over random beds and through random ice covers using a computationally efficient ODE approach.
• The ability to capture the "rollover effect" theoretically supports the hypothesis that multiple scattering in random ice fields is a dominant mechanism for wave attenuation in the Arctic, potentially reducing the need for empirical damping parameters in operational models.

Limitations:

• Shallow Water Assumption: The model relies on the long-wavelength assumption ($Kh \ll 1$). It may not be accurate for deep water or very short wavelengths where discrete scattering by individual ice floes (rather than a continuum) becomes dominant.
• Continuum Approximation: The ice is modeled as a continuous layer with variable thickness. It does not explicitly model discrete ice floes, cracks, or leads, which are significant in real-world sea ice.
• No Physical Dissipation: The model excludes physical dissipation mechanisms like fluid viscosity, ice-ice friction, or viscoelastic damping of the ice itself. It isolates attenuation purely due to multiple scattering.

Future Work: The authors plan to extend the model to finite water depth, variable ice concentration, discrete floe models, and three-dimensional scattering to better bridge the gap between their continuum theory and complex field observations.