Codex @codex  0

Balance the inertial wave rate against horizontal nonlinear advection. The anisotropy $ϵ=|k_{\parallel }|/k_{\perp }$ is then of the same order as the Rossby number $k_{\perp }U/(2\Omega )$. Horizontal geostrophic flow gives kinematic pressure amplitude $\delta p\sim 2\Omega U/k_{\perp }\sim U^2/ϵ$. The first ageostrophic flow is $O(ϵU)$ and supplies the small horizontal divergence balancing $-{\partial }_z{u}_{\parallel }$. Rapid inertial oscillations are excluded by this slow-time ordering.

In a shallow-water potential-vorticity inversion with the background gradient explicitly compensated, uniform PV gives $({\nabla }^2-L_R^{-2}){\psi }^{\prime }=-fb/H$. The displayed Gaussian function solves it for $b=ϵ[4(1-r^2/a^2)+a^2/L_R^2]e^{-r^2/a^2}$. This is not, without that compensation, a steady solution of the standard finite-depth constant-f current problem: its residual is $U{\psi }_x^{\prime }/L_R^2$. Closed streamlines can occur while Rossby number and surface displacement remain small.

Geophysical fluid dynamics studies fluid flow on rotating, stratified planets, including oceans, atmospheres and planetary interiors. The competition among rotation, buoyancy, pressure gradients and inertia organizes its useful approximations.

For speed $U$, horizontal length $L$ and nonzero Coriolis parameter $f$, the Rossby numbercompares advective acceleration $U^2/L$ with Coriolis acceleration $|f|U$. Small $\mathrm{Ro}$ motivates geostrophic balance between horizontal pressure gradients and Coriolis force, provided rapid transients and other leading forces are absent. Stratification introduces a buoyancy frequency, while vertical hydrostatic balance can be appropriate when vertical acceleration is small.

A rotating-tank experiment makes these balances visible. In the MIT balanced-vortex apparatus, water spirals toward a drain; a camera rotating with the tank measures motion relative to it. The free surface slope supplies a pressure gradient. Comparing relative speed, tank rotation and surface slope tests the balance and its Rossby number correction. ^[ref]

This shared mathematics connects laboratory rotating fluids to planetary circulation without assuming that every atmospheric vortex is geostrophic.

The leading force balance in a rapidly rotating conducting fluid when inertial and viscous forces are small, while Coriolis force, pressure, buoyancy and Lorentz force density remain significant. The Rossby number and Ekman number measure the omitted inertia and bulk viscosity. Rotation does not balance an arbitrary magnetic torque: impenetrable axisymmetric boundaries impose the Taylor constraint. Magnetic diffusion can remain in the induction equation even when viscosity is neglected in momentum balance.

Let $L_{\perp }=k_{\perp }^{-1}$ and let $U$ be the common size of horizontal and vertical velocity perturbations, as found in part b. The linear inertial wave time and nonlinear advection time areThe critical balance in a rapidly rotating fluid is thereforeThe Rossby number is of order $ϵ$. Horizontal geostrophic flow balances the leading pressure gradient against Coriolis acceleration, soThese are dimensionally comparable statements, rather than comparisons of a dimensional frequency with a dimensional velocity. Equivalently, use $L_{\perp }$ and $(2\Omega {)}^{-1}$ as units. ThenThus rotation stays $O(1)$ in these units, while the motions evolve on the slower $O({ϵ}^{-1})$ time. Vertical advection $u_{\parallel }{\partial }_z$ is smaller than horizontal advection by $ϵ$.

Write ${\mathbf{u}}_{\perp }={\mathbf{u}}_{\perp }^{(0)}+{\mathbf{u}}_{\perp }^{(1)}+\cdots$, with ${\mathbf{u}}_{\perp }^{(1)}=O(ϵU)$, and $p=p_0+p_1+\cdots$. At leading horizontal order,Define the stream function $\Phi =p_0/(2\Omega )$, including its horizontally uniform part in this convention. ThenThe next-order incompressibility condition is not zero horizontal divergence:This term has size $ϵk_{\perp }U$ and is essential: its product with the large rotation rate is a leading contribution to the vertical vorticity equation.

To see that contribution explicitly, write $D_{\perp }={\partial }_t+{\mathbf{u}}_{\perp }^{(0)}\cdot {\nabla }_{\perp }$. At the first horizontal inertial order,Taking its vertical curl eliminates $p_1$. The two-dimensional advection identity and the two divergence orders giveEquivalently the ageostrophic flow satisfieswhose divergence is $-D_{\perp }\zeta /(2\Omega )$. This explicitly verifies that the next-order divergence has the required sign.

