Start Submission Become a Reviewer

Reading: A new approximation for the dynamics of topographic Rossby waves


A- A+
Alt. Display

Original Research Papers

A new approximation for the dynamics of topographic Rossby waves


Yosef Ashkenazy ,

Solar Energy and Environmental Physics, BIDR, Ben-Gurion University, Midreshet Ben-Gurion 84990, IL
X close

Nathan Paldor,

Institute of Earth Sciences, The Hebrew University of Jerusalem, Jerusalem 91904, IL
X close

Yair Zarmi

Solar Energy and Environmental Physics, BIDR, Ben-Gurion University, Midreshet Ben-Gurion 84990, IL
X close


A new theory of non-harmonic topographic Rossby waves over a slowly varying bottom depth of arbitrary, 1-D, profile is developed based on the linearised shallow water equations on the f-plane. The theory yields explicit approximate expressions for the phase speed and non-harmonic cross-slope structure of waves. Analytical expressions are derived in both Cartesian and Polar coordinates by letting the frequency vary in the cross-shelf direction and are verified by comparing them with the numerical results obtained by running an ocean general circulation model (the MITgcm). The proposed approximation may be suitable for studying open ocean and coastal shelf wave dynamics.

How to Cite: Ashkenazy, Y., Paldor, N. and Zarmi, Y., 2012. A new approximation for the dynamics of topographic Rossby waves. Tellus A: Dynamic Meteorology and Oceanography, 64(1), p.18160. DOI:
  Published on 01 Dec 2012
 Submitted on 14 Sep 2011

1. Introduction

Variations in bottom topography may result in second-class (low-frequency) waves, similar to planetary Rossby waves (e.g. Rhines, 1969a,b, 1989; Csanady, 1982; Pedlosky, 2003). The similarity between the two wave types becomes evident upon considering the general case in which potential vorticity varies with both the Coriolis parameter and water depth. Due to this similarity, rotating tank experiments with varying bottom topography are often used in the study of many aspects of the two types of (Rossby) waves. The potential vorticity considerations make it clear that in the Northern Hemisphere, shallower depths correspond to higher latitudes and thus topographic Rossby waves propagate with the shallower depths on their right, similar to the westward propagation of planetary Rossby waves under the influence of varying Coriolis parameter.

Topographic Rossby waves are of great importance in coastal ocean and lake dynamics. There are numerous observational indications for their existence (e.g. Csanady, 1973, 1976; Saylor et al., 1980; Raudsepp et al., 2003), and their unique dynamics have provided explanation of observed features throughout the world ocean (e.g. Ball, 1963, 1965; Shilo et al., 2007, 2008). Most theories of topographic Rossby waves have concentrated on a particular depth profile (mostly linear and rarely quadratic or exponential) that are amenable to analytical solutions, while general depth profiles received little attention. Here, we propose a simple approximation for topographic wave dynamics in an arbitrary bottom profile in one direction, assuming only that the bottom topography is slowly varying. With this approximation it is possible to estimate the evolution of an initial perturbation, even when the eigenstates of the system are not known. The approximation may be used to study more general coastal, and open ocean, topographic wave dynamics.

The paper is organised as follows: in Section 2.1, we develop the new approximation in Cartesian coordinates and then compare the analytical approximation to numerical simulation results obtained by the MITgcm (Section 2.2). In Section 3, we propose a similar approximation for a polar (cylindrical) coordinate system and a similar comparison to the numerical simulation. We conclude the study with a short discussion in Section 4.

2. A new approximation for topographic Rossby waves in a channel

2.1. Analytical approximation

2.1.1. Zero-order approximation.

We start from the linearised shallow water equations (LSWEs) in a channel of varying depth (see Fig. 1) on the f-plane (see Gill, 1982):

(1 )
(2 )
(3 )
where u and v are the zonal and meridional velocities, η is the free surface, f0 is the Coriolis parameter, g is the gravitation constant, H0 is the unperturbed mean water level and b(y) is the (arbitrary) bottom topography profile. We will see below that, in order to be handled by our proposed perturbation expansion, b(y) has to be slowly varying. We assume periodic boundary conditions in the zonal direction and vanishing meridional velocity at the zonally aligned channel walls. We note that an alternative setting is a channel that is bounded by a wall along the y direction, where the zonal velocity has to vanish. However, such a setting may result in reflecting topographic waves and fast Kelvin waves that propagate along one of the walls of the channel, which complicates the analytical and numerical treatment.

