Nonlinear evolution of unstable solar inertial modes: The case of viscous modes on a differentially rotating sphere

This paper investigates the nonlinear evolution of the Sun's most prominent high-latitude inertial mode ($m=1$) on a differentially rotating sphere, demonstrating through direct numerical simulations that it undergoes a supercritical Hopf bifurcation where Reynolds stresses smooth the differential rotation to saturate the instability at amplitudes comparable to solar observations.

The Gist

Here is an explanation of the paper, translated into everyday language with some creative analogies.

The Big Picture: The Sun's "Wobbly" Dance

Imagine the Sun not as a static ball of fire, but as a giant, spinning, fluid balloon. Because it's made of gas and liquid, it doesn't spin like a solid rock (where the top and bottom move at the same speed). Instead, the Sun's equator spins faster than its poles. This is called differential rotation.

In this fluid dance, there are invisible waves called inertial modes. Think of these like the ripples you see when you shake a bowl of Jell-O. One specific ripple, called the $m=1$ mode, is the loudest and most energetic one we see on the Sun. It's a high-latitude wave that moves back and forth near the poles.

The Problem: Why doesn't the wave just get bigger and bigger?

According to simple physics (linear theory), this specific wave should be unstable. Because the Sun spins faster at the equator and slower at the poles, it's like a shear force that keeps feeding energy into this wave.

If you push a swing at just the right time, it goes higher and higher. In the Sun's case, the differential rotation is constantly "pushing" this wave. Theoretically, the wave should grow until it tears the Sun apart or creates a massive storm. But in reality, the wave stays at a steady, manageable size (about 10 meters per second).

The Question: What stops the wave from growing forever? What is the "brake" that keeps it in check?

The Experiment: A Digital Sandbox

The authors built a computer model to simulate this. They simplified the Sun down to a 2D spinning sphere (ignoring the complex 3D heat and magnetic fields for now) and watched what happened when they turned on the "instability."

They found that the wave does indeed grow, but then it hits a limit and settles down. This is called saturation.

The Analogy: The Car on a Hill

Think of the wave like a car trying to drive up a hill.

1. The Engine (Instability): The differential rotation is the engine, pushing the car up the hill.
2. The Hill (Saturation): As the car goes faster, it starts to drag its feet on the ground (friction).
3. The Result: Eventually, the engine's push equals the friction's drag. The car stops accelerating and cruises at a steady speed.

In the Sun, the "friction" isn't just simple rubbing; it's a clever feedback loop. As the wave gets stronger, it creates a Reynolds stress. Imagine the wave as a giant hand reaching out and smoothing out the speed differences between the equator and the poles. By smoothing out the speed difference, the wave removes its own fuel source. It essentially "eats" the very thing that makes it grow, until it reaches a perfect balance.

The Math: The "Landau" Equation

The paper uses a famous mathematical tool called the Stuart-Landau equation to describe this.

