Codex @codex  0

Absolute vorticity is the vorticity measured in an inertial frame. In a frame rotating with angular velocity $\mathbf{\Omega }$, it is $\mathbf{\omega }+2\mathbf{\Omega }$.

Ertel potential vorticity is the scalar product of absolute vorticity with the gradient of a materially conserved stratifying scalar. In Boussinesq flow with buoyancy $b$,It is materially conserved by inviscid, adiabatic motion.

An inviscid fluid with conservative body force, a materially conserved stratifying scalar $\alpha$, and an equation of state making $\nabla \rho ,\nabla p,\nabla \alpha$ coplanar conserves its Ertel potential vorticity. The stretching of absolute vorticity per unit density cancels the deformation of the scalar gradient; the remaining baroclinic triple product vanishes. In an incompressible density-stratified ideal fluid, choosing the advected density as $\alpha$ gives that cancellation directly.

For dry atmospheric Ertel potential vorticity, use absolute vorticity ${\mathbf{\omega }}_a=\nabla \times \mathbf{u}+2\mathbf{\Omega }$ and potential temperature $\theta =T(p_0/p{)}^{R_d/c_p}$:The conservation statement is Ertel's theorem for smooth, inviscid, adiabatic flow with conservative body forces and a thermodynamic equation of state $p=p(\rho ,\theta )$. Stretching of ${\mathbf{\omega }}_a/\rho$ cancels deformation of $\nabla \theta$ in their material derivative; the remaining baroclinic triple product vanishes because $\nabla p$ lies in the span of $\nabla \rho ,\nabla \theta$. The atmospheric form includes $1/\rho$; the commonly used Boussinesq approximation version has a different normalization.

In a predominantly vertical stratification, $q_E\simeq (f+\zeta ){\partial }_z\theta /\rho$. Thus rotation and stratification jointly label air masses. Material conservation explains how parcels deform and stretch while preserving a useful scalar, supports Rossby waves, and helps diagnose sharp fronts, jets and stratospheric intrusions.

Under an appropriate balance model, potential-vorticity inversion recovers associated wind and temperature fields from a PV distribution together with stratification and boundary data. PV alone does not specify every component of an arbitrary atmospheric flow. Potential-vorticity generation by heating and friction limits conservative parcel reasoning near strong heating, surface stress or mixing.^[ref]

Use eastward $x$, northward $y$, constant layer depth $H$, constant reference density ${\rho }_0$, and a depth-independent horizontal velocity field $(u,v)$. Assume a flat bottom and a rigid upper boundary, giving the two-dimensional incompressibility condition $u_x+v_y=0$. The traditional approximation retains only the vertical component of planetary rotation in the horizontal equations. Take the Coriolis parameter $f=f(y)$, a constant linear drag rate $r$, and no other torques. The component Cauchy momentum equations areDefine relative vorticity $\xi =v_x-u_y$. Apply ${\partial }_x$ to the second equation and subtract ${\partial }_y$ of the first. The pressure gradient force has zero curl because ${\rho }_0$ is constant and mixed partial derivatives commute. Direct expansion of the material-acceleration terms givesThe curl of the rotation terms isThe imposed zonal wind stress contributes $-{\tau }_{s,y}/({\rho }_{0}H)$, and constant linear drag contributes $-r\xi$. HenceSince $Df=v{f}_y$, defining absolute vorticity $\zeta =\xi +f$ givesWithout the incompressibility condition, the left-hand side would contain the stretching term $\zeta {\nabla }_h\cdot \mathbf{u}$. Variable layer depth, variable density or spatially varying linear drag would introduce further terms. No small-Rossby number assumption was needed to obtain this constant-depth vorticity equation.

For a zonal background $U(y)$, its relative vorticity is $\overline{\xi }=-U_y$. The meridional absolute vorticity gradient is thereforeAt $30^{\circ }\mathrm{N}$,A localized jet stream of velocity scale $U=134\mathrm{m}{\mathrm{s}}^{-1}$ and width $L$ hasFor $L=1000$--$2000\mathrm{km}$ this is about $1.7$--$6.8$ times the planetary contribution, so wind curvature cannot generally be neglected. The coefficient in this scaling depends on the chosen jet profile and width convention.

