Codex @codex  0

An oil spill under sea ice first reaches the underside through buoyancy. Oil droplets and gas bubbles rise relative to the surrounding water; their collective buoyancy drives an entraining buoyant plume, while ocean currents bend the plume and carry small droplets away. Droplet size matters: small droplets rise slowly and can be dispersed or partly dissolved before arriving, whereas coalescence produces faster-rising drops. Gas expansion, dissolution and the ambient density stratification change the plume on its ascent. Consequently the arrival footprint can be displaced from the well, and the oil arrival rate at the ice need not equal the discharge rate at the seabed.

The combinations of cover and motion give four distinct collection geometries. Each must then be combined with the four age-and-relief cases below; this accounts for all sixteen combinations without treating age as a substitute for geometry.

For each of these four geometries, the age-and-relief combinations have the following effects.

Thus deformation increases local storage in either age class; motion changes residence time and export in either geometry; and gaps permit direct surfacing in either stationary or moving ice. These distinctions remain important even when two covers have the same mean thickness.

During growth, new sea ice can form around the edges of a lens or at available water-contacting surfaces and enclose oil between successive layers. Oil interrupts direct ice–water contact and changes thermal resistance, so incorporation is not simply freezing a pure oil layer. Oil droplets can be caught in a growing skeletal layer, and connected brine channels allow upward migration when buoyancy overcomes capillary entry resistance. In spring, warming increases the liquid fraction and opens pathways; meltwater and brine circulation can redistribute oil internally. Oil can then reach the surface and melt ponds, darken the ice, or be released to open water when the enclosing ice disappears. Thick ridge pockets may survive longer than flat lenses. The coupling between brine drainage, warming and oil migration is documented by field observations of oil entrainment in first-year ice.

A shorter freezing season gives less time for encapsulation and long-term storage. Less multi-year sea ice, greater open-water exposure and earlier breakup shift more of the spill into surface slicks and mobile fragments. Earlier warming can release stored oil sooner, while reduced cover permits more wave mixing and evaporation. This does not guarantee rapid removal: a plume may still encounter fast ice or deformed reservoirs, and oil-induced darkening promotes ice-albedo feedback. These changes alter exposure, transport and recovery opportunities rather than making every warming effect favourable.

For the oil lens capillary balance, distinguish the three interfacial tensions ${\gamma }_{io}$, ${\gamma }_{ow}$ and ${\gamma }_{iw}$, where the subscripts denote ice, oil and water. Replacing an ice–water area by an ice–oil and oil–water pair costs the effective interfacial energyFor a broad, shallow, non-wetting lens of volume $\mathcal{V}=Ah$, neglect the small edge region. The interfacial energy is $\Gamma A$. Because oil is lighter than water, moving it downward through a depth $h$ costs buoyancy potential energy $\frac{1}{2}({\rho }_w-{\rho }_o)gA{h}^2$. At fixed $\mathcal{V}$ their sum isSetting $E^{\prime }(h)=0$ givesThe second derivative is positive when $\Gamma >0$. Equivalently, if $\vartheta$ is the contact angle measured through the oil, the contact-line balance gives ${\gamma }_{iw}-{\gamma }_{io}={\gamma }_{ow}\cos \vartheta$, so $\Gamma ={\gamma }_{ow}(1-\cos \vartheta )$. These are the interfacial quantities required, rather than ${\gamma }_{io}$ alone. If $\Gamma \le 0$, complete wetting replaces this finite-lens balance; if the ice is rough, pocket geometry replaces the flat-lens approximation. oil–ice–water interfacial experiments give the corresponding contact-angle description.

The printed numerical premise is defective. A crude-oil density near $0.85\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$ means $850\mathrm{kg}{\mathrm{m}}^{-3}$, rather than the printed $0.85\mathrm{kg}{\mathrm{m}}^{-3}$. Moreover, the specified oil–ice tension does not specify $\Gamma$. Even if $1.5\mathrm{N}{\mathrm{m}}^{-1}$ were assumed to be $\Gamma$, taking ${\rho }_w=1025\mathrm{kg}{\mathrm{m}}^{-3}$ would give $h=1.73\mathrm{cm}$ with the literal density and $h=4.18\mathrm{cm}$ with the plausible density repair. Neither calculation proves the claimed bound. With ${\rho }_o=850\mathrm{kg}{\mathrm{m}}^{-3}$, the exact condition isFor an explicitly defined, illustrative effective tension $\Gamma =0.050\mathrm{N}{\mathrm{m}}^{-1}$, the same formula gives $h=7.63\mathrm{mm}$. This is a consistent sub-centimetre example, not a uniquely recoverable correction of the printed interfacial data.

