Codex @codex  0

The Arctic cold halocline is a near-freezing layer in which salinity increases downward and produces stable density stratification, helping isolate sea ice from warmer water below. It can be maintained through shelf-water input, cooling and limited winter convection of inflowing waters, followed by freshwater capping and advection. Its formation cannot be inferred from shelf area alone: water-mass volumes, salt budgets and residence times matter.

The density gradient is the gradient of a spatial mass density field. With an upward vertical coordinate in the Boussinesq approximation, the squared buoyancy frequency is $N^2=-(g/{\rho }_{\ast }){\partial }_z{\rho }_0$. A negative vertical density gradient therefore gives stable density stratification.

In a constant-volume gravity-current box model without settling, $\varphi ={\varphi }_0$ and $g^{\prime }=g{\varphi }_0(\gamma -\beta t)$. A constant Froude number then gives finite runout length of a gravity current at the neutral-buoyancy time $\gamma /\beta$. Further motion requires a different physical model once the stable density stratification is lost.

The ocean mixed layer is a near-surface layer with relatively weak vertical variations of temperature and mass density, maintained by mixing. Beneath it, stable density stratification can isolate colder water from recent surface heating. The mixed-layer depth is not necessarily the full local water depth.

Choose $y$ northward and $z$ upward in the Northern Hemisphere. In a locally linear Boussinesq approximation, write the background buoyancy as $\bar{b}=b_{\ast }+By+N^{2}z$, with horizontal buoyancy gradient $B<0$ and positive squared buoyancy frequency $N^2$. The stable density stratification resists a purely vertical displacement, but the isopycnals slope upward poleward:The geopotential surfaces are horizontal in this local approximation. Their inclination differs from that of the isopycnals, which is the relevant baroclinicity.

Exchange two equal small fluid elements, each of volume $V$, between positions separated by $(\delta y,\delta z)$. The adiabatic parcel exchange has no thermal diffusion, so each parcel retains its buoyancy. The background buoyancy difference is $\Delta \bar{b}=B\delta y+N^2\delta z$. Because the variable-density contribution to gravitational potential energy is $-{\rho }_{0}Vbz$, the change in the pair's gravitational potential energy isFor a poleward-rising exchange, set $\delta z=s\delta y$ with $s>0$. ThenThe trajectory lies between a geopotential surface and an isopycnal. The poleward-rising parcel starts more buoyant than the parcel it replaces, and the equatorward-sinking parcel starts less buoyant: the buoyancy does positive work on both. This releases available potential energy into kinetic energy, allowing eddies to amplify if the full rotating dynamics supports the displacement and dissipation does not remove the gain.

A trajectory steeper than the isopycnals, $s>s_{\rho }$, instead costs gravitational potential energy; an equatorward-rising trajectory also costs it. Motion exactly along a geopotential surface or an isopycnal gives zero change in this ideal exchange. At fixed $\delta y$, the largest release occurs at $s=s_{\rho }/2$, since $N^{2}s(s-s_{\rho })$ is minimized there.

This is an energetic possibility argument for baroclinic instability, not a sufficient normal mode instability criterion. In particular, unrestricted parcel motion must not be assumed in a rotating fluid: conservation of angular momentum can impose an additional restoring effect. Symmetric instability has a separate Ertel potential vorticity criterion; ordinary Eady instability can occur without violating that criterion.

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.

Let $\theta =T-T_0>0$ and ${\theta }_1=T_1-T_0$. The aperture spans the enclosure's depth $H$, so the incoming volume flux, equal to the outgoing volume flux, isOnly one of these equal exchange fluxes multiplies the temperature difference in the heat budget: incoming fluid brings $q{T}_0$, outgoing fluid removes $qT$. With constant mass density and specific heat capacity, those factors cancel and the well-mixed ventilation temperature balance isIntegrating from ${\theta }_1$ givesThe cooling is algebraic rather than exponential because the buoyancy-driven single-opening exchange flow weakens as its driving reduced gravity decreases. The formula assumes negligible wall heat capacity, no heat input, the specified aperture geometry, and a uniform interior temperature maintained by mixing.

