## 1. Introduction

Study of observed subinertial ocean current variability shows a tendency for vertical structures not to fit a simple single baroclinic modal structure. For example, Wunsch (1997) showed that in many cases, both the barotropic and first baroclinic modes are present and that they are sometimes phase locked in the sense that enhances surface currents and thus weakens near-bottom currents (and sometimes in the opposite sense as well). Sanchez de la Lama et al. (2016, hereinafter SLF16) more recently summarized observed vertical structures using empirical orthogonal functions that very often have weak near-bottom velocities. Further, they showed that the variability is often well-represented by nontraditional (“rough bottom”) baroclinic modes that have zero horizontal velocity at the bottom. Indeed, Samelson (1992) shows that a corrugated bottom can lead to weakened near-bottom currents, and SLF16 point out that a uniformly sloping bottom (e.g., Rhines 1970) can, depending on propagation direction and wave scale, also lead to vertical structures with weakened deep currents—hence the expression “rough-bottom modes.” However, SLF16 show that there is not a clear statistical relationship between observed bottom roughness and modified modal structure, a finding that LaCasce (2017) rationalizes in terms of the ubiquity of bottom slopes.

It seems timely to ask whether bottom friction might play a role in weakening deep ocean currents. For example, simple two-layer linear wave problems show that friction decreases deep currents at low frequency (e.g., Allen 1984), and various numerical calculations of nonlinear flows demonstrate that bottom friction also affects the vertical structure of eddy processes (e.g., LaCasce and Brink 2000; Arbic and Flierl 2004; Trossman et al. 2017; Brink 2017). However, it is not completely clear how bottom friction will affect traditional baroclinic modes per se. The present study thus addresses the topic of how strongly bottom friction affects baroclinic modal structures. The question is approached in the context of linear Rossby waves in a continuously stratified ocean.

## 2. Methodology and a calculation

### a. Formulation

*p*is pressure,

*f*is the Coriolis parameter,

*β*is the northward gradient of

*f*,

*N*is the buoyancy frequency, and

*Q*represents an unspecified potential vorticity source. The (

*x*,

*y*,

*z*) coordinates are eastward, northward, and vertical, while

*t*is time. Subscripted independent variables represent partial differentiation. The equation is to be solved subject to a rigid lid and (assuming infinitesimally thin turbulent boundary layers) an Ekman compatibility surface condition:

*τ*

_{0}

^{x},

*τ*

_{0}

^{y}) is the wind stress vector. At the sloping bottom, there is also an Ekman compatibility condition:

*α*

^{x}and

*α*

^{y}are the bottom slopes in the

*x*and

*y*directions, respectively (i.e.,

*h*=

*h*

_{0}+

*α*

^{x}

*x*+

*α*

^{y}

*y*, with

*|α*

^{x}

*x|*,

*|α*

^{y}

*y|*≪

*h*

_{0}), and

*r*is the effective bottom resistance coefficient. In the following, the bottom will often be taken to be flat so that

*α*

^{x}=

*α*

^{y}= 0.

### b. Vertical modal solutions

*n*is taken to have the form

*F*

_{nz}) nonzero at the top of the infinitesimal bottom boundary layer (

*z*= −

*h*

_{0}) as long as Γ ≠ 0, for example, in order to accommodate Ekman pumping when

*r*≠ 0. Thus, the symmetry of the upper and bottom boundary conditions is broken. Further, when bottom friction is nonzero, Γ,

*ζ*

_{n},

*ω*

_{n}, and

*F*

_{n}are complex and wave phase is no longer uniform in the vertical.

### c. A nondimensional parameter

*ω*

_{0R}(the largest frequency, or shortest time scale, available for any wavenumber pairing). Thus, a nondimensional parameter quantifying the frictional effect is

*N*

_{B}is the buoyancy frequency at the bottom (

*z*= −

*h*

_{0}). This expression can be thought of as consisting of two ratios. The first,

*r*(

*h*

_{0}

*ω*

_{0R})

^{−1}, is the ratio of the wave time scale to the barotropic spindown time, a measure of the overall importance of bottom friction. The remaining factors in (7) constitute a Burger number for the wave, that is, the ratio of the natural vertical scale in a stratified system (the deformation scale) to the water depth. If the deformation scale is large relative to the water depth, frictional effects are distributed throughout the water column. If the deformation scale is less than the water depth, frictional effects will be trapped near the bottom. This interpretation emphasizes that frictionally adjusted modal structures depend on both the relative strength of damping and the relative importance of the stratification at a given length scale. Thus, for

*R*≪ 1, modal structures are unperturbed by bottom friction, but for

*R*≥

