Codex @codex  0

For a nonzero-frequency Fourier mode proportional to $e^{i(kx+mz-\omega t)}$, the time-independent source has no oscillatory part. The pressure equation givesTake $k,m\ne 0$ and ${\omega }^2\ne f^2,N^2$ for the following pressure-amplitude formulas. With $P_0$ the amplitude of $p/{\rho }_0$, the momentum and buoyancy equations implyConsequentlybecause the numerator is $(k^2+m^2){\omega }^2-N^2{k}^2-f^2{m}^2$. Propagating linear inertia-gravity waves carry zero perturbation potential vorticity. In degenerate directions the division formulas must be replaced by the original linear equations, or interpreted by a regular limit; the zero-PV property of the nonzero-frequency wave component still follows from $q_t=0$. A steady balanced mode need not have zero potential vorticity.

Put $P=p/{\rho }_0$ and use the downward buoyancy anomaly $b=g{\rho }^{\prime }/{\rho }_0$. This sign convention is required by the printed potential vorticity. The linear rotating Boussinesq approximation givesLet $D=u_x+v_y$ and $\zeta =v_x-u_y$. Horizontal divergence and curl give $D_t-f\zeta =-{\nabla }_H^{2}P$ and ${\zeta }_t=-fD=f{w}_z$. HenceThus the linear potential vorticity is a time-independent field fixed by the initial data. Differentiate the pressure identity twice. Using ${\zeta }_{tt}=-f(P_{zz}+b_z)$ and $b_{ztt}=-N^2(P_{zz}+b_z)$ gives ${\nabla }^2{P}_{tt}=(N^2-f^2)(P_{zz}+b_z)$. Combining terms,The linear pressure equation for rotating stratified flow has source $Q={\rho }_{0}f{N}^{2}q$. “Arbitrary” means that different initial potential-vorticity fields give different time-independent sources; it is not independent forcing that may vary in time. If buoyancy were instead defined upward as $-g{\rho }^{\prime }/{\rho }_0$, both appearances of its sign would change consistently. The pressure source represents the balanced component alongside propagating inertia-gravity waves.

Use the quasi-geostrophic streamfunction convention $u_n=-{\psi }_{n,y}$, $v_n={\psi }_{n,x}$, and write $F=k_I^2$, $f=f_0+\beta y$. In a steady, large-scale basin interior, neglect the material derivatives of relative vorticity and interfacial stretching compared with advection of planetary vorticity. With small Rossby number, weak nonlinear eddy terms, no significant interior friction or topographic forcing, and the specified forcing confined to layer 1, Sverdrup balance isFor equal depths $h$, the depth-integrated meridional transport is $h({\bar{v}}_1+{\bar{v}}_2)=hW/\beta$. Here $W$ is the normalized potential vorticity source appearing in the evolution equation. If the dimensional wind stress curl is used, its usual layer forcing is $W=({\nabla }_h\times \mathbf{\tau }{)}_z/({\rho }_{0}h)$, so $\beta h({\bar{v}}_1+{\bar{v}}_2)=({\nabla }_h\times \mathbf{\tau }{)}_z/{\rho }_0$.

For $\beta >0$ and $W<0$, the upper-layer interior transport is southward. In a closed subtropical basin, negative wind stress curl also corresponds to downwelling Ekman pumping for $f_0>0$ and an anticyclonic gyre. A northward return transport is needed to close the circulation; its narrow western boundary current requires processes outside the frictionless Sverdrup balance. The local interior equations by themselves do not specify the detailed boundary-current structure.

Expanding the layer equation shows the physical budget:The first term is the change of planetary vorticity as a fluid parcel moves north or south. The second is the change of relative vorticity. The last is vortex stretching in layered quasi-geostrophic flow: displacement of the interface changes layer thickness and thus the stretching contribution to potential vorticity. Its opposite signs in the two layer definitions express their thickness changes in opposite directions. The wind stress curl supplies or removes upper-layer potential vorticity. Each $D_n/Dt$ follows its own layer velocity, rather than a common velocity for both layers.

For uniform $W$, choose the local two-layer Sverdrup interiorUnforced background zonal currents and arbitrary additive interface offsets have been set to zero. This choice is also an exact uniform-flow solution of the stated forced equations, not just a leading balance: ${\bar{q}}_1=f_0+\beta y-FVx$, ${\bar{q}}_2=f_0+\beta y+FVx$, and ${\bar{v}}_1{\bar{q}}_{1,y}=V\beta =W$. The local interface slope can be nonzero even though its stretching contribution is constant along each basic-state trajectory. A streamfunction linear in $x$ is an interior-patch description, not a complete globally bounded basin solution.