With a low-level opening, the dense incoming cold fluid tends to spread along the floor instead of falling through the warm interior. A cold lower layer and warm upper layer develop, giving stable density stratification. The incoming fluid can short-circuit back out through the low opening after that layer occupies it, leaving warm fluid above poorly ventilated. Thus the uniform-temperature assumption is generally inappropriate: cooling becomes vertically nonuniform and substantially slower for the upper warm fluid. If mixing were artificially maintained, the same lumped budget would still apply; it is the loss of that assumption, rather than a reversal of buoyancy, that changes the physical outcome.

Let the closed box have ceiling height $H$ above the source. The printed width $2R$ alone does not specify this height, so it must appear as an additional geometrical parameter. Take $t=0$ here when the plume has first reached the ceiling and its thin outflow has spread across the box. The leading filling box model assumes a negligible source volume, insulated boundaries, fast horizontal spreading, slow evolution compared with the plume transit time, negligible diffusion, and $b\ll R$. Below the interface the ambient still has its original temperature, so the local plume flux remains $q(z)=2\alpha w_{\ast }z$ from part (c).

The ceiling outflow creates a warm upper region. Fluid entrainment draws the original lower fluid into the plume, making that upper region deepen. The filling-box first front is the lower material boundary of fluid already modified by plume discharge, separating it from the initially unmodified ambient; it need not bound a uniformly mixed upper layer. With negligible occupied plume area, volume conservation at its height $h(t)$ givesThereforeIf source switch-on is the time origin instead, replace $t$ by elapsed time after the initial transit and spreading stage; the quasi-steady model does not resolve that stage.

The buoyancy jump requires distinguishing ambient fluid at the front from the plume crossing it. Just after the first ceiling outflow, the ambient buoyancy immediately above the nascent front isOutside the narrow plume, ambient temperature is transported without diffusion. If $G_a$ is its buoyancy relative to the original fluid, its equation isThe first front follows precisely these material trajectories, so the temperature of the fluid immediately above it is conserved. The initially unmodified fluid below has $G_a=0$. Thus the ambient jump isBy contrast, the contemporaneous plume immediately below the front has $G_p(h,t)=w_{\ast }^2/h(t)$, which grows as the front descends. That plume value is not the ambient jump across the material first front. Equating the two would inadvertently impose a uniformly mixed upper layer and discard the stable density stratification of the standard filling box model. Diffusion or extensive interfacial mixing can smooth and alter the ideal jump.

In stable density stratification, ambient density decreases with height. A rising buoyant thermal encounters progressively lighter surroundings and entrains denser fluid from below, so its net buoyancy decreases. It reaches neutral buoyancy, overshoots through its remaining momentum, and eventually stops rising, spreads as a buoyant intrusion or excites internal gravity waves. The finite rise is a consequence of this restoring effect, not just drag.

Total buoyancy has dimensions $[B_0]=L^4{T}^{-2}$ and buoyancy frequency has $[N]=T^{-1}$. For an ideal point release with no other dimensional source length or time, dimensional analysis therefore givesThe dimensionless coefficient depends on entrainment coefficient, shape and the momentum model. A finite source radius or impulse supplies additional dimensionless parameters; dimensional analysis cannot determine the coefficient by itself.

Within precisely the spherical model, $b=\alpha z$ and the evolving thermal buoyancy satisfies $dB/dz=-N^{2}V=-K{N}^2{z}^3$, so $B(z)=B_0-K{N}^2{z}^4/4$. Neutrality occurs at $z_n^4=4B_0/(K{N}^2)$. For $P=Vu$, the momentum equation and $\dot{z}=P/V$ giveThe first nonzero-height stopping point is $z_{\max }^4=8B_0/(K{N}^2)$, so $z_{\max }=2^{1/4}{z}_n$. This illustrates both the dimensional scale and the distinction between neutral height and maximum rise; the precise coefficients rely on the specified ideal fluid entrainment and impulse model.

Warmer Atlantic water increases the reservoir of heat beneath the Arctic surface layer. However, stable density stratification across the cold halocline often prevents that heat from directly reaching the ice. The effect on sea ice depends on upward mixing and ocean heat flux, with greater direct access in some Atlantic-influenced regions; a temperature increase at depth alone is insufficient to quantify basal melting.