*O*(1), frictional effects lead to substantial modal adjustments.

*F*

_{n}and integrating over depth. The result is

*O*(1), where

*N*

_{B}is extremely small, as is the case with an exponential stratification.

At first, it may be unintuitive that the frictional effect increases with stratification. After all, stratification inhibits the vertical velocity that is the messenger of spindown in the water column. This is indeed the case, but stronger stratification also leads to an increasing tendency to trap any spindown processes closer to the bottom. Further, stronger stratification, by trapping spindown closer to the bottom, can inhibit near-bottom horizontal flow and hence decrease bottom stresses. Thus, since the bottom boundary condition [(5b)] deals only with conditions very near the bottom, there is no contradiction. Moreover, the frictional modifications to the vertical modal structure often serve to decrease wave damping, as will be seen below.

### d. Asymptotics: Barotropic and bottom-trapped modes

*N*

^{2}=

*F*that are either barotropic or bottom intensified. [Choosing this hyperbolic cosine form in (9a)—as opposed to a simple cosine form—is somewhat arbitrary, since

*γ*is complex, but treating the forms separately is thought to improve clarity. This is equivalent to choosing Real(

*μ*

_{n}

^{2}) < 0 in (4b).] Substitution of (9a) into the bottom boundary condition in (5b) yields

*ω*falls near the inviscid barotropic Rossby wave resonance, |

*γh*

_{0}| ≪ 1, tanh(

*γh*

_{0}) ≈

*γh*

_{0}, and (10) simply reduces to

*γ*)

*h*

_{0}≫ 1, bottom trapping occurs and tanh(

*γh*

_{0}) → 1. Thus, (10) becomes

*R*becomes large, the physically realizable (decaying) solution to (13) is then

*f*[

*N*

_{0}(

*k*

^{2}+

*l*

^{2})

^{1/2}]

^{−1}.

### e. Asymptotics: Baroclinic modes

*μ*

_{n}

^{2}) > 0 in (4b). The parameter

*η*takes on discrete values in accord with the bottom boundary condition, and so

*R*, an approximate solution is found by seeking a perturbation around the

*r*= 0 values:

*O*(

*R*) corrections, the expansions

*ω*

^{(0)}is the frictionless frequency) are used. Thus, (5b) yields

*η*

^{(0)}

*h*

_{0}) = 0,

*η*

^{2}≈

*η*

^{(0)2}+ 2

*η*

^{(0)}

*η*

^{(1)}leads to the solution

*R*(i.e., long waves or weak friction).

*R*, a similar expansion is carried out for small

*R*

^{−1}. In this case, we anticipate that interior velocity near the bottom becomes weak (see section 2f below) so that

*υ*→ 0 at

*z*= −

*h*

_{0}. Thus,

*R*.

### f. Computed results

When there is no bottom friction, Γ is real and it is straightforward to obtain the infinity of solutions for *F*_{n} by solving (6c) for constant stratification. From a practical standpoint, however, solutions of the form (6a) with *r* ≠ 0 raise computational difficulties because solving (6c) involves searches in the complex plane, such that there is sometimes a question as to whether all of the desired roots have been found. For example, it is usually easy to find the solution having the simplest vertical structure, but then some higher-mode solutions can be harder to isolate. A messier, but more algorithmically certain, approach is to expand *F*_{n} in terms of the complete set of inviscid, flat-bottom baroclinic modes and thus replace (6) with an algebraic eigenvalue problem that is straightforward to solve with readily available software.

*N*

^{2}=

*N*

_{0}

^{2}, the inviscid, flat-bottom (where Γ = 0) modes are

*J*

_{0},

*J*

_{1},

*Y*

_{0}, and

*Y*

_{1}are Bessel functions,

*a*

_{m}is a constant,

*q*

_{mB}is

*q*

_{m}evaluated at

*z*= −

*h*

_{0}, and the multiple values of

*λ*

_{m}are found by solving

*F*

_{n}anticipates the multiple solutions to the problem. Multiplication of (4a) by

*G*

_{m}, followed by vertical integration and application of (5a) and (5b) then yields an algebraic eigenvalue problem for

*ω*

_{n}:

*n*= 1, 2, 3, …, Λ), where

*c*

_{m}is the long gravity wave speed associated with inviscid mode

*m.*For example, for constant stratification,

*ω*

_{n}are saved. The first three solutions are identified as

*ω*

_{0},

*ω*

_{1}, and

*ω*

_{2}, and each is associated with a set

*b*

_{0m},

*b*

_{1m}, and

*b*

_{2m}. Using these, the frictional modal structures

*F*

_{0},

*F*

_{1}, and

*F*