• Simple Analogy: Imagine a formula that says:
  • Growth Rate = (How much you want to grow) - (How much you've already grown).
• If you are small, you grow fast.
• If you get too big, the "braking" term kicks in hard, and you stop growing.

The authors proved that this specific type of "stop-and-go" behavior is a Supercritical Bifurcation.

• Supercritical (The Good Kind): If you nudge the system slightly, it gently grows to a new, stable size. It's like a spring that stretches to a new length and stays there.
• Subcritical (The Bad Kind): If you nudge it, it might explode into chaos or jump to a completely different state.

The Sun's waves are "Supercritical." They are well-behaved. They grow, hit a limit, and settle down smoothly.

The "Echoes": Harmonics

When the main wave ($m=1$) gets strong, it doesn't just stay alone. It starts to create "echoes" or harmonics ($m=2, m=3$, etc.).

• Analogy: Think of a guitar string. When you pluck the main note (the fundamental), you also hear higher-pitched overtones.
• The paper found that these overtones are much smaller than the main wave. The second echo is about 8 times smaller, and the third is 25 times smaller.
• Crucially, the authors showed that these "echoes" are not random noise; they are perfectly organized patterns that the main wave forces into existence.

The "Real World" Connection

The team ran their simulation with a viscosity (stickiness) similar to what we think exists on the Sun's surface (caused by supergranules, which are like giant bubbles of rising gas).

The Result: The wave settled at a speed of 28 meters per second.
This is remarkably close to what astronomers actually observe on the Sun (which is around 10–20 m/s). This suggests that the simple physics of "waves smoothing out the spin" is a major reason why the Sun's waves don't go crazy.

The Catch (The "But...")

The authors are careful to say: "Don't overinterpret this."
Their model is 2D (flat). The real Sun is 3D (a sphere with depth). In 3D, the physics is more complex because heat and entropy play a huge role. The 2D model is a simplified "toy model" that captures the essence of the behavior, but the real Sun is a more complicated dance.

Summary in One Sentence

The Sun's biggest waves grow until they "smooth out" the speed differences that feed them, creating a perfect, stable balance that can be predicted by a simple mathematical rule, much like a car finding its top speed on a flat road.

Technical Summary

Here is a detailed technical summary of the paper "Nonlinear evolution of unstable solar inertial modes: The case of viscous modes on a differentially rotating sphere."

1. Problem Statement

The paper investigates the nonlinear evolution and saturation of the high-latitude inertial mode with longitudinal wavenumber $m=1$ on the Sun.

• Context: Observations show this mode has the largest amplitude on the Sun (RMS velocity $>10$ m/s). Linear theory predicts it is unstable due to shear instability caused by latitudinal differential rotation (fast equator, slow poles).
• Gap: While linear stability analysis identifies the instability, it cannot explain the observed finite amplitudes. Fully nonlinear 3D simulations of solar convection are computationally expensive and difficult to interpret regarding specific mode interactions.
• Objective: To numerically and theoretically characterize the saturation mechanism of this instability in a simplified 2D toroidal model, determining the equilibrium amplitude, the role of higher harmonics, and the validity of weakly nonlinear (WNL) theory.

2. Methodology

The authors employ a combination of Direct Numerical Simulations (DNS) and two perturbative Weakly Nonlinear (WNL) theories.

A. Physical Model & Governing Equations

• Domain: A 2D differentially rotating sphere (surface layer).
• Flow: Purely toroidal motions (incompressible viscous flow).
• Equations: The radial vorticity equation is derived from the momentum equation, incorporating:
  • Coriolis force.
  • Latitudinal differential rotation ($\Omega_0(\theta)$) inferred from HMI helioseismic data.
  • Viscosity (modeled via an Ekman number $E$).
  • The $\Lambda$-effect (non-diffusive angular momentum transport) to maintain the initial differential rotation profile.
• Control Parameter: The Ekman number $E = \nu / (r^2 \Omega_{ref})$. The study focuses on the range $10^{-3} \lesssim E < E_c \approx 1.5 \times 10^{-3}$.

B. Numerical Simulations (DNS)

• Solver: Direct Numerical Simulation in the time domain using a 3rd-order Adams–Bashforth temporal scheme and spherical harmonic decomposition (up to degree $L_{max}=200$) for spatial discretization.
• Initialization: Solar-like differential rotation profile + small random noise.
• Process: Simulations run until the system reaches a steady equilibrium state (saturation).

C. Weakly Nonlinear (WNL) Theories

Two classical perturbation methods are applied to derive the Stuart-Landau amplitude equation:
$$\frac{d|A|}{dt} = \sigma_I |A| + \beta_I |A|^3$$

1. Multiscale Expansion: Uses a small parameter $\epsilon$ representing the distance from the critical Ekman number ($E_c - E$). Time is split into fast (oscillation) and slow (amplitude evolution) scales.
2. Amplitude Expansion: Uses the mode amplitude $A$ itself as the small parameter.

• Both methods aim to calculate the linear growth rate ($\sigma_I$) and the nonlinear saturation coefficient ($\beta_I$).

3. Key Contributions

• Validation of Supercritical Hopf Bifurcation: The study confirms that the $m=1$ mode undergoes a supercritical Hopf bifurcation. The system transitions smoothly from a stable state to a finite-amplitude limit cycle as $E$ drops below $E_c$.
• Derivation of Saturation Scaling: The authors derive that the saturated amplitude scales as $|A| \propto \sigma_I^{1/2}$ (or $(E_c - E)^{1/2}$), consistent with the Landau equation where $\beta_I < 0$.
• Characterization of Harmonics: The paper demonstrates that the nonlinear interaction of the fundamental $m=1$ mode generates higher harmonics ($m=2, 3, \dots$). These are identified as asymptotic Koopman modes.
  • The $m=2$ harmonic is North-South (NS) antisymmetric.
  • The $m=3$ harmonic is NS symmetric (same as the fundamental).
  • Their amplitudes scale as $|A|^m$.
• Reynolds Stress Mechanism: The saturation mechanism is identified as the Reynolds stress generated by the mode feeding back onto the differential rotation, smoothing the shear profile until the linear growth rate becomes zero.

4. Key Results

A. Numerical Findings (DNS)

• Bifurcation Type: For $E < E_c$, the system converges to a unique equilibrium state regardless of initial perturbation strength, confirming a supercritical bifurcation.
• Amplitude: At $E = 4 \times 10^{-4}$ (consistent with solar supergranular viscosity), the saturated RMS velocity reaches 28 m/s, which is comparable to observed solar values ($\sim 10$ m/s average, up to 20 m/s).
• Harmonics: In the saturation regime, the $m=1$ mode dominates. However, significant energy is transferred to $m=2$ and $m=3$.
  • At $E = 4 \times 10^{-4}$, the $m=2$ amplitude is $\sim 8\times$ smaller than $m=1$, and $m=3$ is $\sim 25\times$ smaller.
  • As $E$ decreases further (lower viscosity), multiple modes become unstable, and the simple single-mode WNL theory breaks down.

B. Comparison with WNL Theory

• Agreement: The WNL theories (both Multiscale and Amplitude expansion) accurately predict the saturation amplitude and the scaling laws near the critical point ($E \approx E_c$).
• Coefficients: The calculated Landau coefficients ($\sigma_I, \beta_I$) from DNS match the theoretical predictions.
• Spatial Structure: The WNL theory successfully reproduces the spatial eigenfunctions of the fundamental mode and the first two harmonics. The discrepancy between WNL and DNS increases as $E$ moves further from $E_c$.
• Differential Rotation: The theory correctly predicts the modification of the background rotation profile ($\Delta \Omega$) caused by the Reynolds stress, showing good agreement with DNS.

5. Significance and Implications

• Physical Insight: The study provides a clear hydrodynamic explanation for the saturation of solar inertial modes without requiring full 3D convection simulations. It isolates the shear instability and Reynolds stress feedback as the primary drivers.
• Observational Relevance:
  • The predicted scaling ($|A| \propto \sigma_I^{1/2}$) offers a quantitative framework to interpret observed mode amplitudes.
  • The existence of harmonics suggests that peaks observed in solar power spectra at frequencies $2\times$and$3\times$the fundamental$m=1$ frequency (previously reported by Gizon et al.) may be nonlinear harmonics rather than distinct linear modes.
• Limitations & Future Work:
  • The 2D model lacks baroclinic effects and entropy transport, which are crucial in 3D models (where the instability is baroclinic).
  • The model does not yet explain the solar cycle modulation of mode amplitudes.
  • Future work should apply these tools to 3D simulations or use Dynamic Mode Decomposition (DMD) for cases with multiple unstable modes.

Conclusion

The paper establishes that the nonlinear saturation of the solar $m=1$ inertial mode is governed by a supercritical Hopf bifurcation well-described by weakly nonlinear theory. The saturation amplitude is determined by the balance between linear growth and Reynolds-stress-induced smoothing of differential rotation. This 2D framework successfully bridges linear theory and observed solar amplitudes, offering a predictive model for mode harmonics and their scaling behavior.