Equatorial wave

Rossby-gravity waves, first observed in the stratosphere by M. Yanai,[2] always carry energy eastward. But, oddly, their 'crests' and 'troughs' may propagate westward if their periods are long enough. The eastward speed of propagation of these waves can be derived for an inviscid slowly moving layer of fluid of uniform depth H.[3] Because the Coriolis parameter (ƒ = 2Ω sin(θ) where Ω is the angular velocity of the earth, 7.2921 ${\displaystyle \times }$ 10⁻⁵ rad/s, and θ is latitude) vanishes at 0 degrees latitude (equator), the “equatorial beta plane” approximation must be made. This approximation states that “f” is approximately equal to βy, where “y” is the distance from the equator and "β" is the variation of the coriolis parameter with latitude, ${\displaystyle {\frac {\partial f}{\partial y}}=\beta }$.[1] With the inclusion of this approximation, the governing equations become (neglecting friction):

• the continuity equation (accounting for the effects of horizontal convergence and divergence and written with geopotential height):

${\displaystyle {\frac {\partial \phi }{\partial t}}+c^{2}\left({\frac {\partial v}{\partial y}}+{\frac {\partial u}{\partial x}}\right)=0}$

• the U-momentum equation (zonal wind component):

${\displaystyle {\frac {\partial u}{\partial t}}-v\beta y=-{\frac {\partial \phi }{\partial x}}}$

• the V-momentum equation (meridional wind component):

${\displaystyle {\frac {\partial v}{\partial t}}+u\beta y=-{\frac {\partial \phi }{\partial y}}}$.[3]

We may seek travelling-wave solutions of the form[4]

${\displaystyle {\begin{Bmatrix}u,v,\phi \end{Bmatrix}}={\begin{Bmatrix}{\hat {u}}(y),{\hat {v}}(y),{\hat {\phi }}(y)\end{Bmatrix}}e^{i(kx-\omega t)}}$.

Substituting this exponential form into the three equations above, and eliminating ${\displaystyle u,}$ and ${\displaystyle \phi }$ leaves us with an eigenvalue equation

${\displaystyle -{\frac {\partial ^{2}{\hat {v}}}{\partial y^{2}}}+\left({\frac {\beta ^{2}}{c^{2}}}\right)y^{2}\,{\hat {v}}=\left({\frac {\omega ^{2}}{c^{2}}}-k^{2}-{\frac {\beta k}{\omega }}\right){\hat {v}}.}$

for ${\displaystyle {\hat {v}}(y)}$. Recognizing this as the Schrödinger equation for a quantum harmonic oscillator of frequency ${\displaystyle \Omega =\beta /c}$, we know that we must have

${\displaystyle \left({\frac {\omega ^{2}}{c^{2}}}-k^{2}-{\frac {\beta k}{\omega }}\right)={\frac {\beta }{c}}(2n+1),\quad n\geq 0}$

for the solutions to tend to zero away from the equator. For each integer ${\displaystyle n}$ therefore, this last equation provides a dispersion relation linking the wavenumber ${\displaystyle k}$ to the angular frequency ${\displaystyle \omega }$.

In the special case ${\displaystyle n=0}$ the dispersion equation reduces to

${\displaystyle (\omega +ck)(\omega ^{2}-ck\omega -c\beta )=0,}$

but the root ${\displaystyle \omega =-ck}$ has to be discarded because we had to divide by this factor in eliminating ${\displaystyle u}$, ${\displaystyle \phi }$. The remaining pair of roots correspond to the Yanai or mixed Rossby-gravity mode whose group velocity is always to the east [1] and interpolates between two types of ${\displaystyle n>0}$ modes: the higher frequency Poincaré gravity waves whose group velocity can be to the east or to the west, and the low-frequency equatorial Rossby waves whose dispersion relation can be approximated as

${\displaystyle \omega ={\frac {-\beta k}{k^{2}+\beta (2n+1)/c}}.}$

.

Dispersion relations for equatorial waves with different values of ${\displaystyle n}$: The dense narrow band of low-frequency Rossby waves and the higher frequency Poincaré gravity waves are in blue. The topologically protected Kelvin and Yanai modes are highlighted in magenta

The Yanai modes, together with the Kelvin waves described in the next section, are rather special in that they are topologically protected. Their existence is guaranteed by the fact that the band of positive frequency Poincaré modes in the f-plane form a non-trivial bundle over the two-sphere ${\displaystyle {\sqrt {{\bf {k}}^{2}+f^{2}}}=1}$. This bundle is characterized by Chern number ${\displaystyle c_{1}=2}$. The Rossby waves have ${\displaystyle c_{1}=0}$, and the negative frequency Poincaré modes have ${\displaystyle c_{1}=-2.}$ Through the bulk-boundary connection[5] this necessitates the existence of two modes (Kelvin and Yanai) that cross the frequency gaps between the Poincaré and Rossby bands and are localized near the equator where ${\displaystyle f=\beta y}$ changes sign.[6][7]

Discovered by Lord Kelvin, coastal Kelvin waves are trapped close to coasts and propagate along coasts in the Northern Hemisphere such that the coast is to the right of the alongshore direction of propagation (and to the left in the Southern Hemisphere). Equatorial Kelvin waves behave somewhat as if there were a wall at the equator – so that the equator is to the right of the direction of along-equator propagation in the Northern Hemisphere and to the left of the direction of propagation in the Southern Hemisphere, both of which are consistent with eastward propagation along the equator.[1] The governing equations for these equatorial waves are similar to those presented above, except that there is no meridional velocity component ${\displaystyle v(y)}$ (that is, no flow in the north–south direction).

• the continuity equation (accounting for the effects of horizontal convergence and divergence):