_{2}are readily constructed.

A sample calculation with constant *N*^{2}, a flat bottom, and representative parameters illustrates the frictional effects. Specifically, ^{−6} s^{−2}, *h*_{0} = 4500 m, *l* = 2 × 10^{−6} m^{−1}, *r* = 1 × 10^{−4} m s^{−1}, *f* = 0.73 × 10^{−4} s^{−1}, *β* = 2 × 10^{−11} (s m)^{−1}, and *k* is varied over a wide range. The *c*_{1} = 2.3 m s^{−1} in (29a), while *f* and *β* are representative of 30° latitude. Because *r* depends on the strength of total bottom currents (including tides; e.g., Wright and Thompson 1983), it is not obvious what an appropriate deep-sea value might be, but *r* = 5 × 10^{−4} m s^{−1} is commonly used over the continental shelf (e.g., Chapman 1987), and one might expect abyssal currents in most of the ocean to be weaker than representative shelf conditions. Hence the choice of *r* = 1 × 10^{−4} m s^{−1} for these calculations.

Results of the calculation are summarized in the dispersion curves of Fig. 1. For small *k*, the right-hand side of (2b) is small (i.e., the curl of the bottom stress is small), so that frictional effects are weak, and the modal structures (Fig. 2) are essentially *G*_{n}, that is, the forms found in the complete absence of bottom friction. The smaller imaginary part of the modal structure gives rise to phase shifts as a function of the vertical. The real part of frequency *ω*_{R} is essentially the inviscid result, while the imaginary part of frequency *ω*_{I} (Fig. 1b) behaves as one would expect from a perturbation expansion for small *r* (see sections 2d and 2e). However, as *k* increases, *ω*_{I} increases for two of the modes, and then *ω*_{I} for two modes decreases for larger *k*. At the same time, *ω*_{I} for the other mode abruptly begins to increase dramatically for *k* > 0.7 × 10^{−4} m^{−1}. What is happening? Examination of the modal structures for large *k* (Fig. 3) clarifies the situation. Two of the modes adjust so that there is a node in horizontal velocity near the bottom; that is, they can be thought of as becoming modes *n* = 1/2 and *n* = 3/2. The boundary condition [(2b) or (5b)] is being met by having the bottom stress become small even though *r*(*k*^{2} + *l*^{2}) is growing. This sort of behavior, where a linear wave mode structure adjusts so that the effect of bottom friction is minimized, is not unusual in oceanographic problems (e.g., Allen 1984; Power et al. 1989; Brink 2006). On the other hand, the wave mode that has damping increasing with *k* has reached a state, for large *k*, that is strongly bottom trapped (Fig. 3). With a large bottom velocity, there is nothing to mitigate the growing *r*(*k*^{2} + *l*^{2}), and so the wave damping grows strongly as *k* increases and bottom trapping becomes more pronounced. In this case, stratification increases the wave damping because it leads to intensified (rather than weakened) near-bottom currents. It is worth pointing out that this bottom-trapping happens with a flat bottom.

Long-wave vertical modal structures for conditions as in Fig. 1 and *k* = 2.1 × 10^{−7} m^{−1}. The wave frequencies are as follows: *ω*_{0} (solid curves) = 1.04 × 10 ^{−6} +*i*2.22 × 10^{−8} s^{−1}, *ω*_{1} (dashed) = 4.15 × 10^{−9} + *i*1.78 × 10^{−10} s^{−1}, and *ω*_{2} (dash–dotted) = 1.04 × 10^{−9} + *i*4.49 × 10^{−11} s^{−1}. Blue lines indicate the real part of the modal structure, and red lines indicate the accompanying imaginary parts.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

Long-wave vertical modal structures for conditions as in Fig. 1 and *k* = 2.1 × 10^{−7} m^{−1}. The wave frequencies are as follows: *ω*_{0} (solid curves) = 1.04 × 10 ^{−6} +*i*2.22 × 10^{−8} s^{−1}, *ω*_{1} (dashed) = 4.15 × 10^{−9} + *i*1.78 × 10^{−10} s^{−1}, and *ω*_{2} (dash–dotted) = 1.04 × 10^{−9} + *i*4.49 × 10^{−11} s^{−1}. Blue lines indicate the real part of the modal structure, and red lines indicate the accompanying imaginary parts.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

Long-wave vertical modal structures for conditions as in Fig. 1 and *k* = 2.1 × 10^{−7} m^{−1}. The wave frequencies are as follows: *ω*_{0} (solid curves) = 1.04 × 10 ^{−6} +*i*2.22 × 10^{−8} s^{−1}, *ω*_{1} (dashed) = 4.15 × 10^{−9} + *i*1.78 × 10^{−10} s^{−1}, and *ω*_{2} (dash–dotted) = 1.04 × 10^{−9} + *i*4.49 × 10^{−11} s^{−1}. Blue lines indicate the real part of the modal structure, and red lines indicate the accompanying imaginary parts.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

