Internal Waves

While we are all familiar with the ripples on the surface of a lake, a far more pervasive and powerful type of wave moves unseen through the interior of our oceans, atmosphere, and even distant stars. These are internal waves, born from the simple fact that most fluids in nature are not uniform but are stratified, with density changing with depth. These hidden undulations are fundamental transporters of energy and momentum, shaping environments on both terrestrial and cosmic scales. Understanding them reveals a hidden layer of dynamics that governs everything from ocean currents to the life cycle of a star.

This article provides a comprehensive exploration of these fascinating phenomena. First, in the "Principles and Mechanisms" chapter, we will uncover the fundamental physics of internal waves, from the restoring force of buoyancy to their bizarre propagation rules where energy flows at a right angle to the wave's crests. Then, in the "Applications and Interdisciplinary Connections" chapter, we will journey from the deep ocean to the heart of stars, witnessing how these waves create drag on submarines, drive global ocean mixing, and even determine the ultimate fate of planets orbiting too close to their sun.

Imagine a calm lake. If you throw a stone in, you see ripples spread out on the surface. These are surface gravity waves, and they are familiar to all of us. The restoring force that makes them possible is gravity pulling the heavier water down to level the lighter air above it. But what if the density difference wasn't at the surface? What if the water itself wasn't uniform? The deep ocean, the Earth's atmosphere, and even the interiors of stars are not uniform fluids; they are ​​stratified​​, meaning their density changes with height or depth. This stratification sets the stage for a hidden, and far more pervasive, type of wave: the internal wave.

Let's begin with the simplest possible picture of stratification: a two-layer fluid, like oil floating on water, confined between two plates. The oil has a lower density $\rho_1$ and the water a higher density $\rho_2$. If we were to nudge the interface between them, creating a small bump, gravity would work to restore the flat equilibrium. The heavier water in the bump would be pulled down, and the lighter oil in the trough would be pushed up. But, like a pendulum overshooting its lowest point, the interface would oscillate, sending waves propagating along it.

By analyzing the fluid motion, one can derive the "rulebook" for these waves, known as the ​​dispersion relation​​. This mathematical expression tells us the wave's frequency $\omega$ for a given wavenumber $k$ (where the wavelength is $2\pi/k$). For our simple two-layer system, the frequency is given by a wonderfully descriptive formula:

$\omega^2 = \frac{(\rho_2-\rho_1) g k}{\rho_1\coth(kH_1)+\rho_2\coth(kH_2)}$

Here, $g$ is the acceleration due to gravity, and $H_1$ and $H_2$ are the depths of the two layers. Notice the term $(\rho_2 - \rho_1)$. This is the heart of the matter. The restoring force, and thus the wave itself, is driven by the density difference. If the densities were equal, $\omega$ would be zero, and no wave could exist. Also, because the density difference $(\rho_2 - \rho_1)$ inside the fluid is typically much smaller than the difference between water and air, the restoring force is much weaker. This means internal waves generally have much lower frequencies and larger amplitudes than the surface waves we are used to.

While the two-layer model is instructive, most real-world environments like oceans and atmospheres have density that changes smoothly with height. For such a ​​continuously stratified​​ fluid, we need a new way to quantify the "springiness" of the stratification. This is captured by a crucial parameter called the ​​Brunt-Väisälä frequency​​, denoted by $N$.

Imagine you are a scuba diver in a stably stratified patch of ocean. If you were to grab a small parcel of water and pull it downwards, it would become denser than its new surroundings and buoyancy would push it back up. It would overshoot its original position, become lighter than its surroundings, and be pulled back down. It would oscillate! The Brunt-Väisälä frequency $N$ is precisely the natural frequency of this vertical oscillation. It's a measure of the stability of the fluid: a higher $N$ means stronger stratification and a more rapid oscillation. For an internal wave to exist, the fluid must be stably stratified, meaning $N^2$ must be positive.

With the concept of the Brunt-Väisälä frequency in hand, we can now look at the dispersion relation for internal waves in a continuously stratified fluid. What we find is nothing short of bizarre, and it reveals the unique character of these waves. The relation is disarmingly simple:

$\omega = N \cos(\theta)$

Here, $\omega$ is the wave's frequency, $N$ is the local Brunt-Väisälä frequency, and $\theta$ is the angle that the wave's vector of constant phase, $\vec{k}$, makes with the vertical direction.

Let this sink in. For ordinary waves on a string or the surface of water, the frequency depends on the wavelength (the magnitude of $\vec{k}$). Short waves have high frequencies, and long waves have low frequencies. But for internal waves, the frequency depends only on the ​​angle of propagation​​ of the wave's phase. All waves of a given frequency $\omega$, regardless of their wavelength, must travel at the exact same angle to the vertical, $\theta = \arccos(\omega/N)$!. For example, if a wave's frequency is exactly half the Brunt-Väisälä frequency ($\omega = N/2$), its phase fronts must be traveling at an angle of $\pi/3$ radians ($60^\circ$) to the vertical.

