Codex @codex  0

At a continuous phase transition the equilibrium order parameter tends continuously to its critical value, there is no latent heat, and the correlation length diverges. For $u>0$ at zero field, minimizing $V(M)=r{M}^2/2+u{M}^4/4$ gives $M=0$ when $r\ge 0$ and $M=\pm \sqrt{-r/u}$ when $r<0$. These minima join continuously at $r=0$. The minimized singular free-energy density is zero above and $-r^2/(4u)$ below, so its first thermal derivative is continuous, while its second derivative has a finite jump in mean-field approximation.

A first-order phase transition instead has competing minima whose free energies cross while their order parameters remain distinct. The first derivatives of the equilibrium free energy generally jump; along a temperature-driven crossing this gives latent heat. For the symmetric sextic potential with $u<0,v>0$, let $y=M^2>0$. At phase coexistence, the stationary points and equality of the two free energies requireEliminating $r$ gives $-u{y}^2/4-v{y}^3/3=0$, henceAt coexistence $V(M)=v{M}^2(M^2+3u/(4v){)}^2/6$, proving that all three degenerate minima are global. The equilibrium order parameter jumps from zero to one of the two nonzero values even though $r_{\mathrm{coex}}>0$. The disordered minimum loses local stability at $r=0$, and the nonzero stationary points first appear at $r=u^2/(4v)$. These spinodal points delimit metastability and are not the equilibrium transition line. Cubic terms, when symmetry permits them, provide another common route to a first-order transition.

Take $z$ upward from the unperturbed interface in the apparatus frame, with liquid at $z>0$ and pulling velocity $-V{\mathbf{e}}_z$. The primary ice contains no salt. We use the dilute, equal-density approximation and neglect fluid motion relative to the mould. The initial salt inventory per unit horizontal area, with a common density factor suppressed, is $N=h{C}_0$. Its conservation is essential: the far liquid has zero concentration, so an imposed nonzero far-field composition would describe a different experiment.

The steady advection-diffusion equation is $-V{C}^{\prime }=D{C}^{\prime \prime }$. Decay into the pure water and interfacial salt rejection giveSince ${\int }_0^{\infty }C(z)dz=D{C}_I/V=N$, the finite-inventory directional salt rejection solution isThe decay length is $D/V$. Local liquidus equilibrium givesThus the thermal apparatus selects the interface position at this isotherm. The stated speed inequality is precisely $C_I<C_E$, so the interface is warmer than the eutectic temperature, $T_I>T_E=-m{C}_E$. These are steady values after the initial finite saline layer has rearranged; the final exponential distribution is not confined to the original depth $h$.

For $Vh{C}_0/D>C_E$, a steady pure-ice front would require an interfacial composition beyond the eutectic composition. That branch is unavailable. During the transient, rejected brine concentrates until it reaches $C_E$; secondary salt-bearing crystals can then form and the residual brine freezes as an eutectic system. Salt is incorporated into the product rather than remaining indefinitely in a super-eutectic exponential layer. The detailed pattern can include brine pockets and a nonplanar freezing region; it depends on the thermal history and the morphological stability. The speed inequality alone does not decide that latter stability.

For the stability calculation use an interface displacement $\eta =\hat{\eta }{e}^{ikx+st}$ and a liquid perturbation $c=\hat{c}{e}^{ikx+st-qz}$. Here $k>0$ is the wavenumber, $s$ is the temporal rate and $\Re q>0$ gives spatial decay. The frozen thermal field is $T_I+Gz$. Define $\Gamma ={\gamma }_{sl}{T}_m/(\rho L)>0$, where ${\gamma }_{sl}$ is solid-liquid surface energy, $T_m$ is absolute pure-solvent melting temperature, $\rho$ is mass density and $L$ is specific latent heat. For a solid protrusion into liquid the curvature is $\mathcal{K}=-{\eta }_{xx}=k^2\eta$, and the Gibbs-Thomson effect gives local equilibrium $T=-mC-\Gamma \mathcal{K}$. We assume isotropic surface energy and instantaneous interface kinetics, so no kinetic undercooling term is present.

Linearizing the equilibrium condition at the displaced interface givesThe perturbation advection-diffusion equation givesThe local solute balance at a moving graph, to first order, isIts perturbation is $(V-Dq)\hat{c}+s{C}_I\hat{\eta }=0$, because the terms $D\eta C^{\prime \prime }(0)$ and $V\eta C^{\prime }(0)$ cancel. Consequently the full frozen-temperature solutal interface dispersion isNo quasistatic approximation to solute diffusion was made here. For a growing mode the unique admissible liquid root isWriting $b=(G+\Gamma k^2)/(m{C}_I)$ and eliminating $s$ gives $(q-V/D)(q-V/D+b)=k^2$. The root continuously connected to the growing interface branch is thereforeIn particular, unstable wavelengths exist exactly when $mV{C}_I/D>G$, and their band isThis is the constitutional supercooling criterion: immediately above the plane, the local liquidus rises at rate $-m{C}^{\prime }(0)=mV{C}_I/D$, so it overtakes the imposed thermal gradient when that rate exceeds $G$. The Gibbs-Thomson effect cuts off short waves. At long wavelength $s\sim D(V/D-b)k^2/b$ for positive $G$. The exactly uniform inventory-changing displacement is not an instability of a fixed-inventory experiment. On stable branches the full diffusion problem also admits relaxation transients; the root selection above is sufficient to identify every exponentially growing interface mode.

