Topographic Rossby waves explained

Topographic Rossby waves are geophysical waves that form due to bottom irregularities. For ocean dynamics, the bottom irregularities are on the ocean floor such as the mid-ocean ridge. For atmospheric dynamics, the other primary branch of geophysical fluid dynamics, the bottom irregularities are found on land, for example in the form of mountains. Topographic Rossby waves are one of two types of geophysical waves named after the meteorologist Carl-Gustaf Rossby. The other type of Rossby waves are called planetary Rossby waves and have a different physical origin. Planetary Rossby waves form due to the changing Coriolis parameter over the earth. Rossby waves are quasi-geostrophic, dispersive waves. This means that not only the Coriolis force and the pressure-gradient force influence the flow, as in geostrophic flow, but also inertia.

Physical derivation

This section describes the mathematically simplest situation where topographic Rossby waves form: a uniform bottom slope.

Shallow water equations

A coordinate system is defined with x in eastward direction, y in northward direction and z as the distance from the earth's surface. The coordinates are measured from a certain reference coordinate on the earth's surface with a reference latitude

\varphi _{0} and a mean reference layer thickness H _{0} . The derivation begins with the shallow water equations:\begin{align} {\partial u\over\partial t}& + u{\partial u\over\partial x} + v{\partial u\over\partial y} - f _{0v} = - g{\partial η \over\partial x}\\[ 3 pt] {\partial v\over\partial t}& + u{\partial v\over\partial x} + v{\partial v\over\partial y} + f _{0u} = - g{\partial η \over\partial y}\\[ 3 pt] {\partial η \over\partial t}& + {\partial \over\partial x}hu + {\partial \over\partial y}hv =0 , \end{align}

where

u	is the velocity in the x direction, or zonal velocity
v	is the velocity in the y direction, or meridional velocity
h	is the local and instantaneous fluid layer thickness
η	is the height deviation of the fluid from its mean height
g	is the acceleration due to gravity
f _{0}	is the Coriolis parameter at the reference coordinate with f _{0} = 2 \Omega \sin(\varphi _{0)} , where \Omega is the angular frequency of the Earth and \varphi _{0} is the reference latitude

In the equation above, friction (viscous drag and kinematic viscosity) is neglected. Furthermore, a constant Coriolis parameter is assumed ("f-plane approximation"). The first and the second equation of the shallow water equations are respectively called the zonal and meridional momentum equations, and the third equation is the continuity equation. The shallow water equations assume a homogeneous and barotropic fluid.

Linearization

For simplicity, the system is limited by means of a weak and uniform bottom slope that is aligned with the y-axis, which in turn enables a better comparison to the results with planetary Rossby waves. The mean layer thickness

H for an undisturbed fluid is then defined asH = H _{0} + \alpha _{0y} with \alpha = {\left\vert \alpha _{0} \right\vert L \over H _{0}\ll1,}

where

\alpha _{0} is the slope of bottom, \alpha the topographic parameter and L the horizontal length scale of the motion. The restriction on the topographic parameter guarantees that there is a weak bottom irregularity. The local and instantaneous fluid thickness h can be written ash(x,y,t) = H _{0+\alpha} _{0} y +η (x,y,t).

Utilizing this expression in the continuity equation of the shallow water equations yields

{\partial η \over\partial t} + \left(u{\partial η \over\partial x} + v{\partial η \over\partial y}\right) +η \left({\partial u \over\partial x} + {\partial v \over\partial y}\right) + (H _{0+\alpha} _{0y)\left({\partial} u \over\partial x} + {\partial v \over\partial y}\right) + \alpha _{0v=0.}

The set of equations is made linear to obtain a set of equations that is easier to solve analytically. This is done by assuming a Rossby number Ro (= advection / Coriolis force), which is much smaller than the temporal Rossby number Ro_T (= inertia / Coriolis force). Furthermore, the length scale of

η \Delta H is assumed to be much smaller than the thickness of the fluid H . Finally, the condition on the topographic parameter is used and the following set of linear equations is obtained:\begin{align} {\partial u\over\partial t}& - f _{0v} = - g{\partial η \over\partial x}\\[ 3 pt] {\partial v\over\partial t}& + f _{0u} = - g{\partial η \over\partial y}\\[ 3 pt] {\partial η \over\partial t}& + H _{0\left({\partial} u \over\partial x} + {\partial v \over\partial y}\right) + \alpha _{0v=0. \end{align}}

Quasi-geostrophic approximation