For ice-floe oil painting fraction, let $Q$ be the oil volume arriving per day and assume complete retention in a layer of specified thickness $h$ on the passing underside. In time $\Delta t$, a plume of width $d$ sweeps a strip of area $dV\Delta t$, apart from the initial end correction. Coating that strip needs $hdV\Delta t$ of oil. The circular plume area cannot replace the swept width: $AVh$ has dimensions of volume times length per time. ThereforeDeposition becomes discontinuous when $V>V_c$. For $d=100\mathrm{m}$ and $V=10^4\mathrm{m}\mathrm{da}{\mathrm{y}}^{-1}$, the four results areThese percentages refer to the swept, $d$-wide corridor. The percentage of a whole floe of unspecified width is not determined. They are volume-equivalent coverage fractions at the stipulated local thickness; real deposition may instead give a thinner continuous smear, and losses reduce the retained fraction. Rough pockets store more oil per painted area, explaining their smaller coverage at fixed $Q$. The oil is exported as a patchy record of the floe's passage, so later melting can release many dispersed patches rather than a single stationary slick. Low areal coverage is compatible with substantial oil mass and locally thick contamination.

Preparedness should follow those transport stages. Before release, install well-control and subsea capping capability, winter-operable access, warmed pumping equipment, oil storage, under-ice survey capability and drifting position beacons; map currents, sea ice motion and underside traps. Once a blowout starts, stopping or collecting oil near the source limits all later exposure. Survey the actual plume landfall with subsea vehicles and suitable acoustic or optical sensors, track both water and ice trajectories, and sample to distinguish an oil layer from ordinary ice roughness.

Under stable fast ice, targeted access holes and suction collectors can recover accessible reservoirs before encapsulation, provided the underside survey identifies their position and connectivity. Under moving ice floes, collectors must follow or intercept the contaminated floes; fixed surface booms cannot recover oil sealed beneath an intact slab. In broken ice and leads, use ice-capable skimmers and containment where relative ice and water motion allow it. Controlled burning is an option only for an accessible surface layer of sufficient thickness; dispersants concern exposed or subsea droplets with adequate mixing, not a sealed lens. During spring release, return to tracked pockets and floes, collect oil reaching leads and melt ponds, and prepare coastal interception along the forecast export route.

The research questions are specific: detection accuracy against natural underside roughness; droplet-size and plume-landfall prediction; capillary entry into saline pores; trapping capacities and remobilization thresholds of ridges; oil-dependent growth, viscosity and pumping rates at low temperature; timing of spring channel opening; and recovery, burning or dispersion efficiency and environmental effects under realistic ice and mixing conditions. A plan must quantify recoverable mass and operating windows rather than assume any one technique works in all ice states. Arctic response-gap research explicitly treats weather, ice, access and technique limitations.

A convective chimney is a localized, deeply mixed region within a stratified ocean, formed when surface buoyancy loss allows convection to penetrate far below the surrounding mixed layer. The chimney describes the mixed water column and its circulation; instantaneous downward transfer occurs through many smaller plumes, accompanied by compensating motion and lateral exchange, rather than through one permanently uniform downward pipe. Open-ocean convection occurs in the Labrador and Irminger Seas, the Greenland Sea, the northwestern Mediterranean, and episodically in the open Weddell Sea under a large polynya. Dense-water formation on a continental shelf is related but geometrically distinct from deep open-ocean chimney convection.

Observed Greenland Sea convective chimneys include compact, approximately circular, weakly stratified cores with temperature and salinity much more uniform vertically than in the ambient ocean. A mapped example was about $20\mathrm{km}$ across, extended to about $2500\mathrm{m}$ in winter, and rotated anticyclonically. In summer a fresh, low-density cap and shallower lateral intrusions could isolate the deep core, while winter buoyancy loss removed the cap and renewed deep contact. Such cores could persist through several seasons, so presence of a deep homogeneous core is not proof that vigorous convection is occurring continuously. These particular dimensions and seasonal structure come from successive hydrographic surveys of Greenland Sea chimneys.

Formation is favoured by the cyclonic Greenland gyre, which domes density surfaces and can reduce the depth of the stratified barrier. Cold, dry wind outbreaks produce strong sensible, latent and longwave heat loss. Where frazil ice and pancake ice form in the Odden ice tongue, ice growth removes relatively fresh material while brine rejection increases the density of the remaining water. Ice advection can keep exposing water for further freezing and separate local brine release from later freshwater return. This combination of preconditioning, cooling and sustained salt input can make a formerly buoyant surface layer dense enough to convect.

Conversely, a fresh surface cap supplied by Arctic freshwater and ice melt strengthens density stratification. Warmer atmospheric outbreaks reduce surface heat loss and ice production; circulation changes affecting the cold northwesterly winds and the Odden can reduce the brine source. Atlantic advection supplies both heat and salt, with competing effects on density, so temperature alone is not an adequate diagnostic. The historical decline in Odden ice formation and changes in wind regimes offer a mechanism for less favourable chimney production around the period considered, but ice-driven formation is a physical hypothesis, not evidence that every observed chimney requires ice or that every local change follows a single climate index. The Greenland Sea salt-flux model explicitly includes ice production and advection.