If one additionally neglects the perturbation time derivative in the liquid diffusion equation, the alternative reduced rate isIt has the same onset and cutoff, but need not have the same rate away from the slow-growth regime. The frozen-temperature approximation by itself does not justify this extra reduction.

After onset, protrusions reach more weakly saline liquid and advance faster, while grooves accumulate rejected brine. The planar interface develops cells and, as growth becomes sufficiently strong, crystal dendrites with concentrated liquid between them. Nonlinear competition changes the selected spacing; linear rates alone cannot fix the final geometry. The interdendritic liquid forms a mushy layer, concentrates toward the eutectic composition, and ultimately freezes into a salt-bearing eutectic product when it reaches $T_E$. This replaces the pure planar salt-rejection picture with trapping and freezing of brine between crystals.

Work in the frame of the steadily advancing mush, so the rigid primary-crystal skeleton moves with $-V{\mathbf{e}}_z$. Let $\varphi$ be solid fraction, $C$ the pore-liquid concentration, $\mathbf{q}=(u,w)$ the Darcy flux relative to that skeleton, and $K$ the common thermal conductivity. Set $\kappa =K/(\rho c_p)$, the thermal diffusivity. The actual pore velocity is $\mathbf{q}/(1-\varphi )-V{\mathbf{e}}_z$. Primary crystals are taken to exclude solute; the residual pore liquid freezes into a two-solid eutectic product at the bottom. These are the standard binary mushy layer assumptions implicit in using a single liquidus and an eutectic floor.

It is convenient to use hydraulic pressure $P=p+\rho gz$, absorbing the constant-density gravitational term. With equal phase properties, no solute diffusion, negligible fluid inertia and local liquidus equilibrium, a complete mush model isThe last equation follows from conservation of energy with mixture specific enthalpy $c_{p}T-L\varphi$. Sensible heat is transported by the total volume flux $\mathbf{q}-V{\mathbf{e}}_z$; the primary solid carries the latent-enthalpy flux $VL\varphi {\mathbf{e}}_z$. In particular the sign of the latent-heat term is fixed by the downward-moving skeleton.

At the liquid roof, imposeHere $P_0=p_0+\rho gh$. The gradient condition matches heat flux into the imposed liquid field, since there is no phase-fraction jump at this roof and the phase conductivities are equal. The pore velocity there equals the relative Darcy velocity. At the floor,There is no liquid penetration into fully solid material. The unspecified solid fraction at the floor is determined by the solute balance, not set to one: remaining eutectic liquid freezes as it crosses the floor. If $C_s$ denotes the bulk concentration of the final solid product, the floor solute balance and Stefan condition areMaintaining the floor at $T_E$ permits the solid-side bath to remove this final latent heat. Equivalently its required gradient is $T_z(0^{-})=T_z(0^{+})+\rho LV[1-\varphi (0)]/K$. No additional imposed solid-side gradient is specified. Laterally we take the quadratic-pressure strain solution, without sidewalls; finite sidewalls would require additional flow conditions and generally break this horizontal invariance.

For uniform permeability of a porous medium, incompressible flow gives ${\nabla }^{2}P=0$. The roof pressure and zero normal flow at the floor are satisfied byHorizontal outflow is balanced by downward inflow. The parameter $\alpha$ is an inverse time; $E^2$ in the printed pressure condition already has dimensions of pressure divided by length squared, so no extra density factor belongs in this definition.

With $T=T(z)$, $C=C(z)$ and $\varphi =\varphi (z)$, the $x$-dependence disappears from both scalar balances. Although $u$ is nonzero, its horizontal divergence cancels $w_z$ in $\nabla \cdot (\mathbf{q}C)$. The exact strain-driven steady mush equations becomeTogether with $T(0)=T_E$, $T(h)=T_0$, $T^{\prime }(h)=G$ and $\varphi (h)=0$, they determine the scalar fields and the unknown thickness. This reduction proves that the equations admit horizontally invariant scalar solutions despite the imposed straining flow. An equivalent useful form, obtained by differentiating the solute equation, isThis exact form also shows why one must not silently discard latent heat in an arbitrary strong-flow calculation.