Fig. 1.   

Linear (solid) and quadratic (dashed) topographies used for the simulation of topographic Rossby waves. Both profiles span the same height variation over the same latitudinal extent. Here, H0=170 m.

By differentiating eqs. (1) and (2) with respect to t, it is possible to express the zonal and meridional velocities, u, v, as functions of η:

(4 )
(5 )
where the operator is defined by:
(6 )
It is possible to obtain a single equation for η by multiplying eq. (3) by and using eqs. (4) and (5):
(7 )
where . This equation is equivalent to eqs. (1)−(3) and is similar to the one derived in Pedlosky (2003). Gill (1982) developed an equation similar to eq. 7, but this equation is not applicable in our case, since the underlying assumption in its development was the rigid lid approximation, in which η is time-independent.

We non-dimensionalise the equation by scaling the variables as follows: , , and , where the ‘hat’ indicates non-dimensional variables and is the Rossby radius of deformation. To complete the scaling, we also use , although this scaling is not necessary for the mathematical treatment presented below. α is a parameter that characterises the bottom slope and can be chosen to be the maximum of ∣b′(y)∣. Note that the proposed scaling is invalid in the absence of bottom slope, in which case one expects to find only Kelvin and Poincaré waves. α is a small parameter and thus the above scaling of y and t implies ‘shrinking’ of these variables; this indicates that the underlying assumption of this scaling is long meridional wavelength (either compared with the zonal wavelength or compared with the Rossby radius of deformation) and slow wave frequency (compared with f0). The idea behind the different scaling of the zonal and meridional coordinates is the asymmetry of the system, in which the bottom topography is constant in x but varies with y. As it is well known and consistent with the numerical results presented below, topographic waves propagate mainly along isobaths and hardly across bathymetry lines.

For convenience, we now drop the ‘hat’ sign and, unless indicated, the variables are non-dimensional. Under the above scaling, the non-dimensional form of eq. (7) is:

(8 )
(9 )
The high-order time derivative (fast waves) and the derivatives with respect to y (short meridional waves) are multiplied by ɛ, in accordance with the rational behind the above scaling as explained earlier. If one also assumes, in addition, that the zonal wavelength is also long, then the second term in eq. (8) is also smaller than the first and third terms by order ɛ, which leads at zeroth-order approximation, to an equation equivalent to the planetary geostrophy equation (e.g. Primeau, 2002; Vallis, 2006).

Next, we assume that the solution to eq. (8) is of the form:

(10 )
Equation (8) then becomes:
(11 )
where is defined as:
(12 )
We will shortly show that ω is associated with the zeroth-order approximation of the frequency of the waves.

To allow for a perturbation series expansion, is assumed to be of the form:

(13 )
[An expansion of itself instead of the expansion in the exponential yields equivalent results.] The zeroth-order approximation (in ɛ) is obtained by using eq. (13) in eq. (11), leading to the following differential equation for :
(14 )
so that:
(15 )
It is thus clear that ω(y) may be regarded as the zeroth-order frequency. The zeroth-order approximation for the free surface is thus:
(16 )
where is an arbitrary function that should be selected so as to guarantee that the meridional velocity vanishes at the channel walls. We note, however, that since the underlying assumption is that the y-derivatives are small, is a slowly varying function. Below, we present the first-order approximation for .

To evaluate the velocities we first apply the above scaling in eqs. (4) and (5) and scale the velocity components on , the speed of gravity waves, that is and . The non-dimensional counterparts of eqs. (4) and (5) are:

(17 )
(18 )
where [the non-dimensional counterpart of the dimensional operator defined in eq. (6)] is given by:
(19 )
Thus, from eq. (18) it is clear that the zeroth-order approximation for the meridional velocity v0 is determined by the geostrophic balance, unlike the zeroth-order approximation for the zonal velocity u0 [see eq. (17)]. In addition, u0 also contains a secular term (i.e. a term that grows with time), which results from the appearance of y-derivatives on the right-hand side of eq. (17) and since ω is a function of y; see Ashkenazy et al. (2011) for further discussion regarding the appearance of secular terms in the approximation of planetary Rossby waves. Thus, the approximation is valid for limited times such that . It is possible to show that the energy of the system is preserved for each order of approximation, despite the secular terms, and is actually expected owing to the fact that the LSWEs conserve energy (e.g. Vallis, 2006) and from the nature of perturbation theory. We note that secular terms also arise from the classical planetary geostrophic approximation Ashkenazy et al., (2011).