Next, the quasi-geostrophic approximation Ro, Ro_T

\ll 1 is made, such that\begin{align} u = \bar{u} + \tilde{u} & with \bar{u} =- {g \over f _{0}{\partial} η \over \partial y}\\[ 3 pt] v = \bar{v} + \tilde{v} & with \bar{v} = {g \over f _{0}{\partial} η \over \partial x},\\[ 3 pt] \end{align}

where

\bar{u} and \bar{v} are the geostrophic flow components and \tilde{u} and \tilde{v} are the ageostrophic flow components with \tilde{u}\ll\bar{u} and \tilde{v}\ll\bar{v} . Substituting these expressions for u and v in the previously acquired set of equations, yields:\begin{align} & - {g \over f ^{2} _{0 }{\partial} η \over \partial y \partial t} + {\partial \tilde{u} \over\partial t} - f _{0\tilde{v}} = 0 \\[ 3 pt] &{g \over f ^{2} _{0 }{\partial} η \over \partial x \partial t} + {\partial \tilde{v} \over\partial t} + f _{0\tilde{u}} = 0 \\[ 3 pt] &{\partial η \over\partial t} + H _{0\left({\partial} \tilde{u} \over\partial x} + {\partial \tilde{v} \over\partial y}\right) + \alpha _{0} {g \over f _{0}{\partial} η \over \partial x} + \alpha _{0\tilde{v}=0. \end{align}}

Neglecting terms where small component terms (

\tilde{u}, \tilde{v}, {\partial\over\partial t} and \alpha _{0} ) are multiplied, the expressions obtained are:\begin{align} &\tilde{v} = - {g \over f ^{2 }{\partial} _{0} ^{2} η \over \partial y \partial t}\\[ 3 pt] &\tilde{u} = - {g \over f ^{2 }{\partial} _{0} ^{2} η \over \partial x \partial t}\\[ 3 pt] &{\partial η \over\partial t} + H _{0\left({\partial} \tilde{u} \over\partial x} + {\partial \tilde{v} \over\partial y}\right) + \alpha _{0} {g \over f _{0}{\partial} η \over \partial x} =0 . \end{align}

Substituting the components of the ageostrophic velocity in the continuity equation the following result is obtained:

{\partial η \over\partial t} - R ^{2{\partial\over\partial} t}\nabla ^{2} η + \alpha _{0{g} \over f _{0}{\partial} η \over \partial x} =0 ,

in which R, the Rossby radius of deformation, is defined as

R = {\sqrt{gH _{0}\over} f _{0}.}

Dispersion relation

Taking for

η a plane monochromatic wave of the formη = A \cos(k _{xx+k} _{yy-\omega} t + \phi),

with

A the amplitude, k _{x} and k _{y} the wavenumber in x- and y- direction respectively, \omega the angular frequency of the wave, and \phi a phase factor, the following dispersion relation for topographic Rossby waves is obtained:\omega = {\alpha _{0} g \over f _{0}{k} _{x} \over 1+ R ^{2 )}.} _{y}

If there is no bottom slope (

\alpha _{0} =0 ), the expression above yields no waves, but a steady and geostrophic flow. This is the reason why these waves are called topographic Rossby waves.

The maximum frequency of the topographic Rossby waves is

\left\vert \omega \right\vert _{max} = { \left\vert \alpha _{0} \right\vert g \over 2 \left\vert f _{0} \right\vert R },

which is attained for

k _{x} = R ^{-1} and k _{y} = 0 . If the forcing creates waves with frequencies above this threshold, no Rossby waves are generated. This situation rarely happens, unless \alpha _{0} is very small. In all other cases \left\vert \omega \right\vert _{max} exceeds \left\vert f _{0} \right\vert and the theory breaks down. The reason for this is that the assumed conditions: \alpha \ll 1 and Ro_T \ll 1 are no longer valid. The shallow water equations used as a starting point also allow for other types of waves such as Kelvin waves and inertia-gravity waves (Poincaré waves). However, these do not appear in the obtained results because of the quasi-geostrophic assumption which is used to obtain this result. In wave dynamics this is called filtering.

Phase speed

The phase speed of the waves along the isobaths (lines of equal depth, here the x-direction) is

c _{x} = {\omega \over k _{x}={\alpha} _{0g\over} f _{0}{1} \over 1 + R ^{2 )},} _{y}

which means that on the northern hemisphere the waves propagate with the shallow side at their right and on the southern hemisphere with the shallow side at their left. The equation of