Linearizing the vorticity equation now gives the Rossby-wave equation for a sheared zonal currentFor a qualitative estimate at the jet core, freeze the local coefficients. At a smooth maximum $U_{yy}<0$, so $Q_y>\beta$: the restoring absolute vorticity gradient is stronger. The local Rossby-wave dispersion relation in a zonal jet gives, for the same local $U$ and zonal wavenumber $k$ with $l=0$,soThe envelope in the diagram would advance eastward faster, while individual phase contours would advance more slowly eastward, become stationary if $Q_y=U{k}^2$, or move westward if $Q_y>U{k}^2$. The separation between envelope and phase speeds increases. On a plot with time downward, compare horizontal displacement at fixed time: the envelope moves farther right, while a given phase moves less far right, or left after reversal. If the frequency changes sign for a fixed signed $k$, conjugate the mode to restore the convention $\omega \ge 0$; the physical propagation is unchanged.

This conclusion concerns the concave-down jet core. On a flank with $U_{yy}>0$ the correction has the opposite sign; the value at the edge of the Hadley cell alone does not specify the curvature everywhere. Indeed the angular-momentum-conserving profile in 1(a) is a rising inner branch, not a complete smooth jet maximum. A global sheared-current normal mode requires its meridional structure and boundary conditions; the local calculation is the qualitative estimate requested.

Choose a fixed vertical datum so the physical free surface is at $H+\zeta$, with $H$ constant; its depth above the bottom is $h=H+\zeta -b$. The gradient of the constant datum does not enter the pressure force. With $D_H={\partial }_t+u{\partial }_x+v{\partial }_y$, the rotating shallow water equations areThe fluid is hydrostatic with depth-independent horizontal velocity at this order. The pressure-gradient force uses the free-surface slope, not the layer-depth slope: $\nabla \zeta =\nabla h+\nabla b$.

Let $q=v_x-u_y$ be relative vertical vorticity, $q_a=f+q$ the absolute vorticity, and $\delta =u_x+v_y$. Taking the horizontal curl of the two momentum equations and expanding the advective derivatives givesMass conservation is $D_{H}h=-h\delta$. The two stretching factors cancel, proving the exact shallow-water potential vorticity conservation lawThe denominator is $h$, the actual column thickness that changes under horizontal convergence; an arbitrary surface elevation $\zeta$ does not obey that column mass equation. There is no $w{\partial }_z$ term because the reduced layer quantities have no vertical dependence. The vertical motion has already supplied the layer-thickness equation, rather than remaining an independent coordinate of this model.

Now take $H=h_{00}$ and assume $|q|\ll |f|$, $|\zeta |,|b|\ll H$. To first order,For a small Rossby number, leading geostrophic flow satisfies $fv=g{\zeta }_x$ and $fu=-g{\zeta }_y$. ThusFor $f\ne 0$, the barotropic deformation radius is $L_R=\sqrt{gH}/|f|$. Multiplication of the first-order PV by $H$ givesThis is the required shallow-water quasi-geostrophic potential vorticity, including the bottom contribution.

Write $S=fϵ/H$ and $L=L_R$. The initial state has $\psi =0$, so ${\mathcal{Q}}_i=f+(f/H)b$. The specified homogenization sets ${\mathcal{Q}}_f=f$ in $|y|<a$ and leaves it unchanged outside. Hence the requested change isFor $f>0$ this is a descending straight segment from $+Sa$ to $-Sa$, with zero exterior values and jumps at the two edges. Endpoint values do not affect the inversion. It is the anomaly, not the full final PV profile, that appears as the source below.

For the zonal mean, $v=0$, $u=-{\psi }^{\prime }(y)$, and potential-vorticity inversion requireswith decaying disturbance at both infinities. The odd source makes the unique decaying solution odd. A particular interior solution is $S{L}^{2}y$, so writeBoth $\psi$ and ${\psi }^{\prime }$ must be continuous at $y=\pm a$: a jump would introduce a delta function or its derivative into a source that has only finite jumps. Matching at $a$ givesand thereforeDifferentiate the interior streamfunction to obtainThe exterior profiles needed for the sketches are $u=(D/L)e^{-(|y|-a)/L}$ and $\zeta =(f/g)\mathrm{sgn}(y)D{e}^{-(|y|-a)/L}$.