It follows from the fact that ω depends on y that the phase speed, , varies with y, depending on the local slope, b′(y), and on the water depth, ; this is true regardless of the specific bottom topography profile, b(y). Thus, different parts of an initial perturbation will propagate with different phase speeds. This behaviour will be demonstrated and explained below.

It is possible to find the specific bottom topography for which the (zeroth-order) phase speed is constant, that is . Consistent with previous studies (Buchwald and Adams, 1968; Gill, 1982), this specific profile is an exponential profile.

For future reference we note that the dimensional frequency associated with the dispersion relation (12) is given by:

(20 )
and the zonal phase speed is found by dividing ω by k. The dimensional-free surface η0 can be easily calculated by multiplying eq. (16) by H0; the zeroth-order velocities u0, v0, may be found using eqs. (17) and (18), after neglecting terms of order ɛ and using the scaling factors given above (i.e. multiplying these relations by and , respectively).

The LSWEs (1)–(3) can be rescaled using the above scaling. Then, the small parameter multiplies only the time derivative of the zonal momentum equation. Thus, the proposed scaling is relevant when the time derivative in eq. (1) is small compared with the other term in this equation. Such a situation occurs when the wave frequency is small. According to eq. (20), the wave frequency is small when the bottom topography is slowly varying; we thus associate the small parameter α with the bottom slope. In addition, when both the zonal and meridional wavelengths are long compared with the Rossby radius of deformation, or when the water depth approaches zero (i.e. ), the dispersion relation can be approximated as . In this limit, the wave is non-dispersive in the zonal direction and long zonal wavelength (small k) results in a smaller frequency and, hence, more accurate perturbation expansion. This limit is the analogue of the extratropical planetary Rossby waves frequency under the long-wave approximation, for which the wave frequency is . Lastly, when the Rossby radius of deformation is much greater than the zonal and meridional wavelengths and the water depth is sufficiently large, the dispersion relation [eq. (20)] can be approximated as . In this limit, short zonal waves (large k) result in smaller wave frequency and hence better perturbation expansion. This limit is the analogue of equatorial planetary Rossby wave frequency, where .

2.1.2. First- and high-order approximations.

Using eqs. (11) and (13) and writing a differential equation for each order of approximation, we obtain the following general form for :

(21 )
where the coefficients are real and may be found based on using eq. (11); see Adem (1956) for a similar expansion. It is clear that , so that starting from the zeroth-order approximation it is possible to find the first-order coefficients, which, in turn, yield the second-order approximation coefficients, and so on. The general solution to eq. (21) is trivial:
(22 )
The first-order approximation coefficients are:
(23 )
where ω(y) is given in eq. (12). For the special case of , we get and takes a simpler form.

2.1.3. Some remarks.

A few comments regarding the approximations presented above are worth noting.

First, the different Ωn's contain both real and imaginary terms. The imaginary terms ( where j is odd) will lead to exponential growth or decay with time, depending on the sign of . Even if the infinite sum in eq. (13) converges, each order of approximation may contain secular terms and, thus, the approximation is valid for times such that the highest order term is smaller than its preceding term, that is .

Second, the coefficients depend on , b and ω that are functions of y. In addition, in the derivatives of are divided by , such that singular points may occur when decays to zero slower than exponential (while for an exponential decay, the ratio between the derivative of and is constant). This limitation does not affect the solution outside the vicinity of the singularity. Such singular points can be seen in the numerical example discussed below, yet, as expected, other regions are not affected by these singularities.