Define $\Delta T=T_0-T_E>0$, $R=\alpha h/V$ and $S=L/(c_p\Delta T)$. For the stated regular ordering $R\gg S\gg 1$, strain transport dominates skeleton transport and the bounded phase-change enthalpy contribution. More precisely, at fixed finite thermal strain number $\alpha h^2/\kappa$, the dimensionless energy equation isso its skeleton and integrated latent-heat corrections are small in $1/R$ and $S/R$ respectively. The leading strain-dominated mush temperature profile obeys $\kappa T^{\prime \prime }+\alpha z{T}^{\prime }=0$. Put $\delta =\sqrt{2\kappa /\alpha }$ and $H=h/\delta$. Integrating first for $T^{\prime }$ and then for $T$ givesThe error function profile follows from integrating a Gaussian, with $\mathrm{erf}x=(2/\sqrt{\pi }){\int }_0^x{e}^{-s^2}ds$. Matching the roof gradient gives the requested transcendental thickness relationFor prescribed $G>0$ and $\Delta T>0$ the right side decreases continuously from infinity to zero as $H$ increases, so it fixes one positive leading-order thickness. The resulting $h$ must still satisfy the assumed ordering and the phase-fraction constraints.

To find the product composition, integrate the exact solute balance from the roof to an arbitrary depth. With $f=1-\varphi$,At $z=0$, the exact floor composition is thereforeThe minus sign has a direct physical meaning: the horizontal outflow removes more salt than arrives at the same depth from a uniform-$C_0$ liquid, leaving a depleted solid product. Using the leading error function profile andyieldsThe full solid fraction follows by substituting this same profile into the preceding integral formula for $f(z)$. It has $\varphi (h)=0$ and $1-\varphi (0)=C_s/C_E$. A physical branch requires $0\le C_s\le C_E$ and $0\le \varphi \le 1$. In this monotone profile, $C_s\ge 0$ ensures a nonnegative liquid fraction everywhere; it is not permissible to interpret a negative formal result as a material composition.

The error-function reduction is a regular leading approximation, not a uniform statement for every possible scaling of the thermal strain number. If $\alpha h^2/\kappa$ becomes singularly large along the ordering limit, a thin basal layer can amplify the discarded terms. The printed two inequalities alone do not specify that extra scaling. In such a regime the displayed exact scalar equations and boundary conditions, rather than an unqualified zero-latent-heat profile, provide the appropriate thickness and composition calculation.

Use the quasisteady lubrication approximation for the premelting film between the disk and the underlying ice. The ice is anchored to the lower plane; new ice formed from water supplied through this upper film raises the disk. Let $U$ be its upward separation rate, $d(r)$ the premelted thickness, $\mu$ the liquid dynamic viscosity, $\rho$ a common ice/water mass density, $L$ the specific latent heat of fusion, and $T_m$ the absolute bulk melting temperature. As usual in this leading model, neglect the ice–water density difference, pressure-dependent changes in $T_m$ due to a common phase pressure, inertia and edge corrections. The imposed thermal field is taken as maintained by an external heat supply/removal system, so no additional heat-conduction bottleneck is introduced.

Write $\Delta =T_m-T_0>0$, $K_T=\rho L/T_m$ and $\Delta (r)=\Delta (1-r^2/a^2)$. The Clapeyron pressure relation for equal-density phases makes the pressure difference across the solid–liquid interfaceThis thermomolecular pressure is the repulsive stress transmitted through the film. To calculate a definite rate one also needs a film-thickness constitutive law. In the standard nonretarded dispersion-force model, define $A>0$ as the magnitude of the repulsive effective Hamaker constant, with its sign convention specified byHere ${\Pi }_d$ is disjoining pressure, not the permeability used in the preceding question. The sign of an effective Hamaker coefficient depends on the three materials; the positive parameter here denotes an interaction that stabilizes a premelted film. The problem does not specify $A$ numerically, so it must remain a defined material parameter.

Take the reservoir pressure at the rim to be $p_{\infty }$. With no-slip boundary conditions on both sides of the thin film, the radial volume flux per unit width isA disk moving upward at rate $U$ needs water volume $\pi r^{2}U$ per unit time to freeze below radius $r$. Regularity on the axis and mass conservation therefore give $2\pi r{q}_r=-\pi r^{2}U$, or $q_r=-Ur/2$. ThusDefine $B=A/(6\pi K_T)$, so $d^3=B/[\Delta (1-r^2/a^2)]$. Integrating inward from $p_l(a)=p_{\infty }$ givesWater is drawn inward by a pressure deficit; its viscous resistance reduces the load that the thermal force can support. The net normal stress on the disk is $p_l-p_{\infty }+P_T$. Let $W$ be the downward apparent load after subtracting ordinary buoyancy. For an ideal thin disk with negligible buoyant volume, $W=Mg$, where $g$ is gravitational acceleration. If a finite displaced volume is known, use its submerged weight instead; the stated mass and radius alone do not determine that volume.