${\displaystyle {\frac {\partial \phi }{\partial t}}+c^{2}{\frac {\partial u}{\partial x}}=0}$

• the u-momentum equation (zonal wind component):

${\displaystyle {\frac {\partial u}{\partial t}}=-{\frac {\partial \phi }{\partial x}}}$

• the v-momentum equation (meridional wind component):

${\displaystyle u\beta y=-{\frac {\partial \phi }{\partial y}}.}$[1]

The solution to these equations yields the following phase speed: ${\displaystyle c^{2}=gH}$; this result is the same speed as for shallow-water gravity waves without the effect of Earth's rotation.[1] Therefore, these waves are non-dispersive (because the phase speed is not a function of the zonal wavenumber). Also, these Kelvin waves only propagate towards the east (because as Φ approaches zero, y approaches infinity).[3]

Kelvin waves have been connected to El Niño (beginning in the Northern Hemisphere winter months) in recent years in terms of precursors to this atmospheric and oceanic phenomenon. Many scientists have utilized coupled atmosphere–ocean models to simulate an El Niño Southern Oscillation (ENSO) event and have stated that the Madden–Julian oscillation (MJO) can trigger oceanic Kelvin waves throughout its 30- to 60-day cycle or the latent heat of condensation can be released (from intense convection) resulting in Kelvin waves as well; this process can then signal the onset of an El Niño event.[8] The weak low pressure in the Indian Ocean (due to the MJO) typically propagates eastward into the North Pacific Ocean and can produce easterly winds.[8] These easterly winds can transfer West Pacific warm water toward the east, thereby exciting a Kelvin wave, which in this sense can be thought of as a warm-water anomaly that travels under the ocean's surface at a depth of about 150 meters.[8] This wave can be observed at the surface by a slight rise in sea surface height of about 8 cm (associated with a depression of the thermocline) and an SST increase that covers hundreds of square kilometres across the surface of the ocean.[8]

If the Kelvin wave hits the South American coast (specifically Ecuador), its warm water gets transferred upward, which creates a large warm pool at the surface.[8] That warm water also starts to flow southward along the coast of Peru and north towards Central America and Mexico, and may reach parts of Northern California; the wave can then be tracked primarily using an array of 70 buoys anchored along the entire width of equatorial Pacific Ocean, from Papua New Guinea to the Ecuador coast.[8] Temperature sensors are placed at different depths along the buoys' anchor-lines and are then able to record sub-surface water temperature.[8] The sensors send their data in real-time via satellite to a central processing facility. These temperature measurements are then compared and contrasted to historically and seasonally adjusted average water temperatures for each buoy location. Some results indicate deviations from the 'normal' expected temperatures. Such deviations are referred to as anomalies and can be thought of as either warmer-than-normal (El Niño) or cooler-than-normal (La Niña) conditions.[8]

The overall ENSO cycle can be explained as follows (in terms of the wave propagation throughout the Pacific Ocean): ENSO begins with a warm pool traveling from the western Pacific to the eastern Pacific in the form of Kelvin waves (the waves carry the warm SSTs) that resulted from the MJO.[9] After approximately 3 to 4 months of propagation across the Pacific (along the equatorial region), the Kelvin waves reach the western coast of South America and interact (merge/mix) with the cooler Peru current system.[9] This causes a rise in sea levels and sea level temperatures in the general region. Upon reaching the coast, the water turns to the north and south and results in El Niño conditions to the south.[9] Because of the changes in sea-level and sea-temperature due to the Kelvin waves, an infinite number of Rossby waves are generated and move back over the Pacific.[9] Rossby waves then enter the equation and, as previously stated, move at lower velocities than the Kelvin waves and can take anywhere from nine months to four years to fully cross the Pacific Ocean basin (from boundary to boundary).[9] And because these waves are equatorial in nature, they decay rapidly as distance from the equator increases; thus, as they move away from the equator, their speed decreases as well, resulting in a wave delay.[9] When the Rossby waves reach the western Pacific they ricochet off the coast and become Kelvin waves and then propagate back across the Pacific in the direction of the South America coast.[9] Upon return, however, the waves decrease the sea-level (reducing the depression in the thermocline) and sea surface temperature, thereby returning the area to normal or sometimes La Niña conditions.[9]

In terms of climate modeling and upon coupling the atmosphere and the ocean, an ENSO model typically contains the following dynamical equations:

• 3 primitive equations for the atmosphere (as mentioned above) with the inclusion of frictional parameterizations: 1) u-momentum equation, 2) v-momentum equation, and 3) continuity equation
• 4 primitive equations for the ocean (as stated below) with the inclusion of frictional parameterizations:
  • u-momentum,

${\displaystyle {\frac {\partial u}{\partial t}}-v\beta y={\frac {\tau _{x}}{\rho h}},}$

• • v-momentum,

${\displaystyle {\frac {\partial v}{\partial t}}-u\beta y={\frac {\tau _{y}}{\rho h}},}$

• • continuity,

${\displaystyle {\frac {\partial h}{\partial t}}+h\left({\frac {\partial u}{\partial x}}+{\frac {\partial v}{\partial y}}\right)-K_{E}T=0,}$

• • thermodynamic energy,

${\displaystyle {\frac {\partial T}{\partial t}}+u{\frac {\partial T}{\partial x}}-K_{T}h=0.}$.[10]

Note that h is the depth of the fluid (similar to the equivalent depth and analogous to H in the primitive equations listed above for Rossby-gravity and Kelvin waves), K_T is temperature diffusion, K_E is eddy diffusivity, and τ is the wind stress in either the x or y directions.

• Primitive equations
• Beta plane
• Rossby wave
• El Niño
• Equatorial Rossby wave

• Dispersion Relation Diagram for Atmospheric/Oceanic Waves