c _{x} shows that the phase speed varies with wavenumber so the waves are dispersive. The maximum of c _{x} is\left\vert c _{x} \right\vert _{max} = {\alpha _{0g} \over f _{0},}

which is the speed of very long waves (

k ^{2 →} _{y} 0 ). The phase speed in the y-direction isc _{y} = {\omega \over k _{y}={k} _{x} \over k _{y}c} _{x,}

which means that

c _{y} can have any sign. The phase speed is given byc = {\omega \over k} = {k _{x} \over k}c _{x,}

from which it can be seen that

\left\vert c \right\vert \leq \left\vert c _{x} \right\vert as \left\vert k \right\vert = \sqrt{k ^{2 }} _{y} \geq \left\vert k _{x} \right\vert . This implies that the maximum of \left\vert c _{x} \right\vert is the maximum of \left\vert c \right\vert .

Analogy between topographic and planetary Rossby waves

Planetary and topographic Rossby waves are the same in the sense that, if the term

{\alpha _{0} g/f _{0}} is exchanged for - \beta ^{2} _{0 R} in the expressions above, where \beta _{0} is the beta-parameter or Rossby parameter, the expression of planetary Rossby waves is obtained. The reason for this similarity is that for the nonlinear shallow water equations for a frictionless, homogeneous flow the potential vorticity q is conserved:{dq \over dt} =0 with q = {\zeta + f \over h},

with

\zeta being the relative vorticity, which is twice the rotation speed of fluid elements about the z-axis, and is mathematically defined as\zeta = {\partial v \over \partial x} - {\partial u \over \partial y},

with

\zeta > 0 an anticlockwise rotation about the z-axis. On a beta-plane and for a linearly sloping bottom in the meridional direction, the potential vorticity becomesq = {f _{0+\beta} _{0} y + \zeta \over H _{0} + \alpha _{0} y + η }. .

In the derivations above it was assumed that

\begin{align} \beta _{0L} \ll \left\vert f _{0} \right\vert & :smallplanetarynumber \beta \\[ 3 pt] [\zeta] \ll \left\vert f _{0} \right\vert & :smallRossbynumberRo \\[ 3 pt] \alpha _{0L} \ll H _{0&:smalltopographicparameter\alpha} \\[ 3 pt] \Delta H \ll H & :lineardynamics , \end{align}

so

q = { f _{0} \left( 1+ {\beta _{0y} \over f _{0}} + {\zeta \over f _{0}} \right)\over H _{0} \left( 1+ {\alpha _{0y} \over H _{0}+} { η \over H _{0}\right)}={f} _{0\over} H _{0}} \left( 1+ {\beta _{0y} \over f _{0}} + {\zeta \over f _{0}} \right)\left( 1- {\alpha _{0y} \over H _{0}-} { η \over H _{0}} + \ldots\right),

where a Taylor expansion was used on the denominator and the dots indicate higher order terms. Only keeping the largest terms and neglecting the rest, the following result is obtained:

q ≈ {f _{0} \over H _{0}\left(1+{\beta} _{0y} \over f _{0}} + {\zeta \over f _{0}-{\alpha} _{0y} \over H _{0}} - { η \over H _{0}\right).}

Consequently, the analogy that appears in potential vorticity is that

\beta _{0/f} _{0} and - \alpha _{0/H} _{0} play the same role in the potential vorticity equation. Rewriting these terms a bit differently, this boils down to the earlier seen - \beta ^{2} _{0 R} and {\alpha _{0} g/f _{0}} , which demonstrates the similarity between planetary and topographic Rossby waves. The equation for potential vorticity shows that planetary and topographic Rossby waves exist because of a background gradient in potential vorticity.The analogy between planetary and topographic Rossby waves is exploited in laboratory experiments that study geophysical flows to include the beta effect which is the change of the Coriolis parameter over the earth. The water vessels used in those experiments are far too small for the Coriolis parameter to vary significantly. The beta effect can be mimicked to a certain degree in these experiments by using a tank with a sloping bottom. The substitution of the beta effect by a sloping bottom is only valid for a gentle slope, slow fluid motions and in the absence of stratification.

Conceptual explanation

As shown in the last section, Rossby waves are formed because potential vorticity must be conserved. When the surface has a slope, the thickness of the fluid layer