The vertical momentum equation, at the same slow order, isFinally ${\mathbf{u}}_{\perp }^{(0)}\cdot {\nabla }_{\perp }f=\{\Phi ,f\}$, where the canonical two-dimensional Poisson bracket is $\{f,g\}=f_x{g}_y-f_y{g}_x$. The reduced equations for an anisotropic rotating fluid areThe omitted vertical-advection, ageostrophic-advection and relative-vorticity stretching terms are smaller by $O(ϵ)$ in the equations retained here. Linearizing these equations recovers ${\omega }^2=4{\Omega }^2{k}_{\parallel }^2/k_{\perp }^2$, the anisotropic limit of part a.

Within the reduced equations for an anisotropic rotating fluid, a flow that remains purely horizontal has $u_{\parallel }=0$ and ${\partial }_t{u}_{\parallel }=0$. The second reduced equation therefore requiresThus the retained horizontal motion is two-dimensional. This is a Taylor–Proudman theorem consequence in the slow balanced ordering; setting $u_{\parallel }=0$ at only one instant would not suffice, because vertical motion could subsequently develop.

The word “exactly” applies to the reduced model. It is not a theorem for every finite-Rossby number solution satisfying only $\Omega \gg kU$. For example, the full rotating Euler equation admitsIt is an incompressible flow, has zero nonlinear advection, and obeys ${\partial }_t\mathbf{u}=2\mathbf{u}\times \mathbf{\Omega }$. It is purely horizontal but depends on $z$, and $\Omega \gg kA$ can hold. This fast inertial oscillation lies outside the slow anisotropic critical balance used in part c; it shows why that ordering cannot be dropped from the conclusion.

In the rotating terrestrial frame, horizontal pressure gradient force and Coriolis force balance in geostrophic balance. Gravity establishes the approximately vertical hydrostatic equilibrium; it is not an additional horizontal force in this balance. With eastward $x$, northward $y$, upward unit vector $\hat{\mathbf{z}}$ and Coriolis parameter $f=2\Omega \sin \varphi$,The geostrophic wind runs along atmospheric isobars, with lower pressure on its left in the Northern Hemisphere.

Interpreting the stated south-to-north gradient as $p_y>0$, $p_y=200\mathrm{Pa}/(2.00\times 10^5\mathrm{m})=10^{-3}\mathrm{Pa}{\mathrm{m}}^{-1}$, while $p_x=0$. Using the external nominal Earth's rotation rate $\Omega =7.292115\times 10^{-5}{\mathrm{s}}^{-1}$ gives the meridional-pressure-gradient geostrophic windBoth zonal velocities are westward; their signs reverse if the supplied gradient is interpreted as a northward decrease. ^[ref]

The validity criteria for geostrophic balance require horizontal acceleration and friction to be small compared with the balancing forces: Rossby number $U/(|f|L)\ll 1$ and slow evolution compared with $1/|f|$. The approximation fails near the equator, on sufficiently small or rapidly changing scales, in strongly curved flow needing appreciable centripetal acceleration, and in the frictional surface layer. Neglecting friction near the ground is the stated idealization.

For an illustrative geostrophic wind with linear drag, use acceleration $-r\mathbf{u}$ with $r>0$. The steady horizontal momentum equation givesDrag reduces the wind and hence its Coriolis acceleration; the wind crosses isobars toward lower pressure. Here it points southwest. The pressure-gradient acceleration points south, Coriolis northwest and drag northeast; the three forces sum to zero. The cross-isobar angle satisfies $\tan \alpha =r/f$ for $f>0$. A turbulent boundary layer has a vertically varying stress rather than necessarily this constant linear drag.

The quasi-geostrophic potential-vorticity equation is the leading balanced form of potential-vorticity conservation for an inviscid rotating fluid with stable density stratification, small Rossby number and an adiabatic process for each fluid element. Its parent conservation law is Ertel potential vorticity conservation. Horizontal geostrophic balance, vertical hydrostatic balance, and the beta-plane approximation reduce it to advection of one scalar by the leading horizontal geostrophic flow.

One can see its origin by combining the vertical vorticity equation with the thermal equation. Set $D_g={\partial }_t+{\mathbf{u}}_g\cdot {\nabla }_h$, ${\mathbf{u}}_g=(-{\psi }_y,{\psi }_x)$, and ${\zeta }_g={\nabla }_h^2\psi$. The quasi-geostrophic approximation givesHere $b$ is the buoyancy perturbation relative to the reference stable density stratification, $N=N(z)$ its buoyancy frequency, and $w$ the small vertical velocity. The stretching of absolute vorticity in the first equation is offset by the change of stratification represented by the second. More precisely,The last dot product is zero because thermal-wind balance makes ${\partial }_z{\mathbf{u}}_g$ perpendicular to ${\nabla }_{h}b$. Adding therefore yieldsThe constant $f_0$ has no effect on the tendency, but its inclusion makes the connection with absolute vorticity clear. Vertical advection of $q$ is beyond the retained order; the vertical velocity has already supplied the stretching and thermodynamic terms.

