🔭 astrophysics

Effects of density stratification on Rossby waves in deep atmospheres

This paper derives a singularity-free radial wave equation for Rossby waves in a general stratified deep atmosphere using Lagrangian pressure fluctuations, revealing two distinct wave cavities in the Sun and suggesting that radiative interior modes may be observable at the surface due to their nearly constant vorticity in the convection zone.

Original authors: Catherine C. Blume, Bradley W. Hindman

Published 2026-04-24

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Catherine C. Blume, Bradley W. Hindman

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: Finding the Sun's "Hidden Ripples"

Imagine the Sun as a giant, spinning ball of hot gas. We know it has waves moving across its surface, kind of like ripples on a pond. Scientists have recently spotted a specific type of wave called a Rossby wave (or "r-mode"). On Earth, these are the massive weather patterns that swirl around the planet. On the Sun, they are huge, slow-moving spirals of gas.

For a long time, scientists have been trying to figure out how deep these waves go. Do they just skim the surface? Do they dive deep into the Sun's core? Or do they get stuck in the middle?

This paper is like a new map that finally shows us where these waves can travel inside the Sun.

The Problem: The "Mathematical Singularity"

Before this paper, trying to calculate how these waves move deep inside the Sun was like trying to drive a car through a wall.

The Old Way: Previous scientists tried to solve the math using variables that caused the equations to "break" or go to infinity (a singularity) right at the boundary between the Sun's stable inner core and its churning outer layer. It was like hitting a brick wall in your math homework; you couldn't get past it.
The New Trick: The authors, Catherine Blume and Bradley Hindman, decided to change the "language" of the math. Instead of measuring the movement of the gas directly, they measured the pressure changes (specifically, "Lagrangian pressure fluctuation").
The Analogy: Imagine trying to describe a crowd of people moving.
- Old method: You try to track every single person's footstep. If the crowd gets too dense, you get confused and the math breaks.
- New method: You just measure the "squeeze" or pressure of the crowd. It's much smoother, and you can see the whole picture without getting stuck.

The Discovery: Two Separate "Swimming Pools"

Using this new, smoother math, the authors discovered that the Sun isn't just one big room for these waves. It's actually like a house with two separate swimming pools that the waves can swim in, but they can't easily mix.

Pool A: The Radiative Interior (The Deep Core)
- This is the deep, stable center of the Sun.
- Here, the waves behave like they are in a deep, calm ocean. They can travel all the way down to the center.
- The Surprise: The authors found that even though these waves live deep down, their "shadow" (their vorticity or spin) stretches all the way up to the surface. It's like a deep-sea diver whose shadow is cast so clearly on the beach that you can see them from the shore. This means we might actually be able to see these deep waves on the surface!
Pool B: The Convection Zone (The Churning Top)
- This is the outer layer of the Sun where gas boils and churns like a pot of water.
- Here, the waves are trapped in a very shallow "puddle" right near the surface (the top 4 million kilometers). They can't go deeper.
- These waves are much more active and easier to spot, but they are confined to the top layer.

Why This Matters: The "Frequency Shift"

The paper also explains how the depth of the water changes the "pitch" of the wave (its frequency).

The Analogy: Think of a guitar string. If you press down on the string (making it shorter/tighter), the note goes higher. If you loosen it, the note goes lower.
In the Sun:
- Waves in the deep, stable core (Pool A) have a pitch that is almost exactly what we expect from simple 2D math. The difference is so tiny it's almost invisible.
- Waves in the churning top layer (Pool B) have a pitch that is significantly lower (more negative) than expected. The "churning" nature of the gas pulls the frequency down.

The Mystery of the "Ghost" Waves

The authors noticed something weird about the waves from the deep core (Pool A). Even though they live deep down, their "spin" (vorticity) stays almost perfectly constant as they travel up through the churning top layer.

The Metaphor: Imagine a whisper from the bottom of a canyon. Usually, the wind and noise of the canyon would scramble the whisper. But in this case, the whisper travels up the canyon without changing its tone at all.
The Implication: This suggests that the waves we are currently seeing on the Sun's surface might actually be a mix of two things: the shallow waves from the top layer and the deep waves from the core. We might be seeing a "duet" of two different types of waves that look like one.

The Bottom Line

This paper is a breakthrough because:

It fixed the math: By using a new variable (pressure), they removed the "brick walls" that stopped previous scientists.
It found two zones: It clearly showed there are two distinct places where these waves live: the deep core and the shallow surface.
It offers hope: It suggests that the deep waves from the Sun's core might be visible on the surface, which would be a huge deal. It would allow us to use these waves as a "seismic tool" to take an X-ray of the Sun's interior and measure things we can't currently see, like how hot or unstable the gas is deep inside.

In short, the authors built a better pair of glasses that lets us see the hidden, deep rhythms of our Sun.

1. Problem Statement

Rossby waves (r-modes) have recently been observed on the solar surface, offering a potential seismic diagnostic for the poorly constrained properties of the solar interior (e.g., the superadiabatic gradient). However, a fundamental gap exists in understanding their radial eigenfunctions (vertical structure).