Put $S=M^2/f_0$. The thermal wind is $\bar{U}(z)=-Sz$, and the basic potential vorticity is spatially constant: $\bar{q}=2f_0$ when the full reference buoyancy is included, or $f_0$ when that reference contribution has been subtracted. Either way ${\nabla }_h\bar{q}=0$.

For a smooth normal mode, the linear quasi-geostrophic potential-vorticity equation becomesTake $k\ne 0$ and $K>0$, as required for a nontrivial wave with the printed phase-speed parameterization. Except possibly at an isolated critical level of a shear-flow wave, this implies $\hat{q}=0$; regularity then extends that condition through the isolated level. HenceDecay at infinity excludes the growing exponential, giving the decaying vertical structure of a semi-infinite Eady edge waveThe amplitude $C$ is arbitrary. The semi-infinite Eady model supports a wave trapped at its lower boundary with penetration depth $m^{-1}$. The zero-horizontal-wavenumber case has no nonzero decaying solution of this homogeneous vertical equation. Singular neutral interior potential vorticity sheets belong to a different continuous-spectrum class; they are not the smooth decaying edge-wave eigenfunction requested here.

The PDF contains the request for this vertical structure and uses the phase factor $e^{i(kx+ly-kct)}$. The TeX aid omits that request and also inserts a duplicate factor without $c$. The solution uses the original PDF form throughout.

For an exponentially growing normal mode with $k\ne 0$, $\mathrm{Im}(kc)>0$. Such a complex $c$ cannot equal the real basic velocity anywhere, so the interior equation forces $\hat{q}=0$ at every height. Decay then fixes the unique vertical structure $C{e}^{-mz}$, and the lower-boundary condition fixes $c=-S/m$, which is real. This contradiction proves the absence of exponential instability in the semi-infinite Eady model:For $k=0$, there is no advective frequency in the linear interior equation, so an exponentially growing mode would again have $q^{\prime }=0$. The wall buoyancy equation then forces ${\psi }_z^{\prime }(0)=0$, which the nonzero decaying exponential cannot satisfy. Thus these modes do not provide a missing growing branch. This conclusion rules out exponential eigenmode growth, not every possible transient response or singular neutral potential vorticity disturbance.

With a second rigid boundary at $z=D$, the same uniform-shear flow is the Eady model. Both boundaries carry buoyancy perturbations and support Boundary Rossby waves. Their intrinsic propagation directions oppose one another relative to their respective local basic flows. At suitable wavelengths their interacting fields can phase-lock and release basic-state potential energy, producing baroclinic instability. The single-boundary system lacks that second interacting boundary wave.

For a direct comparison, write $\Lambda ={\bar{U}}_z=-S$ and $m=NK/|f_0|$. The vertical solution between the two boundaries is $A\cosh (mz)+B\sinh (mz)$, and their buoyancy conditions areEliminating $A,B$ gives the finite-depth Eady dispersion relationThe bracket is negative for $0<mD<2.399\dots$, so those modes with $k\ne 0$ have a growing and a decaying member. At short wavelengths the two edge waves interact weakly and the frequencies are real. In the semi-infinite limit, the lower-wave root tends to $c=\Lambda /m=-S/m$, recovering the neutral result above.

Let $Q_b=f/H$. At sufficiently small Rossby number, retain advection of this full background shallow-water potential vorticity while dropping the higher-order advection of $\zeta /H$. For uniform-depth Sverdrup balance, this also requires relative-vorticity advection to be small against planetary-vorticity advection; small Rossby number alone does not guarantee this when $\beta$ is very small. Multiplying the resulting steady potential vorticity equation by $H$ givesDefine the topographic potential-vorticity pseudovelocityIts sign is important: it is opposite to the coefficient vector on the left of the preceding equation. Substituting the ocean transport streamfunction expression for relative vorticity yieldsThe pseudovelocity follows contours of $f/H$ and is not the fluid velocity; its units are inverse length squared per time because $\psi$ is a transport streamfunction. If $r=F=0$, the ocean transport streamfunction is constant along each connected background potential vorticity contour. Locally, at regular points,with arbitrary differentiable $\Phi$. Different disconnected components of a level set may carry different functions until boundary conditions identify them. If $f/H$ is constant over an open region, the leading equation gives no restriction there; critical points similarly require regularity and global matching rather than division by a vanishing gradient.