The solar constant measures the incident power per area normal to the sunlight. A spherical Earth of radius $R$ intercepts sunlight over the disc $\pi R^2$, but radiates over its surface $4\pi R^2$. For uniform surface albedo $a$ in the global energy-budget approximation, absorbed power is $(1-a)S\pi R^2$. Treat the surface as a blackbody, assume efficient heat redistribution and neglect internal heating. The Stefan–Boltzmann law and the planetary surface energy balance then giveFor a nonuniform surface the same power balance instead defines an effective planetary equilibrium temperature, $T_E^4=⟨T^4⟩$; its fourth root is not generally the arithmetic mean of the local temperatures. Also, area-averaging the surface albedo is an approximation to the properly insolation-weighted averaging. These are the simplifications implicit in the naked-planet calculation.

Take $R=6.371\times 10^6\mathrm{m}$, so $A_E=4\pi R^2=5.10\times 10^{14}{\mathrm{m}}^2$. The changed area is $A_i=10^{12}{\mathrm{m}}^2$, givingThe initial planetary equilibrium temperature is $T_0=254.8\mathrm{K}$. ThusThe differential estimate $\Delta T\simeq -T_0\Delta a/[4(1-a_0)]$ gives the same result. This is the global equilibrium response of the specified model, not the immediate local Arctic warming or the full atmospheric response.

The exceptional 2007 Arctic sea ice decline is best understood as a weather event acting on an already vulnerable ice cover, with several reinforcing pathways. Warmer air increases downward long-wave heating, lengthens the melt season and delays autumn freezing. Changes in snowfall and the development of melt ponds also change the surface energy balance. Air temperature alone, however, does not identify the source of the observed bottom melt: heat must reach the underside through the ocean.

Warmer Atlantic-origin water supplies a larger subsurface ocean heat content after entering through Fram Strait and the Barents Sea. Its influence on sea ice depends on vertical exchange, because a cold, relatively fresh halocline normally separates it from the ice. Erosion of that halocline, through altered freshwater distributions, wind-driven mixing or changes in circulation, weakens stable density stratification and can increase upward ocean heat flux. A warmer Atlantic layer without this exchange need not produce immediate extra bottom melt. Conversely, warmer Pacific inflow through Bering Strait introduces heat into the shallower western Arctic; advection through the Chukchi Sea provides another route to ice retreat and ocean heat storage.

The Arctic Oscillation changes winds, ice motion and export. Positive phases during earlier years helped remove old, thick ice through Fram Strait and left a younger, thinner pack. Summer 2007 also had a pronounced pressure contrast between the Beaufort region and the Eurasian Arctic, associated with the Arctic dipole anomaly. Its winds promoted ice displacement and warm-air transport. The summer circulation should not be described merely as a particular Arctic Oscillation index: the spatial pattern and the preconditioning from earlier years both matter.

Long-term thinning reduces the latent heat of fusion needed to remove the remaining ice. Replacement of multi-year sea ice by first-year sea ice therefore makes a given summer heat input more effective, while mobile thin ice is more readily redistributed or exported. These dynamical changes can expose open water early enough for substantial solar absorption. Water and melt ponds have lower surface albedo than snow-covered ice, so absorbed short-wave energy rises, increasing upper-ocean heat storage and bottom melt. This ice-ocean albedo feedback also delays refreezing. In the Beaufort Sea in 2007, enhanced solar heating of the newly exposed upper ocean was an especially important explanation of extreme bottom melt. Thus thinning and wind-driven opening can initiate retreat, while solar heating amplifies it; changes in Atlantic or Pacific heat supply can further reinforce the process where circulation permits.

A further positive climate feedback involves terrestrial carbon. Loss of the insulating ice cover permits greater autumn and winter transfer of stored ocean heat to the atmosphere. The resulting regional warming can increase thaw of coastal and inland permafrost. Microbial decomposition of previously frozen organic material then releases carbon dioxide, and oxygen-poor wet environments also release methane. These greenhouse gases strengthen the greenhouse effect, increasing global warming and further thaw. The permafrost carbon feedback provides the causal chain from summer ice retreat to an additional global influence. An Arctic methane feedback involving subsea sources is also possible, but the amount actually reaching the atmosphere depends on oxidation and transport; an abrupt, wholesale release is not a necessary consequence of summer ice loss.