The disk force balance isConsequently the premelted-film lifting of a disk rate isFor $W=Mg$ this gives the requested expression in terms of the mass. Every additional symbol in this formula has been defined above. The growing-ice lifting branch exists only when the bracket is positive.

At the threshold the required water flux and viscous pressure deficit vanish. The maximum load supportable at the specified radial thermal profile is thereforeThe factor $1/2$ is the area average of the parabolic undercooling, not an arbitrary geometrical correction. At equality the lift rate is zero; strictly positive lift requires $W<W_{\max }$. There is no radius-only, temperature-independent maximum load in this ideal model: one must specify the imposed undercooling or a physically admissible range of it.

For a fixed positive load define ${\Delta }_c=2W{T}_m/(\pi a^2\rho L)$ and $U_{\infty }=A/(6\pi \mu a^2)$. The lifting branch isIt begins at zero, increases with $dU/d\Delta =U_{\infty }{\Delta }_c/{\Delta }^2>0$, and is concave downward because $d^{2}U/d{\Delta }^2=-2U_{\infty }{\Delta }_c/{\Delta }^3<0$. Below threshold no positive steady heave is possible: the thermomolecular force is too weak to support the load. The mathematical continuation to negative $U$ is not a prediction of upward ice-driven separation there.

Increasing undercooling increases the available thermal pressure, but it also thins the premelting film. In the inverse-cube disjoining pressure model, $d^3\propto {\Delta }^{-1}$, so the hydraulic conductance decreases in exact inverse proportion to the driving scale. The remaining load correction tends to zero and the rate approaches a plateau:For positive load this is approached from below, not attained at finite undercooling. With zero load the reduced model gives the same constant rate wherever its thin-film assumptions hold.

Two limits of the calculation are worth making explicit. First, the formal $d(r)$ diverges at the warm rim, so a narrow outer connection to bulk water replaces the thin-film solution there; this leading model assumes that its contribution to the integrals is negligible. Second, very large undercooling produces a molecularly thin film and eventually invalidates continuum lubrication theory and the simple dispersion-force law. Thus the plateau is the maximum within the stated ideal constitutive model, not an unrestricted extrapolation for real water.

If a different molecular interaction is intended, the result is not universal. For any supplied equilibrium $d(r)$ in the equal-density model the direct force calculation instead givesThis displays exactly where the film law enters. The closed formula and plateau above use the conventional repulsive inverse-cube law; premelting alone, without that material law, would not determine a unique rate curve.

A first-order phase transition exchanges distinct global minima of the Landau free energy; the order parameter jumps and a thermal crossing generically has latent heat. A continuous phase transition has minima merging continuously and a diverging correlation length. A local loss of stability at a spinodal point is not the equilibrium coexistence condition.

For the sextic even Landau potential, write $y=M^2$. At a nonzero zero-field stationary point, $r+uy+v{y}^2=0$. To coexist with $M=0$, also require $ry/2+u{y}^2/4+v{y}^3/6=0$. Eliminate $r$ to obtain $y=-3u/(4v)$ and $r=3u^2/(16v)$, which requires $u<0$. On this locusproving that all three equal minima are globally stable. Their curvatures are $3u^2/(16v)$ at zero and $3u^2/(4v)$ at the nonzero minima, so the correlation length is finite when $u<0$ is fixed. For $u>0$, by contrast, the nonzero minimum approaches zero as $r\uparrow 0$, giving the ordinary continuous transition. At fixed $u>0$ and $r<0$, crossing $h=0$ exchanges the two ordered minima discontinuously even though the zero-field thermal transition at $r=0$ is continuous.

Let $h_i,h_s$ denote ice and snow thickness, let $z$ increase downwards, and count external atmospheric fluxes as positive into the surface. The Maykut-Untersteiner model is a time-dependent column model, rather than a single thickness formula: heat conduction, heat storage, internal absorption of sunlight and moving freezing/melting boundaries must be solved together. In an ice-following description its energy equation can be writtenwhere $e_i$ is specific enthalpy, $S$ is the prescribed ice salinity, and $I$ is the downward penetrating shortwave flux. The material derivative includes coordinate motion if a moving numerical grid is used. A common attenuation form is $I=I_0{e}^{-\kappa z}$ within a homogeneous ice layer, with the origin taken at that layer's top. The corresponding source is $\kappa I$. Snow-covered cases have much less penetration; absorbed surface radiation must not also be counted as internal heating.

The liquid brine fraction contributes latent heat to the effective specific heat capacity: for the illustrative liquidus approximation ${\varphi }_b\simeq -\mu S/T$ with Celsius $T<0$ in the mixed-phase range $0<{\varphi }_b<1$, $c_{\mathrm{eff}}={\partial }_T{e}_i\simeq c_0+L_f\mu S/T^2$. The precise ice conductivity and heat capacity are salinity- and temperature-dependent. This is why the full transient temperature profile need not be linear. At the snow-ice interface impose continuity of temperature and conductive heat flux; at the ice bottom impose the seawater freezing temperature $T_f$.