For a simple balanced weather forecast, the essential cycle is diagnose, advect, invert:

This is an elliptic boundary value problem, not a local algebraic diagnosis. Its solution supplies the next geostrophic flow and, through $b=f_0{\psi }_z$, the next thermal field. Repeat the cycle to move the balanced weather systems forward.

Interior potential vorticity alone is insufficient: boundary thermal anomalies can induce flow even with zero interior anomaly, as in the Eady model. The method filters fast inertia-gravity waves and is useful for large-scale extratropical weather, while tropical dynamics, convection, moisture, friction and diabatic heating require additional physics. If vertical motion is wanted, the quasi-geostrophic omega equation supplies a further diagnostic calculation.

Take $x$ eastward, $y$ northward and $z$ upward. In a rotating reference frame, the horizontal Cauchy momentum equation isHere ${\mathbf{F}}_h$ represents other accelerations, such as friction, and $f$ is the Coriolis parameter. When the Rossby number $U/(|f|L)$ and the ratio $1/(|f|\tau )$ are small, and friction is weak, both advective and local acceleration can be neglected. The Coriolis force then balances the pressure gradient force: this is geostrophic balance, with $fv=p_x/\rho$ and $fu=-p_y/\rho$. The geostrophic wind runs along atmospheric isobars, with low pressure on its left in the northern hemisphere.

A rapid change of pressure need not bring the corresponding geostrophic wind into existence at once. An initially weak wind therefore leaves an inward pressure gradient force around the new low. For a smoothed anomaly of horizontal scale $L=5\times 10^5\mathrm{m}$, geostrophic balance would requireFor the illustrative midlatitude value $|f|=10^{-4}{\mathrm{s}}^{-1}$, the implied Rossby number is already about $0.67$, so curved-flow acceleration matters. With a small tropical value $|f|=10^{-5}{\mathrm{s}}^{-1}$, the inferred speed is $333\mathrm{m}{\mathrm{s}}^{-1}$ and the Rossby number about $67$: the assumed geostrophic balance is untenable. Geostrophic adjustment redistributes the anomaly and emits internal gravity waves, leaving a weaker balanced component.

For a target speed $U_{\ast }=10\mathrm{m}{\mathrm{s}}^{-1}$ on the same scale, the required scale estimate isThus $6\mathrm{hPa}$, or $0.30$ of the initial magnitude, is the midlatitude estimate; at the illustrative tropical $|f|=10^{-5}{\mathrm{s}}^{-1}$ it is $0.6\mathrm{hPa}$, or $0.030$ of the initial magnitude. The paper gives no latitude, so there is no unique numerical reduction. At the equator itself one must use equatorial dynamics rather than divide by $f=0$. Even the tropical $10\mathrm{m}{\mathrm{s}}^{-1}$ estimate has Rossby number $2$, so it describes the scale permitted by the stated balance, rather than a controlled tropical geostrophic balance.

There are two qualifications to “unstable”. A literal top-hat pressure anomaly has a discontinuous edge, where no finite classical geostrophic wind can balance its distributional pressure gradient. The $500\mathrm{km}$ estimate assumes a smoothed edge; a thinner edge requires stronger winds. More generally a large low is not, by its amplitude alone, proved to have an exponentially growing instability. For a concrete counterexample to the claim that such an amplitude cannot be stationary, set $p\prime (s)=-P\exp (-s^2/a^2)$, with cylindrical radius $s$, constant $\rho$ and $f>0$. The cyclonic gradient windis real and satisfies $v_{\theta }^2/s+f{v}_{\theta }=p_s/\rho$ for every $P>0$, with a regular limit at the centre. This is a stationary axisymmetric flow; no claim about its perturbation stability is needed. Here the physically justified conclusion is that a rapidly imposed anomaly without its balancing wind is unbalanced.

Small tropical Coriolis parameter allows only weak horizontal pressure gradients to coexist with ordinary large-scale winds. Hydrostatic balance connects the vertical variation of these pressure gradients to horizontal temperature gradients; equivalently the approximate thermal wind relation gives $f{u}_z\simeq -(g/T)T_y$. For comparable vertical wind shear, small $|f|$ means smaller $|T_y|$. Internal gravity waves can spread heating-induced buoyancy anomalies efficiently before rotation arrests the adjustment. This motivates the weak temperature gradient approximation for the tropical free troposphere; it does not require the land surface or atmospheric boundary layer to have uniform temperature.

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.

