A reduced model for surface wave-current interactions without spatial scale separation

This paper presents a reduced asymptotic model for the bidirectional interaction between weakly nonlinear surface gravity waves and slowly evolving currents in rotating fluids that eliminates the need for spatial scale separation by coupling a wave amplitude equation with the Craik-Leibovich momentum framework to capture current-induced advection, refraction, and scattering while conserving wave action and energy.

The Gist

Imagine the ocean surface as a busy dance floor. On one side, you have the waves, which are fast, energetic, and constantly moving up and down. On the other side, you have the currents, which are slower, deeper flows that drift lazily across the floor.

For a long time, scientists used a popular rulebook (called the Craik–Leibovich theory) to predict how these two interact. But this old rulebook had a major flaw: it treated the waves like a fixed, unchangeable background. It was as if the dancers (waves) were just a painted backdrop on the wall, and the slow walkers (currents) could push against them, but the dancers couldn't push back. The waves were "prescribed"—meaning scientists just guessed how strong they were, rather than calculating how they actually moved.

The New Model: A Two-Way Conversation
In this paper, Onuki and Fujiwara propose a new, upgraded model. They want to treat the waves and currents as equal partners in a conversation.

Here is the core idea in simple terms:

1. The Waves Push Back: In their new model, the waves aren't just a static backdrop. They are dynamic. As the slow currents push the waves, the waves change shape and speed. Crucially, because the waves change, they generate a force (called the Stokes drift) that pushes back on the currents. It's a true two-way street.
2. No "Big vs. Small" Rule: Usually, scientists simplify math by assuming waves are tiny compared to currents, or vice versa. This new model breaks that rule. It allows the waves and currents to be the same size. This means it can accurately describe complex situations where a current swirls right next to a wave, causing the wave to bend, scatter, or speed up in complicated ways.
3. The "Narrow Band" Trick: To make the math solvable without a supercomputer, they make one specific assumption: the waves are all roughly the same "pitch" (frequency), like a choir singing the same note, even if they are singing in different directions. This allows them to track the "volume" (amplitude) of the wave field without tracking every single water molecule.

The "Energy Bank" Analogy
One of the most important claims of this paper is about conservation.

Think of the ocean system as a bank account.

• Old Model: The waves were like a gift card you couldn't spend. You could use them to move the currents, but you couldn't take energy out of the waves to change the currents, nor could you put energy back into the waves.
• New Model: The waves and currents share a single, closed bank account. If the currents slow down, the waves might speed up, and vice versa. The total amount of "energy money" in the system stays exactly the same. The authors prove mathematically that their new equations respect this rule perfectly. They also show that the "momentum" (the push) is conserved, meaning the system doesn't magically create or lose movement.

Why This Matters (According to the Paper)
The paper suggests that in the real ocean, waves aren't just passive passengers. They are active participants. When currents get turbulent (creating what scientists call "Langmuir circulation"—those long, parallel lines of foam you see on the ocean), the waves might actually be helping to drive that turbulence, not just reacting to it.

By using this new model, scientists can finally simulate a scenario where the waves and currents evolve together, feeding off each other's energy, without needing to separate them into "big" and "small" categories. It's a more honest, balanced, and energy-consistent way to look at the churning ocean surface.

In a Nutshell
The authors have built a mathematical "bridge" that connects the fast world of surface waves with the slow world of deep currents. Unlike previous models that treated waves as a fixed script, this new model lets the waves improvise and react, ensuring that energy and momentum are always accounted for, just like a perfectly balanced ledger.

Technical Summary

Technical Summary: A Reduced Model for Surface Wave–Current Interactions Without Spatial Scale Separation

Problem Statement
The interaction between surface gravity waves and mean currents is a fundamental driver of turbulent mixing in the ocean surface boundary layer, particularly through the generation of Langmuir circulations. The standard theoretical framework for this phenomenon is the classical Craik–Leibovich (CL) theory, which incorporates wave effects into the wave-averaged momentum equation via the Stokes drift and the associated vortex force. However, in many existing applications, the Stokes drift is prescribed externally rather than being solved for dynamically. This "prescribed-wave" limit treats the surface wave field as a static catalytic agent rather than an explicit dynamical energy reservoir. Consequently, standard CL models fail to capture the mutual feedback where spatially inhomogeneous and time-dependent currents advect, refract, scatter, and modulate wave amplitudes. Furthermore, existing two-way coupled theories often rely on asymptotic scale separation (e.g., WKB-type ray descriptions), which limits their ability to describe multidirectional scattering and spectral redistribution induced by currents with horizontal scales comparable to the wavelength.