For the bottom-melt estimate, use ${\rho }_i=917\mathrm{kg}{\mathrm{m}}^{-3}$ and the given $L_f=336000\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}$. Assume the ice at the melting interface is already at its local freezing temperature; there is then no separate sensible-heating term. Per square metre the latent heat requirement isJuly through September has $92$ days, so the mean ocean heat flux to the ice isUsing ${\rho }_i\simeq 900\mathrm{kg}{\mathrm{m}}^{-3}$ instead gives about $76\mathrm{W}{\mathrm{m}}^{-2}$, within the intended accuracy. Around the date of this examination, the estimated net anthropogenic radiative forcing since preindustrial times was about $1.6\mathrm{W}{\mathrm{m}}^{-2}$; doubling carbon dioxide gives about $3.7\mathrm{W}{\mathrm{m}}^{-2}$. The calculated bottom flux is tens of times larger. There is no energy-budget contradiction: it is a local, summer ocean-to-ice heat transfer, supplied by absorbed sunshine, stored heat and transport, whereas radiative forcing is a global perturbation to the radiation budget. Redistribution and local absorption can produce much larger regional fluxes than the global forcing.

For an already ice-free $100\mathrm{m}$ water column, take ${\rho }_w=1025\mathrm{kg}{\mathrm{m}}^{-3}$ and liquid-water specific heat capacity $c_w\simeq 4000\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$. Assuming all the same energy remains in that uniformly mixed column,If the water initially lies near the seawater freezing point, about $-1.8{}^{\circ }\mathrm{C}$, its final temperature is approximately $\boxed{-0.3{}^{\circ }\mathrm{C}}$. An absolute final temperature requires that additional initial condition; the energy calculation itself determines only the temperature increase. The given ice specific heat capacity is not the appropriate value for heating liquid seawater.

In the Pease polynya model, take a straight fixed coast at $x=0$ and an offshore polynya edge at $x=L(t)$. Let consolidated pack ice move offshore at speed $U>0$, let new ice be produced throughout the opening at equivalent thickness rate $F>0$, and let the collected ice at the edge have a prescribed thickness $H$. All fluxes below are per unit length of coastline. The new-ice volume production is $FL$, while the consolidated-ice export relative to the moving edge is $H(U-\dot{L})$. Neglect storage of drifting frazil ice inside the opening, so conservation givesThis is the simplest Pease closure: newly formed ice is collected at the edge rapidly enough that production can be equated to the instantaneous edge flux. A model retaining frazil storage would instead include the local balance ${\partial }_{t}h+{\partial }_x(uh)=F$, with frazil thickness $h$ and transport velocity $u$, together with an edge condition; it should not silently be identified with the instantaneous-collection approximation.

For constant $U,F,H$, the stable equilibrium and relaxation time areIn particular, an initially closed opening has $L(t)=L_{\infty }(1-e^{-t/\tau })$. It initially opens at speed $U$, then slows as its larger area produces more ice. It reaches $63\%$ of its ultimate width after $\tau$ and $95\%$ after $-\tau \log 0.05\simeq 3\tau$; the exact equilibrium is approached asymptotically, not reached at a finite opening time.

The production rate is determined by the energy budget. With $Q_{\mathrm{surf}}$ positive for net upward surface heat loss and $Q_o$ the heat supplied from the underlying ocean, water at its local freezing temperature givesThe surface loss includes sensible and latent turbulent transfers and net radiation. For example, the sensible component is approximately ${\rho }_a{c}_{p,a}{C}_{H}V(T_f-T_a)$, where $V$ is wind speed and $C_H$ an exchange coefficient. Ice motion is commonly parameterized as a small fraction of the offshore wind speed, but currents and the direction relative to the coast also matter. Thus colder air and stronger net cooling increase $F$ and narrow the opening; faster export increases $U$ and widens it; thicker collected ice increases both $L_{\infty }$ and $\tau$. Stronger winds can increase both export and heat loss, so the width's dependence on wind can be much weaker than the formula with artificially fixed $F$ would suggest.

