Codex @codex  0

A convective chimney is a localized deeply mixed region formed when surface buoyancy loss overcomes ocean density stratification. Small plumes, compensating flow, rotation and lateral exchange maintain its structure; it need not be a permanently uniform descending column. Brine rejection can increase its surface density forcing. A fresh summer cap can isolate a deep homogeneous core until winter forcing renews contact with the surface, so a surviving core need not imply continuous active convection.

For equal heat extraction, the density increase caused by brine rejection can exceed that caused by cooling near the freezing point because the thermal density response is weak. Using conservation of salt during freezing, the ratio is $R=(\partial \rho /\partial S)S_0{f}_i/[(1-f_i)|\partial \rho /\partial T|\Delta T]$, with $f_i=c_{p,w}\Delta T/L_f$. This is a parcel comparison, not a guarantee that freezing overcomes any ambient density stratification.

A layer with a pronounced vertical change in salinity is a halocline. When relatively fresh surface water lies over saltier water, its contribution to density stratification inhibits vertical exchange. A cold Arctic halocline can separate sea ice from warmer Atlantic water; increased heat at depth does not imply an equally increased ocean heat flux at the ice unless the intervening exchange also changes.

Two deep layers of incompressible flow with upper mass density ${\rho }_1$, lower mass density ${\rho }_2$, and tangential velocities $U_1,U_2$ have this dispersion relation when surface tension is absent. Decaying velocity potentials and the linearized kinematic boundary condition determine their amplitudes; the dynamic boundary condition for an inviscid interface gives the displayed relation. For $k>0$ the squared departure of the phase velocity $c$ from the density-weighted mean velocity is $g({\rho }_2-{\rho }_1)/[k({\rho }_1+{\rho }_2)]-{\rho }_1{\rho }_2(U_1-U_2{)}^2/({\rho }_1+{\rho }_2{)}^2$. A negative value gives Kelvin-Helmholtz instability despite stable density stratification.

Michael E. McIntyre studies geophysical fluid dynamics, including atmospheric and oceanic motion. His graduate notes develop the connection between rotation, density stratification and large-scale motion, providing a physical foundation for interpreting waves and circulation rather than treating weather maps as purely descriptive patterns.^[ref]

That connection is useful because a local parcel force balance and the evolution of a global circulation concern different scales. Conservation laws and carefully stated approximations determine when a reduced dynamical model can represent the observed flow.

In the limiting model of infinitely strong stable density stratification, small Froude number and fixed vertical scales suppress vertical displacement. Each horizontal level evolves by two-dimensional vortex dynamics; $z$ labels independent copies of the horizontal Poisson equation. The leading normalized Ertel potential vorticity is $Q/N^2=\zeta$. This is a constrained asymptotic model: finite buoyancy frequency restores vertical motion and coupling, and arbitrarily short vertical scales invalidate the independent-layer limit. McIntyre's discussion of potential-vorticity inversion explains the limiting independent-layer picture.

Ocean anoxia is a condition in which dissolved oxygen is absent or nearly exhausted in part of the ocean. Oxygen is supplied by exchange with the atmosphere, biological production and transport of ventilated water. Respiration and oxidation consume it, so weak renewal combined with a large organic-matter supply can deplete a water body.

Density stratification can restrict exchange between oxygenated surface water and deeper water. Measurements of dissolved oxygen and circulation therefore help distinguish low supply from high consumption. Euxinia additionally requires dissolved sulfide; anoxia alone does not specify which reduced chemical species are present.

Open-ocean convection is vertical mixing driven when near-surface seawater becomes dense enough to overturn the existing density stratification. Surface heat loss can increase mass density, while salinity changes also matter: freshening can stabilize the upper layer, whereas salt enrichment can promote sinking. Deep mixing therefore depends on the prior state of the water column as well as the current surface forcing.

Profiles of temperature and salinity reveal the mixed layer and its renewal. Labrador Sea observations show why comparable winter cooling need not yield the same mixing depth when the upper ocean has different freshwater content.^[ref] A convective chimney describes a localized deeply mixed region, rather than requiring a permanently uniform column sinking straight downward.

Write the background mass density as $\bar{\rho }(z)={\rho }_0+{\rho }_0^{\prime }z$ and its hydrostatic pressure as $\bar{p}(z)$, whereApply an Eulerian perturbation of a fluid variable to the background fields:Thus $\pi$ is the Eulerian pressure perturbation, in pascals; it is not the background hydrostatic pressure.