Short-wave vertical modal structures for conditions as in Fig. 1 and *k* = 1.5 × 10^{−4} m^{−1}. The wave frequencies are: *ω*_{0} (solid curves) = 7.01 × 10^{−8} + *i*3.033 × 10^{−7} s^{−1}, *ω*_{1} (dashed) = 1.32 × 10^{−7} + *i*7.93 × 10^{−11} s^{−1}, and *ω*_{2} (dash-dot) = 1.21 × 10^{−7} + *i*5.54 × 10^{−10} s^{−1}. Blue lines indicate the real part of the modal structure, and red lines indicate the accompanying imaginary parts.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

Short-wave vertical modal structures for conditions as in Fig. 1 and *k* = 1.5 × 10^{−4} m^{−1}. The wave frequencies are: *ω*_{0} (solid curves) = 7.01 × 10^{−8} + *i*3.033 × 10^{−7} s^{−1}, *ω*_{1} (dashed) = 1.32 × 10^{−7} + *i*7.93 × 10^{−11} s^{−1}, and *ω*_{2} (dash-dot) = 1.21 × 10^{−7} + *i*5.54 × 10^{−10} s^{−1}. Blue lines indicate the real part of the modal structure, and red lines indicate the accompanying imaginary parts.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

Short-wave vertical modal structures for conditions as in Fig. 1 and *k* = 1.5 × 10^{−4} m^{−1}. The wave frequencies are: *ω*_{0} (solid curves) = 7.01 × 10^{−8} + *i*3.033 × 10^{−7} s^{−1}, *ω*_{1} (dashed) = 1.32 × 10^{−7} + *i*7.93 × 10^{−11} s^{−1}, and *ω*_{2} (dash-dot) = 1.21 × 10^{−7} + *i*5.54 × 10^{−10} s^{−1}. Blue lines indicate the real part of the modal structure, and red lines indicate the accompanying imaginary parts.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

The asymptotic expressions for the imaginary part of frequency for barotropic or bottom-trapped modes [(11) and (14)] are overplotted as blue broken lines in Fig. 1b. The weak-friction limit in (11) is seen to replicate the calculations quite well for *k* < 0.4 × 10^{−4} m^{−1}, while the strong-friction limit (the blue dashed line for *ω*_{I} > 1.5 × 10^{−7} s^{−1}) captures only the right magnitude and trend in this example.

The asymptotic results for baroclinic modes [(21) and (23)] are overplotted as green broken lines in Fig. 1b. They both provide excellent agreement with the direct calculations in the appropriate *R* range. Note that three of the expressions for wave damping rates—(11), (14), and (21)—do not depend on wave orientation, although the large *R*, higher-mode expression in (23) does depend on orientation through *ω*^{(0)}.

Up to this point, all results have been for the case of constant *N*^{2}. One might ask whether the results change substantially when one uses a more realistic, surface-intensified exponential stratification. The inviscid baroclinic modes in this case obey (25), and we choose *z*_{S} = 350 m (following SLF16) and ^{−5} s^{−2} in order to obtain a first internal mode gravity wave speed of 2.3 m s^{−1} as in the case with constant *N*^{2}. The resulting modal structures (such as the dashed line in the left panel of Fig. 5) have larger amplitudes near the surface and relatively constant values in the lower portions of the water column where *N*^{2} is small. After solving (28), Fig. 4 is obtained for the imaginary part of the frequency. We emphasize that the only difference between this and the calculations leading to the lower panel of Fig. 1 is the use of exponential stratification. Two points stand out. First, the long-wave damping is much weaker (note the changed vertical scale) with exponential stratification than in the case for constant stratification. This is not surprising in light of the relatively weak (compared to the upper water column) bottom expression of the baroclinic modes with weak deep stratification. Second, the transition between inviscid modal structures (for longer waves) and the bottom- or surface-intensified large-friction (short wave) modes occurs at roughly the same wavelengths. For example, the long-wave barotropic wave’s flat *ω*_{I} curve begins to bend at about *k* = 3 × 10^{−5} m^{−1} with constant *N*^{2} (Fig. 1) and around *k* = 4 × 10^{−5} m^{−1} with the exponential stratification (Fig. 4). Similarly, the first baroclinic mode reaches its maximum *ω*_{I} at *k* = 6.4 × 10^{−5} m^{−1} for constant *N*^{2} (Fig. 1) and at 6.2 × 10^{−5} m^{−1} for exponential stratification (Fig. 4). We conclude from this insensitivity that the vertical structure of the density stratification does not substantially affect either 1) our conclusions about the importance of bottom friction for changing modal structures or 2) the utility of the *R*_{I} parameter [(8)]. [It is evidently more appropriate to use a depth-averaged *N*^{2} than the actual near-bottom value (7).]