Third, the higher-order approximations result in a non-typical wave propagation, that is a different evolution than that of a monochromatic wave propagation where ω is constant. Thus, when computing the phase speed using, for example, longitude-time (Hovemöller) diagram, the phase speed changes with time (i.e. the crest or trough lines of the longitude-time diagram are not linear). We observe such changes of phase speed with time in the numerical simulation–this support the proposed analytical result. The reason for this cross-bathymetry phase speed dependence might be due to the set of eigenstates that make up the initial perturbation, each covering its own spatial range, that propagate in its corresponding phase speed. When starting from a true eigenstate of the system (which is hard to find for non-trivial bottom topography), there will be no change in wave frequency neither with time nor in space. The advantage of the proposed approximation is that there is no need to know the eigenstates of the system (corresponding to the bottom topography) in order to evaluate the short time evolution of an initial perturbation.

2.2. Comparing the analytical approximations with numerical simulations

We use the Massachusetts Institute of Technology general circulation model (MITgcm, see Adcroft et al., 2002, 2004) to validate the analytical approximation. We consider a channel on the f-plane (f0=10−4 s−1) covering 120 km in the zonal and meridional directions, with periodic boundary conditions in the zonal direction and walls in the meridional boundaries. Free-slip boundary conditions and vanishing meridional velocities at the channel walls are assumed. The lateral resolution is 2 km. There is only one layer and, since the MITgcm can handle partial cells, one can consider depth variations even in a single layer ocean. The maximal depth in the domain is H0=170 m, at the southern wall, and it decreases monotonically to a depth of 100 m at the northern edge, such that ɛ ≈ 3.4×10−7; linear and quadratic depth profiles are considered (Fig. 1). We use these two profiles in the numerical verification of the analytical results to demonstrate the dependence of the zonal phase speed on both the depth gradient and the actual depth [eq. (20)], so the wave characteristics in the two depth profiles are different even though the overall bottom difference across the 120 km wide channel is identical. The gravitation constant used is g=9.81 m s−2. The simulations were run for 1 week, a typical time of topographic Rossby waves evolution.

The initial conditions for the free surface of the two simulations are sinusoidal functions in both directions (two cycles in the zonal direction and half a cycle in the meridional direction), as shown in Fig. 2a; the initial velocities were calculated based on the geostrophic approximation (Fig. 2b, c), which is not an eigensolution of the governing equations. As expected from the analytical considerations presented above, the wave propagates to the West, such that the shallower depths are on the right relative to the direction of propagation. However, there is a noticeable difference between the two cases of linear and quadratic bottom topography: in the linear topography, the northward parts of the wave propagated faster than the southern parts, whereas the opposite occurred in the quadratic topography (Fig. 3). The agreement between the analytical approximation and the numerical results is good, both for the zeroth-order (Fig. 3c, d) and first-order (Fig. 3e, f) approximations. Similar comparisons between the numerical results and the analytical approximations for the zonal and meridional velocities are shown in Figs. 4 and 5, respectively. Here, the first-order approximation is closer to the numerical results. Singular behaviour may occur close to the channel's walls (Fig. 3e, 4e, f, 5e), as expected from eq. (23).

Fig. 2.   

The initial (a) free surface η, (b) zonal velocity u and (c) meridional velocity v for the linear bottom topography case. The corresponding fields of the quadratic bottom topography are very similar to the linear topography ones. Note that the proportion of the figure does not reflect the same zonal and meridional extent. Also note that the meridional wavelength is four times larger than the zonal one.

Fig. 3.   

The numerical (upper panels), zeroth-order approximation (middle panels) and first-order approximation (lower panels) free surface η after 7 d of simulation, for the linear (left panels) and quadratic (right panels) bottom topographies. The agreement between the numerical and analytical results is good. In all cases, the wave propagates westward (such that the shallower depths are on the right), where for the linear bathymetry case, faster wave propagation is observed at higher latitudes.

Fig. 4.   

As Fig. 3 for the zonal velocity u. Note the large amplitude close to the boundary of the channel (small and large y). The first-order approximation (lower panels) is more similar to the numerical results (upper panels) compared to the zero-order approximation (middle panels).

Fig. 5.   

As Fig. 3 for the meridional velocity v.