For the brine rejection mechanism, let ${\dot{h}}_i>0$ be ice growth and let $s_w,s_i$ be water and ice salt mass fractions. Newly formed ice removes mass ${\rho }_i{\dot{h}}_i$ per area per time. Compared with removing water at concentration $s_w$, it leaves an excess salt flux in the liquid ofThe units are $\mathrm{kg}{\mathrm{m}}^{-2}{\mathrm{s}}^{-1}$. With the exam's approximation $s=10^{-3}S$ for salinity $S$ in psu, the corresponding mixed-layer salinity tendency is, to first order,If ${\beta }_S=(1/{\rho }_w)(\partial \rho /\partial S)>0$, the surface buoyancy input is negative, with scale $-g{\beta }_S({\rho }_i/{\rho }_w){\dot{h}}_i(S_w-S_i)$. The enhanced density produces downward saline plumes that deepen mixing if this buoyancy loss overcomes the available density stratification. This explains why exporting new ice while retaining its rejected salt can strengthen overturning, whereas later local melt reverses the freshwater effect.

For the numerical comparison, the cover gives the specific heat capacity of ice, not of seawater. A seawater value must therefore be supplied explicitly. Take $c_{p,w}=4000\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$ and use the stipulated ${\rho }_w=1025\mathrm{kg}{\mathrm{m}}^{-3}$. A $100\mathrm{m}$ column has initial water mass per areaCooling it by $1.8\mathrm{K}$ removesThe specified thermal density relation givesThis assumes the whole column mixes and neglects entrainment and advected heat. Using the supplied ice heat capacity $2100$ for liquid seawater would change $Q$ incorrectly.

At the freezing point the same heat extraction produces a pure-ice mass per areaHere $f_i$ is the fraction of the original water mass removed into ice. With the explicit ice-density approximation ${\rho }_i=900\mathrm{kg}{\mathrm{m}}^{-3}$, the resulting ice thickness isThe water-equivalent thickness uses a reference water density: with fresh water at $1000\mathrm{kg}{\mathrm{m}}^{-3}$, it is $m_i/1000=\boxed{2.196\mathrm{m}\simeq 2.2\mathrm{m}}$. The depth of the original seawater removed is instead $m_i/1025=\boxed{2.143\mathrm{m}}$. These three quantities are different; $2.2\mathrm{m}$ is not the physical ice thickness for the stated density approximation.

Let the initial salinity be $S_0$. Complete exclusion of salt from the ice leaves all the initial salt in a liquid mass $M_1=M_0-m_i$. Conservation of salt during freezing givesThe initial salinity is not printed, so a unique numerical salinity increase requires another explicit assumption. For $S_0=35\mathrm{psu}$,The total density increase after cooling followed by freezing is then about $0.9464\mathrm{kg}{\mathrm{m}}^{-3}$. The factor $1/(1-f_i)$ is required because salt concentration refers to the remaining liquid mass. Using ice thickness divided by the original water depth as the removed mass fraction would violate this accounting.

For this assumed salinity, the freezing-to-cooling density efficiency shows that equal heat extraction at freezing produces about $0.7664/0.18=\boxed{4.26}$ times the density increase of the initial cooling. At low temperature the thermal expansion response is weak, whereas brine rejection changes salinity efficiently. This supports strong freezing-driven convection where water is already near its freezing point. It does not mean freezing always creates deep chimneys: the initial density stratification, freshwater supply, mixing depth and later melting remain decisive. The freezing point itself decreases as salinity rises; keeping it at $-1.8^{\circ }\mathrm{C}$ is the exam's approximation. Allowing that small additional cooling and a temperature-dependent equation of state would require a refined heat budget.

Important other sites of freezing-enhanced downward transfer include the Ross and Weddell Antarctic shelves and their coastal polynyas, Adélie Land and Cape Darnley, Arctic shelves such as the Barents and Kara Seas, and the Sea of Okhotsk. On shelves the saline dense water can cascade down a continental slope and feed deeper waters; the same salt-removal balance operates without requiring an isolated open-ocean chimney. A later direct observational example is bottom-water formation associated with Cape Darnley sea-ice production.

Atmospheric warming of the Greenland Sea can reduce surface cooling and ice formation, while increased freshwater input can inhibit convection. If these changes weaken the Atlantic meridional overturning circulation and its northward heat transport, their ocean-mediated tendency is relative cooling, or reduced warming, of northwestern Europe compared with otherwise identical atmospheric warming. This is not a prediction of net European cooling under all emissions scenarios, nor does one chimney determine the whole overturning circulation: other convection regions, winds, freshwater pathways and coupled ocean adjustment also matter. The physical conclusion is the sign of the heat-transport response, not a specified century-scale temperature decrement.