Write the parabolic depth as $H(x)=ax(L_x-x)$ with $a>0$, and use $f=f_0+\beta y$, $\beta >0$. The topographic potential-vorticity pseudovelocity isIt is westward everywhere in the interior. In the Northern Hemisphere it is also southward on the western slope and northward on the eastern slope. Its characteristics are the background potential vorticity contourswhich bend towards the south on approaching either shallow side. Along these characteristics forcing accumulates into a transport response and drag spreads it across neighboring characteristics. The actual current follows contours of $\psi$, not contours of this pseudovelocity in a forced region.

The sketch uses an explicit illustrative positive drag, a single connected impermeable boundary with $\psi =0$, $H=4x(1-x)$, and $f=1+0.25y$ in normalized coordinates, with small forcing amplitude $F_0=10^{-3}$; linearity makes its value affect amplitudes rather than contour shapes. Solid and dashed contours distinguish the two signs of $\psi$; arrows indicate the pseudovelocity, not the real current. The positive lower-half wind curl tends to drive cyclonic circulation and the negative upper-half wind curl anticyclonic circulation; topographic steering bends the gyres, so their dividing contour need not coincide with the forcing's zero line. The plot is an example, not a uniquely specified circulation: the question supplies neither a drag magnitude nor full boundary data.

The westward propagation of transport information explains the connection with western boundary currents: a broad wind-driven interior generally needs a narrow western return region to satisfy the impermeable boundary condition. For constant depth the characteristic direction is purely westward and the Stommel boundary layer width is $r/\beta$ for this drag normalization. Here strong depth gradients also steer the return along slopes and can make topographic boundary currents important; a conventional flat-bottom western-current profile does not follow unchanged. The equation becomes singular where $H=0$. The drawn shore limits impose zero transport; a physically resolved shoreline would need a positive-depth cutoff or a separate near-shore model, and the small-Rossby number approximation need not remain uniform there.

Write $\zeta =(\nabla \times \mathbf{u})\cdot {\mathbf{e}}_z$ and $q=\zeta -f\eta /h_0$, so the PDF's vector potential vorticity is $\mathbf{Q}=q{\mathbf{e}}_z$. The local TeX loses some of this boldface. Taking the vertical curl of the momentum equation gives ${\zeta }_t=-f\nabla \cdot \mathbf{u}$, while mass conservation gives ${\eta }_t=-h_0\nabla \cdot \mathbf{u}$. ThusThis is the linear perturbation of shallow-water potential vorticity, scaled by $h_0$; it should not be confused with the full nonlinear conserved ratio $(f+\zeta )/h$.

Since $\nabla \cdot ({\mathbf{e}}_z\times \mathbf{u})=-\zeta$, taking the divergence of momentum gives $(\nabla \cdot \mathbf{u}{)}_t=f\zeta -g{\nabla }^2\eta$. Differentiate mass conservation and use $\zeta =q+f\eta /h_0$ to obtainFor $\mathbf{Q}=0$, substitution of the stated plane wave gives the linear rotating shallow-water dispersion relation, with $c=\sqrt{g{h}_0}$:The two even branches approach $\pm c|k|$ at large $|k|$ and have values $\pm |f|$ at $k=0$. The original dispersion diagram below also compares the speed magnitudes.

For $k\ne 0$, the phase velocity and group velocity magnitudes areThe printed convention $e^{i(kx+\omega t)}$ means the signed phase velocity is $-\omega /k$ and the signed group velocity is $-d\omega /dk$. The speed magnitudes above are independent of that sign convention.

For $f\ne 0$, longer wavelengths have greater phase velocity but smaller group velocity: individual crests and wave packets therefore answer “faster or slower” differently. At long wavelength, $|k|\ll |f|/c$, the frequency is nearly the inertial frequency $|f|$, $|c_p|\sim |f|/|k|$, and $|c_g|\sim c^2|k|/|f|$. Rotation has its largest relative effect here, and the waves are strongly dispersive. At short wavelength, $|k|\gg |f|/c$, both speeds tend to $c$ and rotation gives only a small, weakly dispersive correction. If $f=0$, the gravity waves are nondispersive, with speed $c$; at $k=0$ the phase velocity formula is undefined.

Adopt the Boussinesq approximation, hydrostatic approximation, traditional approximation, and inviscid quasi-geostrophic approximation. Assume stable constant $N_0^2>0$, nonzero Coriolis parameter $f_0$, a beta plane, small Rossby number, $\beta L/|f_0|$ of that small order, and small displacement of the background stratification. Vertical advection of the buoyancy perturbation is then higher order. Define the quasi-geostrophic streamfunction byThe buoyancy and vertical vorticity equations at the retained order areDifferentiate the first in $z$. Since $J({\psi }_z,{\psi }_z)=0$, its material-derivative term becomes $f_0{D}_g{\psi }_{zz}$. Elimination of $w_z$ gives diabatically forced quasi-geostrophic potential vorticity:Thus it is the height gradient of heating, rather than heating alone, that produces potential vorticity at this order.

