Final states of two-dimensional turbulence above large-scale topography: stationary vortex solutions and barotropic stability

This paper elucidates the final states of two-dimensional topographic turbulence by demonstrating that localized vortices locked to topographic extrema follow Gaussian profiles with a "sinh"-like potential vorticity relation, and by revealing that the stability of these vortex-topography configurations depends critically on the background flow energy.

The Gist

Imagine the ocean as a giant, swirling bathtub. Now, imagine that the bottom of this bathtub isn't flat; it has a big hill in the middle and deep valleys in the corners. This is what scientists call topography.

When you stir this bathtub (adding energy), the water doesn't just swirl randomly forever. Eventually, it settles down into a specific, calm pattern. This paper is about figuring out exactly what that final, settled pattern looks like and why it happens.

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

1. The Big Picture: The "Background" vs. The "Stars"

When the ocean turbulence settles, it creates two things:

• The Background Flow: A gentle, large-scale current that covers the whole ocean floor. Think of this as the fog or the haze in the air. It follows simple, predictable rules (like a straight line).
• The Localized Vortices: These are the big, swirling storms (cyclones) and anti-storms (anticyclones). Think of these as bright stars popping out of the fog.

For a long time, scientists could describe the "fog" perfectly, but they were clueless about the "stars." They didn't know exactly how the storms formed or why they stayed in certain spots.

2. The Discovery: The "Sinh" Shape and the "Gaussian" Cookie

The authors looked at high-resolution computer simulations of this swirling ocean. They found two surprising things about the storms:

• The "Sinh" Connection: If you plot the speed of the water against the pressure, the storms don't follow a straight line. Instead, they follow a specific curve that looks like the mathematical function "sinh" (pronounced sinch). It's like saying the storms have a very specific "personality" that is different from the background fog.
• The Gaussian Shape: The storms themselves aren't messy blobs. They are shaped like perfect bell curves (or Gaussian shapes). Imagine a perfectly round, fluffy cloud or a muffin top. They are smooth, round, and centered exactly on the hills and valleys of the ocean floor.

3. The Recipe: A New Model

The authors created a new "recipe" to describe the final state of the ocean. They said:

"Take the gentle background fog (the linear part) and drop two perfect, fluffy Gaussian storms (one on the hill, one in the valley) right on top of it."

When they tested this recipe against their computer simulations, it worked perfectly. It was like finding the exact blueprint for how the ocean settles down after a storm.

4. The Mystery: Why Do Storms Sit Where They Do?

Here is the most interesting part. The storms don't just sit anywhere; they have a relationship with the ocean floor:

• Low Energy (Calm Ocean): When the water isn't moving too fast, a Cyclone (spinning one way) likes to sit on top of a Hill, and an Anticyclone (spinning the other way) likes to sit in a Valley. They are "locked" in place.
• High Energy (Wild Ocean): When the water is moving very fast, the rules flip! The Anticyclone moves to the Hill, and the Cyclone moves to the Valley.

Why does this happen?
The authors used a stability test (like a physics stress test) to figure this out.

• Think of the storm as a ball and the ocean floor as a landscape.
• In a calm ocean, a ball is stable if it sits in a bowl (valley) or on a hill, depending on how heavy it is. The "Cyclone on Hill" setup is a stable, comfortable spot.
• In a wild ocean, the landscape changes. The "Cyclone on Hill" becomes a wobbly, unstable spot where the ball wants to roll away. The "Anticyclone on Hill" becomes the new comfortable, stable spot.

5. Why Should We Care?

This isn't just about math games.

• Real Oceans: In the real ocean, we see giant, permanent storms sitting in deep ocean basins (like the Lofoten Basin in the Atlantic). This paper explains why they are there and why they stay there.
• Predicting the Future: By understanding these "final states," scientists can better predict how ocean currents will behave, which affects climate, marine life, and even how heat is transported around the planet.

The Takeaway

This paper is like solving a puzzle where the pieces are swirling water. The authors found that the final picture is made of a simple background layer plus two perfect, round storms. They also discovered that the "comfort zone" for these storms changes depending on how much energy is in the water. It's a beautiful example of how chaos (turbulence) eventually finds a simple, stable order.

Technical Summary

1. Problem Statement

The paper addresses the poorly understood nature of the final states of freely decaying two-dimensional (2D) quasi-geostrophic (QG) turbulence over topography.

• Background: Classical theories (minimum-enstrophy principle) predict a linear relationship between potential vorticity (PV) and streamfunction, leading to flows anti-correlated with topography (cyclones over depressions, anticyclones over elevations). However, this contradicts oceanic observations of long-lived anticyclones trapped in topographic depressions.
• Gap: Recent high-resolution simulations show that final states consist of a background flow and localized, long-lived vortices. While the background flow follows a linear PV-streamfunction relation, the specific structure, functional form, and stability mechanisms of the localized vortices remain unclear.
• Objective: To characterize these localized vortices, develop an empirical model for quasi-stationary final states, and explain the observed vortex-topography correlations (locking) through linear stability analysis.

2. Methodology