Irrotational flow permits $\mathbf{u}=\nabla \varphi$, and incompressible flow imposes $\nabla \cdot \mathbf{u}=0$, hence $\Delta \varphi =0$. In the laboratory frame, the no-penetration condition at the translating sphere is ${\partial }_r\varphi (a,\theta )=U\cos \theta$.

For separation of variables, take the angular factor $g(\theta )=\cos \theta$, which satisfies $\frac{1}{\sin \theta }(\sin \theta g^{\prime }{)}^{\prime }=-2g$. The radial equation iswith solutions $r$ and $r^{-2}$. Rest at infinity excludes the $r$ term. The normal boundary condition fixes the remaining coefficient, giving the potential flow around a translating sphereThere is no azimuthal velocity. With fluid density ${\rho }_f$, the kinetic energy isThe angular integral is $2\pi$, and the radial integral is $1/(3a^3)$, so $K_f={\rho }_f\pi a^3{U}^2/3$. Since $M=(4\pi a^3/3){\rho }_s$ and $k={\rho }_f/{\rho }_s$, $K_f=kM{U}^2/4$. Thus the added mass of a sphere is $kM/2$.

During a fall through $h$, gravity does work $Mgh$ and buoyancy does work $-kMgh$. Starting from rest, conservation of energy givesConsequently $U=\sqrt{4(1-k)gh/(2+k)}$. The mass cancels: the reduced speed comes from both buoyancy and the fluid's added mass.

Let $K$ be intrinsic permeability of a porous medium, $\varphi$ the mobile pore fraction, and let $h$ measure thickness normal to the sloping cap. Define cold and reservoir-temperature properties byThe imposed inequality makes ${\delta }_c>0$, so the cold fluid is still buoyant. Assume $K>0$, $\sin \theta >0$, positive densities and viscosities, and the usual positive thermal-expansion and viscosity coefficients. Away from localized fronts, neglect the along-current thickness gradient. Hydrostatic water pressure and Darcy's law then give the upslope Darcy velocitiesTypically $u_h>u_c$: warming increases buoyancy and decreases dynamic viscosity.

Current depths and the moving moving thermal front in a porous current. For this slope-driven porous gravity current, the injected mass flux per unit well length fixes the cold plateau thickness:The moving thermal front in a porous current moves at $V_T=\Gamma u_c$. Since $u_c$ is explicitly a Darcy velocity, the fluid mass stored per unit horizontal area is $\varphi \rho h$. Apply mass conservation across the moving moving thermal front in a porous current:Writing $\gamma =\varphi \Gamma$ and $r=u_c/u_h$, this gives the distal warm depthFor a retarded moving thermal front in a porous current take $0\le \gamma <1$; otherwise this two-plateau geometry needs reconsideration. The factor $(1-\gamma )/(1-\gamma r)$ is important: the distal mass flux is not generally $Q$, because mass is being stored as the cold region replaces the warm region. Setting the two fluxes equal would silently assume a stationary moving thermal front in a porous current.

If the model absorbs porosity into storage and uses pore velocity as its transport speed, the same formulas use $\gamma =\Gamma$. The explicit Darcy convention in the PDF instead gives $\gamma =\varphi \Gamma$; porosity must be specified or absorbed consistently. In the stationary-front limit $\gamma =0$, the two steady depths reduce to $Q/({\rho }_c{u}_c)$ and $Q/({\rho }_h{u}_h)$. In the Boussinesq approximation, use a common reference density in the mass factors while retaining ${\delta }_c,{\delta }_h$ in the driving force.

Leakage thresholds. A normal thickness $h_i$ generates cap overpressure ${\Pi }_i={\delta }_{i}g{h}_i\cos \theta$. Before leakage modifies the current, the two plateau values areFor the caprock leakage threshold $\Delta p>0$, the candidate critical mass-injection rates are thereforeCold-only leakage occurs for $Q_c^{\ast }\le Q<Q_h^{\ast }$, provided $Q_c^{\ast }<Q_h^{\ast }$. Both plateaus can leak once $Q\ge \max (Q_c^{\ast },Q_h^{\ast })$. In the usual common-density approximation, ${\mu }_c>{\mu }_h$ and $r<1$ ensure $Q_c^{\ast }<Q_h^{\ast }$, giving the expected sequence: no leakage, cold-only leakage, then leakage from both cold and warm regions. Without that approximation the ordering must be checked; the given inequality $\Delta \rho >\alpha \Delta T$ alone does not establish it. These are onset criteria computed on the nonleaking current, not a post-leakage mass budget.