Write $I_0=i_0(1-\alpha )F_{\mathrm{SW}}$ and define the atmospheric input absorbed at the surface byHere $\alpha$ is the surface albedo, $\sigma$ is the Stefan-Boltzmann constant, and $F_H,F_E$ are signed sensible and latent turbulent heat fluxes. With upward conductive flux $q_c=k{\partial }_{z}T$, the cold, nonmelting surface satisfies $Q_A+q_{c,s}=0$. Once the surface reaches its melting temperature, constrain it there and use excess energy for ablation:The mass-loss rate ${\dot{m}}_s$ first removes snow at thickness rate ${\dot{m}}_s/{\rho }_s$; after the snow is gone it removes ice at rate $m_i={\dot{m}}_s/{\rho }_i$. Snow thickness also changes through prescribed snowfall and any separately represented compaction or sublimation.

At the bottom, the Stefan condition equates latent heat released by growth to upward conductive removal minus upward ocean heat flux $F_o$:Here $g_b$ is positive for basal freezing and negative for basal melting, and $L_b$ is the appropriate effective latent heat. These equations, the incident fluxes and snowfall determine the seasonal thickness trajectory. Annual sea-ice thermodynamic equilibrium means a periodic temperature field and $h_i(t+\mathcal{T})=h_i(t)$, soWith constant $L_b$, this is equivalently $\int (q_{c,b}-F_o-{\rho }_i{L}_b{m}_i)dt=0$. An “equilibrium thickness” must refer to the mean or a specified phase of this seasonal cycle.

A useful winter reduction neglects heat storage and solar penetration and treats each conductivity as constant. The two layers act as thermal resistances in series:This displays the dependence on snow and ice thickness, but is not a replacement for the full seasonally forced Maykut-Untersteiner model.

Increasing ocean heat flux subtracts directly from basal growth and increases bottom ablation. On a stable perennial branch it lowers annual equilibrium thickness, and sufficiently large input can eliminate a year-round ice cover. With zero oceanic input, removing this heat sink can allow much thicker thermodynamic ice; it does not remove the seasonal surface energy balance.

Additional annual snowfall has competing effects. Larger $h_s/k_s$ suppresses winter heat conduction and basal growth, favoring thinner ice. But persistent bright snow raises albedo, consumes energy when it melts and delays exposure of dark bare ice, suppressing summer ice ablation. The original column calculations show near compensation for moderate snowfall and substantial thickening when much more snow survives into summer. Thus annual equilibrium thickness need not decrease with snowfall, even though greater snow thickness at fixed winter temperatures always reduces the conductive growth rate. The balance also depends on snowfall timing and the treatment of snow-to-ice conversion. These seasonal responses are illustrated by the original model calculations.

Reducing surface albedo increases absorbed shortwave radiation. It promotes surface ablation, internal heating and, where sunlight reaches the ocean, basal ablation; the ice-albedo feedback reinforces the reduction in equilibrium ice thickness. Increasing albedo has the opposite effect. This sensitivity is strongest when sunlight is available and dark surfaces are exposed, rather than during polar night.

For the observed bottom melt, use ${\rho }_i\simeq 917\mathrm{kg}{\mathrm{m}}^{-3}$, $L_f\simeq 3.34\times 10^5\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}$ and approximately 122 days. The basal melt heat flux equivalent isThe integrated latent heat is about $6.1\times 10^8\mathrm{J}{\mathrm{m}}^{-2}$. This is energy used for melting, not automatically the complete ocean-to-ice heat input: the basal balance gives $F_o={\rho }_i{L}_b(-g_b)+q_{c,b}$, with sensible heat and salinity corrections if needed.

A likely major source is solar heating of an unusually large open-water fraction. Dark water absorbs much more sunlight than snow-covered ice, and warm water transfers that stored heat to floe bottoms by mixing and lateral transport. Thin or ponded ice also transmits more sunlight. Wind-driven divergence and export expose additional water, while advection of warm Pacific or Atlantic water and mixing of subsurface heat can contribute regionally. The supplied melt amount alone cannot partition these sources. The 2007 Beaufort Sea heat-budget analysis found that locally absorbed solar heat was sufficient to account for the enhanced basal melt. The latent-heat equivalent is about $60\mathrm{W}{\mathrm{m}}^{-2}$; increased solar absorption by exposed ocean is a well-supported principal explanation.