Imaginary part of frequency for conditions as in Fig. 1, but the stratification is exponential (^{−5} s^{−2}; *z*_{S} = 350 m). Note the change in vertical scale relative to Fig. 1b.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

Imaginary part of frequency for conditions as in Fig. 1, but the stratification is exponential (^{−5} s^{−2}; *z*_{S} = 350 m). Note the change in vertical scale relative to Fig. 1b.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

Imaginary part of frequency for conditions as in Fig. 1, but the stratification is exponential (^{−5} s^{−2}; *z*_{S} = 350 m). Note the change in vertical scale relative to Fig. 1b.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

The results presented here are fairly representative of many calculations (not detailed here), executed with both constant and exponential *N*^{2}, in that, in all cases, most wave modes evolve toward a state with weak near-bottom velocities and with decreasing wave damping as *r*(*k*^{2} + *l*^{2}) increases. However, in every calculation, there is always one wave mode that is increasingly bottom intensified and where the damping becomes large. In some calculations, the bottom-trapped wave mode is continuous with the inviscid barotropic mode as *k* varies (e.g., Fig. 4), but in some cases (such as Fig. 1), the bottom-trapped mode is continuous with one of the small-*r* baroclinic modes. Also, for a given set of parameters, the transition from the nearly inviscid modal structure toward the strongly frictional structures tends to occur at roughly the same wavenumber range (e.g., where *k* ≈ 0.7 × 10^{−4} m ^{−1} and *R*_{I} = 3.7 for Fig. 1b) for each of the three gravest modes. This transition occurs where the deformation scale [≈2*πf*/(*N*_{0}*k*) = 4100 m here] is comparable to the ocean’s 4500-m depth. Finally, even with a sloping bottom (*α*^{x} and/or *α*^{y} ≠ 0), the results do not change qualitatively: for large *r*(*k*^{2} + *l*^{2}), there is a single, strongly damped, bottom-trapped wave mode while all other modes adjust so that near-bottom velocity is small.

## 3. Forced solutions

*N*

^{2}as (4a) and (5) so that the (generally complex) free mode frequency for mode

*n*is given by

*δ*

_{nm}is the Kronecker delta (=1 for

*n = m*and =0 for

*n*≠

*m*). Bear in mind that the function Γ is complex when

*r*≠ 0, and even when

*r*= 0, it can be either positive or negative. Note that the normalization implied by this condition reduces to a more traditional form [(26)] when Γ = 0, that is, when the bottom is flat and there is no dissipation.

*F*

_{n}and (4a) by

*f*

^{2}

*P*

_{t}and subtract and integrate over depth. After applying boundary conditions [(5a) and (5b) and (31b) and (31c)], the following expression is obtained:

*P*in terms of the modes

*F*

_{m}:

