Wave-induced drift in third-order deep-water theory

This paper investigates third-order deep-water wave-induced drift using a reduced Hamiltonian formulation, revealing that while classical Stokes drift slightly misestimates particle motion, incorporating difference harmonic terms significantly improves agreement with nonlinear trajectory calculations across various sea states.

The Gist

Imagine the ocean as a giant, invisible conveyor belt. When you see a wave rolling toward the shore, it looks like the water is moving forward. But if you drop a cork into the water, you'll notice something strange: the cork bobs up and down, and it moves forward a little bit, but it doesn't travel all the way with the wave. The wave passes through the water, but the water itself mostly stays put, just dancing in circles.

This paper is about figuring out exactly how much that cork (or a plastic bottle, or a tiny fish) actually gets pushed forward by the wave, and how our old math rules for calculating that push are slightly off.

Here is the story of the research, broken down into simple concepts:

1. The Old Rule: The "Perfect Circle" Myth

For a long time, scientists used a simple formula (called Stokes Drift) to guess how far things drift. They imagined the water particles moving in perfect circles, like a Ferris wheel.

• The Problem: This formula is based on "linear" math, which assumes waves are gentle and small. It's like assuming a car always drives at a steady 30 mph.
• The Reality: Real waves are messy. They have steep crests and flat troughs. When waves get bigger or interact with each other, the "Ferris wheel" isn't a perfect circle anymore; it's a lopsided loop. The old formula was slightly underestimating the push at the surface and overestimating it deep down.

2. The New Tool: The "Hamiltonian" Recipe

The author, Raphael Stuhlmeier, used a sophisticated mathematical tool called the Hamiltonian formulation. Think of this as a high-tech recipe book for waves.

• Instead of just looking at the surface, this recipe accounts for the "hidden ingredients" inside the wave.
• It breaks waves down into layers of complexity:
  • Level 1 (Linear): The basic wave.
  • Level 2 (Quadratic): When two waves meet, they create "ghost waves" (called bound harmonics) that stick to the main wave.
  • Level 3 (Cubic): When three waves interact, they tweak the speed and shape of the whole group.

3. The Big Discovery: The "Difference" Waves

The most exciting part of the paper is about what happens when you have two different waves traveling together (like a big wave and a small wave).

• The Analogy: Imagine two runners on a track. One is fast, one is slow. If they run side-by-side, you see a pattern of "beats" or pulsing where they overlap.
• The Science: When these waves interact, they create a "difference wave." This isn't a new wave traveling forward; it's a slow, deep pulse that moves underneath the main action.
• The Result: The old math ignored this deep pulse. The new math includes it. The author found that this "difference wave" acts like a hidden current deep underwater.
  • Near the surface: The old math was okay, but slightly too low.
  • Deep underwater: The old math was way off. The new math shows that the "difference wave" actually pushes particles forward more than we thought, even deep down where the main waves don't reach.

4. Testing the Theory: From Single Waves to Random Chaos

The author didn't just do this for one perfect wave. They tested it in three scenarios:

1. The Solo Wave (Monochromatic): A single, perfect wave. The new math confirmed that the old formula is a bit too small at the surface and too big deep down.
2. The Duo (Bichromatic): Two waves interacting. Here, the "difference wave" effect became very important. It created a "return flow" deep down (water moving backward) before the net drift took over.
3. The Chaos (Random Waves): Real ocean conditions with many waves of different sizes. Even in this mess, the new formula (which includes the "difference" terms) gave a much better prediction of how much drift happens at different depths.

5. Why Should You Care?

You might think, "So what? It's just a tiny bit of drift." But this matters for:

• Pollution: If you spill oil or plastic in the ocean, where does it go? The old math might say it stays near the surface; the new math suggests it might drift deeper or further than we thought.
• Climate & Ecology: Tiny plankton and bacteria rely on these currents to move around. Getting the drift calculation wrong means we might misunderstand how nutrients or bacteria spread through the ocean.
• Safety: Ships and offshore platforms need to know exactly how much the water is moving to stay safe.

The Bottom Line

The ocean is more complex than a simple Ferris wheel. Waves are like a complex dance where partners interact, creating hidden currents deep below the surface.

This paper says: "Stop using the simple, old map. We have a new, more detailed map that includes the hidden 'difference waves.' It shows that the ocean's conveyor belt is slightly stronger at the surface and surprisingly active deep down than we previously believed."

By using this new, more accurate math, we can better predict where floating trash goes, how marine life travels, and how the ocean really moves.

Technical Summary

Here is a detailed technical summary of the paper "Wave-induced drift in third-order deep–water theory" by Raphael Stuhlmeier.

1. Problem Statement

The paper investigates the kinematics of fluid particles beneath unidirectional, deep-water waves, specifically focusing on wave-induced drift (Stokes drift) up to the third order of nonlinearity.

• The Core Issue: While the classical Stokes drift formulation (derived from linear theory) is widely used in oceanography and engineering, it relies on simplifications that may not accurately capture particle transport in steep or complex wave fields.
• The Gap: There is a lack of quantitative comparisons between the classical Stokes drift approximation and the actual Lagrangian drift computed directly from particle trajectory mappings using higher-order nonlinear wave theories (2nd and 3rd order).
• The Challenge: Comparing Eulerian (fixed coordinate) and Lagrangian (particle-following) descriptions is difficult, particularly for unsteady, irregular wave fields where particle paths are not closed loops. Furthermore, standard Eulerian theories often obscure the distinction between free waves and bound harmonics, which are crucial for drift calculations.