In Fig. 6, we present the numerically calculated latitudinal variation of the zonal phase speed in the two cases and the corresponding analytical estimates. The phase speed is negative, indicating that indeed the wave travels to the West (such that the shallower depths are on the right). The observation of faster (i.e. more negative) zonal phase speeds at higher y for the linear topography, and the vice versa for the quadratic topography is clearly evident. The agreement between the numerically computed zonal phase speed and the analytical speed is good, where, as expected, the first-order approximation yields a closer agreement to the numerical results. We note that the agreement between the numerical and analytical approximations becomes less accurate with time due to the presence of the secular term in the analytical solution; see Section 2.1.2. As predicted by the analytical approximation, the numerically calculated phase speed depends on time, and the one shown in Fig. 6 is based on the mean phase speed of the first 14 d of the simulations.

Fig. 6.   

Zonal phase speed as a function of latitude for the numerical (solid), zero-order (long dashed) and first-order (short dashed) results, for the (a) linear and (b) quadratic bottom topographies. The wave propagates westward as the phase speed is negative but for the linear topography the wave propagation speed increases northwards and the opposite occurs for the quadratic topography. The phase speed calculated based on the first-order approximation is closer to the numerical one compared to the phase speed calculated based on the zero-order approximation. The phase speed was calculated based on the mean of the first 14 d of simulation.

The initial states used here are eigenmodes (harmonic functions) of the quasi-geostrophic approximation. These do not preserve their structure with time, as the perturbation propagates with different speeds for different latitudes; thus, the quasi-geostrophic solution provides an inaccurate approximation for the set-up studied here. This inaccuracy is most probably due to the more significant bottom slope, which cannot be considered to be small as required in the quasi-geostrophic approximation.

The reason for the difference between the zonal phase speeds of the two runs is clear from eq. (20). For the linear topography case, only the local water depth affects the zonal phase speed and thus, as the depth decreases, the denominator becomes smaller and hence the phase speed increases. In contrast, for the quadratic topography, the local bottom slope appears in the numerator and as it decreases towards higher y, so does the phase speed. In this case, the effect of the local bottom slope in reducing the zonal phase speed towards larger y outweighs the effect of local depth that tends to increasing the zonal phase speed at larger y.

It is possible to approximate the eigenstates of LSWEs in the case of a linear depth profile using Airy Functions and Parabolic Cylinder Functions Cohen et al. [(2010); see also the planetary counterpart in Paldor and Sigalov (2008)]. These functions oscillate at large depths and decay (faster than exponential) at shallower depths. Higher eigenmodes extend farther into shallower depths, where each of them propagates with its own phase speed. These straightforward considerations can explain the different phase speeds for different y's by viewing an arbitrary initial perturbation as a collection of eigenstates, each propagating with its own speed in a different subregime of the domain, leading to different phase speeds at different latitudes.

Although we consider here the case of topographic wave within a channel, our approximation is also relevant for coastal shelf and open ocean topographic wave dynamics. The reason is that the problem that we solve is essentially an initial value problem in the meridional direction. Thus, the channel width can be as large as desired—it is possible to approximate the dynamics of initially meridionally concentrated perturbation within the channel. The perturbation can be located far from the channel walls, to mimic open ocean dynamics, or close to the shallow boundary to mimic coastal shelf dynamics. We have verified these claims by performing additional numerical simulations (results not shown) with a channel that is four times wider than the one considered in the numerical simulations discussed above. We found similar agreement between the analytical approximation and the numerical results, with better performance for the perturbation that is located far from the channel walls. As expected, the first-order approximation was more accurate than the zeroth-order one.

We performed an additional numerical simulation with a sinusoidal (two waves in the y-direction) bottom topography. Here, only the zeroth-order approximation exhibited a reasonable agreement with the numerical simulations. We speculate that the secular terms of the approximation diverge faster for rough topography profile.

It is unlikely to find exact analytical solutions of the LSWEs for arbitrary bottom topography. The analytical approximation presented here provides the wave structure and zonal phase speed for bottom topography that is uniform in the zonal direction (i.e. the long-shelf direction), but is arbitrary in the meridional (cross-shelf) direction, and the approximation is valid for small ɛ values. Note that the time of validity of the approximate solution becomes shorter as ɛ increases. In addition, the details of the wave's meridional structure are not required for calculating the (zeroth-order) dispersion relation.

3. An approximation for topographic Rossby waves in circular basins

3.1. Analytical approximation

3.1.1. Zero-order approximation.