For $L\gg a$, put $A=a/L$, $Y=y/L$, so $|Y|\le A$. Uniform Taylor expansion givesConsequentlyFor $f>0$ the leading interior current is westward and approximately parabolic.

For $L\ll a$ and away from the edges, $u\simeq -S{L}^2$ and $\zeta =(f/g)\psi \simeq ϵy$: the surface follows the bottom slope centrally. Near $|y|=a$, the exponentially thin corrections give positive edge jets of width $L$ and peak scale $SLa/2$ when $f>0$. Those jets decay in the exterior. The elevation is odd, has opposite-sign extrema just inside the edges, and decays to zero outside; at the positive edge $\zeta (a)\simeq ϵ(a-L)/2$. The following sketch uses the full matched solution with $L/a=0.08$. For negative $f$, the current and PV-change signs reverse, while this elevation profile does not.

Let ${\mathbf{\omega }}_a=2\mathbf{\Omega }+\nabla \times \mathbf{u}$ be absolute vorticity, where $\mathbf{\Omega }=(0,0,f/2)$ is constant. Let $\rho$ be density and $\alpha$ a materially conserved stratifying scalar, so $D\alpha /Dt=0$. For the incompressible density-stratified ideal fluid one may choose the advected density itself; in a thermodynamic description one may use an entropy-like scalar with $\rho =\rho (p,\alpha )$. The Rossby–Ertel potential vorticity isErtel's theorem states that it is conserved following a parcel under inviscid, adiabatic motion with conservative body force and the stated stratifying-scalar condition.

Start from the rotating momentum equationHere $\Phi$ includes gravity and any conservative centrifugal potential. Taking its curl givesMass conservation is $D\rho /Dt=-\rho \nabla \cdot \mathbf{u}$, specializing to $D\rho /Dt=0$ here. Combining these equations yieldsDifferentiating the material conservation of $\alpha$ gives $D({\partial }_i\alpha )/Dt=-({\partial }_i{u}_j){\partial }_j\alpha$. Dotting these two equations, the velocity-gradient terms cancel after relabelling their indices. ThereforeIf $\alpha =\rho$, the triple product vanishes immediately. More generally, an equation of state $\rho (p,\alpha )$ makes $\nabla \rho$ a linear combination of $\nabla p$ and $\nabla \alpha$, so it also vanishes. HenceThis proves the conservation law and specifies why $\alpha$ must be a suitable stratifying scalar: an arbitrary unrelated passive tracer need not annihilate the baroclinic triple product. No viscosity or scalar diffusion is permitted in this theorem.

Let $\zeta$ denote free-surface height and $h=\zeta -b$ the water depth. The inviscid rotating shallow water equations with constant Coriolis parameter $f$ areIt is $\zeta$, not $h$, that supplies the horizontal pressure gradient. For a fixed bottom, mass conservation equivalently gives $D_{H}h=-h(u_x+v_y)$.

Write $\xi =v_x-u_y$ for relative vorticity. Differentiate the second momentum equation in $x$ and the first in $y$, then subtract. The pressure terms cancel. Expanding the differentiated advection terms givesConsequently the absolute vorticity $q_a=\xi +f$ satisfies $D_H{q}_a=-q_a{\nabla }_H\cdot {\mathbf{u}}_H$. Combining this with mass conservation yields exact shallow-water potential vorticity conservation:A sloping bottom changes the depth and hence the stretching, but contributes no curl to the gravitational pressure force.