A numerical example illustrates the sensitivity of the forced problem to the imposed scales. We consider a case with exponential stratification [(25a)] with ^{−5} and *z*_{S} = 350 m. In addition, *r* = 1 × 10^{−4} m s^{−1}, *f* = 7.3 × 10^{−5} s^{−1}, *β* = 2 × 10^{−11} (m s)^{−1}, and *h*_{0} = 4500 m. With scales representative of weather systems, *k = l* = 3.14 × 10^{−6} m^{−1} (Fig. 5, left), modal structures *F*_{n} are indistinguishable from the inviscid case *G*_{n}, and *b*_{0} = −1.49 × 10^{−2} + *i*5.21 × 10^{−5} m^{−1/2} and *b*_{1} = −5.41 × 10^{−2} + *i*4.32 × 10^{−5} m^{−1/2}. With scales representative of midocean eddies at 30°N (Chelton et al. (2007), results are hardly changed: for *k = l* = 2.1 × 10^{−5} m^{−1} (not shown), modal structures *F*_{n} are virtually indistinguishable from the inviscid case and *b*_{0} = −1.49 × 10^{−2} + *i*3.46 × 10^{−4} m^{−1/2} and *b*_{1} = −5.40 × 10^{−2} + *i*2.82 × 10^{−4} m^{−1/2}. These largely unperturbed modal structures argue strongly that, for the ocean mesoscale, bottom friction cannot account for weakened near-bottom velocities as found by SLF16 (a similar conclusion would be reached even for *r* = 5 × 10^{−4} m s^{−1}). For short scales, approaching the submesoscale, frictional effects become pronounced (Fig. 5, right): when *k = l* = 1.05 × 10^{−4} m^{−1}, *b*_{0} = −5.87 × 10^{−4} + *i*7.42 × 10^{−5} m^{−1/2} and *b*_{1} = −5.36 × 10^{−2} − *i*2.81 × 10^{−5} m^{−1/2}. For these short scales, the frictional modal modification is large [in accord with (8b)], and the formerly barotropic mode is bottom intensified so that *b*_{0} becomes two orders of magnitude smaller than for the larger-scale examples. On the other hand, the short-scale baroclinic modes are adjusting toward zero bottom velocity even though the surface amplitude hardly changes; hence, Real (*b*_{1}) is largely unchanged. Further, the largest damping of these three examples occurs at intermediate length scales (*k* = *l* = 2.1 × 10^{−5} m^{−1}), where *ω*_{I1} = *i*2.34 × 10^{−9} s^{−1}. When the horizontal length scale is smaller (*k* = *l* = 1.05 × 10^{−4} m^{−1}), velocities near the bottom for baroclinic modes are becoming smaller so that damping is smaller (*ω*_{I1} = *i*9.84 × 10^{−11} s^{−1}). It thus seems likely that, in the open ocean, modal properties are not strongly affected by bottom friction until very short (less than around 100 km) wavelengths are reached.

Modal structures for the two gravest modes with *r* = 1 × 10^{−4} m s^{−1}, *f* = 7.3 × 10^{−5} s^{−1}, *β* = 2 × 10^{−11} (m s)^{−1}, *h*_{0} = 4500 m, and exponential stratification (^{−5} s^{−2}; *z*_{S} = 350 m). The solid lines are the gravest mode, and the dashed lines are the next (first baroclinic) mode. Blue is the real part, and red is the imaginary part. The modes are normalized according to (32). (left) Here, *k = l* = 3.14 × 10^{−6} m^{−1}. The red lines are indistinguishable from the vertical axis. (right) Here, *k = l* = 1.05 × 10^{−4} m^{−1}.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

Modal structures for the two gravest modes with *r* = 1 × 10^{−4} m s^{−1}, *f* = 7.3 × 10^{−5} s^{−1}, *β* = 2 × 10^{−11} (m s)^{−1}, *h*_{0} = 4500 m, and exponential stratification (^{−5} s^{−2}; *z*_{S} = 350 m). The solid lines are the gravest mode, and the dashed lines are the next (first baroclinic) mode. Blue is the real part, and red is the imaginary part. The modes are normalized according to (32). (left) Here, *k = l* = 3.14 × 10^{−6} m^{−1}. The red lines are indistinguishable from the vertical axis. (right) Here, *k = l* = 1.05 × 10^{−4} m^{−1}.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

Modal structures for the two gravest modes with *r* = 1 × 10^{−4} m s^{−1}, *f* = 7.3 × 10^{−5} s^{−1}, *β* = 2 × 10^{−11} (m s)^{−1}, *h*_{0} = 4500 m, and exponential stratification (^{−5} s^{−2}; *z*_{S} = 350 m). The solid lines are the gravest mode, and the dashed lines are the next (first baroclinic) mode. Blue is the real part, and red is the imaginary part. The modes are normalized according to (32). (left) Here, *k = l* = 3.14 × 10^{−6} m^{−1}. The red lines are indistinguishable from the vertical axis. (right) Here, *k = l* = 1.05 × 10^{−4} m^{−1}.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

*r*clearly affect response amplitude.

Pressure response for wind forcing at a 48-day period. Blue curves are the real part of pressure, and red curves are the imaginary part (although the red curves are almost indistinguishable from the vertical axis in both cases). Computed for *r* = 1 × 10^{−4} m s^{−1}, *f* = 7.3 × 10^{−5} s^{−1}, *β* = 2 × 10^{−11} (m s)^{−1}, *h*_{0} = 4500 m, and exponential stratification (^{−5} s^{−2}; *z*_{S} = 350 m). Shown are (left) *k = l* = 3.14 × 10^{−6} m^{−1}, and (right) *k = l* = 1.05 × 10^{−4} m^{−1}.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

Pressure response for wind forcing at a 48-day period. Blue curves are the real part of pressure, and red curves are the imaginary part (although the red curves are almost indistinguishable from the vertical axis in both cases). Computed for *r* = 1 × 10^{−4} m s^{−1}, *f* = 7.3 × 10^{−5} s^{−1}, *β* = 2 × 10^{−11} (m s)^{−1}, *h*_{0} = 4500 m, and exponential stratification (^{−5} s^{−2}; *z*_{S} = 350 m). Shown are (left) *k = l* = 3.14 × 10^{−6} m^{−1}, and (right) *k = l* = 1.05 × 10^{−4} m^{−1}.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

Pressure response for wind forcing at a 48-day period. Blue curves are the real part of pressure, and red curves are the imaginary part (although the red curves are almost indistinguishable from the vertical axis in both cases). Computed for *r* = 1 × 10^{−4} m s^{−1}, *f* = 7.3 × 10^{−5} s^{−1}, *β* = 2 × 10^{−11} (m s)^{−1}, *h*_{0} = 4500 m, and exponential stratification (^{−5} s^{−2}; *z*_{S} = 350 m). Shown are (left) *k = l* = 3.14 × 10^{−6} m^{−1}, and (right) *k = l* = 1.05 × 10^{−4} m^{−1}.

Citation: Journal of Physical Oceanography 48, 9; 10.1175/JPO-D-18-0070.1

One might ask whether there is really any benefit to be had by expanding the solution to (31) in terms of the frictional vertical modes *F*_{n}(*z*) [(4) and (5)]. After all, one could solve the system by expanding in terms of the flat-bottom, inviscid modes *G*_{n}(*z*). Comparing solutions using the two approaches, we find that convergence requires the same number of vertical modes using either approach. The advantage of the frictional/sloping bottom modes lies in (36a), where one simply solves a sequence of uncoupled ordinary differential equations. When the solution is expanded in terms of the inviscid modes, the equivalent of (36a) includes frictional coupling among all the modes. Solving such a system of coupled differential equations for anything except the simplest functional forms of wind stress would be tedious indeed.

## 4. Discussion

One might ask to what extent the above results are relevant to actual oceanographic conditions. Specifically, can frictional effects realistically account for the tendency for many observed subinertial ocean current variations to have an *n* = 1/2 structure (SLF16)? After all, bottom friction, at some level, is ubiquitous in the ocean, and it would often not entail a sensitivity to wave propagation angle relative to bottom slope. It is, of course, a stretch to apply the present linear results to ocean currents that are known to be nonlinear, but the attempt nonetheless is made. The results in the previous section argue strongly that, for the ocean mesoscale, bottom friction alone cannot account for the SLF16 modal modification, even if a presumably too-large friction coefficient of 5 × 10^{−4} m s^{−1} is applied. For substantially shorter wavelengths (i.e., the ocean submesoscale), the modal structures do adjust to have near-zero bottom velocities for the baroclinic modes, as expected (Fig. 5). It thus seems that bottom friction by itself is unlikely to distort wave modal structures in the open ocean for the most energetic (mesoscale) eddy length scales. However, this conclusion should be treated as tentative, since it is based entirely on linear physics, while it is likely that nonlinearity plays a substantial role in the actual ocean (e.g., Chelton et al. 2007). Certainly, it has been demonstrated that linearizing about the mean vorticity field associated with a steady flow leads to a modification of modal structures that can either enhance or weaken near-bottom velocities (Killworth et al. 1997). Further, other studies (e.g., Trossman et al. 2017; LaCasce 2017) show that the presence of a sloping or irregular bottom, or of bottom friction, can lead to surface-intensified flow. In contrast to midoceanic conditions, frictional effects in shallower water are relatively more important [as measured by *r*/(*ωh*_{0})], and stratification is generally stronger in the upper few hundred meters than at abyssal depths. Our analysis suggests that, for water depths shallower than a few hundred meters, bottom friction is clearly expected to affect both linear wave (e.g., Brink 2006) and nonlinear eddy (e.g., Brink 2017) structures.

Perhaps the more unexpected result here is the development of the strongly bottom-trapped, strongly damped mode at short horizontal length scales. A bottom-intensified inviscid mode is of course also found in the presence of a sloping bottom (Rhines 1970), but bottom trapping with a flat bottom was not anticipated. Thus, the near-bottom surface quasigeostrophic (SQG; e.g., Held et al. 1995) phenomenology has been broadened. When the mode is strongly bottom trapped, it is unlikely to be stimulated by a surface wind stress (see section 3), but it could still be readily excited in an initial value problem, in cases with a body forcing, or in problems involving flow over varying bottom topography. Since the bottom-trapped mode is strongly damped, it is most likely to be found close to where it is excited.

## Acknowledgments

We appreciate thoughtful comments from Brian Arbic and Joe LaCasce. Partial funding for this article is provided by the National Science Foundation Physical Oceanography section through Award OCE-1433953.

## REFERENCES

Allen, J. S., 1984: A simple model for stratified shelf flow fields with bottom friction.

,*J. Phys. Oceanogr.***14**, 1200–1214, https://doi.org/10.1175/1520-0485(1984)014<1200:ASMFSS>2.0.CO;2.Arbic, B. K., and G. R. Flierl, 2004: Baroclinically unstable geostrophic turbulence in the limits of strong and weak bottom Ekman friction: Application to midocean eddies.

,*J. Phys. Oceanogr.***34**, 2257–2273, https://doi.org/10.1175/1520-0485(2004)034<2257:BUGTIT>2.0.CO;2.Brink, K. H., 1989: Evidence for wind-driven current fluctuations in the western North Atlantic.

,*J. Geophys. Res.***94**, 2029–2044, https://doi.org/10.1029/JC094iC02p02029.Brink, K. H., 2006: Coastal-trapped waves with finite bottom friction.

,*Dyn. Atmos. Oceans***41**, 172–190, https://doi.org/10.1016/j.dynatmoce.2006.05.001.Brink, K. H., 2017: Surface cooling, winds, and eddies over the continental shelf.

,*J. Phys. Oceanogr.***47**, 879–894, https://doi.org/10.1175/JPO-D-16-0196.1.Chapman, D. C., 1987: Application of coastal-trapped wave theory along the California coast.

,*J. Geophys. Res.***92**, 1798–1816, https://doi.org/10.1029/JC092iC02p01798.Chelton, D. B., M. G. Schlax, R. M. Samelson, and R. A. de Szoeke, 2007: Global observations of large oceanic eddies.

,*Geophys. Res. Lett.***34**, L15606, https://doi.org/10.1029/2007GL030812.Held, I. M., R. T. Pierrehumbert, S. T. Garner, and K. L. Swanson, 1995: Surface quasigeostrophic dynamics.

,*J. Fluid Mech.***282**, 1–20, https://doi.org/10.1017/S0022112095000012.Killworth, P. D., D. B. Chelton, and R. A. de Szoeke, 1997: The speed of observed and theoretical long extratropical planetary waves.

,*J. Phys. Oceanogr.***27**, 1946–1966, https://doi.org/10.1175/1520-0485(1997)027<1946:TSOOAT>2.0.CO;2.LaCasce, J. H., 2017: The prevalence of oceanic surface modes.

,*Geophys. Res. Lett.***44**, 11 097–11 105, https://doi.org/10.1002/2017GL075430.LaCasce, J. H., and K. H. Brink, 2000: Geostrophic turbulence over a slope.

,*J. Phys. Oceanogr.***30**, 1305–1324, https://doi.org/10.1175/1520-0485(2000)030<1305:GTOAS>2.0.CO;2.Power, S. B., J. H. Middleton, and R. H. J. Grimshaw, 1989: Frictionally modified continental shelf waves and the subinertial response to wind and deep-sea forcing.

,*J. Phys. Oceanogr.***19**, 1486–1506, https://doi.org/10.1175/1520-0485(1989)019<1486:FMCSWA>2.0.CO;2.Rhines, P. B., 1970: Edge-, bottom-, and Rossby waves in a rotating stratified fluid.

,*Geophys. Fluid Dyn.***1**, 273–302, https://doi.org/10.1080/03091927009365776.Samelson, R. M., 1992: Surface-intensified Rossby waves over rough topography.

,*J. Mar. Res.***50**, 367–384, https://doi.org/10.1357/002224092784797593.Sanchez de la Lama, M., J. H. LaCasce, and H. K. Fuhr, 2016: The vertical structure of ocean eddies.

,*Dyn. Stat. Climate Syst.***1**, dzw001, https://doi.org/10.1093/climsys/dzw001.Trossman, D. S., B. K. Arbic, D. N. Straub, J. G. Richman, E. P. Chassignet, A. J. Wallcraft, and X. Xu, 2017: The role of rough topography in mediating impacts of bottom drag in eddying ocean circulation models.

,*J. Phys. Oceanogr.***47**, 1941–1959, https://doi.org/10.1175/JPO-D-16-0229.1.Wright, D. G., and K. R. Thompson, 1983: Time-averaged forms of the nonlinear stress law.

,*J. Phys. Oceanogr.***13**, 341–345, https://doi.org/10.1175/1520-0485(1983)013<0341:TAFOTN>2.0.CO;2.Wunsch, C., 1997: The vertical partition of oceanic horizontal kinetic energy.

,*J. Phys. Oceanogr.***27**, 1770–1810, https://doi.org/10.1175/1520-0485(1997)027<1770:TVPOOH>2.0.CO;2.