🔢 mathematics

Inverse construction of exact atmospheric models and the 3.7--4.4\,mHz acoustic resonances

This paper establishes an inverse method to construct exact atmospheric models that explain the observed 3.7–4.4 mHz acoustic resonances by demonstrating that a finite-thickness mesospheric cut-off barrier, modeled via a Rosen-Morse potential, accounts for the distinct damping behaviors of the fundamental and overtone modes through quantum-like tunnelling and leakage mechanisms.

Original authors: O. V. Kaptsov

Published 2026-10-07

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: O. V. Kaptsov

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.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

When a massive volcano erupts or a great earthquake strikes, the ground does not just shake; it sends a pulse of energy into the air above. This pulse travels upward through the layers of our atmosphere, bouncing and shifting as it encounters changes in temperature and density. Sometimes, this energy gets trapped, creating a low, humming vibration that circles the planet. Scientists have long known that these atmospheric vibrations occur at specific, steady frequencies, much like the notes produced by a musical instrument. For decades, researchers have tried to understand why these notes have the specific pitches they do and why some of them fade away quickly while others ring out for a long time. The answer lies not just in the temperature of the air, but in the invisible shape of the atmosphere itself, which acts as a complex filter for sound waves.

A new study by O. V. Kaptsov takes a fresh look at this problem by treating the atmosphere as a mathematical puzzle that can be solved in reverse. Instead of starting with a model of the air and guessing what sounds it might make, the researcher started with the sounds themselves and worked backward to build a perfect model of the air that would produce them. By using a specific mathematical approach to describe how sound waves move through a gas, the study constructs exact atmospheric profiles that match real-world observations. The goal was to isolate the vertical column of air above a single point and determine exactly how its structure controls the frequencies we hear, setting aside the complications of the Earth's curvature, winds, or the movement of the ground itself to focus purely on the vertical resonator.

The research begins by acknowledging that the atmosphere is not a simple, uniform tube. It has a "cut-off" frequency, a threshold that acts like a gatekeeper for sound waves. Waves with frequencies below this threshold can get trapped near the ground, while those above it can escape into space. Previous models suggested that if the atmosphere were perfectly smooth and transparent to these waves, it would trap the lowest frequency vibration perfectly but would not allow the higher frequency vibrations to exist in the way we observe them. In reality, however, we see both a strong, long-lasting low note and a weaker, quickly fading high note. The study demonstrates that a perfectly transparent atmosphere cannot explain this observation. To match what we actually see, the atmosphere must contain a specific type of barrier—a layer of air that is thick enough to reflect some waves but thin enough to let others leak through.

To solve this, the author built a mathematical model of the atmosphere using a shape known as a Rosen–Morse barrier. This model creates a smooth profile that closely matches the real temperature of the Earth's atmosphere from the ground up to the edge of space. When the researcher tested this model, it produced two distinct resonances that matched the real world almost exactly. The first resonance, at 3.76 millihertz, corresponds to the low-frequency hum observed after major events like the eruption of Mount Pinatubo or the 2022 Hunga Tonga explosion. In this model, this low note sits just below the top of the atmospheric barrier. Because it is trapped below the barrier, it can only escape by a slow process called tunneling, which allows it to ring out for a very long time with a high quality factor of 173.

The second resonance, at 4.37 millihertz, represents the higher-frequency note that is often seen but fades much faster. In the model, this note sits above the top of the barrier. Because it is not trapped below the gate, it bounces off the barrier and the upper atmosphere, losing energy quickly. This results in a much lower quality factor of 20, explaining why this vibration damps out so rapidly compared to the lower one. The study confirms that the difference in how long these two sounds last is not due to the air getting thicker or thinner in a simple way, but because the two frequencies occupy different positions relative to this invisible barrier. One is hidden below the gate, leaking slowly; the other is above it, bouncing away.

The power of this approach was tested against real data from seven different major events, including volcanic eruptions and earthquakes from 1982 to 2022. The researcher used computer simulations to generate the specific atmospheric conditions for the time and place of each event, then calculated the resulting frequencies. In every case where a low-frequency oscillation between 3.7 and 3.8 millihertz was identified by observers, the model successfully reproduced it. The study also found that the exact frequency of this low note is determined by the depth of the lower atmosphere and the height of the barrier, rather than by solar activity or the specific details of the upper atmosphere. This means that even as the sun's activity changes, the fundamental "note" of the atmosphere remains remarkably stable, governed by the structure of the air closest to the ground.

Ultimately, this work clarifies that the atmosphere's ability to hold onto sound is a matter of scattering and structure, not just temperature. The observed spectrum of atmospheric vibrations requires a barrier of finite thickness that connects the lower atmosphere to the thermosphere, a feature that allows the low-frequency mode to be trapped while letting the higher mode escape. By constructing these exact models, the study provides a clear, analytical explanation for why the Earth's atmosphere hums at the specific pitches it does, and why some of those hums linger while others vanish. It shows that the key to understanding these global vibrations lies in the precise shape of the atmospheric cut-off, a discovery that turns a complex physical problem into a solvable equation.

Technical Summary: Inverse Construction of Exact Atmospheric Models and the 3.7–4.4 mHz Acoustic Resonances