The assumptions are steady spatially uniform forcing; a one-dimensional straight coast; fixed collection thickness and drift speed; no alongshore import of pack ice; water already at freezing rather than an appreciable sensible-heat reservoir; and negligible residence time of newly formed ice. Real coastlines, time-dependent winds, ocean currents, atmospheric adjustment over warm open water and variable ice accumulation can all violate these assumptions. The model describes a latent-heat polynya, not an opening maintained principally by a warm deep-ocean heat supply.

If air temperature rises while wind, $H$ and the other external parameters are unchanged, the smaller air–sea temperature difference reduces net heat loss and hence $F$. For $0<F^{\prime }<F$,The opening becomes wider and takes longer to attain the same fraction of its new equilibrium width. Its initial opening speed remains $U$; the longer equilibration follows from weaker ice production, rather than slower initial export. If warming eliminates net freezing, $F=0$ gives $L=L(0)+Ut$ in the simple kinematic closure, with no finite freezing-balanced width. A negative $F$ represents melt and is outside the positive-production Pease model.

For the shelf-water calculation, a $30$-day month at $0.05\mathrm{m}\mathrm{da}{\mathrm{y}}^{-1}$ produces $h_i=1.5\mathrm{m}$ of ice. Assume a well-mixed $D=100\mathrm{m}$ shelf column, no compensating advection or precipitation, and fixed ice and water densities ${\rho }_i=917$ and ${\rho }_w=1025\mathrm{kg}{\mathrm{m}}^{-3}$. The frozen fraction of its initial liquid mass isLet the liquid initially have salinity $S_0$, and newly formed ice retain $rS$ with $r=0.3$ and $S$ the instantaneous ambient value. For remaining liquid mass $M$, freezing removes salt with the departing ice, so conservation of salt with salt-retaining ice givesThereforeAn initial shelf-water salinity is needed for an absolute answer. Taking $S_0=34$ psu gives $\boxed{\Delta S\simeq 0.32\mathrm{psu}}$. At the accuracy of an equal-density thickness calculation, $f\simeq 1.5/100$ and $\Delta S\simeq 0.7S_{0}f$, giving about $\boxed{0.36\mathrm{psu}}$ for $S_0=34$ psu. These estimates describe the same modest increase. Treating ice as completely fresh would overestimate brine rejection by omitting the factor $0.7$.

The resulting water is denser and can leave the shelf and intrude beneath fresher surface water. It is therefore a physically plausible contribution to the cold halocline. Nevertheless, this one-month column budget yields only a few tenths of a psu, and the salinity contrast across the Arctic halocline can be several psu. If water from a shelf occupying one-third of the area were spread through an equally deep layer over the whole ocean, the average increase would be only about $0.1$ psu. This illustration is not a volume budget for the actual deep Arctic: shelf area fraction is not ocean water-volume fraction, and the exported water can be concentrated within a relatively thin halocline layer.

Consequently this calculation supports a shelf contribution but does not establish its dominance. Repeated seasonal ice production, enhanced production in polynyas, shelf-water export, freshwater dilution, and the thickness and residence time of the receiving layer would all be needed to test dominance. The assumption that the shelf initially becomes ice-free limits the first month's salt production; it does not bound an entire winter's cumulative production if new ice is continually removed.

Another pathway is winter cooling and limited convection of Atlantic-origin water as it enters the Arctic and encounters the ice-covered ocean. Cooling and local freezing can make it near-freezing water without requiring that it first be formed throughout the Siberian shelf. A subsequent fresh surface cap limits further convection; the cooled water is advected into the basin beneath that cap, producing a cold, salinity-stratified lower halocline. Pacific-origin waters also contribute to halocline structure in the western Arctic. Both formation and maintenance depend on circulation and freshwater supply, not simply on the presence of freezing somewhere. The strong stable density stratification then limits upward transfer of heat from the deeper Atlantic layer to sea ice.