For the steady response linearized about rest, assume $\beta \ne 0$. The quasi-geostrophic potential-vorticity equation reduces to $\beta {\psi }_x=f_0{R}_z/N_0^2$. Choosing the horizontally uniform pressure gauge so that $\psi \to 0$ in the east gives the steady quasi-geostrophic response to localized heatingSetUsing the error function to integrate the Gaussian function givesThese velocities vanish as $x\to +\infty$. A depth-dependent but horizontally uniform addition to $\psi$ is dynamically irrelevant for the requested horizontal velocity. At $z=0$, streamlines are contours of $e^{-y^2/L^2}I(x)$. For $f_0,\beta ,R_0>0$, flow approaches the heating region from the west on its northern side, turns southward, and returns westward on its southern side. The contours are open, not closed gyres: $I$ tends to a nonzero constant in the west.

For $-H<z<0$, $R_z<0$ and the forcing decreases potential vorticity; steady balance requires $\beta v_g=f_0{R}_z/N_0^2<0$, giving southward flow. Below $z=-H$, $R_z>0$ and the meridional flow reverses. A constant heating rate independent of $z$ would have no such interior potential vorticity source.

For the nonlinear diagnostic quasi-geostrophic omega equation, write $\zeta ={\nabla }_h^2\psi$. Apply ${\nabla }_h^2$ to buoyancy evolution and $f_0{\partial }_z$ to vorticity evolution:Subtracting eliminates the pressure tendency and yieldsThe beta term is linear and must remain when the nonlinear terms are neglected. For the steady weak-forcing solution, $f_0\beta {\psi }_{xz}=f_0^2{R}_{zz}/N_0^2$, soThis also follows immediately from the steady linear buoyancy equation. It satisfies $w=0$ at $z=0$, decay at depth and horizontal infinity, and the linear quasi-geostrophic omega equation; with these homogeneous conditions the elliptic problem has no extra decaying homogeneous solution. The rigid-lid interpretation is an additional boundary approximation for this diagnostic response.

At $x=0$, the displayed $v_g$ and $w$ provide the requested meridional/vertical arrows. For positive heating there is upwelling concentrated near $y=0$, maximal in depth at $z=-H$, with northward flow below that depth and southward flow above. The vertical circulation tends to zero far from the forcing. A western zonal wake persists, but its meridional and vertical components vanish as $x\to -\infty$. Horizontal ageostrophic flow supplies the divergence associated with the upwelling; the projected $(v_g,w)$ slice is not itself a two-dimensional incompressible flow and need not have closed streamlines.

The steady calculation fixes an interior circulation, not its attainability from every possible initial and boundary state. For example, a rigid lid with initially uniform boundary buoyancy conserves that buoyancy at linear order because $R(0)=w(0)=0$. The formal steady profile instead generally has nonuniform $b^{\prime }(0)=f_0{\psi }_z(0)$. Such additional initial boundary data require a time-dependent adjustment and cannot simply be imposed on this particular steady solution. No boundary buoyancy data are supplied for the requested steady problem.

In an adiabatic quasi-geostrophic approximation, hydrostatic approximation and geostrophic balance give $b^{\prime }=f_0{\psi }_z$. The leading buoyancy equation is $D_g{b}^{\prime }+N^{2}w=0$; hence $w=-D_g(f_0{\psi }_z/N^2)$. At the rigid surface $z=\alpha y$, the no-normal-flow condition is $D_g(z-\alpha y)=0$, so $w=\alpha v_g=\alpha {\psi }_x$. The small boundary slope allows the disturbance condition to be evaluated at the flattened boundary $z=0$ to the retained order.

For thermal wind $U=\Lambda z$, choose the basic quasi-geostrophic streamfunction $\overline{\psi }=-\Lambda zy$. Its buoyancy anomaly is $\overline{b}=-f_0\Lambda y$, and its interior quasi-geostrophic potential vorticity is constant. Therefore the linear disturbance equations areThe term $-\Lambda {\psi }_x^{\prime }$ is the advection of the basic meridional buoyancy gradient. Since $q^{\prime }$ is simply advected at each fixed height, $q^{\prime }=0$ initially implies $q^{\prime }=0$ forever. At the lower boundary define $a=\alpha N^2/f_0$; thenThe PDF uses $\alpha$ for the lower-boundary slope throughout; the TeX's isolated $z=cy$ is a transcription error.