The derivations of the previous section can be straightforwardly extended to basins other than the channel studied above. Since circular basins seem to be closer to the geometry of a typical lake geometry, we find it useful to discuss this geometry in details.

The LSWEs in polar coordinates on the f-plane and with varying depth H0 − b(r) are given by:

(24 )

where vr and vθ are the radial and azimuthal velocities, η is the free surface displacement, f0 is the constant Coriolis parameter, g is the gravitation constant, H0 is the unperturbed water depth and b(r) is the (radial) bottom topography. The independent variables are the radial, azimuthal and time coordinates (r, θ and t, respectively). Note that we assume here, similar to the previous section, that the depth changes only in the radial direction, r.

The velocities may be expressed in terms of the free surface:

(25 )
(26 )
where the time operator is given by eq. (6). From the set eq. (24) it is possible to obtain a single equation for the free surface by applying the operator to the continuity (i.e. the third) equation and using eqs. (25) and (26) to eliminate the velocity components from this equation. The result is:
(27 )
Note that the equation is singular at r=0.

We now perform the following scaling to switch to non-dimensional (or scaled) variables: , , , and , where the ‘hat’ indicates non-dimensional (scaled) variable. α is a parameter that characterises the mean bottom slope. We now omit the ‘hat’ sign and, unless otherwise indicated, the different variables designate non-dimensional variables.

Using the above scaling, eq. (27) transforms to:

(28 )
where . Using eq. (25), the boundary condition of vanishing radial velocity at the basin periphery is:
(29 )
The solution to eq. (28) is assumed to be of the form:
(30 )
and substituting this form in eq. (28) leads to:
(31 )
where ω(r) is:
(32 )
ω(r) is associated below with the zeroth-order approximation for the frequency.

We next assume that is of the form:

(33 )
To zeroth-order approximation in ɛ one obtains:
(34 )
such that may be regarded as the zeroth-order term in the power series of the frequency. is then given by:
(35 )
where should satisfy the boundary conditions of zero radial velocity at the basin periphery. Since the radial derivatives are assumed to be small (following the scaling given above), should be a slowly varying function. Below we present the first-order approximation for η (i.e. η1).

It is possible to find the zero-order velocities , using the scaling and and using eqs. (25) and (26).

Since ω is a function of r, the phase speed varies as a function of r, depending on the local slope b′(r) and local water depth . There exists a bottom topography profile for which the (zero-order) phase speed is constant, that is . This particular profile is given by:

(36 )
where C1 is an arbitrary constant. ω must be negative since for positive ω the constant C1 must vanish (or else b(r) tends infinity at r=0), and for such positive ω, b(r) > 1 implying the impossible scenario where the depth profile is larger than the total thickness of water. Thus, to zeroth-order approximation, the constant rotation of the topographic wave is clockwise (anti-cyclonic). When C1=0 (or when ), the profile is parabolic with the maximum/minimum b(r) located at the centre of the basin.

The dimensional zeroth-order frequency is given by:

(37 )
and the dimensional depth profile that yields a constant ω is:
(38 )
The zeroth-order dimensional velocities and free surface may be found straightforwardly using the above scaling, and by neglecting terms of order ɛ.

3.1.2. First- and high-order approximations.

In polar coordinates, the higher-order approximation is essentially the same as that in Cartesian coordinates given in eqs. (21) and (22), with the trivial replacement of y by r, and using eqs. (31) and (33).

The coefficients of the first-order approximations for Ω are:

(39 )
where ω(r) is given in eq. (32). Also here, for constant ω the coefficients a1,1 and a1,2 are zero while a1,0 takes a simpler form.

3.2. Comparison between the analytical approximation and numerical simulations

As in Section 2.2, we used the MITgcm for the numerical verification of our analytical approximations. The spatial dimensions of the square used in the simulation are 120×120 m with a resolution of 2 km. The Coriolis parameter is f0=10−4 s−1, the gravitation constant is g=9.81 m s−2 and H0=90 m. The (dimensional) b(r) is shown in Fig. 7, and it follows eq. (38) so that the zeroth-order frequency is constant and corresponds to a topographic wave rotation of one cycle per 10 d.

Fig. 7.   