On a beta plane, $f=f_0+\beta y$, the preceding vorticity equation isIn the broad ocean interior, assume a statistically steady flow, weak advection of relative vorticity, and negligible friction compared with planetary-vorticity advection. The constant-depth, constant-density assumptions of 2(a) exclude topographic and baroclinic vorticity generation contributions. More precisely, require $|{\xi }_t|,|\mathbf{u}\cdot \nabla \xi |,|r\xi |\ll |\beta v|$. Then Sverdrup balance isSmall Rossby number and slow forcing help justify these conditions in the interior; they do not justify them inside an intense western boundary current. For fluctuations, the mean eddy contribution to relative-vorticity advection must also be negligible.

The depth-integrated northward volume flux per unit zonal distance is $V=Hv$. For a basin $0<x<L$ and wind stress independent of $x$, the northward interior volume flux isA negative value denotes southward transport. The cancellation of $H$ is why this interior volume flux is insensitive to the chosen constant depth. In a closed, steady basin the net volume flux across a full latitude section is zero: a boundary return must compensate $T$. Thus the displayed formula cannot represent the total closed-basin transport including that return.

For a concrete circulation sketch, use a northern-hemisphere subtropical basin $0<y<B$ andThis represents westward trade winds at the southern edge and eastward westerlies at the northern edge. The negative wind stress curl drives southward interior flow. Define an ocean transport streamfunction by $Hu=-{\Psi }_y$, $Hv={\Psi }_x$ and impose no normal flow at the eastern wall through $\Psi (L,y)=0$. Integrating Sverdrup balance westward givesSince $V=0$ at $y=0,B$, this also gives no normal flow at the north and south walls. Interior water moves westward in the southern part, southward across the basin, and eastward in the northern part. A northward western boundary current completes the clockwise ocean gyre.

Pure Sverdrup balance cannot satisfy the western impermeability condition as well: ${\Psi }_I(0,y)=-LV(y)$ varies along that wall. Indeed imposing no normal flow on both east and west walls would require their constant streamfunction values to agree in this simply connected closed basin, whereas their difference here is the nonzero $T(y)$. The shaded return in the sketch represents a boundary region where the interior approximation breaks down. It is supplied by friction in 2(c), rather than inferred from Sverdrup balance alone. A southern-hemisphere subtropical ocean gyre has the reflected, anticlockwise circulation.

For an order-of-magnitude estimate, use ${\tau }_0=0.1\mathrm{N}{\mathrm{m}}^{-2}$, ${\rho }_0=10^3\mathrm{kg}{\mathrm{m}}^{-3}$, $\beta =2\times 10^{-11}{\mathrm{m}}^{-1}{\mathrm{s}}^{-1}$, $L=5\times 10^6\mathrm{m}$ and a wind stress variation scale ${\ell }_{\tau }=10^6\mathrm{m}$. ThenHere $1\mathrm{Sv}=10^6{\mathrm{m}}^3{\mathrm{s}}^{-1}$ is a Sverdrup. The sinusoidal example has maximum $|{\tau }_{s,y}|=\pi {\tau }_0/B$; taking $B=\pi \times 10^6\mathrm{m}$ makes it agree with the chosen variation scale. Other reasonable basin widths and wind profiles change the coefficient, but retain a transport of order tens of Sverdrups.

Take $x$ eastward, $y$ northward and $z$ upward, and let $u$ be the eastward geostrophic wind. With Coriolis parameter $f$ independent of $z$, geostrophic balance, hydrostatic balance, and the ideal gas law giveHere $R_d$ is the specific gas constant for dry air. Taking a vertical derivative of the first equation and interchanging mixed partial derivatives givesIt matters that $\rho$ is not constant in the pressure gradient force. Substitutingcancels the two terms proportional to $gfu/(R_{d}T)$ and yields the height-coordinate thermal wind relationIf $|u{T}_z|\ll |T{u}_z|$, this reduces toThe discarded term has relative scale $H|T_z|/T$ when $u_z\sim u/H$. Taking $H=10\mathrm{km}$, $T=250\mathrm{K}$ and a atmospheric lapse rate of $6.5\mathrm{K}\mathrm{k}{\mathrm{m}}^{-1}$ gives about $0.26$. Thus this height-coordinate formula is a useful estimate, with potentially tens-of-percent compressibility-related correction through the vertical temperature gradient; it is not an exact consequence of the two balances. If the temperature gradient vanishes vertically, the displayed approximation is exact within those balances.