2. Methodology

The author employs a weakly nonlinear approach using the Hamiltonian formulation of the water wave problem (Zakharov and Krasitskii) to derive velocity fields and free surface elevations up to the third order.

• Hamiltonian Expansion: The study utilizes the reduced Hamiltonian formulation to separate weak nonlinear effects, specifically:
  1. Dispersion corrections (frequency shifts).
  2. Changes to wave shape (bound harmonics).
  3. Energy exchange between modes (neglected for short-time scales).
• Constant Magnitude Approximation: The author assumes that complex wave amplitudes have constant magnitudes but time-dependent phases. This allows the neglect of the slow evolution of the wave envelope (governed by the Zakharov equation) over the timescale of a few wave periods, which is sufficient for calculating particle drift.
• Particle Trajectory Mapping: Instead of relying solely on analytical approximations, the study numerically integrates the ordinary differential equations (ODEs) for particle trajectories:
  $$\mathbf{x}'(t) = \mathbf{u}(\mathbf{x}, t)$$
  where $\mathbf{u}$ is the velocity field derived from the 1st, 2nd, and 3rd order potential expansions.
• Wave Configurations: The study systematically analyzes:
  1. Monochromatic waves: Steady, periodic Stokes waves.
  2. Bichromatic waves: Unsteady waves with two harmonics (introducing difference harmonics).
  3. Multichromatic waves: Focused wave groups and random irregular seas.
• Drift Calculation: The Lagrangian mean velocity ($u_L$) is computed by averaging particle displacements over a wave period (or group period). The Stokes drift ($u_S$) is defined as the difference between Lagrangian and Eulerian mean velocities ($u_L - u_E$).

3. Key Contributions

• Derivation of a Modified Stokes Drift Formula: The paper derives a new approximation for Stokes drift in bichromatic and multichromatic waves that includes difference harmonic terms (quartic terms in amplitude). This extends the classical Kenyon (1969) formulation.
• Quantitative Validation: The study provides a rigorous numerical comparison between the classical Stokes drift, the new modified formula, and direct numerical integration of particle trajectories using 1st, 2nd, and 3rd order Eulerian velocity fields.
• Clarification of Depth Dependence: It elucidates how different nonlinear effects (bound harmonics vs. dispersion corrections) dominate at different depths.
• Spectral Analysis: The modified drift formulation is applied to parametric wave spectra (Phillips, Pierson-Moskowitz, Ochi-Hubble) to quantify the impact on total Stokes transport.

4. Key Results

A. Monochromatic Waves

• Surface Drift: The classical Stokes drift formula (based on linear theory) slightly underestimates the drift at the surface for steep waves.
• Deep Drift: The classical formula slightly overestimates the drift at depth.
• Higher-Order Effects: The inclusion of 3rd-order dispersion corrections (frequency shifts) reduces the drift velocity slightly compared to the linear prediction, bringing it closer to exact solutions (e.g., Longuet-Higgins).

B. Bichromatic and Multichromatic Waves

• Role of Difference Harmonics: In unsteady waves, the interaction between modes creates "difference harmonics" ($k_1 - k_2$). These terms decay much slower with depth than the primary free waves.
• Modified Drift Formula: The inclusion of a term representing the drift of these difference harmonics (Term II in Eq. 40/44) significantly improves agreement with numerical trajectory integration, particularly at greater depths.
• Depth Filtering:
  • Near the surface: 1st-order free waves dominate the drift.
  • At depth: 2nd-order difference harmonics dominate, creating a "return flow" (negative drift) in localized regions under wave groups, though the net drift over a full period remains positive (in the direction of propagation).
• Phase Shifts: Third-order theories introduce asymmetric frequency corrections that cause phase shifts between wave modes, complicating direct comparisons between Eulerian and Lagrangian formulations.

C. Random and Irregular Waves

• Averaging: For random seas, local "return flows" are averaged out, resulting in a net forward drift throughout the water column.
• Spectral Impact: When applied to energy spectra (e.g., Pierson-Moskowitz), including difference harmonic terms increases the surface drift by >40% for certain conditions and increases the total Stokes Transport by approximately 3.3%.
• High-Frequency Sensitivity: The results are highly sensitive to the high-frequency cut-off of the spectrum, confirming previous findings that resolving high frequencies is critical for accurate drift estimation.

5. Significance and Implications

• Operational Modeling: The findings suggest that operational ocean models (like WAVEWATCH III), which currently use the simplest infinite-depth Stokes drift formulation, may be underestimating surface drift and total transport, particularly for steep or complex sea states.
• Environmental Transport: Improved accuracy in Stokes drift calculations is crucial for modeling the transport of marine pollutants (microplastics), oil spills, and biological organisms.
• Theoretical Clarity: The work bridges the gap between Eulerian wave theories and Lagrangian particle kinematics, demonstrating that while the classical Stokes drift is a robust first approximation, incorporating difference harmonics is essential for accurate depth-dependent predictions in nonlinear wave fields.
• Limitations: The study assumes infinite depth and neglects viscosity and directional spreading. The author notes that finite depth effects are significant for bound harmonics and will be addressed in future work.

In conclusion, Stuhlmeier demonstrates that while the classical Stokes drift is remarkably accurate for linear waves, third-order nonlinear effects, particularly those arising from difference harmonics, are necessary to accurately predict particle drift at depth and total transport in realistic, unsteady sea states.