This relation immediately sets a fundamental speed limit on the system. Since $\cos(\theta)$ cannot be greater than 1, the wave frequency $\omega$ can never exceed the Brunt-Väisälä frequency $N$. $N$ acts as a natural cutoff frequency for the internal wave field.

The weirdness does not stop there. If you watch the crests and troughs of an internal wave—its planes of constant phase—moving at the angle $\theta$, you might naturally assume that the wave's energy is also flowing in that direction. You would be wrong. The energy of the wave packet travels with the ​​group velocity​​, $\vec{v}_g$, and for internal waves, this direction is wildly different from the phase velocity.

In one of the most elegant and surprising results in all of fluid dynamics, it turns out that the group velocity vector is ​​perpendicular​​ to the phase velocity vector.

$\vec{v}_g \cdot \vec{k} = 0$

This means that the energy of the wave flows at a right angle to the direction the crests are moving. Imagine a series of parallel lines (the wave crests) moving diagonally upwards and to the right. The energy of this wave would be flowing diagonally downwards and to the right, exactly along the lines of the crests. This counter-intuitive behavior is a direct consequence of the strange dispersion relation. The group velocity, which tells us how energy propagates, is the gradient of the frequency in wavenumber space ($\vec{v}_g = \nabla_{\vec{k}} \omega$). A careful calculation based on the dispersion relation confirms this perpendicular relationship. The phase propagates at an angle $\theta$ to the vertical, while the energy propagates at an angle $\theta$ to the horizontal. This orthogonal relationship between phase and energy propagation is the single most important property for understanding how internal waves transport energy and momentum through the ocean and atmosphere.

Armed with these fundamental principles, we can now follow the life of an internal wave packet on its journey through a realistic, non-uniform fluid.

• ​​Birth and Character:​​ A wave is born, perhaps from the tide flowing over an undersea mountain. Its motion is a mix of horizontal and vertical sloshing. How is its energy distributed? The ratio of horizontal to vertical kinetic energy is not fixed; it depends on the wave's geometry. For a wave confined in a channel, this ratio is $(k_z/k_x)^2$, where $k_z$ and $k_x$ are the vertical and horizontal wavenumbers. This means waves that are long and flat (small $k_x$ compared to $k_z$) have motion that is overwhelmingly horizontal.

• ​​Carrying Momentum:​​ Like all waves, internal waves carry not just energy but also momentum. This isn't the familiar momentum of a billiard ball, but a more subtle quantity called ​​pseudo-momentum​​. When a wave propagates through a fluid, the organized motion of the fluid particles represents a net momentum that can be deposited into the background flow when the wave breaks or dissipates. The vertical flux of horizontal pseudo-momentum turns out to be proportional to $-k_x k_z$. This flux is how internal waves, generated near the Earth's surface, can propagate upwards into the stratosphere and deposit their momentum, driving vast, slow-moving wind patterns like the Quasi-Biennial Oscillation.

• ​​The Journey Through a Changing World:​​ The ocean and atmosphere are not uniform. The stratification, and thus the Brunt-Väisälä frequency $N(z)$, changes with depth. As a wave packet travels into a region with different $N$, it must adapt. A fundamental principle called ​​wave action conservation​​ dictates its evolution. In a slowly varying medium, this principle implies that as a wave travels into a region of stronger stratification (higher $N$), its energy density must increase. This process, known as shoaling, causes wave amplitudes to grow as they propagate, pushing them towards instability and breaking.

• ​​Reflection and Turning Points:​​ What happens if a wave with frequency $\omega$ propagates into a region where the stratification becomes so weak that the local $N(z)$ drops to the value of $\omega$? From our dispersion relation, $\omega = N(z) \cos(\theta)$, this implies $\cos(\theta)=1$, or $\theta=0^\circ$. The phase propagation becomes purely vertical. At this depth, known as a ​​turning point​​, the horizontal wavenumber $k_x$ goes to zero, and the vertical component of the group velocity also vanishes. The wave can no longer propagate vertically and is reflected, its energy sent back in another direction. This phenomenon is responsible for trapping wave energy in specific layers of the ocean and atmosphere, creating regions of intense mixing.