For the quasi-geostrophic approximation, choose constant reference depth $H_0$, speed $U$, and horizontal length $L$. Require Rossby number $\mathrm{Ro}=U/(|f|L)\ll 1$, slow evolution on the advective time $L/U$, smooth horizontal scales, and $L/L_D=O(1)$ or smaller, where $L_D=\sqrt{g{H}_0}/|f|$ is the Rossby deformation radius. Write $\zeta ={\zeta }_0+\eta$, $b=b_0+\delta b$, and $h=H_0+\eta -\delta b$. The leading momentum balance is geostrophic balance:It gives $\eta /H_0=O(\mathrm{Ro}(L/L_D{)}^2)$, while we require $\delta b/H_0=O(\mathrm{Ro})$ and slopes on the chosen horizontal scale. Relative vorticity divided by $f$ is also $O(\mathrm{Ro})$. ThusThe derivative of the constant first term vanishes. In advecting the first-order anomaly, the ageostrophic flow contributes only at the next order. Therefore the shallow-water quasi-geostrophic potential vorticity equation isThis retains the balanced vortical branch, including Rossby waves driven here by the topographic potential-vorticity gradient. It filters the freely oscillating inertia–gravity branches of the full shallow water equations, and does not retain independent boundary Kelvin wave dynamics. Small Rossby number alone is insufficient for the reduction: fast waves must be absent at leading order, as in well-prepared balanced initial data.

In the channel take $H_0=h_{00}$ and measure the bottom anomaly from mid-channel: $\delta b=\epsilon (y-L/2)$. Define $\beta =\epsilon f/H_0$. The resting potential vorticity is $\beta (y-L/2)$, and the linear Rossby wave equation isFor nonzero zonal wavenumber $k$, wall conditions $\psi =0$ selectThe $k=0$ modes are stationary. Write $K^2=k^2+(n\pi /L{)}^2+L_D^{-2}$ and take real $A$ for the following sketch. The total potential vorticity is $q=\beta (y-L/2)-K^2\psi$. For a contour whose undisturbed position is $y_0$, its small meridional displacement isThus potential-vorticity conservation says a northward displacement in a positive background gradient carries a negative anomaly, $q^{\prime }=-\beta \delta y$. Potential-vorticity inversion turns this into positive $\psi$ and anticyclonic motion around the bulge. At a positive cosine crest, $v={\psi }_x$ is northward on its western flank and southward on its eastern flank. Advection raises the contour on the west and lowers it on the east, so the crest moves west. Algebraically, $\delta y_t=v$ and $\delta y=(K^2/\beta )\psi$ give $\omega /k=-\beta /K^2$, the same phase velocity as above. For $\beta <0$ all propagation directions reverse.

For a square basin use a nonzero-frequency normal mode $\psi =\mathrm{Re}\{\Psi (x,y)e^{-i\omega t}\}$, with $\Psi =0$ on all walls. Its equation isSet $\alpha =\beta /(2\omega )$ and $\hat{\Psi }=e^{i\alpha x}\Psi$. Substitution cancels the first derivative and givesThe square's Dirichlet Laplacian eigenfunctions are $\hat{\Psi }=A\sin (m\pi x/L)\sin (n\pi y/L)$, with $m,n$ positive integers. Consequently the square-basin topographic Rossby mode hasFor $\beta >0$, the positive-frequency branch therefore tends to $\beta L/[2\pi \sqrt{m^2+n^2}]$ as $L/L_D\to 0$, as required. The negative-frequency solution is its complex-conjugate representation for a real disturbance; the physical frequency magnitude is $|\beta |/[2\sqrt{{\pi }^2(m^2+n^2)/L^2+L_D^{-2}}]$.

Let $\mathbf{U}=(u,v)$, $\zeta =v_x-u_y$ and $D/Dt={\partial }_t+\mathbf{U}\cdot \nabla$. The inviscid shallow water equations areTaking the vertical curl of the momentum equation, with constant $f$, gives $D\zeta /Dt=-(f+\zeta )\nabla \cdot \mathbf{U}$. Combining this with the depth equation provesThus shallow-water potential vorticity is materially conserved: changes in absolute vorticity are exactly compensated by stretching or compression of the layer depth.

The relation to a Bernoulli function presupposes steady flow. In that case define $B=g(h+H)+\frac{1}{2}|\mathbf{U}{|}^2$. The vector identity $(\mathbf{U}\cdot \nabla )\mathbf{U}=\nabla (\frac{1}{2}|\mathbf{U}{|}^2)+\zeta \hat{\mathbf{z}}\times \mathbf{U}$ givesSteady volume conservation allows a transport stream function with $h\mathbf{U}=({\psi }_y,-{\psi }_x)$. Hence $\hat{\mathbf{z}}\times (h\mathbf{U})=\nabla \psi$, soIn particular, $\mathbf{U}\cdot \nabla B=0$: the Bernoulli function is constant along each streamline. Where the flow permits a single-valued dependence on the stream function, the displayed derivative describes how that constant changes between neighbouring streamlines. For an unsteady flow, shallow-water potential vorticity remains conserved but the stated steady Bernoulli function relation does not generally hold.