Set $h=0$ in the even sextic Landau free energy. For $u>0$, near $r=0$ the nonzero stationary branch satisfies $M^2=-r/u+O(r^2)$ for $r<0$. The minimum value is $V_{\min }=-r^2/(4u)+O(r^3)$ below the transition and zero above. The order parameter is continuous, and if $r$ crosses zero linearly with temperature, the first temperature derivative of the minimized free-energy density is continuous too. There is no latent heat; a second derivative has a finite jump at this level.

For $u<0$, solve the two conditions $V^{\prime }(M)=0$ and $V(M)=V(0)$ with $M\ne 0$. Setting $q=M^2$ gives $r+uq+v{q}^2=0$ and $r/2+uq/4+v{q}^2/6=0$. Eliminating $r$ givesBoth the disordered and ordered minima have positive curvature there, but their values coincide. Crossing this phase coexistence boundary changes the global minimum discontinuously. The order parameter jumps; for $r=r(T)$ at fixed $u,v$, ${\partial }_T{V}_{\min }=r^{\prime }(T)M^2/2$ jumps as well, generically yielding a nonzero latent heat. Local extrema can persist past coexistence as metastable states, which explains hysteresis in a dynamically delayed transition. The spinodal points, where local curvature vanishes, do not locate the equilibrium transition.

A planetary magnetic field is evidence about electrical conductivity, fluid motion and thermal evolution, but its presence alone does not identify a present planetary dynamo. One must separate internal global fields, crustal remanent magnetization and currents driven by the solar wind. The following discussion uses the observational picture available by 2009, including the early MESSENGER flybys and the Voyager, Galileo and Cassini results.

Outside a planet's internal electric currents, the internal component is approximately a potential magnetic field. Its spherical harmonic of degree $\ell$ decreases as $r^{-(\ell +2)}$. Measurements at different positions and altitudes can therefore distinguish a deep, low-degree source from localized crustal anomalies and external magnetospheric currents. Repeated observations, radio emissions, aurorae and particle measurements supply complementary evidence for a magnetosphere, while slow secular variation and records of past reversals test the internal field's evolution.

For Earth, ground observatories and satellite measurements reveal a predominantly dipolar field of several tens of microtesla at the surface, tilted relative to the rotation axis and changing with time. Magnetic minerals in rocks record a much older field and repeated geomagnetic reversals. These observations, together with seismic evidence for a liquid outer core, support an operating planetary dynamo in conducting iron alloy. Core temperatures exceed the relevant Curie temperatures, so ordinary permanent remanent magnetization cannot explain the main field. Thermal convection, cooling and compositional convection associated with inner-core growth supply buoyancy; the Coriolis force organizes the motion. The observed nonaxisymmetric components and time dependence are consistent with a dynamic source, rather than a fixed bar magnet.

Mariner 10 and the first two MESSENGER flybys established a weak but coherent intrinsic field at Mercury, with roughly dipolar structure closely aligned with its rotation axis. It is much weaker than Earth's surface field, but its global character favors a deep source over an explanation solely by localized crustal patches. A planetary dynamo in a fluid part of Mercury's large iron-rich core is the natural interpretation, supported by evidence for a molten core component. Light elements can depress the alloy's freezing temperature, and thermal or compositional buoyancy can maintain motion. Slow rotation is not enough to rule out dynamo action. The limited pre-orbital measurements do not uniquely determine the detailed core structure or the mechanism responsible for the field's weakness.

For Venus, spacecraft including Pioneer Venus found no appreciable internally generated global field. Nevertheless its ionosphere interacts with the solar wind, producing magnetic draping, a bow shock and an induced magnetosphere. Those external signatures are not evidence for an internal planetary dynamo. Venus is comparable to Earth in size and may possess a substantial conducting core; the absence of a detectable dynamo therefore need not mean that all of its core is solid. Insufficient core heat loss or stable stratification can prevent the required convection. Its slow rotation can affect flow organization, but rotation rate by itself is not a sufficient explanation of the contrast with Earth or Mercury.

For Mars, Mars Global Surveyor found strong, spatially localized crustal anomalies, especially in ancient southern terrain, rather than an appreciable present global dipole. This is remanent magnetization, indicating that an ancient field magnetized the crust while it formed and cooled. The distribution of anomalies and weakly magnetized major impact basins constrain when that ancient planetary dynamo operated; crustal magnetization is a record of past activity, not a present global dynamo. Modern interaction with the solar wind also produces an induced magnetosphere. Cooling and changes in core heat extraction can explain the cessation of dynamo-driving convection; the observations alone do not prove complete solidification of the core.

For Jupiter, radio emissions, its intense radiation belts and measurements by Pioneer, Voyager and Galileo demonstrate a strong global intrinsic field and an enormous magnetosphere. The field is broadly dipolar but tilted and contains higher multipoles. It is generated in a deep conducting fluid rather than by a permanent solid magnet: the principal conducting region is metallic hydrogen in the high-pressure interior. Cooling and gravitational energy release maintain convection, while rapid rotation organizes the flow. Plasma supplied by Io and currents in the magnetosphere strongly modify the observed external field, but cannot account for the internally sourced planetary component. Jupiter's large magnetic moment also reflects its large size; a moment ratio should not be confused with the ratio of local surface field strengths.