• ​​The Inevitable End:​​ Our discussion so far has neglected a key aspect of the real world: friction. Real fluids have viscosity ($\nu$), and real oceans and atmospheres have thermal diffusivity ($\kappa$). These effects act to damp the wave's energy, converting its organized motion into heat. The damping rate for an internal wave is found to be approximately $\frac{1}{2}(\nu+\kappa)k^2$, where $k$ is the magnitude of the total wavevector. The $k^2$ dependence is crucial: it means that small-scale waves (large $k$) are damped out extremely quickly, while large-scale waves (small $k$) can travel for thousands of kilometers across entire ocean basins with little energy loss. This is why the dominant internal waves in the ocean are the great internal tides, with wavelengths of hundreds of kilometers, while the small ripples they generate break and dissipate locally, driving much of the mixing in the deep ocean.

Finally, let us step back and view the ocean not as a stage for a single wave, but as a concert hall for a grand, chaotic symphony of them. Tides, winds, and currents are constantly generating internal waves of all frequencies, wavelengths, and directions. What does the resulting "internal wave climate" look like?

It turns out that this chaos has a profound underlying order. As waves propagate, they can amplify and steepen until they become unstable and break, much like a surface wave breaking on a beach. This wave breaking creates turbulence and mixing. This process acts as a universal cap on wave energy. The internal wave field is said to be ​​saturated​​: on average, it holds as much energy as it can before the waves break and dissipate. This saturation hypothesis leads to a remarkable prediction for the spectrum of kinetic energy in the ocean. The kinetic energy spectrum as a function of vertical wavenumber $m$, denoted $E_K(m)$, is predicted to follow a universal power law:

$E_K(m) \propto N^2 m^{-3}$

This is a simplified version of the celebrated Garrett-Munk spectrum, which has been found to describe the internal wave field in much of the world's oceans with surprising accuracy. It tells us that the energy is not distributed evenly; there is far more energy in large vertical scales (small $m$) than in small vertical scales. It is a statistical fingerprint of the ocean's hidden turmoil, a beautiful and simple law that emerges from the complex life cycles of countless individual waves, each following the strange and wonderful rules we have just explored.

Having grasped the fundamental principles of internal waves—their curious sideways propagation and their ability to carry energy and momentum—we are now equipped to go on a journey. It is a journey that will take us from the familiar dark depths of our own oceans to the fiery, turbulent hearts of distant stars. Along the way, we will discover that these unseen undulations are not mere curiosities of fluid dynamics; they are a central actor in the grand theatre of the cosmos. They are the silent, invisible orchestra conducting the flow of energy and momentum, shaping environments on scales from a few meters to millions of kilometers. What we have learned in theory, we shall now see in magnificent practice.

Let us begin at home, in the vast, stratified waters of the Earth's oceans. The sea is not a uniform bathtub; it is layered like a cake, with density changing with depth due to variations in temperature and salinity. This stratification is the perfect stage for the ballet of internal waves.

Imagine an Autonomous Underwater Vehicle (AUV) gliding silently through the pycnocline, a region where the ocean's density changes rapidly. Just as a boat leaves a V-shaped wake on the surface, this submersible vehicle carves an invisible wake in the density layers beneath. This wake is a train of internal waves, primarily "lee waves" that trail behind the vehicle. The amplitude of these waves is not negligible; a vehicle just a meter in diameter moving at a slow crawl can cause vertical displacements of the water column of dozens of centimeters.

Now, we must remember a fundamental principle of physics: there is no such thing as a free lunch. To create these waves, the AUV must expend energy. This manifests as a unique form of resistance known as ​​internal wave drag​​. For any object moving through a stratified fluid—be it a submarine, a whale, or our AUV—part of its propulsive energy is constantly being drained away to generate these waves that radiate into the deep. This wave drag is entirely separate from the familiar friction of the water against the hull. At certain speeds, especially slower ones, the energy required to push aside the internal layers of the ocean can become equal to, or even greater than, the conventional form drag. This is a critical, and often counter-intuitive, consideration for engineers designing efficient underwater vehicles. The quietest, most energy-efficient path is not always the fastest, as moving too slowly can paradoxically maximize the energy lost to the silent song of internal waves.

The dance is not limited to moving objects. The very floor of the ocean is a participant. When a steady ocean current, itself a layered river, flows over a submerged mountain range or ridge, the effect is profound. The topography forces the fluid layers up and over the obstacle, generating immense internal lee waves that can propagate vertically for kilometers and horizontally for hundreds of kilometers. This process exerts a drag force on the topography itself, effectively trying to slow the Earth's rotation, albeit by an infinitesimal amount! More importantly, this wave generation is a colossal energy transfer mechanism. The waves carry energy and momentum from the mean ocean current and deposit it elsewhere, often by breaking, much like surface waves on a beach. This "topographic wave drag" acts as a gargantuan, slow-motion egg beater, driving a significant fraction of the mixing that occurs in the deep ocean. This mixing is essential for the global ocean circulation—the great conveyor belt that transports heat around the planet and regulates our climate.