Bottom topography b(r) as a function of the radius r of the circular basin. The dashed line indicates the unperturbed free surface, H0=90 m. This depth profile yields a constant frequency ω for the zeroth-order approximation [eq. (36)]. We choose C1=0.

In Fig. 8, we present a comparison between the numerical simulation and the zeroth-order approximation of the free surface. The simulation is initiated with a Gaussian profile in the radial direction, multiplied by a sinusoidal function in the azimuthal direction, with a wave number, k=2 (see Fig. 8a). The initial velocities were chosen to be in geostrophic balance with the free surface and the initial perturbation rotated clockwise. The numerically calculated free surface for t=3 and 8 d is shown in Fig. 8b, c. It is evident that the initial perturbation rotates with constant angular velocity for different values of r, confirming the analytical prediction of constant ω for the profile given in eq. (38). In Fig. 8d, the zeroth-order approximation free surface for t=8 d is depicted. The orientation of the numerical and analytical free surface for t=8 d (Fig. 8c, d) is similar; yet, while the zeroth-order approximation exhibits only rotation of the initial perturbation, the numerically calculated free surface exhibits also radial propagation. This difference may be attributed to viscosity and non-linear terms that are included in the numerical model. We note that the first-order approximation exhibits dependence of wave frequency on the radial coordinate; this difference between the analytical first-order approximation and the numerical results may be linked to higher wave frequency close to the periphery of the basin, making the perturbation expansion less accurate.

Fig. 8.   

(a) Initial free surface η(θ, r). (b) The free surface of the numerical simulations for t=3 d. (c) As (b) for t=8 d. (d) As (c) but for the zero-order approximation. Note that the initial perturbation rotates clockwise and completes one full cycle in approximately 10 d. Also note that angular frequency ω(r) is almost constant as a function of radius r, as predicted by the theory.

We have also compared numerical simulations (results not shown) to the zeroth-order approximation for other bottom topography profiles. As expected, the numerical simulations of these cases indicated different angular rotations at different values of r, that is ω does vary with r. The remarks listed in Section 2.1.3 for a periodic zonal channel are also relevant to the current case of a circular basin.

4. Discussion and summary

We have studied topographic Rossby waves dynamics on the f-plane using the LSWEs. A slowly varying arbitrary bottom topography profile in one of the directions was assumed. We propose a new approximation, according to which the time and cross-bathymetry directions are considered to be slow and long, that is changes in the cross-bathymetry direction are small. The small non-dimensional parameter is the square of a number representing the slope of the bottom topography (which we choose to be the maximum of the absolute value of the slope). It then follows that the cross-bathymetry coordinate appears as a parameter, that is the underlying differential equation, for each order of approximation, does not contain derivatives with respect to the cross-bathymetry coordinate. This equation leads to a wave frequency (and hence phase speed) that explicitly depends on the cross-bathymetry coordinate where the wave propagation depends on the local slope of the bottom topography and on the local water depth. [We note that it is possible to apply an alternative expansion in which the wave number depends on y (i.e. k(y)) instead of the ω(y) dependence. We however find the latter easier to interpret and implement as an initial value problem.] The analytical approximation is valid for slowly varying arbitrary initial conditions in the cross-bathymetry direction that satisfy the boundary conditions; the analytical approximation exhibits close similarity with the numerical simulations carried out using the MITgcm. The proposed approximation may be relevant to coastal shelf and open ocean wave dynamics since the channel can be very wide and the evolution of a meridionally localised perturbation is very well approximated by the analytical theory. Both linear and quadratic bottom topographies were used in the simulations, both spanning the same overall height difference. Yet, the linear bottom topography yielded faster wave propagation towards shallower depths, while the quadratic bottom topography exhibited the opposite behaviour, due to the smaller local slope at shallower depths.

We have developed a similar approximation for a circular basin using the LSWEs in polar coordinates on the f-plane and a general bottom topography profile in the radial direction. The approximation may be found useful in the study of lakes with (nearly) circular symmetry. This case shares similar characteristics with a channel of periodic boundary conditions. We have found a specific form of the bottom topography that yields the same angular rotation as a function of the radial dimension; constant (zeroth-order) topographic wave rotation may be just clockwise rotation.