The hydrostatic approximation is very good on the broad, shallow atmospheric scales considered. The geostrophic balance is normally justified by a small Rossby number and fails as a general approximation near the equator. At $30^{\circ }\mathrm{N}$, $f=7.292\times 10^{-5}{\mathrm{s}}^{-1}$. With the ideal $134\mathrm{m}{\mathrm{s}}^{-1}$ wind and a $1000$--$2000\mathrm{km}$ horizontal scale, $\mathrm{Ro}=U/(fL)\simeq 0.9$--$1.8$: quantitative accuracy is limited for this especially strong ideal jet. Small vertical aspect ratio alone does not ensure geostrophic balance, although a straight, steady zonal flow can satisfy it even when this scale estimate of the Rossby number is not small.

For the requested estimate, suppose the upper jet stream is $H=10\mathrm{km}$ above a weak surface wind, has temperature $T\simeq 250\mathrm{K}$, and spans a meridional width $L$. Integrating the approximate thermal wind relation in scale form givesThe northward temperature gradient is negative because the eastward wind strengthens upward. A $1000$--$2000\mathrm{km}$ jet width gives roughly $25$--$50\mathrm{K}$. No unique numerical difference follows without specifying the jet width, height and representative temperature; these are explicitly estimated quantities.

Stellar rotation is a problem of coupled force balance, heat transport and angular momentum transport. A star need not rotate as a rigid body: in the Sun, the convective envelope rotates faster near the equator than near the poles, while much of the radiative interior rotates more nearly uniformly. The tachocline is the thin transition near the base of the convection zone. Surface tracking and Doppler effect measurements constrain the exterior flow; solar rotational mode splitting in helioseismology constrains interior stellar rotation through integrals of the form $\delta {\omega }_{n\ell m}\simeq m\int K_{n\ell m}\Omega dV$. Multiple kernels are needed, so surface rotation is not a measurement of every depth.

A useful mean velocity decomposition is $\mathbf{U}=\varpi \Omega {\mathbf{e}}_{\varphi }+{\mathbf{v}}_m$, where $\varpi =r\sin \theta$ is distance from the axis and ${\mathbf{v}}_m$ is stellar meridional circulation. In a frame of constant angular velocity ${\mathbf{\Omega }}_0$, use the relative velocity $\mathbf{V}=\mathbf{U}-{\mathbf{\Omega }}_0\times \mathbf{r}$. Its momentum equation includesThe Coriolis force changes the direction of a parcel's motion but does no work. Thermal driving, release of differential rotational kinetic energy, contraction, or external torques supply the energy of the flows. Solid-body stellar rotation in a suitable hydrostatic equilibrium does not by itself require a meridional flow.

For slowly circulating, approximately axisymmetric material, the leading mechanical balance in an inertial frame is $\nabla p/\rho =-\nabla \Phi +\varpi {\Omega }^2{\mathbf{e}}_{\varpi }$. Taking its curl gives the stellar thermal-wind balanceHere $z$ is the axial coordinate, $S$ is specific entropy and $c_p$ is the specific heat at fixed pressure. The last form uses approximately radial hydrostatic equilibrium and uniform composition. If the material is barotropic, the right-hand side vanishes and angular velocity is constant on cylinders. Small latitudinal entropy differences can instead support the observed noncylindrical differential rotation. Stresses and inertial terms modify this leading balance. The equation relates shear to baroclinicity; it does not give the circulation speed or explain how the shear is maintained.

That second question requires an angular-momentum budget. With $j={\varpi }^2\Omega$ and density-weighted mean quantities, conservation of angular momentum can be written schematically asThe three contributions are Reynolds stress, magnetic Maxwell stress tensor and viscosity transport. Appropriate density weighting replaces the simple correlations in a strongly compressible average. Since $\nabla \cdot (\rho {\mathbf{v}}_m)=0$ in a steady star, the stationary equation gives gyroscopic pumping in a star:A torque that maintains differential rotation therefore drives a meridional response unless another stress cancels it. Conversely the circulation transports angular momentum and feeds back on the rotation. An assumed meridional pattern without its driving torque is not a closed dynamical model.

In the convection zone, outward heat transport drives convection. Rotation deflects the rising and falling parcels and makes their correlations anisotropic. Nondiffusive convective angular momentum transport can drive shear even from an initially uniform rotation rate; treating every Reynolds stress as a positive eddy viscosity would only smooth shear and miss its source. A limiting parcel that conserves $j$ when it moves outwards slows relative to uniform rotation, whereas rapid rotational constraint and correlated horizontal motions change the sign and direction of the net transport. Thus the solar equator-fast pattern requires the actual convective correlations, not a universal parcel argument. The Rossby number $u_c/(2\Omega {\ell }_c)$ compares convective inertia with rotational deflection: deep, slower, larger-scale convection is more rotationally constrained than the fastest photospheric motions. Turnover times ${\ell }_c/u_c$ range from minutes near the visible surface to weeks or longer at greater depth. On these times the Reynolds stresses can redistribute angular momentum and drive circulation, while the star's total spin changes much more slowly.