If we thought the ocean was a grand stage, let us now look inside a star, where the same play unfolds on a truly cosmic scale. It is a testament to the profound unity of physics that the same equations governing waves in a puddle also describe the inner workings of a sun. The stable "radiative zones" inside stars, where energy is transported by photons, are stratified fluids of hot plasma. They are, therefore, perfect environments for internal waves. In astrophysics, these are known as ​​internal gravity waves​​, or simply ​​g-modes​​.

By analyzing the physics, we find that these waves have a peculiar character dictated by their dispersion relation. This relation, the fundamental "rulebook" for wave motion, shows that the frequency of a g-mode depends critically on the angle of its propagation. Unlike sound waves, which travel at the same speed in all directions, internal gravity waves travel most readily in a direction that is nearly horizontal. Their frequency, $\omega$, is always less than the local buoyancy frequency, $N$. This property makes them ideal vehicles for transporting energy and momentum over vast distances laterally within a star.

But where do these stellar waves come from? The answer lies in the "boiling" convective zones found in many stars (including our Sun). These zones are regions of chaotic, churning turbulence. At the boundary where a turbulent convective zone meets a calm, stratified radiative zone, the turmoil acts like a relentless storm on the ocean surface. The convective plumes, like giant hands, constantly push and pull on the stable layers, launching a continuous spectrum of internal waves that propagate away into the star's interior. In a beautiful transformation, the chaotic kinetic energy of convection is converted into the organized, propagating energy of waves.

Of course, these waves do not travel forever. The stellar interior is not perfectly transparent to them. As a wave propagates, the periodic compressions and rarefactions of the gas cause it to emit and absorb thermal radiation. This "photon viscosity" acts as a drag on the wave, a process neatly modeled as ​​Newtonian cooling​​. It damps the wave's amplitude, causing it to deposit its energy as heat along its path. The wave's energy, born in the convective furnace, is thus gradually returned to the star as thermal energy in the quiet radiative layers.

We now arrive at the most dramatic consequences of this hidden music. Internal waves are not merely passive carriers of energy; they are active sculptors of the very stars and planetary systems through which they move.

One of their most powerful roles is the transport of angular momentum. Stars are not solid bodies; they spin at different rates at different depths. Imagine an internal wave launched from the core, propagating outwards into layers that are rotating at a different speed. As the wave travels, it approaches a location where the fluid's own rotation speed perfectly matches the horizontal speed of the wave pattern. This is a ​​critical layer​​. At this point, a remarkable thing happens: the wave is completely absorbed, dumping its entire cargo of angular momentum into the background flow. It is a perfect, one-way transfer of momentum. This mechanism is believed to be one of the most important processes for redistributing angular momentum inside stars, coupling the spin of the core to that of the envelope and explaining why the cores of some stars rotate much more slowly than we would otherwise expect.

Furthermore, as a wave travels from the dense stellar core into the more tenuous outer layers, its amplitude grows, often dramatically, to conserve its energy flux. Like the crack of a whip, the displacement of the fluid can become so large that the wave becomes unstable and ​​breaks​​. This wave breaking creates intense, localized turbulence and mixing. This is a revolutionary concept: internal waves can induce mixing in regions of a star that are supposed to be perfectly stable and unmixed. This can bring fresh hydrogen fuel down to a burning shell or dredge up the "ash" of nuclear fusion (like helium and carbon) from the core to the surface, fundamentally altering a star's evolution and observable chemical composition.

Perhaps the most astonishing application of this physics is in the realm of interacting binary stars and exoplanets. Consider a "hot Jupiter"—a gas giant planet in a tight orbit around its star. The planet's gravity raises a tidal bulge on the star. As the star rotates, this bulge is dragged through the star's stratified interior, continuously generating internal waves at the boundary of the radiative core. These tidally excited waves propagate into the core, where they damp out and deposit their energy. This entire process acts as a powerful brake, draining orbital energy from the planet. The internal waves inside the star are the agents of tidal friction, causing the planet's orbit to decay and, in many cases, spelling its eventual doom as it spirals into its parent star.

But in this cosmic dance of death, there is a surprising twist of life. The energy deposited in the star's core by these decaying tidal waves acts as a supplemental heat source. This steady "tidal heating" means the core does not need to burn its nuclear fuel quite as rapidly to maintain the star's overall luminosity. By reducing the rate of nuclear fusion, the planet's presence can actually extend the star's main-sequence lifetime. In a final, beautiful irony, the planet, through the medium of internal waves, helps its star live longer, even as the very same process seals the planet's own fate.

From the drag on a submarine to the lifetime of a star, the story of internal waves is a story of the unity of physics. The simple, elegant principles of buoyancy and pressure in a layered fluid, when applied across the vast canvas of nature, paint a picture of a universe interconnected by an invisible, dynamic, and powerful web of waves.