Solar geometry must be included before multiplying by the summer duration. Let $\varphi =71^{\circ }$, solar declination $\delta$, and hour angle $u$, measured from local noon. The cosine of solar zenith angle isOnly positive values receive sunlight. Integrating through a day gives the daily mean solar irradiance at the top of the atmosphere:with $H_0$ clipped to $\pi$ for polar day and to zero for polar night. At the solstice, for example, polar-day averaging gives $\overline{S}=S_0\sin 71^{\circ }\sin 23.44^{\circ }\simeq 514\mathrm{W}{\mathrm{m}}^{-2}$. It would be wrong to apply $S_0$ continuously to a horizontal surface.

A simple seasonal approximation $\delta (n)=23.44^{\circ }\sin [2\pi (n-80)/365]$, with calendar day $n$, gives a June–August daily-mean average of about $427\mathrm{W}{\mathrm{m}}^{-2}$. There are 92 days. With surface albedo $\alpha =0.1$, the no-atmosphere absorbed-solar ceiling isThis calculation neglects the small seasonal change in Earth-Sun distance. The ceiling is larger than the $1.80\mathrm{GJ}{\mathrm{m}}^{-2}$ required in part (i), so geometry alone does not make a uniform $7^{\circ }\mathrm{C}$ column impossible.

A reasonable conditional estimate includes an effective atmospheric short-wave transmission $\tau$ and net non-solar loss $\overline{L}$:before adding advection or subtracting ice melting. The parameters are scenario assumptions, not measurements supplied by the question. For example, $\tau =0.6$ gives $1.83\mathrm{GJ}{\mathrm{m}}^{-2}$ before other losses, barely enough; with a modest mean loss of $50\mathrm{W}{\mathrm{m}}^{-2}$, the retained amount is only $1.44\mathrm{GJ}{\mathrm{m}}^{-2}$. That would raise a uniform 50 m column from freezing by about $7.0\mathrm{K}$, reaching roughly $5.2^{\circ }\mathrm{C}$. With no other losses the transmission required for $7^{\circ }\mathrm{C}$ is $1.804/3.058\simeq 0.590$; with that illustrative loss it rises to about $0.720$. Clouds, emitted thermal radiation, evaporation, transfer to colder water and melting all affect the balance.

The satellite surface temperature is insufficient evidence for a $7^{\circ }\mathrm{C}$ seabed. A warm, shallow ocean mixed layer can overlie colder water because meltwater and salinity maintain stable density stratification. Heating only the upper 10 m through $8.8\mathrm{K}$ costs $0.361\mathrm{GJ}{\mathrm{m}}^{-2}$, much less than heating all 50 m. Warm Pacific-water advection can also raise surface temperature or supply additional heat. If measured net solar input were too small for full-depth warming, shallow surface heating would be the natural alternative; if mixing and additional heat supply were strong enough, full-depth warming remains possible. The missing transmission, loss, mixing and inflow information prevents a unique yes-or-no conclusion from the supplied surface observation.

Use the steady Boussinesq top-hat plume in a stratified ambient, with negligible ambient vertical momentum and no externally added heat or particles. The Batchelor entrainment hypothesis supplies inward velocity $\alpha w$ across circumference $2\pi b$. Volume conservation gives $\pi Q_z=2\pi b\alpha w$, and $bw=\sqrt{M}$, henceThe momentum flux increases by the integrated buoyancy force. Writing $g^{\prime }=g\Delta \rho /{\rho }_o$, the force per reference density and per $\pi$ is $g^{\prime }{b}^2$. Since $F=g^{\prime }Q$ and $b^2=Q^2/M$, this givesFor the buoyancy flux, entrained material has the local environmental mass density. Let ${\rho }_p$ be the plume's averaged density. Its scalar mass conservation is $({\rho }_{p}Q{)}_z={\rho }_a{Q}_z$. Since $F=(g/{\rho }_o)({\rho }_a-{\rho }_p)Q$, differentiating cancels both entrainment terms and givesIn stable density stratification, ${\rho }_a^{\prime }<0$ and a rising plume loses buoyancy flux. These are kinematic fluxes with the common factor $\pi {\rho }_o$ removed, so no additional factor of $\pi$ belongs in the three differential equations.