For small slopes $\cos \theta \simeq 1$ and $\tan \theta \simeq \sin \theta \simeq \theta$. If another thickness convention is used, its hydrostatic column and projected flux must be changed consistently; one should not mix a normal thickness with a vertical-pressure formula lacking the cosine.

Cross-current heat conduction. Heat transfer from the warm formation makes temperature vary across the carbon dioxide depth and introduces a warming time controlled by thickness, thermal diffusivity and the surrounding rock's heat capacity. It smooths the sharp thermal adjustment and causes gradual changes in buoyancy, viscosity and the velocity profile. The cold fluid generally warms, becomes more mobile and requires less depth to carry a prescribed mass flux; the enhanced cold-region overpressure and its distinct leakage zone tend to shrink. Warming may occur before a parcel reaches the idealized advective moving thermal front in a porous current, especially for a thin current. A single supplied $\Gamma$ no longer describes all heat transport. Quantitative depths and thresholds then require a coupled temperature equation and thermal boundary data; their exact changes cannot be inferred from the linear property laws alone.

Let $g=Gm/r^2>0$, $H_P=-dr/d\log P>0$, $\nabla =d\log T/d\log P$, and ${\nabla }_{\mu }=d\log \mu /d\log P$. Consider a small radial fluid displacement $\xi$ whose sound-crossing time is short enough to keep its pressure equal to its surroundings. Its composition is frozen and its heat exchange negligible. This is the local stellar convective stability test.

For gas-pressure-dominated ideal gas matter, $\rho \propto \mu P/T$. The ambient gradient is thereforeThe displaced element conserves specific entropy and mean molecular weight, giving $(d\log \rho /dr{)}_{\mathrm{parcel}}=-(1-{\nabla }_{\mathrm{ad}})/H_P$, where ${\nabla }_{\mathrm{ad}}=({\Gamma }_2-1)/{\Gamma }_2$. The parcel-minus-environment mass density difference isIts buoyancy acceleration is minus $g$ times this difference. ThusPositive stellar buoyancy frequency squared gives a restoring force. The gas-pressure-dominated Ledoux criterion is thereforeEquality is marginal in this ideal adiabatic test. For monatomic gas ${\nabla }_{\mathrm{ad}}=2/5$. An inward increase of mean molecular weight has ${\nabla }_{\mu }>0$ and stabilizes the layer; uniform composition recovers the Schwarzschild criterion. For a radiative layer one substitutes the stellar radiative temperature gradient, ${\nabla }_{\mathrm{rad}}=3\kappa LP/(16\pi acGm{T}^4)$.

More generally define $\delta =-(\partial \log \rho /\partial \log T{)}_{P,\mu }$ and $\phi =(\partial \log \rho /\partial \log \mu {)}_{P,T}$. The same displacement argument gives $N^2=(g/H_P)[\delta ({\nabla }_{\mathrm{ad}}-\nabla )+\phi {\nabla }_{\mu }]$ and stability for $\nabla <{\nabla }_{\mathrm{ad}}+(\phi /\delta ){\nabla }_{\mu }$. Both coefficients equal one in the gas-dominated ideal-gas limit used above. This is local dynamical stability against convection; it is distinct from global radial collapse and from instabilities requiring heat or composition diffusion.

Take a static, nonrotating, nonmagnetic spherical equilibrium with $\nabla P=-\rho \nabla \Phi$ and ${\nabla }^2\Phi =4\pi G\rho$. Use a fluid displacement $\mathbf{\xi }(\mathbf{r})e^{-i\omega t}$. Write $p_1,{\rho }_1,{\varphi }_1$ for Eulerian fluid perturbations. The Lagrangian pressure perturbation and corresponding mass density change obey $\Delta P=p_1+\mathbf{\xi }\cdot \nabla P$, $\Delta \rho ={\rho }_1+\mathbf{\xi }\cdot \nabla \rho$. Linearizing mass conservation, the Euler equations for an inviscid fluid and Poisson equation for Newtonian gravity givesThere is no equilibrium acceleration to multiply a perturbed mass density. The adiabatic process condition, with composition carried by the parcel, closes the system:These are the complete linear adiabatic stellar oscillation equations, including the perturbation of self-gravity. Neglect of heat exchange is appropriate when oscillation periods are short compared with relevant thermal relaxation times; it does not determine nonadiabatic excitation or damping.

For the pressure and buoyancy modes, separate angular dependence using spherical harmonics:The horizontal amplitude ${\xi }_h$ has dimensions of length. Angular differentiation gives $\nabla \cdot \mathbf{\xi }=[r^{-2}(r^2{\xi }_r{)}^{\prime }-\ell (\ell +1){\xi }_h/r]Y_{\ell }^m$, and tangential momentum gives ${\omega }^{2}r{\xi }_h=p_{\ell }/\rho +{\varphi }_{\ell }$. Define the adiabatic sound speed, stellar buoyancy frequency and Lamb frequency byThe adiabatic mass density relation becomes ${\rho }_{\ell }/\rho =p_{\ell }/(\rho c_s^2)+(N^2/g){\xi }_r$. Substitution produces a radial form of the full oscillation equations for nonzero $\omega$:The apparent $N^2/g$ factor is evaluated through its defining gradient at the centre rather than by dividing two zeros.

