A consistent phase-averaged model of the interactions between surface gravity waves and currents

This paper presents a consistent phase-averaged model of two-way interactions between surface gravity waves and ocean currents, derived from a variational structure that ensures conservation of momentum and energy, and applies it to the generation of inertial oscillations by surface waves.

The Gist

Imagine the ocean as a giant, bustling dance floor. On this floor, there are two main groups of dancers: the Currents (the slow, steady flow of water moving across the ocean) and the Waves (the fast, rhythmic bobbing of the surface).

For a long time, scientists studied these two groups separately. They would say, "Okay, let's assume the currents are moving this way, and see how the waves react," or "Let's assume the waves are doing this, and see how they push the currents." It was like watching a dance where one partner is frozen in place while the other moves, or where they take turns leading but never actually dance together.

This paper, written by Jacques Vanneste and William Young, introduces a new way to watch the dance. They built a consistent model where the waves and currents are constantly talking to each other, pushing and pulling in real-time, while strictly obeying the fundamental laws of physics (specifically, the conservation of energy and momentum).

Here is the breakdown of their discovery using simple analogies:

1. The Problem: The "Broken" Dance

In the past, when scientists tried to simulate how waves and currents interact, they often used a "one-way street" approach.

• The Old Way: Imagine a wind blowing on a river. The wind creates waves. The old models would calculate how the river moves, then calculate how the waves move on top of that river. But they often forgot that the waves actually push back on the river!
• The Result: This created a "leaky" system. Energy would mysteriously appear or disappear in the computer simulations. It was like a dance where the partners occasionally pass through each other or lose their balance, breaking the rules of physics.

2. The Solution: The "Consistent Wave-Current Model" (CWCM)

The authors created a new set of rules (a mathematical model) where the waves and currents are locked in a perfect, two-way conversation.

• The Doppler Shift (The Moving Walkway):
  Imagine you are running on a moving walkway at an airport. If the walkway moves with you, you go faster. If it moves against you, you go slower.
  In the ocean, waves "feel" the current moving beneath them. The authors figured out exactly how to calculate this speed. They realized that the "speed" the wave feels isn't just the speed of the water at the surface, but a weighted average of the water speed all the way down to the bottom. It's like the wave is feeling the "mood" of the entire column of water, not just the surface.

• The Pseudomomentum (The Invisible Push):
  Waves carry a hidden "kick" called pseudomomentum. Think of it like a surfer carrying a heavy backpack. Even if the surfer isn't moving forward, the backpack has weight.
  When waves move, this "backpack" pushes against the current. The authors showed that this push is exactly what balances the equation. If the waves push the current one way, the current must push back, ensuring that the total energy of the system stays constant.

3. The Secret Ingredient: The "Variational Principle"

How did they ensure the math was perfect? They didn't just guess the equations; they built them from a "blueprint" called a Variational Principle.

• The Analogy: Imagine you are trying to find the most efficient path for a hiker to cross a mountain. Nature always chooses the path that uses the least amount of "effort" (or energy) to get from point A to point B. This is a fundamental rule of the universe.
• The Application: The authors started with the most basic laws of fluid motion (how water moves) and applied this "least effort" rule. By doing so, they derived their new model. Because they started with this fundamental rule, they guaranteed that their model would never break the laws of conservation. Energy and momentum are mathematically locked in; they can't vanish.

4. A Real-World Example: The "Inertial Oscillation"

To prove their model works, they revisited a famous problem from 1970 by a scientist named Hasselmann.

• The Scenario: Imagine a calm ocean. Suddenly, a strong wind blows. This wind creates waves.
• The Mystery: Hasselmann showed that these waves can actually start a slow, spinning motion in the water (called an inertial oscillation), like a spinning top.
• The Old Confusion: Previous models were confused about where the energy for this spinning came from. Did the waves give up their energy to spin the water? Or did the wind do extra work?
• The New Clarity: Using their new model, the authors showed that the waves do not lose their energy to spin the water. Instead, the wind has to do extra work to spin the water. The waves act as a messenger, telling the water to spin, but the wind pays the energy bill. The math is now crystal clear: the energy balance adds up perfectly.

Why Does This Matter?

This isn't just a math puzzle; it matters for real life.

• Climate and Weather: Ocean currents and waves play a huge role in how heat is moved around the planet. If our computer models of the ocean are "leaky" (losing energy), our weather forecasts and climate predictions will be wrong.
• Safety: Better models mean better predictions for storms, tsunamis, and how ships should navigate rough seas.

The Takeaway

Vanneste and Young have built a perfectly balanced scale for the ocean. They showed that to understand how the ocean moves, you can't look at the waves and the currents separately. You have to treat them as a single, dancing unit where every push and pull is accounted for. By using the "least effort" blueprint of nature, they created a model that is not only accurate but also respects the universe's strict rules on energy and momentum.

Technical Summary

1. Problem Statement

The paper addresses the challenge of modeling the two-way interaction between surface gravity waves and ocean currents. While the effects of currents on waves (advection and refraction via the Doppler shift) and the effects of waves on currents (via the Stokes–Coriolis and vortex forces) are individually well-understood, existing coupled models often suffer from inconsistencies.

• The Issue: Current operational models typically couple wave action transport equations with current models (like the Craik–Leibovich system) in an a posteriori manner. This often involves prescribing the Stokes velocity ($u_S$) or the current velocity ($u_E$) as fixed inputs for the other equation.
• The Consequence: This decoupling violates fundamental conservation laws. Specifically, without a consistent variational structure, the exchange of energy and momentum between waves and currents is not rigorously defined, leading to ambiguities regarding which kinetic energy definition (Eulerian $|u_E|^2/2$ vs. Lagrangian $|u_L|^2/2$) is conserved.