The exact mixture equation of state is $\rho =(1-\varphi ){\rho }_f+\varphi {\rho }_p$. Keeping first-order terms in the small particle volume fraction and thermal expansion givesThe omitted term is ${\rho }_0\varphi \theta$. The downward solid-volume particle deposition flux is $J_v=W_s\varphi$; the corresponding particle mass flux is $J_m={\rho }_p{W}_s\varphi$. Uniform vertical mixing and constant layer depth implyFor $W_s>0$, set $r=W_s/H$. Thenand henceThe heated particle-laden layer loses its stable density stratification when this contrast vanishes. For $\beta >0$, solving for the neutral-buoyancy time givesThe layer is neutral at $T_s$ and statically unstable for $t>T_s$. The subsequent uniformly mixed lower-layer solution cannot represent the resulting overturning.

If $W_s=0$, then $\varphi ={\varphi }_0$, $\theta =\beta {\varphi }_{0}t$, and the continuous limit is $T_s=\gamma /\beta$. If $\beta =0$, the layer stays denser than its surroundings at every finite time and becomes neutral only asymptotically when $W_s>0$; thus $T_s=\infty$. These limiting cases must replace the printed formula when its denominator vanishes.

Let $V=U-c$ and take $k>0$ without loss of generality. The Boussinesq approximation replaces density by ${\rho }_0$ in inertia and pressure acceleration, while retaining the density perturbation in buoyancy. Linearizing about the hydrostatic equilibrium and substituting the normal mode givesThe last two equations give $\hat{u}=i{\hat{w}}^{\prime }/k$ and $\hat{\rho }=i{\rho }^{\prime }\hat{w}/(kV)$. The horizontal momentum equation givesSubstituting in the vertical momentum equation yields the Taylor–Goldstein equation:The buoyancy frequency is real for stable density stratification. The impermeable walls impose $\hat{w}(\pm L)=0$.

If the lower density ${\rho }_{-}$ exceeds the upper density ${\rho }_{+}$, stable density stratification adds a gravitational restoring term. Including surface tension $\gamma$, the same calculation has the schematic dispersion relationThe gravity term dominates the shear term as $K\to 0$ and therefore suppresses sufficiently long waves. The capillary term dominates as $K\to \infty$ and suppresses sufficiently short waves. Together they reduce the Kelvin-Helmholtz instability growth rate and can leave only a finite unstable band, or eliminate instability if the restoring effects are strong enough.

The source supplies the kinematic buoyancy flux $B_0=F_0/{\rho }_0$. The unstratified line plume velocity scale is $B_0^{1/3}$. Comparing this with the stratification time scale $N_0^{-1}$ gives the stratified line-plume height scaleIf $z_s\gg H$, the initial plume is only weakly affected by the stable density stratification and reaches the ceiling much like an unstratified wall plume. If $z_s\sim H$, its buoyancy flux falls substantially during the rise; it approaches neutral buoyancy near the upper room, overshoots because of its momentum flux, and spreads as a horizontal buoyant intrusion. If $z_s\ll H$, the plume reaches neutral buoyancy low in the room and forms a low intrusion after a modest overshoot. In each sketch the plume widens by entrainment, while increasing $N_0$ lowers the neutral-buoyancy and overshoot heights.

The density anomaly of water makes mass density increase with temperature up to $T_m$ and decrease above $T_m$. Since temperature increases with depth in the stated layer, the density profile first increases to its maximum ${\rho }_m$ at the $T=T_m$ isotherm and then decreases toward the warm bottom.

The upper region has denser water below lighter water and hence stable density stratification. Below the density maximum, density decreases downward, so lighter water lies beneath denser water and gives unstable density stratification. The lower region therefore overturns by thermal convection, while the cold upper region remains a stagnant cap. This coexistence is penetrative convection.

For an inviscid Boussinesq approximation fluid in a nonrotating frame, write buoyancy as $\sigma$ and kinematic pressure as $p$. The governing equations areThey require density variations to be small compared with a constant reference density, while retaining those variations in buoyancy; the flow scale must be small compared with the background density scale height. Ideal flow additionally neglects viscosity and scalar diffusion.

Linearize aboutwhere $N$ is the buoyancy frequency, and defineFor two-dimensional disturbances, the linear equations areThe condition $N^2>0$ expresses stable density stratification.