It had been shown previously that when changes in bottom topography are very moderate with respect to the changes in the wavelength, it is possible to replace the mean water depth by the local water depth (e.g. Pedlosky, 2003). In essence, our proposed approximation is based on similar arguments, formulating the approximation rigorously and developing also the first-order approximation. The main advantage of the proposed approximation over existing ones is the ability to find the dispersion relation for arbitrary bottom topography. Extension of the present work to basins with elliptical shape is possible, which may be more relevant to lake dynamics in some cases.


  1. Adcroft, A, Campin, J.-M, Heimbach, P, Hill, C and Marshall, J. 2002. MITgcm Release Manual. MIT/EAPS: CambridgeMA. Online at: 

  2. AdcroftA. CampinJ.-M. HillC. MarshallJ. Implementation of an atmosphere-ocean general circulation model on the expanded spherical cube. Mon. Weather Rev. 2004; 132(12): 2845–2863. 

  3. AdemJ. A series solution for the barotropic vorticity equation and its application in the study of atmospheric vortices. Tellus. 1956; 8(3): 364–372. 

  4. Ashkenazy, Y, Paldor, N and Zarmi, Y. 2011. On the meridional structure of extra-tropical Rossby waves. Tellus A. 63(4), 817–827. 

  5. BallF. K. The effect of rotation on the simpler modes of motion of a liquid in an elliptic paraboloid. J. Fluid Mech. 1963; 22: 529–545. 

  6. BallF. K. Second-class motions of a shallow liquid. J. Fluid Mech. 1965; 23: 545–561. 

  7. BuchwaldV. T. AdamsJ. K. Propagation of continental shelf waves. Proc. R. Soc. Lond. A. 1968; 305(1481): 235–250. 

  8. CohenY. PaldorN. SommeriaJ. Laboratory experiments and a non-harmonic theory for topographic Rossby waves over a linearly sloping bottom on the f-plane. J. Fluid Mech. 2010; 645: 479–496. 

  9. CsanadyG. T. Wind-induced barotropic motions in long lakes. J. Phys. Oceanogr. 1973; 3: 429–438. 

  10. CsanadyG. T. Topographic waves in Lake Ontario. J. Phys. Oceanogr. 1976; 6: 93–103. 

  11. CsanadyG. T. Circulation in the Coastal Ocean. Reidel: Dordrecht, 1982 

  12. GillA. E. Atmosphere–Ocean Dynamics. Academic Press: London, 1982 

  13. PaldorN. SigalovA. Trapped waves on the mid-latitude beta-plane. Tellus. 2008; 60A: 742–748. 

  14. PedloskyJ. Waves in the Ocean and Atmosphere. Springer-Verlag: Berlin Heidelberg, New York, 2003 

  15. PrimeauF. Long Rossby wave basin-crossing time and the resonance of low-frequency basin modes. J. Phys. Oceanogr. 2002; 32: 2652–2665. 

  16. RaudseppU. BeletskyD. SchwabD. Basin-scale topographic waves in the Gulf of Riga. J. Phys. Oceanogr. 2003; 33: 1129–1140. 

  17. RhinesP. B. Slow oscillations in an ocean of varying depth. 1. Abrupt topography. J. Fluid Mech. 1969a; 37: 161–189. 

  18. RhinesP. B. Slow oscillations in an ocean of varying depth. 2. Islands and seamounts. J. Fluid Mech. 1969b; 37: 191–205. 

  19. RhinesP. B. Deep planetary circulation and topography: simple models of midocean flows. J. Phys. Oceanogr. 1989; 19: 1449–1470. 

  20. SaylorJ. H. HuangC. K. ReidR. Vortex modes in southern Lake Michigan. J. Phys. Oceanogr. 1980; 10: 1814–1823. 

  21. Shilo, E, Ashkenazy, Y, Rimmer, A, Assouline, S, Katsafados, P and co-authors. 2007. The effect of wind variability on topographic waves: Lake Kinneret case. J. Geophys. Res. 112, C12024. 

  22. ShiloE. AshkenazyY. RimmerA. AssoulineS. MahrerY. Wind spatial variability and topographic wave frequency. J. Phys. Oceanogr. 2008; 38: 2085–2096. 

  23. VallisG. K. Atmospheric and Oceanic Fluid Dynamics. Cambridge University Press: Cambridge, 2006 

comments powered by Disqus