For incompressible flow, mass conservation becomes $D{\rho }_{\mathrm{tot}}/Dt=0$. Linearizing about rest givesThe Euler equations are ${\rho }_{\mathrm{tot}}D\mathbf{u}/Dt=-\nabla p_{\mathrm{tot}}-{\rho }_{\mathrm{tot}}g\hat{\mathbf{z}}$. The background pressure gradient and weight cancel. Dropping products of perturbations givesUnder $|{\rho }_0^{\prime }z|\ll {\rho }_0$, replace $\bar{\rho }$ by ${\rho }_0$ in inertia but retain the background mass density gradient in the displacement-induced buoyancy. This is the Boussinesq approximation. Its linear Boussinesq equations here areThe absolute value makes explicit that the small-density-variation condition concerns magnitude, including the usual negative ${\rho }_0^{\prime }$ of stable density stratification.

Use plane wave perturbations proportional to $e^{i(k_{x}x+k_{z}z-\omega t)}$, with amplitudes $\hat{u},\hat{w},\hat{\pi },\hat{\varrho }$. The linear equations becomeTo eliminate pressure without dividing by $k_x$, multiply the vertical momentum equation by $k_x$ and subtract $k_z$ times the horizontal momentum equation:Use incompressible flow and the equation for mass density to obtain ${\omega }^2{\rho }_0(k_x^2+k_z^2)\hat{w}=-g{\rho }_0^{\prime }{k}_x^2\hat{w}$. Hence, for a nontrivial internal gravity wave,For stable density stratification, ${\rho }_0^{\prime }<0$ and $N=\sqrt{-g{\rho }_0^{\prime }/{\rho }_0}$ is the buoyancy frequency. Then $0\le {\omega }^2\le N^2$ and the normal modes oscillate rather than grow. For ${\rho }_0^{\prime }>0$, a mode with $k_x\ne 0$ has ${\omega }^2<0$ and grows exponentially: denser fluid lies above lighter fluid. Uniform density, ${\rho }_0^{\prime }=0$, is marginal for this restoring mechanism.

The case $k_x=0$, $k_z\ne 0$, has zero frequency in the displayed relation. Incompressible flow then requires $\hat{w}=0$, so it is a degenerate non-propagating case rather than a finite-frequency vertical oscillation. These equations concern the incompressible model of part (b); its $N^2$ does not contain the adiabatic compression correction needed in part (a).

The central implication is that a helical conducting flow can continually regenerate a large-scale magnetic field, rather than merely stretch an existing component for a finite time. The alpha effect competes with microscopic and turbulent diffusivity, so dynamo action requires a sufficiently strong regenerative response on a sufficiently large scale.

In a planetary dynamo, convection in an electrically conducting interior provides motion. Rotation through the Coriolis force, together with density stratification and boundaries, can give the motion a preferred handedness and nonzero kinetic helicity. This makes the alpha effect plausible, although its sign and magnitude vary spatially. A geodynamo must also satisfy the magnetic matching conditions at the core boundaries and the global excitation threshold. Rotation alone does not automatically establish the required helicity correlations.

For stellar and galactic fields, differential rotation supplies the Omega effect, converting a poloidal magnetic field into a toroidal magnetic field. An alpha effect can regenerate the poloidal component, closing an alpha-Omega dynamo loop. The alpha-squared dynamo provides another possibility when sufficiently helical motions dominate over shear. A large-scale nearly axisymmetric mean field is compatible with the Cowling anti-dynamo theorem because the fluctuating motions and fields that generate the mean-field electromotive force need not be axisymmetric.

The growing Beltrami field calculation is a kinematic onset model, not a prediction of unlimited growth. The growing Lorentz force density modifies the flow, leading to dynamo quenching and saturation. In nearly ideal closed systems, conserved magnetic helicity constrains the simultaneous large- and small-scale field evolution; helicity transport through boundaries and finite resistivity can therefore be important. Realistic planetary and astrophysical applications require nonuniform transport coefficients, actual geometry and boundary conditions, and a nonlinear saturation mechanism.

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.