Material conservation of shallow-water potential vorticity carries $P=0$ from the reservoir into the channel. At leading order in a long, slowly varying channel, $v$ and $v_x$ are negligible, so $\zeta \simeq -u_y$. Zero absolute vorticity therefore implies $u_y=f$. Integration from $y=0$, followed by cross-channel geostrophic balance, gives the zero-potential-vorticity rotating channel flow:In particular, $u_2=u_1+fb$ and $h_2=h_1-(fb/g)(u_1+fb/2)$. The velocity increases towards the shallower wall. Differentiating the depth-based specific energy confirmsThis cross-channel constancy is also the $P=0$ consequence of $\nabla B=-P\nabla \psi$.

For a reservoir of finite depth, exact $P=0$ requires relative vorticity $\zeta =-f$, an anticyclonic rotation cancelling the planetary vorticity. A reservoir at rest in the rotating frame has $\zeta =0$ and $P=f/h$, so it does not have exactly zero shallow-water potential vorticity. A very deep, nearly quiescent reservoir can nevertheless supply the approximation $P\simeq 0$ because $f/h$ is small. The idealization must describe vanishing absolute vorticity or this deep-reservoir limit, rather than ordinary finite-depth rest in the rotating frame.

Assume inviscid, unforced evolution, so shallow-water potential vorticity is conserved. If a column has initial depth $H_1$ and relative vorticity ${\zeta }_1$, moving to $H_2>H_1$ at fixed Coriolis parameter givesIts absolute vorticity is amplified by vortex stretching. A column initially without relative vorticity develops cyclonic spin: positive in the Northern Hemisphere and negative in the Southern Hemisphere. More generally the increase is in the direction of its original absolute vorticity, not necessarily the direction of $f$ if the initial spin is unusually strong.

Conservation of shallow-water potential vorticity with fixed relative vorticity givesFor ordinary planetary-vorticity-dominated columns, moving towards the equator decreases the magnitude of absolute vorticity, so the column must become shallower and spread horizontally. This is the compression counterpart of vortex stretching. For $\zeta =0$, the depth ratio is $f_2/f_1$, in either hemisphere before reaching the equator. The displayed general formula is essential if the fixed relative vorticity is large: the sign of the depth change then depends on the original absolute vorticity, and conservation cannot continue through zero absolute vorticity with positive finite depth and nonzero conserved potential vorticity.

Let $\mathbf{\zeta }=\nabla \times \mathbf{u}$ and $\mathbf{\omega }=\mathbf{\zeta }+2\Omega {\mathbf{e}}_z$ be the relative vorticity and absolute vorticity. The nonlinear acceleration identity rewrites the momentum equation asTaking the curl removes the pressure, kinetic energy gradient and shearing-sheet tidal potential. ThusBoth $\nabla \cdot \mathbf{u}$ and $\nabla \cdot \mathbf{\omega }$ vanish. The vorticity equation is consequentlyFor the absolute-vorticity flux tensor $T_{ij}=u_i{\omega }_j-{\omega }_i{u}_j$, use the flux convention $(\nabla \cdot \mathbf{T}{)}_j={\partial }_i{T}_{ij}$, with the first index denoting transport direction. This givesThe scalar $\theta$ is a buoyancy displacement variable, so dimensional consistency gives it units of length; the printed potential-temperature terminology presupposes this normalization. Its curl source directly forces only the $y,z$ components of absolute vorticity. Hence radial vorticity conservation in a shearing sheet takes the local flux formIts volume integral is constant when the net boundary flux vanishes. However $D_t{\omega }_x=\mathbf{\omega }\cdot \nabla u_x$ still permits vortex stretching; the radial component is not generally a materially conserved scalar in three dimensions. Using divergence on the second index of the stated tensor would reverse the transport sign, so the tensor convention matters.