Regularity at the centre excludes singular solutions. At a free surface the Lagrangian pressure perturbation vanishes, $\Delta P=0$; outside the star the gravitational perturbation decays as $r^{-\ell -1}$. For a model with mass density tending to zero at its surface, continuity of ${\varphi }_{\ell }$ and its radial derivative gives ${\varphi }_{\ell }^{\prime }(R)=-(\ell +1){\varphi }_{\ell }(R)/R$. If the equilibrium mass density jumps to vacuum, include the displaced-surface mass sheet: the outward-minus-inward derivative jump is $4\pi G\rho (R){\xi }_r(R)$, so the interior condition is ${\varphi }_{\ell }^{\prime }(R)=-(\ell +1){\varphi }_{\ell }(R)/R-4\pi G\rho (R){\xi }_r(R)$. An atmospheric boundary condition can replace the ideal free surface.

These conditions make ${\omega }^2$ an eigenvalue, not an arbitrary local sound frequency. With conservative boundary conditions the adiabatic operator is self-adjoint, giving real ${\omega }^2$; negative values describe instability. The frequencies depend on $P(r),\rho (r),{\Gamma }_1(r)$, self-gravity, stratification and boundaries, as well as angular degree $\ell$ and radial order. In a spherical nonrotating star they are degenerate in $m$. For homologous equilibrium structures, their scale isthus the typical oscillation time measures inverse square root of mean mass density, while individual frequencies probe the interior adiabatic sound speed and stellar buoyancy frequency profiles.

For radial modes, write ${\xi }_r=r\eta$. Eliminating the pressure and gravitational perturbations gives the radial stellar pulsation equationThis is a Sturm-Liouville problem. Multiplication by $\eta$ and integration, with vanishing boundary terms, gives its Rayleigh quotientFor constant ${\Gamma }_1>4/3$, $P^{\prime }<0$ makes both numerator contributions nonnegative. At ${\Gamma }_1=4/3$, a homologous displacement is neutral; for constant ${\Gamma }_1<4/3$ the same trial displacement makes the quotient negative. For the uniform-density stellar model, $\eta =$ constant is an exact mode and ${\omega }^2=(3{\Gamma }_1-4)GM/R^3$. With varying ${\Gamma }_1$, the integral criterion, rather than a universal pointwise threshold, controls radial stability.

The Cowling approximation neglects ${\varphi }_{\ell }$ while retaining the equilibrium gravitational field. It is useful for short-wavelength modes, but is not needed for the full derivation above. In a locally slowly varying region, take both remaining amplitudes proportional to $e^{i\int k_{r}dr}$ and retain the leading derivative terms. ThenEliminating either amplitude gives the acoustic-gravity propagation relationPositive $k_r^2$ is oscillatory propagation; negative $k_r^2$ means an evanescent wave. The high-frequency branch, ${\omega }^2>S_{\ell }^2,N^2$, describes stellar acoustic modes, restored chiefly by compressibility and pressure. At frequencies well above $N$, this gives ${\omega }^2\simeq c_s^2[k_r^2+\ell (\ell +1)/r^2]$. The low-frequency propagating branch in stable stratification, ${\omega }^2<N^2,S_{\ell }^2$, instead describes stellar gravity modes, restored by buoyancy. At $\ell =0$ there is no such nonradial gravity-wave cavity.

A stellar acoustic mode is trapped between an inner turning point near $\omega =S_{\ell }$ and an outer reflecting region. Low-$\ell$ modes penetrate deeply; radial modes reach the centre. Higher-degree modes turn farther out. Standing waves require the WKB quantization conditionwhere the phase $\alpha$ depends on the central or turning-point behaviour and surface reflection. For high radial order and small degree, the leading acoustic travel-time result is the large frequency separationThe $\ell /2$ term is the leading central angular phase shift; smaller frequency separations depend on detailed interior gradients. The near-surface mass density stratification sets the acoustic cutoff frequency. In a plane-parallel isothermal atmosphere with mass density scale height $H_{\rho }$, ${\omega }_{\mathrm{ac}}\simeq c_s/(2H_{\rho })$. Modes below this cutoff can reflect and form a cavity; waves sufficiently above it escape and need an outgoing-wave boundary condition. Nonadiabatic damping, driving and rotation alter real-star mode properties, but the adiabatic frequency problem isolates their dependence on the equilibrium structure.