This gives the practical advection-and-inversion description. Start with interior three-dimensional quasi-geostrophic potential vorticity, advect it with ${\mathbf{u}}_g$, and recover $\psi$ fromThe stratified quasi-geostrophic inversion is elliptic for positive $N^2$. It must be supplied with lateral and vertical boundary conditions and any circulation or mean-flow constraints; interior potential vorticity alone is not sufficient. At rigid horizontal lids, $w=0$ makes boundary buoyancy satisfy $D_{g}b=0$, and $b=f_0{\psi }_z$ supplies the vertical derivative data for inversion. Boundary buoyancy perturbations can therefore support flow even when the interior anomaly is zero. Once $\psi$ is found, its horizontal derivatives give velocity, its vertical derivative gives buoyancy, and the thermodynamic equation diagnoses $w$. This balanced evolution excludes the independent fast-wave initial data present in the full equations.

Finally, a gradient in $f$ supplies a restoring mechanism even in an otherwise uniform fluid. Linearizing about rest with constant $N$, the quasi-geostrophic potential-vorticity equation becomesSubstitution of $\psi \propto e^{i(kx+ly+mz-\omega t)}$ gives the Rossby wave dispersionFor $\beta >0$, its zonal phase speed is westward relative to a resting background. A meridionally displaced parcel retains potential vorticity, so it acquires a relative-vorticity anomaly that induces the velocity tending to return the disturbance. Background shear and stratification modify the full potential-vorticity gradient; they can alter propagation or permit instability. On a constant f-plane the planetary contribution vanishes, but gradients of relative vorticity, boundary buoyancy or topography can still support balanced wave motion. Close enough to the equator that $f_0$ vanishes, this mid-latitude quasi-geostrophic approximation and its scaling must be replaced.

The basic quasi-geostrophic streamfunction $\bar{\psi }=-\Lambda yz$ gives $\bar{u}=\Lambda z$, $\bar{v}=0$ and constant quasi-geostrophic potential vorticity $\bar{Q}=f$. DefineBecause the background has no potential-vorticity gradient, the linear potential-vorticity equation is simplyIts characteristics have constant $y,z,x-\Lambda zt$, so the general solution is $q^{\prime }=\mathcal{F}(y,z,x-\Lambda zt)$. In particular, $q^{\prime }=F(y,z)G(x-\Lambda zt)$ is the requested arbitrary separated family; one product is not the most general initial field.

Linearizing the vertical velocity requires retaining advection of ${\bar{\psi }}_z=-\Lambda y$ by $v^{\prime }={\psi }_x^{\prime }$. ThereforeThe impermeable boundary at $z=0$ imposesThis is conservation of boundary buoyancy, not merely a condition that ${\psi }_z^{\prime }$ vanish.

For zero interior potential vorticity, write ${\psi }^{\prime }=\mathrm{Re}\{\varphi (z,t)e^{i(kx+ly)}\}$. With $K_h=\sqrt{k^2+l^2}>0$, stratified quasi-geostrophic inversion givesThe decaying solution isNo time dependence has been assumed to obtain this vertical structure. Applying the boundary condition gives $-\mu \dot{A}-i\Lambda kA=0$, henceThus the Eady edge wave travels eastward with horizontal phase velocity $\omega /k=\Lambda /\mu$ for $k\ne 0$.

The propagation mechanism can be made explicit using boundary parcels. Thermal wind gives the background buoyancy gradient ${\bar{b}}_y=f{\bar{\psi }}_{zy}=-f\Lambda$. A small meridional parcel displacement $d$ produces boundary buoyancy perturbation $b^{\prime }=-d{\bar{b}}_y=f\Lambda d$. Meanwhile the decaying potential-vorticity inversion gives $b^{\prime }=f{\psi }_z^{\prime }=-f\mu {\psi }^{\prime }$ at the boundary. Consequently ${\psi }^{\prime }=-\Lambda d/\mu$ there, and the parcel equation $d_t=v^{\prime }={\psi }_x^{\prime }$ becomesA displaced buoyancy contour therefore induces a velocity field that advances that displacement eastward. No interior potential-vorticity gradient is needed because the boundary buoyancy supplies the invertibility data.

Vortex stretching is what sets the induced flow's vertical reach. In the interior, the zero-potential vorticity condition requiresBoundary isopycnal displacement decreases with height as $e^{-\mu z}$, so fluid columns are stretched or compressed over depth ${\mu }^{-1}=f/(N{K}_h)$. Conservation of potential vorticity makes this stretching balance the relative vorticity. The resulting inversion factor $b^{\prime }=-f\mu {\psi }^{\prime }$ gives the propagation speed: stronger density stratification reduces penetration and speed, while larger $f$ increases them. Zero interior potential vorticity therefore does not imply zero interior flow.