For the Keplerian shearing sheet background $\mathbf{U}=-3\Omega x{\mathbf{e}}_y/2$, the only spatial derivative is ${\partial }_x{U}_y=-3\Omega /2$. Since $U_x=0$ and all fields are independent of $y$, $\mathbf{U}\cdot \nabla \mathbf{U}=0$ and $\nabla \cdot \mathbf{U}=0$.

The Coriolis acceleration exactly balances the tidal term:With constant pressure and $\theta =0$, the buoyancy displacement variable equation also has zero advection, zero source $U_x$ and zero diffusion. Thereforeis a steady solution. Its absolute vorticity is ${\mathbf{\omega }}_0=(\Omega /2){\mathbf{e}}_z$, since the background relative vorticity is $-3\Omega {\mathbf{e}}_z/2$.

The basic state has velocity fieldand buoyancy $B=N^{2}z$. For disturbances independent of $y$, the linearized equations areSubstituting a plane wave proportional to $\exp [i(kx+mz-\omega t)]$ and eliminating $p$, $v$, and $b$ gives the dispersion relationThus the requested coefficients are

An instability exists precisely when some wavenumber pair makes ${\omega }^2<0$. Since the Brunt–Väisälä frequency satisfies $N^2>0$, this is possible exactly whenThe basic relative vorticity is $\Lambda {\mathbf{e}}_z$, so its absolute vorticity is $(f+\Lambda ){\mathbf{e}}_z$. Its Ertel potential vorticity is thereforeThe instability criterion can consequently be written as $fQ<0$: the vertical absolute vorticity has the opposite sign to the planetary vorticity. This is inertial instability, approached most directly by disturbances with $|m/k|\gg 1$.

Letbe the vertical absolute vorticity. Taking the vertical curl of the momentum equation removes both the tidal potential and the isothermal pressure force, because each is a gradient field. The two-dimensional vorticity equation isThe surface-density equation isThe quotient rule now givesSince $\zeta =W/\Sigma$ is the vortensity,

Flow visualization and the idealization. A streamline is everywhere tangent to the instantaneous velocity field at one fixed time. A streakline is the locus, at one observation time, of particles that previously passed through a specified injection point. A marked particle's own trajectory is a pathline. These coincide geometrically in steady flow, but need not do so in an unsteady draining flow.

The intended ideal liquid is approximately incompressible and inviscid: viscosity and associated shear stresses are negligible on the scales and times being considered, the Mach number is small, and the density is effectively constant. This does not require zero vorticity; the initially rotating state is already a counterexample. It also does not justify ignoring viscosity in an arbitrarily thin boundary layer or at a real drain.

Spiralling flow and angular momentum. In an axisymmetric inviscid flow, pressure and gravity exert no torque about the vertical axis. A parcel therefore conserves axial specific angular momentum,The late, shallow remnant in the intended draining picture is supplied by parcels initially near the outer rim, so $r_0\simeq R$. The slow-drain angular-momentum amplification then givesThe assumption about the parcels' origin is necessary: torque conservation by itself does not assign $r_0=R$ to every parcel at every time.

A continuous dye feed near the outer free surface produces inward spirals that wind increasingly rapidly near the drain. The figure uses a late quasisteady model, $u_r=-D/r$ and $u_{\varphi }=C/r$ with positive constants $D,C$, giving $\varphi ={\varphi }_0-(C/D)\ln (r/r_0)$. Its streaklines and streamlines approximately coincide over an observation interval short compared with the overall draining time. Strictly, individually marked particles at the initial surface trace pathlines, not a continuously injected streakline.

The free-surface depression. Let $z=H(r)$ be the surface. Slow radial draining makes radial inertial corrections small compared with the azimuthal centripetal acceleration. The radial Euler equation and vertical hydrostatic equilibrium giveThus $g{H}^{\prime }(r)=u_{\varphi }^2/r=C^2/r^3$, with $C={\omega }_0{R}^2$. Integrating from $r$ to the rim givesFor $r\ll R$ this becomesThe surface rises outward because the inward pressure gradient force must sustain the circular motion. The dimensions are metres, and the $1/r^2$ profile is not valid inside the finite drain or vortex core.