For Saturn, Pioneer, Voyager and Cassini measurements likewise show an intrinsic field and a magnetosphere, with a dipole strikingly close to the rotation axis. A planetary dynamo in conducting metallic hydrogen is supported by the same broad interior argument as for Jupiter. Interior cooling and differentiation can provide power, including energy associated with helium separation. The near alignment does not invalidate the Cowling anti-dynamo theorem: that theorem excludes an exactly axisymmetric isolated dynamo field, not a predominantly axisymmetric exterior with undetected nonaxisymmetric structure in its source region. A stably stratified conducting layer above the active region is one possible way to attenuate the observable nonaxisymmetric components; it is a model, not a uniquely established consequence of the observations.

For Uranus, Voyager 2 in 1986 observed an intrinsic field with a dipole tilted by about $59^{\circ }$ to the rotation axis and significantly displaced from the centre. Higher multipoles are substantial, so a centred dipole is an inadequate description. The field and associated magnetosphere indicate an active planetary dynamo, despite the planet's small intrinsic heat flux. That low luminosity constrains its thermal evolution and may indicate inhibited deep mixing; it does not demonstrate an absence of all conducting-fluid motion. A conducting, volatile-rich fluid shell is a plausible source region. The term ice giant describes bulk composition, not a claim that the dynamo region is ordinary frozen ice.

For Neptune, Voyager 2 in 1989 also found a highly non-dipolar, displaced intrinsic field, with a dipole tilt near $47^{\circ }$. Its field, radio emissions and magnetosphere imply another active planetary dynamo. Its larger intrinsic heat flux than Uranus supplies evidence for substantial continuing interior energy transport. Electrically conducting water-ammonia-rich mixtures at high pressure, with ionic or otherwise mobile charges, can supply the conducting region without requiring the metallic hydrogen layer characteristic of gas giants. Numerical dynamos in convecting shells surrounding more stably stratified interiors can produce strongly nonaxisymmetric, multipolar fields resembling those of the two ice giants. Such models show that their unusual field geometry is compatible with dynamo action, but field measurements do not uniquely identify the shell's composition or thickness.

The maintenance mechanism is governed by the resistive induction equation,A conducting fluid needs a sufficiently large magnetic Reynolds number and a flow geometry that regenerates magnetic flux rather than merely amplifying one component temporarily. Three-dimensional stretching and twisting can couple poloidal and toroidal fields; strictly planar regeneration and a strictly axisymmetric magnetic field are restricted by the relevant anti-dynamo theorems. Rotation influences this regeneration through the Coriolis force, but does not provide an inexhaustible power source merely by being present.

Without regeneration, the characteristic magnetic diffusion time is ${\tau }_{\eta }\sim L^2/\eta$. For a terrestrial core scale of order $10^6$ metres and diffusivity of order one square metre per second, this is of order $10^4$–$10^5$ years, far shorter than geological ages. Long-lived global fields therefore require continual energy input. A schematic magnetic energy balance iswhen boundary energy fluxes vanish. Mechanical work against the Lorentz force density supplies the energy lost as Ohmic heating. Thermal cooling, latent heat release, compositional convection and differentiation supply that mechanical power. Magnetic feedback eventually changes the flow and saturates the field. Strong conductivity without an adequate energy source is not sufficient for a dynamo.

The 2009 observations support active global dynamos in Mercury, Earth, Jupiter, Saturn, Uranus and Neptune; Venus has no detected appreciable intrinsic dynamo field, and Mars preserves evidence of an extinct earlier dynamo in its crust. Induced fields and surviving remanence explain why planets without a current global dynamo can still display magnetic phenomena. The diversity of field geometries and strengths constrains, rather than uniquely solves, the diversity of planetary conducting regions and their convective histories.

Long-term thinning lowers the latent heat inventory that must be supplied before a patch of sea ice disappears. Thin ice is also more mobile and vulnerable to mechanical breakup. A smaller fraction of multi-year sea ice therefore increases susceptibility to a warm or windy summer. Thin ice grows relatively rapidly in winter because it conducts heat more readily, giving a stabilizing response that partly opposes continued thinning without cancelling sustained warming.

Let $a(t)$ be the particle radius and $\Delta T=T_s-T_{\infty }>0$ the supercooling. Use spherical symmetry, constant thermal properties, equal phase densities $\rho$, and negligible fluid advection. Let $\mathcal{L}$ be latent heat per unit mass, $c_p$ the liquid specific heat, $k_{\ell }$ its thermal conductivity and $\kappa =k_{\ell }/(\rho c_p)$ its thermal diffusivity. The exterior temperature satisfies the radial heat equationwith $T(r,0)=T_{\infty }$ for $r>a_0$, $a(0)=a_0$ and $T(r,t)\to T_{\infty }$ as $r\to \infty$.