For imposed boundary vertical velocity, retain the same decaying structure and writeThe boundary relation becomes the forced Eady edge wave equationThe integrating factor $e^{i\omega t}$ gives, for ${\omega }_0\ne \omega$,The arbitrary homogeneous coefficient $B$ sets the initial disturbance. With $A(0)=0$, this is $A=C(e^{-i{\omega }_{0}t}-e^{-i\omega t})/[i(\omega -{\omega }_0)]$.

At resonance, ${\omega }_0=\omega$, integration instead yieldsAn initially undisturbed fluid has ${\psi }^{\prime }=Ct{e}^{-\mu z}\cos (kx+ly-\omega t)$. The imposed pumping reinforces the same boundary normal mode every cycle, with no dissipative loss, so the linear amplitude grows secularly. The divergence describes resonance of the undamped linear model; sufficiently large wave amplitude invalidates the assumed small-disturbance approximation.

Potential-vorticity inversion is the diagnostic recovery of the flow from its advected vortical scalar, using the balance relation and specified boundary conditions. Its usefulness is a separation of tasks: transport the scalar by the current velocity, invert to find the new velocity, and repeat. The inverse is spatially nonlocal even though the transport law is local. The four models differ chiefly in the operator being inverted and the information needed at boundaries.

In two-dimensional vortex dynamics, planar incompressibility permits $\mathbf{u}=(-{\psi }_y,{\psi }_x)$, and vertical vorticity is $\zeta ={\nabla }^2\psi$. The curl of inviscid momentum givesThere is no vortex stretching in strictly planar motion, so vorticity is materially conserved. The two-dimensional potential-vorticity inversion is the Poisson equation ${\nabla }^2\psi =\zeta$. On the plane, with suitable integrability and far-field behavior, its Green's function is $(2\pi {)}^{-1}\log |\mathbf{x}-{\mathbf{x}}^{\prime }|$, givingDifferentiation gives the planar Biot-Savart kernel. Here ${\psi }_h$ is a harmonic function fixed by boundary, circulation and far-field data. A rigid impermeable wall is a streamline, so $\psi$ is constant on it; constants or circulations around distinct holes require specification. The vorticity alone cannot distinguish two flows differing by a permitted harmonic contribution, such as a uniform current on an unbounded plane. A constant added to $\psi$ has no physical effect.

In nonrotating layerwise-two-dimensional vortex dynamics, strong stable density stratification constrains motion to almost horizontal layers. In the ideal independent-layer limit, vertical velocity vanishes at leading order, $(u,v)=(-{\psi }_y,{\psi }_x)$ andhold separately at each fixed $z$. Height is a parameter: the horizontal inverse Laplacian contains no vertical derivative. With background buoyancy $N^{2}z$, the leading Ertel potential vorticity is $Q=\mathbf{\omega }\cdot \nabla b\simeq N^2\zeta$, so the materially advected normalized scalar is $Q/N^2\simeq \zeta$. The same horizontal Poisson equation is inverted layer by layer, using each layer's boundary data. This is an asymptotic constrained model, not arbitrary three-dimensional Euler motion with vertical velocity set to zero: its validity requires small Froude number and controlled vertical scales. At finite buoyancy frequency, the weak vertical motion and tilted isopycnals restore coupling. The inversion may be independent by layer while the reconstructed fields vary with height. The horizontal Euler equations for an inviscid fluid determine the pressure gradient, while hydrostatic balance relates its vertical derivative to buoyancy.

In shallow-water quasi-geostrophic inversion, begin with the exact shallow-water potential vorticity $(f+\zeta )/h$. Linearizing its anomaly about depth $H$ gives, after multiplying by $H$, $\zeta -f_0\eta /H$. For small Rossby number, geostrophic balance gives $\psi =g\eta /f_0$, while $\zeta ={\nabla }_h^2\psi$. Keeping the beta plane background gradient givesThus $q-\beta y$ is inverted with the modified Helmholtz equation $({\nabla }_h^2-L_D^{-2})\psi =q-\beta y$, under the specified boundary and far-field conditions. The velocity follows from $\psi$ and surface elevation from $\eta =f_0\psi /g$. A horizontal Fourier mode of the anomaly satisfiesThe Rossby deformation radius therefore limits the reach of the balanced response. Layer stretching, represented by $-\psi /L_D^2$, is absent from ordinary planar vorticity inversion. Unlike the pure Poisson equation, the finite-deformation operator also responds to the constant part of $\psi$, which specifies mean elevation rather than an arbitrary streamfunction offset. Specified bottom topography adds a known source to $q$; it is subtracted along with the background before inversion.