Use the inviscid, nonrotating Boussinesq approximation, with a stable background mass density $\bar{\rho }(z)$ and reference density ${\rho }_{\ast }$. Write $B=-g{\rho }^{\prime }/{\rho }_{\ast }$ for the perturbation buoyancy and $\pi =p^{\prime }/{\rho }_{\ast }$ for the kinematic pressure. The buoyancy frequency is $N_0^2=-(g/{\rho }_{\ast }){\bar{\rho }}_z>0$. Dropping products of perturbations in the Boussinesq equations gives the Linearized Boussinesq equationsThe last equation expresses incompressible flow; the second follows by advecting the background mass density gradient. Neglecting rotation and viscosity is part of this internal gravity wave model.

For a plane internal gravity wave proportional to $\exp [i({\mathbf{k}}_h\cdot {\mathbf{x}}_h+mz-\omega t)]$, put $\kappa =|{\mathbf{k}}_h|>0$ and $K^2={\kappa }^2+m^2$. Eliminating the horizontal velocity, pressure and buoyancy from the linear equations givesFor example, the horizontal momentum and continuity equations give $\pi =-\omega mw/{\kappa }^2$; substituting $B=-i{N}_0^{2}w/\omega$ into vertical momentum yields the displayed dispersion relation. Thus an IGW has $0<|\omega |\le N_0$ in this model.

On the positive-frequency branch, the phase velocity normal to a constant-phase plane and the group velocity areConsequently $({\mathbf{k}}_h,m)\cdot {\mathbf{c}}_g=0$, and phase and group velocity are perpendicular. Equivalently, the dispersion relation is homogeneous of degree zero in the wave vector, so differentiating with respect to its scale proves the same orthogonality. Energy travels with the group velocity, along the phase planes. At the degenerate limit $m=0$, $\omega =N_0$ and the group velocity vanishes; the orthogonality statement then has this limiting interpretation.

The Boussinesq approximation requires $\delta =(\rho -{\rho }_0)/{\rho }_0=\varphi ({\rho }_p/{\rho }_0-1)\ll 1$. The reference mass density ${\rho }_0$ can then be used in inertia, while the small density excess is retained in buoyancy. Dilute particle volume fraction $\varphi \ll 1$ alone does not suffice if ${\rho }_p/{\rho }_0$ is exceptionally large. Define $G=g({\rho }_p-{\rho }_0)/{\rho }_0$, so the reduced gravity is $g^{\prime }=G\varphi$.

The shallow water equations require depth small compared with the horizontal evolution scale, weak vertical acceleration and approximately hydrostatic pressure. Here the ambient is deep and quiescent, so its leading contribution is a hydrostatic reference pressure; the excess pressure in the layer is ${\rho }_{0}G\varphi (h-z)$. The model also needs nearly uniform streamwise velocity across the section and the stated absence of secondary circulation. The advancing gravity current nose is a separate region where these assumptions can fail.

A nearly uniform particle volume fraction can be maintained by turbulent mixing with eddy diffusivity $K_z$, provided $h^2/K_z\ll h/w$ and the mixing time is short compared with horizontal advection and bulk evolution. Equivalently, $wh/K_z\ll 1$ keeps the bulk settling-induced gradient small. The particles should respond rapidly enough to follow the mixing motions. A thin boundary layer adjacent to the absorbing bed or walls need not be well mixed.

Uniform bulk concentration does not imply zero sedimentation. The relative downward particle velocity still supplies a particle deposition flux of approximately $w\varphi$ to an absorbing boundary. Mixing redistributes the remaining particles and maintains the bulk concentration profile while its overall level decreases; it need not suspend every particle indefinitely. Deposition requires negligible resuspension in the model.

At a fixed height, first determine the stored quantities per unit vertical height and per unit length along the line source. With $\Delta \rho =\rho -{\rho }_0$, the triangular-profile line plume has the useful profile integralThe geometric volume is $\mathcal{A}=2b$; the mass is $\mathcal{D}={\rho }_0{I}_0+\Delta \rho I_1=b(\rho +{\rho }_0)$. The upward momentum is $\mathcal{P}=W({\rho }_0{I}_1+\Delta \rho I_2)$, while the buoyancy force is $\mathcal{B}=-g\Delta \rho I_1$. ThusHere positive buoyancy corresponds to a rising light plume. If instead buoyancy is defined kinematically, divide $\mathcal{B}$ by ${\rho }_0$; the density-weighted convention here is the one matching the printed equations.

Multiply each local density by the vertical velocity to obtain the fluxes through a horizontal section, still per unit source length:In particular, the stored momentum equals the mass flux $Q$, and the buoyancy flux is $F=g({\rho }_{0}V-Q)$. The different triangular-profile factors in storage and flux must not be replaced by top-hat coefficients.