h is not constant. The conservation of the potential vorticity forces the relative vorticity \zeta or the Coriolis parameter f to change. Since the Coriolis parameter is constant at a given latitude, the relative vorticity must change. In the figure a fluid moves to a shallower environment, where h is smaller, causing the fluid to form a crest. When the height is smaller, the relative vorticity must also be smaller. In the figure, this becomes a negative relative vorticity (on the northern hemisphere a clockwise spin) shown with the rounded arrows. On the southern hemisphere this is an anticlockwise spin, because the Coriolis parameter is negative on the southern hemisphere. If a fluid moves to a deeper environment, the opposite is true. The fluid parcel on the original depth is sandwiched between two fluid parcels with one of them having a positive relative vorticity and the other one a negative relative vorticity. This causes a movement of the fluid parcel to the left in the figure. In general, the displacement causes a wave pattern that propagates with the shallower side to the right on the northern hemisphere and to the left on the southern hemisphere.

Measurements of topographic Rossby waves on earth

From 1 January 1965 till 1 January 1968, The Buoy Project at the Woods Hole Oceanographic Institution dropped buoys on the western side of the Northern Atlantic to measure the velocities. The data has several gaps because some of the buoys went missing. Still they managed to measure topographic Rossby waves at 500 meters depth. Several other research projects have confirmed that there are indeed topographic Rossby waves in the Northern Atlantic.^{[1] [2] [3]}

In 1988, barotropic planetary Rossby waves were found in the Northwest Pacific basin.^[4] Further research done in 2017 concluded that the Rossby waves are no planetary Rossby waves, but topographic Rossby waves.^[5]

In 2021, research in the South China Sea confirmed that topographic Rossby waves exist.^{[6] [7]}

In 2016, research in the East Mediterranean showed that topographic Rossby Waves are generated south of Crete due to lateral shifts of a mesoscale circulation structure over the sloping bottom at 4000 m (https://doi.org/10.1016/j.dsr2.2019.07.008).

Notes and References

1. Louis . John P. . Petrie . Brian D. . Smith . Peter C. . January 1982 . Observations of Topographic Rossby Waves on the Continental Margin off Nova Scotia . Journal of Physical Oceanography . 12 . 1 . 47–55 . 10.1175/1520-0485(1982)012<0047:ootrwo>2.0.co;2 . 1982JPO....12...47L . 0022-3670.
2. Oey . L-Y. . Lee . H-C. . December 2002 . Deep Eddy Energy and Topographic Rossby Waves in the Gulf of Mexico . Journal of Physical Oceanography . 32 . 12 . 3499–3527 . 10.1175/1520-0485(2002)032<3499:deeatr>2.0.co;2 . 2002JPO....32.3499O . 0022-3670. free .
3. Hamilton . Peter . July 2009 . Topographic Rossby waves in the Gulf of Mexico . Progress in Oceanography . 82 . 1 . 1–31 . 10.1016/j.pocean.2009.04.019 . 2009PrOce..82....1H . 0079-6611.
4. Schmitz . William J. . March 1988 . Exploration of the Eddy Field in the Midlatitude North Pacific . Journal of Physical Oceanography . 18 . 3 . 459–468 . 10.1175/1520-0485(1988)018<0459:eotefi>2.0.co;2 . 1988JPO....18..459S . 0022-3670. free .
5. Miyamoto . Masatoshi . Oka . Eitarou . Yanagimoto . Daigo . Fujio . Shinzou . Mizuta . Genta . Imawaki . Shiro . Kurogi . Masao . Hasumi . Hiroyasu . May 2017 . Characteristics and mechanism of deep mesoscale variability south of the Kuroshio Extension . Deep Sea Research Part I: Oceanographic Research Papers . 123 . 110–117 . 10.1016/j.dsr.2017.04.003 . 2017DSRI..123..110M . 0967-0637.
6. QUAN . QI . CAI . ZHONGYA . JIN . GUANGZHEN . LIU . ZHIQIANG . 2021-03-22 . Topographic Rossby Waves in the Abyssal South China Sea . Journal of Physical Oceanography . 51 . 6 . 1795 . 10.1175/jpo-d-20-0187.1 . 2021JPO....51.1795Q . 0022-3670. free .
7. Wang . Qiang . Zeng . Lili . Shu . Yeqiang . Li . Jian . Chen . Ju . He . Yunkai . Yao . Jinglong . Wang . Dongxiao . Zhou . Weidong . October 2019 . Energetic Topographic Rossby Waves in the Northern South China Sea . Journal of Physical Oceanography . 49 . 10 . 2697–2714 . 10.1175/jpo-d-18-0247.1 . 2019JPO....49.2697W . 202189386 . 0022-3670. free .