In stratified quasi-geostrophic inversion, continuously varying vertical displacement replaces the shallow layer's single stretching term. For positive $N^2(z)$ and constant nonzero $f_0$, the three-dimensional quasi-geostrophic potential vorticity isGeostrophic balance and hydrostatic balance give velocity $(-{\psi }_y,{\psi }_x)$ and buoyancy perturbation $b=f_0{\psi }_z$. Invert the anisotropic elliptic boundary value problem using interior $q$, prescribed top and bottom buoyancy as Neumann boundary conditions, and the lateral, far-field and global-flow data. Boundary buoyancy evolves by its own material conservation when the boundary is impermeable. It cannot be discarded: the Eady edge wave in the preceding solution has zero interior potential vorticity anomaly but nonzero boundary buoyancy and a nonzero induced flow.

For constant $N,f_0$, the stretched coordinates for quasi-geostrophic inversion $Z=Nz/|f_0|$ turn the inversion operator into the three-dimensional Laplacian. Its infinite-space Green's function is $-1/(4\pi |\mathbf{X}-{\mathbf{X}}^{\prime }|)$ in $(x,y,Z)$ coordinates. Hence a horizontal disturbance of scale $L$ communicates vertically over scale $|f_0|L/N$, explicitly coupling layers. In a bounded domain with homogeneous difference data, multiplying the difference of two inversions by its streamfunction and applying integration by parts givesThis proves uniqueness up to a constant when the relevant boundary terms vanish; a mean or pressure reference fixes the remaining constant. Compatible boundary fluxes are needed for a purely Neumann boundary-value problem.

Across all four models, inversion provides the induced flow, while scalar transport supplies time evolution. The independent horizontal inverse in the first two models becomes a finite-deformation inverse in shallow water and a vertically coupled elliptic inverse in continuous stratification. The principle applies within the stated balance or constrained model. In the unrestricted Boussinesq equations or shallow water equations, gravity waves and divergent motion can share the same potential vorticity, so the vortical scalar alone does not reconstruct every possible flow.

In the Boussinesq approximation, replace density by a constant reference density ${\rho }_0$ in inertia and incompressible flow, retaining the small density departure only in the gravitational force. Write $\rho =\bar{\rho }(z)+{\rho }^{\prime }$ and subtract the resting hydrostatic pressure. The upward buoyancy perturbation acceleration and squared buoyancy frequency areThus stable density stratification has $N^2>0$. Let $p$ henceforth denote the pressure anomaly divided by ${\rho }_0$. The ideal, nonrotating Boussinesq equations areThe buoyancy equation follows by materially conserving the total buoyancy $b_0(z)+\sigma$, where $b_0^{\prime }=N^2$; it does not assume constant $N$.

For an $x$-periodic flow of period $P$, define $\bar{a}=P^{-1}{\int }_0^{P}adx$. An infinite-domain average with vanishing end contributions works as well. The necessary pressure condition is no mean zonal pressure gradient, $\overline{p_x}=0$: periodic pressure satisfies it, whereas a pressure containing a term proportional to $x$ does not. Use incompressible flow to write the zonal momentum equation in conservative form,Averaging its zonal derivative terms yields the exact Reynolds stress relation

For two-dimensional small disturbances about rest, the Linearized Boussinesq equations are $u_t=-p_x$, $w_t=-p_z+\sigma$, ${\sigma }_t=-N^{2}w$, and $u_x+w_z=0$. Differentiate the two momentum equations to eliminate pressure: $(u_z-w_x{)}_t=-{\sigma }_x$. Differentiating in $x$ and using $u_{xz}=-w_{zz}$ gives ${\nabla }^2{w}_t={\sigma }_{xx}$. One further time derivative givesThe equation remains valid for $N=N(z)$. No vertical derivative of $N$ was taken: ${\partial }_x^2(N^{2}w)=N^2{w}_{xx}$.

For constant $N$ and $k\ne 0$, a plane internal gravity wave hasBoth signs give the prescribed lower-boundary velocity. The radiation condition excludes energy arriving from infinity. At positive frequency, the vertical group velocity isFor a source below the fluid, energy must propagate upward, so the outgoing solution has $m<0$, even though its vertical phase velocity $\omega /m$ is downward.