Apply volume conservation, mass conservation and vertical momentum flux balance to a slice of the triangular-profile line plume. Ambient ingestion contributes $2e$ of volume and $2{\rho }_{0}e$ of mass per unit height and source length. The ambient is quiescent, so it supplies no leading vertical momentum. The integrated driving force is the stored buoyancy $\mathcal{B}$. Using the quantities from part (a), the non-Boussinesq triangular-profile line-plume balances areTogether with Batchelor entrainment, these are three equations for $W,\rho ,b$. They neglect viscous boundary stresses and streamwise pressure-force corrections within the integral-plume approximation.

Multiplying volume conservation by $g{\rho }_0$ and subtracting $g$ times mass conservation cancels the ambient sources exactly, givingThis is buoyancy conservation in the homogeneous ambient; it is a consequence of the first two balances, not an additional independent equation.

Now take the Boussinesq approximation in inertia and entrainment, retaining the small density deficit in buoyancy. Then $Q={\rho }_{0}bW$, $M=(2/3){\rho }_{0}b{W}^2$ and $F=(2/3)gbW({\rho }_0-\rho )$. ConsequentlyThe mass source is $2\alpha {\rho }_{0}W=3\alpha {\rho }_{0}M/Q$. Substitution into the mass conservation, momentum and buoyancy balances gives the unsteady Boussinesq triangular-profile line-plume balancesThese use density-weighted mass and momentum fluxes. Removing ${\rho }_0$ from some flux definitions while retaining it in the source would mix incompatible conventions.

Let $z$ measure distance inward from the surface. While $a\ll R$, curvature is negligible and this is a planar two-phase Stefan problem. Use a common reference mass density $\rho$, specific heat capacity $c_p$, thermal conductivity $k$ and thermal diffusivity $\kappa =k/(\rho c_p)$; retain the mass density differences only in the subsequent buoyancy estimates. Define $\Delta T=T_m-T_s>0$, $\Theta =\overline{T}-T_m\ge 0$, and $\vartheta =\Theta /\Delta T$. The paper's parameter $S=L/(c_p\Delta T)$ is the latent-to-sensible heat ratio, reciprocal to the more usual Stefan number.

For superheated two-phase Stefan similarity, solving $T_t=\kappa T_{zz}$ with $T(0)=T_s$, $T(a)=T_m$ and $T(\infty )=\overline{T}$ gives the Neumann solution of the Stefan problem. With $\eta =z/(2\sqrt{\kappa t})$,The solid expression holds for $0<z<a$ and the liquid expression for $z>a$. The interface gradients inserted in the Stefan condition determine the positive constant $\lambda$:The right-hand bracket decreases from infinity and crosses zero, while the equivalent left-hand side $S\sqrt{\pi }\lambda e^{{\lambda }^2}$ increases, proving uniqueness. The liquid-side flux subtracts from the roof cooling: it is heat brought to the interface by the superheated interior.

For large $S$, $\lambda \ll 1$. Keeping the superheat term while using $\mathrm{erf}\lambda \simeq 2\lambda /\sqrt{\pi }$ and $\mathrm{erfc}\lambda \simeq 1$ yields $S{\lambda }^2+(\vartheta /\sqrt{\pi })\lambda \simeq 1/2$. Hence the small-interface-parameter growth with liquid superheat approximation isThe solid temperature is approximately linear across its thin crust, while the liquid thermal boundary layer has thickness of order $\sqrt{\kappa t}$. For fixed superheat and $S\gg 1$, $A\simeq \sqrt{2/S}$; for $\vartheta \gg \sqrt{S}$, $A\simeq \sqrt{\pi }/\vartheta$. This retains the strong suppression of growth by a very hot interior.

To estimate overturn, take local gravitational acceleration $g$ at the roof. The liquid's mass density contrast is $\Delta {\rho }_l={\rho }_f\alpha \Theta$, while the solid contrast is $\Delta {\rho }_s={\rho }_s-{\rho }_f$. A liquid thermal layer of thickness $\delta$ has Rayleigh number $\Delta {\rho }_{l}g{\delta }^3/({\mu }_l\kappa )$. As $\delta \sim \sqrt{\kappa t}$, it reaches its critical value after a time scaling asIn viscous overturn of a dense crust, an unstably stratified viscous crust has a Rayleigh-Taylor instability growth time of order ${\mu }_s/(\Delta {\rho }_{s}ga)$. Setting this equal to its age, with $a=A\sqrt{\kappa t}$, givesThese are formation-to-overturn times. A specified critical liquid Rayleigh number ${\mathrm{Ra}}_c$ multiplies ${\tau }_l$ by ${\mathrm{Ra}}_c^{2/3}$; logarithmic amplification and geometrical constants also affect the crust estimate. They are omitted for the requested order-of-magnitude scaling.