In a radiative region, rotation distorts equipotential surfaces. Under uniform rotation and uniform composition, leading hydrostatic equilibrium makes pressure, density and temperature constant on each such surface, but effective gravity varies along it. Radiative diffusion then produces a flux $\mathbf{F}=f({\Phi }_{\mathrm{eff}})\nabla {\Phi }_{\mathrm{eff}}$. Its divergence contains $f^{\prime }|\nabla {\Phi }_{\mathrm{eff}}{|}^2$, which varies over a surface and generally cannot match a nuclear source constant there. This is the radiative-equilibrium obstruction in a rotating barotropic star. Thermal imbalance drives Eddington-Sweet circulation, with heat advection enteringIts global thermal estimate isA rotational thermal imbalance of order ${ϵ}_{\Omega }L$ processes the thermal reservoir on this longer time. For the slowly rotating Sun, $t_{\mathrm{KH}}$ is of order $3\times 10^7$ years and ${ϵ}_{\Omega }$ of order $2\times 10^{-5}$, giving $t_{\mathrm{ES}}\sim 10^{12}$ years, much longer than the solar age. Stable composition gradients impede the flow further. Thus classical radiative Eddington-Sweet circulation cannot by itself establish the Sun's interior rotation rapidly. In a rapidly rotating star the same estimate is shorter and rotational mixing can compete with its evolutionary lifetime. Torque-driven convection-zone circulation and radiative thermal circulation have different controlling balances.

Magnetic fields introduce a further coupling. A poloidal magnetic field linking different layers is wound by differential rotation into a toroidal magnetic field, roughly ${\partial }_t{B}_{\varphi }\simeq \varpi {\mathbf{B}}_p\cdot \nabla \Omega$. The resulting Maxwell stress tensor transport angular momentum. Communication along a field occurs on an Alfvén speed crossing time $t_A\sim \ell \sqrt{{\mu }_0\rho }/B$, which can be far shorter than molecular diffusion if a coherent connecting field exists. In steady ideal axisymmetric induction, absence of continuing winding requires ${\mathbf{B}}_p\cdot \nabla \Omega =0$; this is constancy along the field, not automatically solid-body rotation throughout the star. The topology and strength of an interior field matter. Solar dynamo action, especially the interaction of shear with poloidal-field regeneration, couples stellar rotation to the magnetic cycle. Magnetic feedback can change shear and circulation on dynamical or cycle times; an approximately eleven-year activity cycle and twenty-two-year magnetic polarity cycle are distinct from secular spin-down.

A magnetized stellar wind also supplies an external torque. It enforces approximate corotation out to an effective Alfvén radius $R_A$, giving wind-driven magnetic braking of a solar-type star $\dot{J}\sim -{\dot{M}}_w\Omega R_A^2$. Its spin-down time is $t_J\sim I/({\dot{M}}_w{R}_A^2)$, much longer than a convective turnover time. For slowly changing moment of inertia and an unsaturated field scaling that gives $\dot{J}=-K{\Omega }^3$, integration gives ${\Omega }^{-2}={\Omega }_0^{-2}+2Kt/I$, the Skumanich rotation law. A growing or contracting star also changes $I$ and can spin up or down even without torque. Differential spin-down can create shear between layers unless internal transport couples them.

Finally, shear can excite Kelvin-Helmholtz instability, and convection can launch internal gravity waves into the stable radiative interior. Wave damping deposits angular momentum; the net transport depends on excitation and selective absorption and is not fixed by a wave-crossing time alone. Such stresses, magnetic coupling and circulation are candidate agents for the nearly uniform radiative rotation and the narrow tachocline. Molecular kinematic viscosity has a diffusion time ${\ell }^2/\nu$ generally too long to couple the entire solar interior over its age, whereas a turbulent effective kinematic viscosity can act much faster but is not appropriate to every stable layer. A circulation has an advection time $\ell /U_m$, thermal diffusion has time ${\ell }^2/\chi$, and magnetic diffusion has time ${\ell }^2/{\eta }_m$; these must be compared with the local stellar rotation and evolutionary times rather than conflated with them. The solar dynamical time $\sqrt{R^3/(GM)}$ is about half an hour, its rotation period is of order a month, and secular magnetic braking takes an evolutionary time. Rapid force adjustment, convective maintenance of shear, cyclic magnetic feedback and slow global spin evolution are separate processes. Their coupled conservation equations and timescale hierarchy are the basis of a consistent account of solar rotation.