The authors employed a combination of high-resolution numerical simulations, statistical analysis of flow fields, empirical modeling, and linear stability theory.

• Numerical Simulations:

  • Governing Equations: Unforced, single-layer QG flow on an $f$-plane over a doubly periodic domain ($[-\pi, \pi) \times [-\pi, \pi)$).
  • Topography: An idealized large-scale sinusoidal topography $\eta(x,y) = \cos x + \cos y$, featuring one elevation (center) and one depression (corners).
  • Initial Conditions: Random relative vorticity fields concentrated around specific centroid wavenumbers ($k_{ini}$) to control initial enstrophy, with fixed energy levels (low energy $E=0.25E^\#$ and high energy $E=2E^\#$).
  • Solver: GPU-accelerated pseudo-spectral code (GeophysicalFlows.jl) with spectral filtering for dissipation.

• Data Analysis & Decomposition:

  • The flow field was decomposed into a background flow and localized vortices.
  • The background flow was assumed to satisfy a linear PV-streamfunction relation ($q_b = \mu_b \psi_b$).
  • The residual field (vortices) was analyzed via scatter plots of PV vs. streamfunction to identify functional relationships.

• Empirical Modeling:

  • Proposed a superposition model: $f_{total} = f_{background} + f_{vortex}$.
  • Background: Modeled using the linear relation derived from the topography.
  • Vortices: Modeled as Gaussian profiles centered at topographic extrema, assuming they follow a "sinh"-like PV-streamfunction relation (analogous to flat-bottom turbulence).
  • Parameters ($\mu_b$, vortex strength $\Gamma_v$, core size $L_v$) were determined by matching integral invariants (Energy, Enstrophy, Quartic Casimir) between the model and simulation data.

• Stability Analysis:

  • Conducted linear stability analysis (Rayleigh's criterion) on the constructed stationary solutions.
  • Analyzed the eigenvalue problem for perturbations around the base flow (superposition of topographic flow and Gaussian vortex) to determine stability conditions based on the signs of the background linear factor ($\mu_b$) and vortex strength ($\Gamma_v$).

3. Key Contributions

1. Identification of "sinh" Relation: Demonstrated that after subtracting the background flow, the residual PV-streamfunction relationship of localized vortices follows a "sinh" functional form ($q \sim \alpha \sinh \beta(\psi - \psi^\#)$), similar to flat-bottom turbulence.
2. Empirical Stationary Solution: Developed the first explicit empirical model for quasi-stationary final states in topographic turbulence. The model successfully reconstructs the global flow field by superimposing a linear background flow with Gaussian vortices locked to topographic extrema.
3. Stability Mechanism for Vortex Locking: Provided a theoretical explanation for the observed vortex-topography correlations via linear stability analysis. The study proves that the stability of a vortex depends on the background energy level (represented by $\mu_b$).
4. Robustness Verification: Confirmed that the "sinh" trend and Gaussian profiles persist even in complex random topography and high-energy regimes, despite the loss of strict stationarity.

4. Key Results

• Low-Energy Regime ($E < E^\#$):

  • Phenomenology: Two quasi-stationary vortices emerge: a cyclone locked to the topographic elevation (center) and an anticyclone locked to the depression (corners).
  • Model Performance: The empirical model (Linear Background + Gaussian Vortices) accurately reproduces the PV, streamfunction, and velocity fields.
  • Stability: Linear stability analysis reveals that the cyclone/elevation and anticyclone/depression configurations are stable when the background energy is low ($\mu_b > 0$). Conversely, the anti-correlated configurations are unstable.

• High-Energy Regime ($E > E^\#$):

  • Phenomenology: Vortices exhibit different behaviors. When initialized with large-scale structures, an anti-correlation emerges: anticyclones lock to elevations and cyclones to depressions (though with more unsteadiness).
  • Stability: Stability analysis shows that for high background energy ($\mu_b < 0$), the stability conditions reverse. The anticyclone/elevation and cyclone/depression configurations become stable, explaining the observed locking in high-energy simulations.

• General Features:

  • The "sinh" relationship and Gaussian vortex shapes are robust features across different topographies (random vs. sinusoidal) and energy levels.
  • The background flow acts to homogenize PV, reducing the linear coefficient $\mu_b$ as vortex strength increases.

5. Significance

• Theoretical Advancement: This work bridges the gap between the minimum-enstrophy principle (which fails to predict vortex locking) and observational reality. It provides a concrete mathematical framework (stationary solutions) for the final states of 2D topographic turbulence.
• Oceanic Relevance: The findings explain the formation and persistence of quasi-permanent anticyclones in oceanic basins (e.g., the Lofoten Basin) and the phenomenon of eddy saturation. The stability analysis offers a mechanism for why certain vortex-topography pairings are preferred over others depending on the energy state of the system.
• Predictive Capability: The empirical model allows for the explicit construction of final states based on integral invariants, offering a tool for predicting flow structures in geophysical fluid dynamics without running full transient simulations.
• Future Directions: The authors suggest extending this barotropic framework to baroclinic (stratified) systems, which is crucial for realistic ocean modeling, noting that stratification may introduce additional complexity to these stability mechanisms.