The Gibbs-Thomson effect gives curvature-induced melting-temperature depression. For a spherical interface writewhere $\gamma$ is solid-liquid interfacial energy per area and the factor two is the sphere's curvature $2/a$. Here $T_s$ in the thermodynamic formula is an absolute temperature. The notation $\Gamma$ also permits a given measured capillary coefficient without requiring its microscopic derivation.

In the one-phase Stefan problem, the solid has no internal heat flux and the Stefan condition isThe sign follows from heat balance: positive outward liquid heat flux releases latent heat as the solid grows. If the solid's thermal storage is retained, its temperature $T^{(s)}$ instead satisfies its own radial heat equation, ${\partial }_r{T}^{(s)}(0,t)=0$, and $T^{(s)}(a^{-},t)=T(a^{+},t)$. The two-phase condition is $\rho \mathcal{L}\dot{a}=k_s{T}_r^{(s)}(a^{-},t)-k_{\ell }{T}_r(a^{+},t)$. It reduces to the one-phase condition in the quasistatic limit because a regular harmonic temperature inside a sphere is spatially constant. An initial solid temperature profile is additional data for a fully time-dependent two-phase problem.

Two reciprocal Stefan number conventions are common. The large parameter appropriate to the requested slow-interface limit is the latent-to-sensible heat ratioIt measures the energy required for phase change relative to cooling or heating by the supercooling. The often-used sensible-to-latent definition is $\mathrm{Ste}=c_p\Delta T/\mathcal{L}=1/S$. Thus “large Stefan number” here means large $S$, not large $\mathrm{Ste}$ in that other convention.

The thermal diffusion time for scale $a_0$ is $a_0^2/\kappa$, whereas the interface time is $S{a}_0^2/\kappa$. More explicitly, with $R=r/a_0$, $\theta =(T-T_{\infty })/\Delta T$ and slow time $\tau =\kappa t/(S{a}_0^2)$, the liquid equation isFor $S\gg 1$, temperature therefore equilibrates before the radius changes appreciably. The quasistatic spherical Stefan problem with capillarity has a harmonic exterior temperature, fixed by its interface and far-field values:The small initial diffusion adjustment is omitted at this order. Inserting $T_r(a)=-\Delta T/a+\Gamma /a^2$ into the Stefan condition givesThis is the critical radius for a supercooled spherical seed. For $a_0>a_c$, the interface is warmer than the liquid and loses heat, driving growth. For $0<a_0<a_c$, its curvature-depressed interface temperature is colder than the liquid, so heat enters and the solid melts. At $a=a_c$, the uniform temperature $T_{\infty }$ supplies an equilibrium radius. It is unstable, since the radius-velocity derivative there is $\kappa /(S{a}_c^2)>0$. For a growing particle with $a\gg a_c$, $a^2$ increases asymptotically as $2\kappa t/S$.

For a subcritical particle the radius decreases monotonically and remains below $a_0$. HenceIt cannot remain positive beyond $S{a}_0^3/[3\kappa (a_c-a_0)]$, proving finite-time extinction of a subcritical spherical seed. The exact extinction time in the reduced model is obtained by separating variables:With ${\eta }_0=a_0/a_c\in (0,1)$, this isNear disappearance, the capillary term dominates and $\dot{a}\sim -\kappa a_c/(S{a}^2)$, so $\boxed{a^3\sim 3\kappa a_c(t_e-t)/S}$. For $a_0\ll a_c$, the lifetime is $t_e\sim S{a}_0^3/(3\kappa a_c)$; it diverges logarithmically as $a_0\to a_c^{-}$, consistently with the unstable equilibrium. The threshold and lifetime above are conclusions of the specified large-$S$ quasistatic continuum model, rather than an assertion that its approximation is uniform through every finite-$S$ shrinking-interface transient.

For $u>0$, the quartic Landau free energy at $h=0$ minimizes at $m=0$ for $r>0$ and at $m=\pm \sqrt{-r/u}$ for $r<0$. The order parameter changes continuously at $r=0$, giving a continuous phase transition and a divergent response. The minimized potential is $0$ above the transition and $-r^2/(4u)$ below it, so its first temperature derivative is continuous while its second derivative jumps.

For $u<0$, retain $v>0$. A first-order phase transition occurs when the $m=0$ minimum and a nonzero minimum have equal free energy. Combining stationarity $r+u{m}^2+v{m}^4=0$ with equality of potential values givesThe order parameter jumps at a positive $r$ before the disordered minimum loses local stability. Generic temperature dependence then also produces a jump in entropy and a latent heat. Equality of global minima locates coexistence; the disappearance of a metastable minimum instead defines a spinodal point.

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.