The terrestrial estimate. Use the vertical component of Earth's rotation as the ideal initial angular speed, ${\omega }_0\simeq {\Omega }_E\sin \lambda$. At Cambridge, taking $\lambda \simeq 52.2^{\circ }$ and ${\Omega }_E=7.292\times 10^{-5}{\mathrm{s}}^{-1}$, this gives ${\omega }_0\simeq 5.76\times 10^{-5}{\mathrm{s}}^{-1}$. Hence, for $R=1\mathrm{m}$, $r=0.01\mathrm{m}$ and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$,Taking the whole terrestrial angular speed rather than its local vertical projection would give about $2.7\mu \mathrm{m}$. The Cambridge estimate uses the projection appropriate to the vertical vortex. This very small head assumes terrestrial rotation supplies the only initial swirl; imposed residual rotation would change it quadratically.

Curl of the azimuthal velocity. For axisymmetric $u_{\varphi }=r\omega (r)$ with no vertical dependence, the rotational part of the vorticity is vertical:The ideal quasisteady radial sink has no curl outside its core. The source's symbol $\Omega$ denotes this vertical vorticity component, not the fluid's angular velocity. In general its magnitude is the absolute value of the bracket.

The two limiting states. For initial uniform rotation, $\omega ={\omega }_0$, whereas the late approximation has $\omega =C/r^2$. ThereforeThe late flow can have large angular velocity and nonzero circulation while being locally irrotational. The drain core is excluded from the domain of the $1/r$ velocity field.

What is actually conserved. Inviscid does not mean that each parcel's vorticity magnitude is fixed in three dimensions. For an incompressible barotropic liquid with conservative body forces, the vorticity equation isVertical compression of a material column changes its vorticity through vortex stretching, even while its volume is conserved by horizontal expansion. In the columnar shallow-layer limit, vorticity dilution in a shallow draining layer givesThus a remnant much shallower than its original column can have much smaller vertical vorticity. At finite depth there is generally a residual value; the exactly irrotational late profile is a leading limiting approximation, not an exact parcelwise conversion of nonzero vorticity to zero in finite time.

Kelvin's circulation theorem conserves circulation around a material loop, not around a fixed spatial circle. A loop whose original radius is near $R$ has circulation about $2\pi {\omega }_0{R}^2$, consistent with the late $u_{\varphi }=C/r$ field. Fluid also leaves through the drain, so the remaining tank contents are not a fixed material system. If the calculation is expressed explicitly in the Earth-fixed rotating frame, the corresponding shallow-water potential vorticity is $(\zeta +f)/H$, with $f=2{\Omega }_E\sin \lambda$; the conserved quantity involves absolute vorticity.

Take a Northern-Hemisphere cyclonic vortex, with positive relative vorticity. For a fixed-depth incompressible layer, conservation of absolute vorticity givesThe anticlockwise flow on the east side moves parcels northward, increasing their planetary vorticity and creating negative relative vorticity anomalies. On the west side it moves parcels southward and creates positive relative vorticity anomalies. Initially,so the positive anomaly grows on the western flank while the eastern flank loses positive vorticity. The cyclonic pattern consequently shifts westward. The induced anomalies themselves drive meridional flow and further relative vorticity anomalies, transmitting a Rossby wave pattern through the surrounding fluid. The restoring process is this exchange between planetary vorticity and relative vorticity, rather than a parcel executing a permanently westward translation.

The sketch shows the initial tendency, with the generated signs separated from the initial positive core. An isolated disturbance contains many wavevectors and disperses; westward phase velocity does not imply that all its wave energy travels westward. A Southern-Hemisphere cyclonic vortex reverses the signs of the initial and generated relative vorticities, while the same beta plane mechanism remains.

For a composition-fixed equation of state, let $\dot{\theta }=D\theta /Dt$ and let $\mathbf{F}$ be nonconservative force per unit mass. The vorticity and scalar-gradient evolution equations give the source terms in the title. Spatially varying heating projected along absolute vorticity and the curl of friction can generate or destroy Ertel potential vorticity. Heating that changes potential temperature does not automatically change PV everywhere: its spatial gradient and orientation matter. Unresolved irreversible mixing can also alter the parcelwise conservation picture.