Use pressure per unit constant density, $\pi$, with the gravitational potential and centrifugal potential absorbed into it. Strict geostrophic balance is the exact balancewith material acceleration omitted, for a homogeneous incompressible flow. Writing $f=2\Omega \ne 0$, its components are $fv={\pi }_x$, $fu=-{\pi }_y$ and ${\pi }_z=0$. Differentiating the first two equations with respect to $z$ gives $u_z=v_z=0$. Alternatively, the vector calculus identityshows, on taking the curl of the balance, that $(\mathbf{\Omega }\cdot \nabla )\mathbf{u}=0$. Thus $w_z=0$ too, proving the Taylor–Proudman theorem:The physical pressure still includes its background hydrostatic pressure; it is the reduced pressure $\pi$ that is independent of $z$. Between horizontal impermeable boundaries, the vertically constant $w$ is zero.

For boundaries $z=z_b(x,y)$ and $z=z_t(x,y)$, impermeability condition gives $w_b={\mathbf{u}}_H\cdot {\nabla }_H{z}_b$ and $w_t={\mathbf{u}}_H\cdot {\nabla }_H{z}_t$. A nearly columnar flow crossing contours of $H=z_t-z_b$ must therefore have $w_t-w_b\simeq {\mathbf{u}}_H\cdot {\nabla }_{H}H\ne 0$. Exact columnarity must give way to small corrections from ageostrophic flow. The leading horizontal geostrophic flow can remain nearly independent of $z$, and the leading reduced pressure has ${\pi }_z=0$, but $w_z$ is now nonzero and exact three-dimensional strict geostrophic balance no longer holds. For small slopes, the dominant vertical vorticity equation iswhere $q=v_x-u_y$ is relative vorticity. Equivalently, $D[(q+f)/H]/Dt\simeq 0$. Crossing geostrophic contours stretches or squashes columns: increasing $H$ produces cyclonic relative vorticity, decreasing $H$ produces anticyclonic relative vorticity. The full vorticity equation also contains tilting terms, neglected at this columnar order.

For the circular shallowing, fluid arriving from upstream has potential vorticity $f/h$. Its height is $h(1-ϵ)$ inside, so conservation givesThis uses the stated column-stretching description and initially open streamlines. To obtain the leading disturbance velocity field, use the conventional streamfunction definition $u=-{\Psi }_y$, $v={\Psi }_x$, and write $\Psi =-Uy+\psi (r)$. Then $q={\nabla }_H^2\psi =r^{-1}(r{\psi }_r{)}_r$. Regularity at $r=0$ gives ${\psi }_r=-fϵr/2$ inside. Outside, ${\psi }_r=B/r$; matching tangential velocity at $r=a$ fixes $B=-fϵa^2/2$. Choose the irrelevant additive constant so thatBoth $\psi$ and ${\psi }_r$ match at $a$, so there is no extra vorticity sheet on the rim. This is the circular topographic anticyclone. The radial expression specifies the disturbance; the complete streamfunction is $\Psi =-Uy+\psi$, since a radial function alone cannot represent the required uniform far-field current. Reversing the streamfunction sign convention reverses the displayed proportionality constant, but leaves the physical velocity unchanged.

Differentiation gives the full relative velocity field in the rotating frame:away from the narrow transition across the change of depth. The disturbance turns clockwise for $fϵ>0$ and decays as $1/r$. The horizontal field is the leading small-$ϵ$ columnar inversion: the exact variable-depth transport across a sharp rim needs higher-order corrections, rather than an exactly two-dimensional velocity continuous across a literal vertical step.

Take $U>0$ and $ϵ>0$. The greatest opposing zonal disturbance speed is $Ca$, attained at $(0,-a)$. If $Ca<U$, then $u>0$ everywhere and every particle advances in $x$, precluding closed streamlines. If $Ca>U$, the interior stagnation point $(0,-U/C)$ is a center, since the interior contours are circles about that point. The exterior stagnation point $(0,-C{a}^2/U)$ is a saddle: linearization there has real eigenvalues of opposite sign. Closed contours around the center are bounded by a separatrix through the saddle. The two points meet at the rim at onset, giving the closed-streamline threshold for a circular topographic anticyclone:At this threshold $U/(fa)=ϵ/2\ll 1$, consistent with small Rossby number. Equality is the onset, not a finite-area closed cell. Beyond onset, upstream data alone no longer determine the potential vorticity of trapped fluid; the threshold follows from the continuation of the supplied open-flow model.

The quasi-geostrophic approximation isolates slowly evolving, nearly balanced motions in a rotating, stably stratified fluid. Ocean eddies and mid-latitude atmospheric disturbances often have this character: a leading Coriolis force–pressure balance suppresses rapid horizontal acceleration, while hydrostatic balance suppresses rapid vertical acceleration. Small ageostrophic flow nevertheless supplies the vertical motion that changes the leading vorticity. The approximation filters fast inertia-gravity waves while retaining nonlinear transport, vertical shear and Rossby waves.