To obtain this selection causally, let $T=\mu t$ and $Z=\mu z$, and replace the constant amplitude by $a(Z,T)$. The method of multiple scales treats $(z,t,Z,T)$ as independent during differentiation: the carrier derivatives hold $Z,T$ fixed, while $a_T$ holds $Z$ fixed and $a_Z$ holds $T$ fixed. Only after differentiating do we restrict to $T=\mu t$, $Z=\mu z$. Acting on the complex wave, the operators becomeWith $K^2=k^2+m^2$, the wave equation has residualThe zeroth-order term is the dispersion relation; the first-order term vanishes precisely whenIf $a_0(T)$ is the switched-on boundary amplitude, with $a_0(T)=0$ for $T<0$, the outgoing envelope is $a(Z,T)=a_0(T-Z/b)$ for $b>0$. Its characteristics carry the boundary signal into initially undisturbed fluid. If $b<0$, characteristics instead arrive from the upper half-space; zero initial disturbances and no incoming signal cannot establish the proposed boundary-forced wavetrain. This is the internal-wave envelope radiation condition.

The leading incompressible flow polarization gives $u^{\prime }=-(m/k)w^{\prime }$. Hence, to second order in amplitude and leading order in the slow modulation,There is no meridional transport here. Envelope transport gives ${\partial }_t|a{|}^2=-b{\partial }_z|a{|}^2$, so an initially vanishing mean develops asThese are the leading expressions for internal-wave momentum deposition; higher slow-modulation corrections to the stress are beyond the displayed approximation. In particular, a spatially uniform fully established wavetrain has zero stress divergence, but the front has already produced its mean flow.

The mean wave energy density per unit mass follows from ${\sigma }^{\prime }=-i{N}^2{w}^{\prime }/\omega$:Thus $\bar{u}=(k/\omega )\mathcal{E}$. In the original frame the mean kinetic energy is ${\bar{u}}^2/2=O(a^4)$. In a frame translating at $c=\omega /k$, the undisturbed velocity is $-c$, and the change in mean kinetic energy is $-c\bar{u}+O(a^4)=O(a^2)$. Indeed $\mathcal{E}-c\bar{u}=0$ at this order behind the front. Mean-flow energy is therefore essential at the same order as wave energy in the translating frame. This follows also from the frame transformation of energy, $E\mapsto E-cP$ up to the fixed background constant, and is why conservation of momentum cannot be omitted from the moving-frame conservation of energy budget.

A convective chimney is a narrow, deeply mixed region within the surrounding stratified ocean, rather than a uniformly descending solid column. Greenland Sea observations around 2001 found features roughly ten kilometres wide with winter mixing extending to about $2.5$ km. Long-lived examples had an anticyclonically rotating core; in summer a lighter surface cap could isolate the surviving deep core, with winter open-ocean convection reconnecting it to the surface. Wadhams's 2005 observations describe this structure. A large mixed depth does not automatically mean renewal of water at that depth: the mixed water must also be dense enough to enter the surrounding deep-water mass.

Formation requires surface buoyancy loss strong enough to overcome density stratification. The regional cyclonic gyre raises density surfaces and preconditions the centre by reducing the stratification that winter mixing must penetrate. Cold, dry winter air removes heat; sea ice formation can provide additional salt forcing. Around the period relevant to this paper, warmer air, altered winds associated with the North Atlantic Oscillation, and less formation of the Odden ice tongue were identified as reducing these forcing mechanisms. Fresh surface water also resists deep mixing. Wadhams and collaborators' 2002 study discusses these possibilities, while emphasizing that its observed 2001 chimney formed through cooling without local ice formation.

When salt-poor ice forms, brine rejection leaves salt in the liquid and increases its salinity and mass density. Exporting the newly formed ice allows further freezing and repeated salt input into the remaining water. Dense parcels then sink, helping deepen the mixed region and supply open-ocean convection.

For the numerical comparison, let the liquid specific heat capacity be $c_w$ and initial salinity be $S_0$. These are not supplied by the question: take the usual approximate values $c_w=4.0\mathrm{kJ}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$ and $S_0=35$ psu for a numerical estimate. The ice heat capacity cannot replace $c_w$ during liquid cooling. Per unit area of a $100$ m column,The stated thermal density response givesExtracting the same $Q$ at the freezing point converts a liquid mass fraction $f_i=c_w(1.8)/L_f\simeq 0.02143$ into salt-free ice. Conservation of salt during freezing givesUsing the stated salinity-to-density conversion,The first-order approximation $\Delta S\simeq S_0{f}_i$ gives $0.75$ psu and a factor $4.2$. Thus freezing-to-cooling density efficiency strongly favours freezing for this cold seawater parcel, where thermal expansion is weak. It does not ensure penetration through an arbitrarily strong fresh surface cap.