Limitations of Prior Work:
- Geophysical models: Typically assume quasi-2D systems (shallow-water approximation) or stable, non-convective atmospheres, which do not apply to the Sun's deep, stratified interior.
- Stellar models: Often rely on asymptotic expansions or assume fully radiative or fully convective stars.
- Mathematical Singularities: Previous derivations using variables like Lagrangian displacement ( $\xi$ ) or Eulerian pressure ( $P_1$ ) result in wave equations containing the buoyancy frequency squared ( $N^2$ ) in the denominator. Since $N^2$ transitions from positive (radiative interior) to negative (convection zone), these equations possess internal singularities at the boundary where $N^2=0$ , making it difficult to solve for modes spanning the entire solar interior.
- Approximations: Many studies utilize the "traditional approximation" (neglecting horizontal Coriolis terms), which the authors argue may be invalid for Rossby waves in the solar radiative interior.

2. Methodology

The authors derive a new radial wave equation valid for a deep atmosphere with general gravitational stratification, specifically applied to the Standard Solar Model (Model S).

Governing Equations: They start with linearized, inviscid fluid equations on a $\beta$ -plane (local Cartesian coordinates) including both vertical and horizontal components of the Coriolis force (non-traditional terms).
Key Variable Transformation: Instead of using displacement or Eulerian pressure, they adopt the Lagrangian pressure fluctuation ( $\delta P$ ) as the primary working variable.
- They define $\delta P$ via the relation $\frac{\partial}{\partial t}\delta P = \frac{\partial P_1}{\partial t} - gW$ (where $W$ is vertical velocity).
- They employ a low-frequency approximation ( $\omega \ll \Omega$ ) suitable for Rossby waves, retaining terms up to first order in frequency to break pure geostrophic balance while maintaining the quasi-hydrostatic approximation.
Derivation: By taking the curl and divergence of the momentum equations and eliminating other variables, they derive a second-order Ordinary Differential Equation (ODE) solely in terms of $\delta P$ .
Numerical Solution: They solve the resulting Sturm-Liouville problem numerically using the shooting method on Model S, applying boundary conditions of $\delta P = 0$ at the surface and $\partial_z \delta P = 0$ at the center.

3. Key Contributions

Singularity-Free Formulation: The use of $\delta P$ yields a wave equation where the coefficients are non-singular even as $N^2$ changes sign. This eliminates the internal singularities that plagued previous attempts to model Rossby waves across the radiative-convection boundary.
Identification of Two Distinct Wave Cavities: The derived equation clearly reveals two separate propagation regions:
1. Radiative Interior: Where $N^2 > 0$ .
2. Convection Zone: Where $N^2 < 0$ .
Rejection of Traditional Approximation: The authors demonstrate that the "traditional approximation" (ignoring horizontal Coriolis terms) is not valid for Rossby waves in the solar radiative interior, particularly for tesseral modes, because the vertical extent of these modes is comparable to the stellar radius.
Sturm-Liouville Structure: The vertical equation is cast in a standard Sturm-Liouville form, facilitating the application of perturbation theory for future helioseismic inversions.

4. Results

Two Families of Modes:
- Radiative Interior Modes ( $n_{ri}$ ): These have positive eigenvalues (effective depth $h$ ). Their radial eigenfunctions are vertically propagating in the radiative interior and evanescent in the convection zone.
- Convection Zone Modes ( $n_{cz}$ ): These have negative eigenvalues. They are confined to a narrow cavity just below the photosphere and are evanescent in the radiative interior.
Radial Eigenfunctions (The "Surprise"):
- For radiative interior modes, the vertical vorticity eigenfunctions are nearly constant throughout the convection zone, maintaining large amplitude all the way to the photosphere. This suggests these deep-seated modes could be observable at the surface.
- For convection zone modes, the eigenfunctions are tightly confined near the surface (within $\sim 4$ Mm).
Frequency Shifts:
- Radiative Interior: The frequency shift relative to the 2D shallow-water case is negligible (order of $10^{-4}$ ), making them difficult to distinguish from the 2D limit based on frequency alone.
- Convection Zone: These modes exhibit significant frequency shifts (up to 10% for low azimuthal orders) due to the unstable stratification.
Spectral Appearance:
- Radiative interior modes have extremely small frequency spacing between radial overtones, causing them to blend into a single broad peak in realistic observations (e.g., 10-year datasets).
- Convection zone modes have larger spacing but still appear as a cluster of unresolved peaks in current observational windows.

5. Significance and Implications

Observability of Deep Modes: The finding that radiative interior modes have constant vorticity amplitudes in the convection zone challenges the assumption that these waves are trapped deep inside. It raises the possibility that current surface observations of Rossby waves may actually be detecting modes originating in the radiative interior.
Reinterpretation of Observations: The paper suggests that previous observational attempts to determine radial behavior (which showed conflicting results) might be seeing an entanglement of both convection zone and radiative interior modes.
Diagnostic Potential: The sensitivity of the eigenfrequencies to the buoyancy frequency ( $N^2$ ) and density implies that, if individual radial overtones could be resolved, solar Rossby waves could be used to measure the superadiabatic gradient in the convection zone, a parameter currently poorly constrained by p-mode helioseismology.
Theoretical Framework: The derivation provides a robust, singularity-free mathematical framework for studying Rossby waves in any deep, stratified astrophysical fluid, moving beyond the limitations of shallow-water or fully convective/radiative approximations.