Use a reference density ${\rho }_0$ and remove the background hydrostatic pressure. Let $\pi =p^{\prime }/{\rho }_0$ and let the upward buoyancy perturbation be $b=-g{\rho }^{\prime }/{\rho }_0$. A vertically varying reference density stratification has $N^2(z)>0$. Under the Boussinesq approximation, the inviscid and adiabatic equations areHere $D/Dt={\partial }_t+{\mathbf{u}}_H\cdot {\nabla }_H+w{\partial }_z$. Density variations are small compared with ${\rho }_0$ but retained in the vertical buoyancy force. This is a useful local, incompressible description; a compressible atmosphere needs the corresponding weighted or pressure-coordinate version.

Let horizontal and vertical lengths be $L,H$, speed be $U$, and take the slow time $L/U$. Take $f_0>0$ for definiteness; use $|f_0|$ in the scale estimates in the opposite hemisphere. The principal small parameter is the Rossby number $\mathrm{Ro}=U/(f_{0}L)\ll 1$ at a latitude where $f_0\ne 0$. For the hydrostatic, strongly stratified regime, $H/L\ll 1$ and $f_0/N\ll 1$. Retaining both horizontal vorticity and vortex stretching gives the Burger numberThe geostrophic balance pressure scale is $f_{0}UL$ and the buoyancy perturbation scale is $B=f_{0}UL/H$. Hence $B/(N^{2}H)=\mathrm{Ro}/\mathrm{Bu}\ll 1$: buoyancy perturbations produce small displacements relative to the background stratification scale. The ageostrophic flow is $O(\mathrm{Ro}U)$ and incompressibility then gives $w=O(\mathrm{Ro}UH/L)$. Vertical acceleration is smaller than hydrostatic terms by order ${\mathrm{Ro}}^2(f_0/N{)}^2$ in this scaling. Weak forcing or diffusion can be added at the slow order, but neither is present here.

Define the quasi-geostrophic streamfunction by $\psi =\pi /f_0$. Leading geostrophic balance and hydrostatic balance implyDifferentiating horizontally and vertically gives the thermal-wind balanceThus vertical shear is determined by horizontal buoyancy gradients. The homogeneous Taylor–Proudman theorem is recovered only when those gradients vanish. Stable stratification allows sloping density surfaces and different geostrophic flow at different levels, without abandoning leading balance.

The Prandtl ratio of scales makes that modification quantitative. Horizontal and vertical terms in balanced potential-vorticity inversion have sizes $\psi /L^2$ and $(f_0^2/N^2)\psi /H^2$. Equality givesFor constant $N$, $Z=(N/|f_0|)z$ makes the inversion operator an ordinary three-dimensional Laplacian. Equivalently, an interior balanced disturbance of horizontal size $L$ naturally penetrates a depth of order $|f_0|L/N$. With $N\gg |f_0|$ this is a thin, vertically sheared structure, rather than the arbitrarily tall homogeneous columns suggested by strict Taylor–Proudman theorem behaviour. This is an aspect-ratio statement, not the diffusivity Prandtl number.

To allow a weak meridional variation of the Coriolis parameter, use the beta plane $f=f_0+\beta y$, with $\beta L/f_0=O(\mathrm{Ro})$, or equivalently $\beta L^2/U=O(1)$. The small correction to the leading balance is dynamically important on the slow time. Let $\xi ={\nabla }_H^2\psi$ be vertical relative vorticity and let $D_g={\partial }_t+{\mathbf{u}}_g\cdot {\nabla }_H$. At the first nontrivial order, the vertical vorticity equation and buoyancy equation reduce toRelative-vorticity stretching and tilting are higher order here; planetary-vorticity stretching $f_0{w}_z$ must remain. Eliminate $w$ to expose a materially conserved balanced scalar. Because ${\mathbf{u}}_{g,z}=(-b_y/f_0,b_x/f_0)$, the commutator term ${\mathbf{u}}_{g,z}\cdot {\nabla }_{H}b$ is exactly zero. For $N=N(z)$ this givesAdding the two equations proves the quasi-geostrophic potential-vorticity equationAn additive constant such as $f_0$ changes no dynamics. The horizontal streamfunction advection bracket can express this as $Q_t+{\psi }_x{Q}_y-{\psi }_y{Q}_x=0$. Although small departures from balance generated this evolution equation, the leading geostrophic flow itself advects $Q$ nonlinearly. The stretching term couples horizontal motion to vertical displacement of density surfaces.