For $10^4$ such bergs per year, the rate of heat absorption north of the Antarctic Convergence isRelative to the specified oceanic polar heat flux,The bergs themselves move equatorward carrying a thermal deficit; their melting removes heat from warmer water. This is equivalent to poleward heat transport by equatorward iceberg drift, explaining the sign of the inferred poleward contribution. The estimate includes the sensible and latent terms from part (b), and does not include the separate effects of released freshwater on density stratification or ocean circulation.

A steric sea level change is caused by a change of water volume at fixed ocean mass. Warming usually lowers seawater mass density, producing thermosteric sea level change; freshening also lowers density and produces a halosteric contribution. Both depend on the temperature, pressure and salinity dependence of the equation of state. Regional freshening or warming can therefore change local sea level even when the global ocean mass has not increased.

An eustatic sea level change concerns the ocean water amount or basin capacity. On the coming-century timescale, the principal positive mass contributions are melting mountain glaciers and ice caps, and loss of grounded ice from Greenland and Antarctica, including iceberg discharge. Terrestrial water storage also matters: groundwater depletion transfers water to the ocean, whereas reservoir filling retains it on land. Snowfall accumulation can partly compensate ice loss. Basin changes associated with geology matter over longer periods, and vertical land motion changes locally measured relative sea level. Melting floating sea ice has already largely displaced its eventual meltwater volume; its direct contribution is small compared with transfer of grounded ice into the ocean. Small density and salinity corrections do not alter that distinction.

For scale, a contemporary budget available in 2010 gave the following approximate contributions for 1993–2003, in $\mathrm{mm}\mathrm{y}{\mathrm{r}}^{-1}$:These are period-specific estimates with appreciable uncertainties, particularly for Antarctica, rather than exact closure of the budget. For comparison, over 1961–2003 the observed rate was about $1.8\mathrm{mm}\mathrm{y}{\mathrm{r}}^{-1}$ and the estimated thermal expansion contribution about $0.4\mathrm{mm}\mathrm{y}{\mathrm{r}}^{-1}$. The shorter interval need not define a permanent acceleration by itself, but the magnitudes show that thermal expansion and land-ice loss are both essential.

With continued warming, increasing ocean heat content sustains thermosteric sea level change, including after surface warming slows because the deep ocean equilibrates gradually. Increased ice ablation generally increases glacier runoff, although small glaciers eventually lose much of their available ice. Surface melt and changes in outlet-glacier flow can increase the ice-sheet contribution; ocean-driven loss of ice shelves can accelerate discharge of grounded ice. Greater snowfall provides some compensation but cannot be assumed to cancel these losses. Thus, under continued substantial greenhouse gas emissions, the global rate is likely to rise above its historical few millimetres per year, with an increasing contribution from the great ice sheets. The precise century-scale increase depends on emissions, ocean heat uptake and ice sheet dynamics; a single fixed rate extrapolated from one decade is inadequate.

For the idealized mixed-layer calculation, the stated density response givesAt fixed column mass, $\rho H=(\rho +\Delta \rho )(H+\Delta \eta )$, so to first orderThis is only the thermosteric sea level change from the prescribed layer and warming. The much larger total observed rise since the nineteenth century does not, by itself, prove deeper warming: part of that rise is eustatic sea level change. After independently allowing for the mass contribution, a larger remaining steric signal would require warming below $200\mathrm{m}$, a different temperature distribution, or corrections to the simplified equation of state. Direct temperature profiles do independently establish subsurface ocean warming. It reaches depth through winter mixing and open-ocean convection, especially in the North Atlantic and parts of the Southern Ocean; through subduction of surface waters into the interior; and through advection and turbulent mixing. Formation and spreading of intermediate and deep waters transmit surface properties far beyond the immediate region of heat uptake. Warm surface water need not simply sink locally: cooling, freshwater effects and density stratification govern the actual pathways.

Using $R=6.371\times 10^6\mathrm{m}$, the ocean area is $A_o=0.70(4\pi R^2)=3.57\times 10^{14}{\mathrm{m}}^2$. A liquid runoff volume $\Delta V=600\mathrm{k}{\mathrm{m}}^3=6.00\times 10^{11}{\mathrm{m}}^3$ givesThis is the contribution from that input, before allowing for other exchanges in the global ocean freshwater budget.