Problem Statement
Large explosive eruptions and earthquakes excite atmospheric oscillations near 3.7 mHz and 4.4 mHz that couple to the solid Earth, producing harmonic Rayleigh waves. While numerical normal-mode computations for spherical Earth-atmosphere models have successfully identified these frequencies, the physical mechanism governing the distinct damping characteristics of these two modes remains to be analytically isolated. Specifically, the fundamental mode (near 3.7 mHz) exhibits weak damping (high quality factor), whereas the first overtone (near 4.4 mHz) is strongly damped. Previous interpretations suggest these are acoustic modes trapped between the ground and the lower thermosphere, but the specific role of the atmospheric cut-off structure in controlling this spectrum requires clarification.

Methodology
The paper isolates the vertical atmospheric resonator, treating the atmosphere as a hydrostatic ideal-gas profile with vertically propagating linear acoustic-gravity waves. Sphericity, winds, variable composition, and coupling to the elastic Earth are excluded to focus on the vertical scattering structure.

Inverse Formulation: The linearized wave equation is transformed via travel time $\tau = \int dz/c$ and a gauge substitution into a Klein–Gordon equation, $\ddot{u} = u_{\tau\tau} - Q(\tau)u$ , where $Q(\tau)$ acts as a potential. The paper establishes a constructive inverse method: for any potential $Q(\tau)$ admitting a positive zero-energy solution $\psi$ , the corresponding atmospheric sound speed profile $c(\tau)$ and height profile $z(\tau)$ can be explicitly constructed by quadrature.
Exactly Solvable Models: The author utilizes specific potentials to generate model atmospheres:
- Darboux–Crum Potentials: Used to construct "reflectionless" (transparent) atmospheres that reproduce the U.S. Standard Atmosphere between 11 and 100 km.
- Rosen–Morse Barrier: A specific potential form $Q(\tau) = q_0 + q_1 \tanh(\tau/T_w) + q_2 \text{sech}^2(\tau/T_w)$ is fitted to the standard atmosphere to create an exactly solvable model with a mesospheric cut-off barrier. This allows for a closed-form hypergeometric resonance condition.
Validation: The exact analytical solutions are compared against full-wave numerical computations for event-specific atmospheres (derived from NRLMSIS 2.1) for seven historical events (including Pinatubo, El Chichón, and Hunga Tonga).

Key Contributions and Results

Explicit Atmospheric Reconstruction: The paper provides a general theorem (Theorem 3.1) showing that for any potential $Q$ with a positive zero-energy solution, the atmospheric profile is explicitly constructible. This extends known reflectionless profiles (where $Q$ is constant) to a broader class of potentials.
Transparency vs. Leakage: A transparent atmosphere constructed via Darboux–Crum transformations reproduces the middle atmosphere temperature profile (within 3 K) and is perfectly transparent to vertical waves above its cut-off frequency (~4.23 mHz). However, such a model supports a perfectly trapped fundamental mode but no damped overtone. This result demonstrates that the observed spectrum (a leaky fundamental and a damped overtone) cannot be explained by a realistic temperature profile alone; it requires a cut-off barrier of finite thickness that connects the lower atmosphere to the thermosphere.
Mechanism of Damping: The Rosen–Morse barrier model yields resonances at 3.76 mHz and 4.37 mHz with leakage quality factors ( $Q_{leak}$ $Q_{l e ak}$ ) of 173 and 20, respectively. The analysis reveals the physical mechanism:
- The fundamental (3.76 mHz) lies below the top of the mesospheric barrier (4.12 mHz). It is a quasi-bound state that leaks energy only via quantum tunnelling into the thermosphere, resulting in weak damping.
- The first overtone (4.37 mHz) lies above the barrier top. It is confined only by partial reflection from the barrier and the thermospheric gradient, leading to strong damping.
Event-Specific Validation: Full-wave computations for event-specific atmospheres reproduce the reported 3.7–3.8 mHz fundamental oscillation in all six events where such a mode was identified. The model confirms that resonance frequencies are set by the lower and middle atmosphere structure and are insensitive to solar activity within the tested range of NRLMSIS parameters (variations in F10.7 and Ap change frequencies by at most 0.01 mHz).
WKB Approximation: A WKB quantization condition is derived for the fundamental mode, relating its frequency to the cut-off of the lower atmosphere and the cavity depth. This recovers the "quarter-wave" picture proposed by Tahira (1995) with high precision.

Significance
The paper claims that the vertical acoustic spectrum is controlled primarily by the scattering structure of the atmospheric cut-off, rather than simply by the realism of the temperature stratification. The observed difference in damping between the 3.7 mHz and 4.4 mHz modes is a direct consequence of these modes occupying different scattering regimes relative to the mesospheric barrier: the fundamental tunnels through the barrier, while the overtone reflects off it.

The study provides an exactly solvable framework that links specific atmospheric potentials to observable resonance frequencies and quality factors, offering a rigorous analytical explanation for the distinct behaviors of atmospheric oscillations excited by major geophysical events. The results are consistent with previous numerical normal-mode computations and observational data, while providing a clearer physical interpretation of the leakage mechanisms.

More like this