For nonzero zonal wavenumber let $\mu =N|k|/|f_0|>0$. The zero-interior-potential vorticity condition gives ${\hat{\psi }}_{zz}-{\mu }^2\hat{\psi }=0$. In the semi-infinite domain choose decay as $z\to \infty$. Writing $B_0(t)={\hat{\psi }}_z(0,t)$, potential-vorticity inversion then givesThis sloping-boundary Eady edge wave has all its time dependence in the boundary buoyancy equation; the interior responds instantaneously through elliptic potential-vorticity inversion. A single boundary wave oscillates without exponential growth. Vertical shear and topography contribute oppositely to its propagation. The contributions cancel when $\alpha =f_0\Lambda /N^2$, the slope of a background isopycnal. For $f_0>0$, positive $\Lambda$ without slope gives eastward propagation, while a positive slope without shear gives westward propagation. The phase velocity varies as $1/|k|$, and the trapping depth is $1/\mu$. Its frequency is independent of $|k|$ on each signed branch, so its zonal group velocity is zero in this ideal semi-infinite model.

With a rigid horizontal upper boundary, put $B_D(t)={\hat{\psi }}_z(D,t)$ and $M=\mu D$. Solving the same second-order inversion problem with both derivative data givesDifferentiation at $0$ and $D$ checks the prescribed values. At the upper boundary $w^{\prime }=0$, so $({\partial }_t+ik\Lambda D)B_D=ik\Lambda \hat{\psi }(D)$. Combining this with the lower-boundary condition givesThus two Boundary Rossby waves interact across the layer. For normal modes proportional to $e^{-ikct}$ their dimensional wave-speed matrix isFor $\Lambda \ne 0$, define $C=c/(\Lambda D)$ and $\tilde{\alpha }=a/\Lambda$. Its characteristic polynomial is the Eady model with a sloping lower boundary dispersion relationThe source writes $\tilde{\mu }=NkD/f_0$; for $k,f_0>0$ this is $M$, and the dispersion relation is unchanged by replacing that signed quantity by its absolute value. Its two roots areIf $\tilde{\alpha }>1$, the second term is positive, so $\Delta >0$ for every $M>0$: all modes are spectrally stable. If $\tilde{\alpha }<1$, the increasing function $M\tanh M$ runs from zero to infinity, so there is a unique $M$ withThere the squared term is zero and $\Delta <0$, giving an exponentially growing member of the complex conjugate pair. By continuity it lies in an unstable band. This proves both the requested instability for $0\le \tilde{\alpha }<1$ and instability for negative slopes in the sloping-boundary Eady model. It does not assert that every wavelength is unstable. Strongly negative slopes move the resonant band to large $M$, where the coupling is exponentially weak.

At $\tilde{\alpha }=1$ the speeds are real, but their crossing at $M=\coth M$ can give a defective neutral mode and algebraic growth; absence of exponential instability is weaker than boundedness of every initial disturbance. The nondimensional slope is undefined when $\Lambda =0$; the dimensional matrix remains valid and then has two real speeds. The instability for $\tilde{\alpha }<1$ is the counterpropagating wave instability mechanism: lower and upper boundary waves can match their laboratory phase speeds and exchange energy with the vertical shear.

The first-order perturbation of the vortensity isTaking the curl of the linearized momentum equation and using the linearized continuity equation yieldsThis is linearized vortensity conservation. It is the perturbative form of material conservation of potential vorticity in an inviscid barotropic fluid; the forcing $-\nabla \Psi$ is a gradient and therefore creates no vorticity.

Put $D(x,y)=H(x,y)-h(x,y)$ and let $\mathbf{u}=(u,v)$ be the depth-independent interior velocity. The kinematic boundary conditions on the sloping upper and lower surfaces areIntegrating mass conservation through the layer givesThe inviscid vertical-vorticity equation isConsequently the shallow-water potential vorticityobeys the forced evolution equationPositive upward Ekman pumping removes layer thickness and raises the potential vorticity of the remaining column.

For a steady small-Rossby number flow, $q\simeq f/D$, soEquivalently,The same result follows from the integrated vortex-stretching balance

When $\nabla D=(c,0)$ with $c>0$,If the upper pumping is absent or weak and the upper surface is locally level, the impermeable-bottom condition is $w_b=-cu$. The stipulated $w_b<0$ then gives $u>0$, and in the Northern HemisphereThe steady interior flow is therefore directed northeastward, along contours of $f/D$ in the unforced limit. As a parcel moves eastward into deeper water, its vortex column stretches; a poleward displacement increases $f$ and preserves potential vorticity. Nonzero Ekman pumping drives motion across the $f/D$ contours. This is topographic potential-vorticity steering and the associated topographic Sverdrup balance.