2. Methodology

The authors derive a Consistent Wave-Current Model (CWCM) using a rigorous variational approach based on Hamilton's principle.

• Starting Point: The derivation begins with the exact Lagrangian for rotating, incompressible Euler equations with a free surface.
• Wave-Mean Decomposition: They introduce a Lagrangian wave-mean decomposition of the flow map, $\phi = \bar{\phi} + \xi$, where $\bar{\phi}$ represents the mean current and $\xi$ represents wave displacements.
• Approximations:
  1. GLM-like Average: They adopt a specific Generalized Lagrangian Mean (GLM) framework where the mean flow map preserves volume and keeps the mean free surface flat (filtering out mean surface waves). This ensures the Lagrangian mean velocity $u_L$ is exactly divergence-free ($\nabla \cdot u_L = 0$).
  2. Linear Wave Ansatz: A WKB ansatz is applied to the wave displacement $\xi$, assuming a separation of scales between rapid phase variations and slow envelope modulations.
  3. Whitham Averaging: The Lagrangian is averaged over the wave phase to obtain an averaged Lagrangian.
• Variational Formulation: The authors apply Hamilton's principle ($\delta \int L dt = 0$) to this averaged Lagrangian. Crucially, they treat the wave spectrum as a continuum of wavetrains, introducing the wave action density $N(x, k)$ as a dynamical variable.

3. Key Contributions and Model Formulation

The resulting CWCM couples the wave action transport equation with the Craik–Leibovich (CL) system through two specific, consistent mechanisms:

A. The Doppler Shift

The absolute frequency $\Omega$ in the wave action equation is defined as:
$$\Omega(x, k) = \sigma(x, k) + k \cdot U(x, k)$$
where $\sigma$ is the intrinsic frequency. The crucial innovation is the definition of the Doppler velocity $U$:
$$U(x, \kappa) = \int_{-d}^{0} u_L(x, z) Q(x, z, \kappa) \, dz$$
Here, $u_L$ is the Lagrangian mean velocity, and $Q$ is a weight function dependent on the wavenumber $\kappa$ and depth. This ensures the Doppler shift accounts for the vertical shear of the Lagrangian current, not just the Eulerian mean.

B. The Pseudomomentum

The coupling to the current equations occurs via the pseudomomentum $p(x, z)$, which acts as a forcing term in the CL system:
$$p(x, z) = \iint Q(x, z, \kappa) N(x, k) k \, dk$$
The CL equation is formulated in terms of $u_L$ and $p$:
$$\partial_t(u_L - p) + (f + \nabla \times (u_L - p)) \times u_L + \nabla \pi = u^\circ$$
This formulation avoids the explicit use of the Stokes velocity $u_S$ or Eulerian velocity $u_E$ as primary variables, instead using $u_L$ and $p$.

4. Key Results

Conservation Laws

The variational derivation guarantees that the CWCM strictly conserves:

• Wave Action: Transported by the phase-space velocity.
• Mass: Ensured by the divergence-free constraint on $u_L$.
• Momentum: The total momentum density is shown to be $\int u_L dz$. The wave contribution appears as a flux (radiation stress) rather than a local momentum density, resolving the "wave momentum myth."
• Energy: The total conserved energy is:
  $$E = \int \frac{1}{2}|u_L|^2 \, dz + \iint \sigma N \, dk$$
  Significance: The conserved kinetic energy is the Lagrangian kinetic energy ($|u_L|^2/2$), not the Eulerian counterpart. The paper demonstrates that the Eulerian kinetic energy budget is complex and non-conservative in isolation, whereas the Lagrangian budget is simple and consistent.

Application: Hasselmann's Problem

The authors apply the model to the generation of inertial oscillations by surface waves (originally studied by Hasselmann, 1970).

• Setup: A uniform wind forcing generates a wave spectrum, which in turn drives inertial currents.
• Finding: The model confirms that the energy for the inertial oscillations does not come from the wave energy itself. Instead, the wind performs work that increases both the wave energy and the Lagrangian kinetic energy of the currents simultaneously.
• Energetics: The analysis clarifies that the "growth of swell" forces the currents, and the energy balance is only consistent when using the Lagrangian kinetic energy definition.

5. Significance and Implications

• Theoretical Consistency: The CWCM provides the first mathematically consistent framework where wave action and current dynamics are coupled without violating conservation of energy or momentum. It resolves long-standing ambiguities about the definition of mean kinetic energy in wave-current interactions.
• Variational Advantage: By deriving the model from a variational principle, the authors ensure that symmetries (and thus conservation laws via Noether's theorem) are preserved automatically, avoiding the complexity of high-order asymptotic expansions.
• Practical Impact: The model suggests that operational ocean models should prioritize the Lagrangian mean velocity $u_L$ and the pseudomomentum $p$ over Eulerian variables when coupling waves and currents. This is particularly relevant for modeling Langmuir circulations, inertial oscillations, and wave-driven mixing.
• Future Directions: The framework is robust enough to be extended to include stratification (treating buoyancy as a passive scalar) and turbulence parameterizations, providing a solid foundation for next-generation ocean circulation models.

In summary, Vanneste and Young present a rigorous, variational-based model that unifies wave and current dynamics, proving that the Lagrangian mean velocity is the natural variable for describing energy-conserving wave-current interactions.