For the sea-ice dilution estimate, the ocean volume is $V_o=A_o(4000\mathrm{m})=1.428\times 10^{18}{\mathrm{m}}^3$. Let its initial mean salinity be $S_0$; the question does not specify it, so use $S_0\simeq 35$ psu for a numerical estimate. Interpret the $600\mathrm{k}{\mathrm{m}}^3$ as liquid freshwater-equivalent melt, with ${\rho }_f\simeq 1000\mathrm{kg}{\mathrm{m}}^{-3}$. Salt mass remains fixed while liquid water mass increases, givingConsequentlyApproximating the densities as equal gives $-1.47\times 10^{-5}$ psu per year. If the quoted volume instead denotes literal solid-ice volume, use ${\rho }_i\Delta V$ for the added mass; with ${\rho }_i=917\mathrm{kg}{\mathrm{m}}^{-3}$ the result is $-1.32\times 10^{-5}$ psu per year. Either convention gives the same order of magnitude. Floating-ice melt can change mean ocean salinity even though its direct sea-level effect is small: displaced volume and dissolved salt concentration are different budgets.

An observed net freshening is a statement about the sum of all freshwater exchanges, not a separate measurement of every positive contribution. Let $M$ be liquid ocean mass, $R$ net freshwater runoff from land, $I$ net global sea-ice melt minus formation, and $P-E$ net precipitation minus evaporation, all expressed as masses per unit time. The global ocean freshwater budget and salt conservation giveThe runoff term must include non-ice river input as well, or those other exchanges must already have been netted out. If the measured freshening equals that predicted from $I$ alone, positive runoff is still possible providedor an equivalent compensating term exists in the complete budget. For instance, excess evaporation can remove water at the same rate at which land-ice runoff supplies it; compensating net ice formation elsewhere must also be included if $I$ was initially only the Arctic contribution. Gross seasonal melt cannot be compared with an annual global trend without subtracting refreezing. Such compensation also means that the runoff-only sea-level increment need not equal the total retained increment. If runoff is additional, sea-ice melt is already net and global, and every other freshwater term is zero, the two contributions must add: under those stronger assumptions the alleged equality would be inconsistent. The apparent paradox is resolved by closing the complete freshwater and salt budgets.

Off northeast Greenland, the Northeast Water Polynya is organized by the interaction of local currents and winds with coastal geometry and relatively fixed ice barriers. A regional northward coastal current encounters the Norske Øer shelf-ice barrier; divergence and shelter on its lee side help create the southern opening. This local northward current should not be replaced in the explanation by assuming uniform southward flow along the entire Greenland coast. In summer, weak winds, the current and the barrier can maintain the opening, while solar heating and the seasonal heat budget allow it to expand northward without continuous new-ice formation.

Winter conditions differ. Strong northerly winds can remove newly formed sea ice from the northern region around Ob Bank, with the Ob Bank ice barrier and Kronprins Christian Land controlling shelter and divergence. Large air–sea temperature differences then sustain rapid freezing and export. Thus the recurring geographical opening has a changing balance of dynamical export and seasonal thermodynamics; long-lived does not mean identical formation conditions throughout the year. Ice barriers themselves can also change between years.

The North Water Polynya, or Pikialasorsuaq, occupies the Smith Sound and northern Baffin Bay region between Greenland and the Canadian Arctic. An ice arch across the narrow upstream channel restricts the arrival of thick drifting Arctic ice. Northerly winds and predominantly southward transport through the channel remove newly formed ice into Baffin Bay. The combination of blocked upstream ice supply and continuing downstream export maintains low ice concentration despite strong winter cooling. The location and stability of the arch therefore matter as well as the local wind. Warm Atlantic-origin water carried along the Greenland side can provide a sensible-heat contribution, with tidal and other mixing helping deliver heat upward; the entire opening should not be attributed either exclusively to warm-water convection or exclusively to wind.

Both regions permit much greater winter ocean-to-atmosphere heat transfer than nearby thick ice, and both can concentrate ice formation and brine rejection. The production of saline shelf water affects local density stratification and exchanges with surrounding waters; it is not automatically equivalent to deep or abyssal water formation. In spring and summer the early availability of light and open water promotes biological production and supports rich marine food webs. These persistent or recurrent openings also provide habitat and access of particular importance to Arctic animals and coastal communities. The key Arctic mechanisms are ice-barrier control of supply, current- and wind-driven divergence, and a seasonally varying contribution from ocean heat.