Methodology
The authors propose a reduced asymptotic model that couples the CL mean-momentum equation to a physical-space amplitude equation for the wave field, eliminating the need for spatial-scale separation between the waves and the vortical flow. The derivation proceeds through the following steps:

1. Scaling and Decomposition: The governing incompressible Euler equations in a rotating frame are non-dimensionalized. The flow is decomposed into a fast oscillatory wave component and a slow vortical mean flow. The wave steepness $\epsilon$ serves as the small parameter. The mean current velocity is scaled as $O(\epsilon)$ relative to the wave orbital velocity, implying a time-scale separation of $O(\epsilon^{-2})$ between wave and mean-flow evolution.
2. Multiple-Time-Scale Expansion: A multiple-time-scale expansion is performed in powers of wave steepness. The wave field is assumed to be narrow-band in intrinsic frequency around a carrier frequency $\omega$, but not localized in propagation direction; its horizontal wavenumber support is concentrated near the circle $|\mathbf{k}| = \kappa$.
3. Quasi-Linear Closure: The derivation adopts a quasi-linear approximation where quadratic products of wave variables are neglected in the wave equation. This isolates current-induced modulation while excluding quartic resonant wave–wave interactions.
4. Solvability and Reconstitution: At second order ($O(\epsilon^2)$), a solvability condition is derived for the wave amplitude. The authors employ a "reconstitution" technique (following Wagner et al., Thomas, and Thomas & Yamada) to improve the linear dispersion relation accuracy.
5. Phenomenological Modifications: To ensure physical consistency, the derived system undergoes specific modifications:
   • Dispersion: A higher-order time-derivative term is added to the wave amplitude equation to recover a more accurate linear dispersion relation.
   • Incompressibility: The incompressibility constraint is applied to the Lagrangian-mean velocity ($\mathbf{U}_L = \mathbf{U} + \mathbf{U}_s$) rather than just the Eulerian mean, ensuring the coupled system satisfies energy conservation.
   • Stokes Correction: The interaction term is modified to include the Stokes drift ($\mathbf{U}_s$) alongside the mean flow ($\mathbf{U}$), effectively approximating the third-order Stokes correction for deep-water waves.

Key Contributions and Results
The resulting reduced model consists of a coupled system of equations:

• A wave-averaged momentum equation (CL type) where the Stokes drift is determined dynamically by the wave amplitude.
• A complex amplitude equation for the wave field that accounts for current-induced advection, refraction, and multidirectional scattering without WKB assumptions.
• An explicit expression for the Stokes drift as a quadratic functional of the wave amplitude.

The paper establishes that this system possesses rigorous conservation properties:

1. Wave Action: The wave equation admits a conserved wave action invariant, derived from the time-scale separation between waves and currents.
2. Energy: The coupled system admits a closed energy budget. The total energy invariant represents the sum of the wave field energy (proportional to wave action) and the kinetic energy of the Lagrangian-mean current. This extends the classical CL result by incorporating energetic feedback from the current growth/decay to the wave field.
3. Momentum: In the absence of rotation, the system conserves the sum of the Eulerian-mean momentum and the wave pseudomomentum.

Significance and Claims
The authors claim that this model provides a "tractable bidirectional extension" of the classical Craik–Leibovich framework. Its primary significance lies in its ability to represent regimes where current-induced wave evolution feeds back significantly on the mean flow, without relying on spatial-scale separation. By solving for the wave amplitude and mean flow simultaneously, the model treats the Stokes drift not as an external forcing but as a dynamic component of the energy reservoir.

The paper notes that while the model is currently limited to inviscid, homogeneous fluids with constant depth and weak nonlinearity, it serves as a foundation for future extensions. These include the restoration of wave–wave interactions (cubic nonlinear terms), the inclusion of viscosity and buoyancy (Boussinesq approximation), and the application to variable bathymetry and wave breaking. The authors suggest that numerical implementation of this model will allow for direct testing of the hypothesis that bidirectional wave evolution materially alters the growth, saturation, and energetics of Langmuir circulations compared to simulations using prescribed Stokes drift.