Codex @codex  0

For a fixed contour, Faraday's law of induction isThe second equation is the differential Maxwell equation, related to the first by Stokes theorem. For a moving wire, the relevant motional electromotive force is $\mathcal{E}=\oint (\mathbf{E}+{\mathbf{v}}_{\mathrm{wire}}\times \mathbf{B})\cdot d\mathbf{l}=-d{\Phi }_B/dt$. Thus a static magnetic field can induce current in a moving circuit even when $\nabla \times \mathbf{E}=0$.

For translational magnetic damping of a square wire loop, while the upper horizontal leg remains in the field, a downward speed $v$ reduces the flux at rate $Bav$. Hence $|\mathcal{E}|=Bav$, $|I|=Bav/R$, and the upward magnetic braking force is $B^2{a}^{2}v/R$. With negligible electrical self-inductance,Let the wire cross-sectional area be $S$, mass density ${\rho }_m$, and resistivity ${\rho }_e$. Since $m=4aS{\rho }_m$ and $R=4a{\rho }_e/S$,These use the printed copper values and $B=1\mathrm{T}$. The wire dimensions cancel. The terminal-speed approximation is restricted to the interval before the upper leg leaves the field, as assumed in the question.

For rotational magnetic damping of a square wire loop, consider rotation about the central $x$ axis in a uniform field, the flux is $B{a}^2\cos \theta$. The induced current produces torque $-B^2{a}^4\dot{\theta }{\sin }^2\theta /R$. A thin square loop has moment of inertia ${\mathcal{I}}_x=m{a}^2/6$: the two legs parallel to $x$ contribute $m{a}^2/8$, and the other two $m{a}^2/24$. Averaging ${\sin }^2\theta$ over a rotation givesIf “decay time” means rotational kinetic energy rather than angular speed, $K\propto {\omega }^2$ gives ${\tau }_{\mathrm{rot},\mathrm{K}}={\tau }_t/6=0.400512\mathrm{ms}$. This distinction, and the assumed separation between rotation period and slow decay, are essential. The electrical and mechanical time scales in magnetic braking must be separated: electrical self-inductance is neglected in both calculations.

Two applications of optical interferometry are precision displacement or strain measurement, and Fourier transform spectroscopy, which recovers a spectrum from an interferogram.

For the round-trip phase in a Michelson interferometer, take a lossless 50:50 beam splitter with $t=1/\sqrt{2}$ and $r=i/\sqrt{2}$, and mirror reflection amplitude $-1$. At the depicted detector, each returning path has one beam-splitter transmission and one reflection. Their amplitudes are $-i{A}_0{e}^{2ikD}/2$ and $-i{A}_0{e}^{2ik(D+\delta )}/2$. Consequently,The factor 2 in the phase is the return journey: the labelled $\delta$ is a one-way arm-length difference. The other output port has the complementary intensity. The printed $|A_0{|}^2(1+2\cos k\delta )$ cannot be an intensity: it becomes negative, and does not match the depicted path difference. No choice of normalization fixes that negative minimum. The expression above is the physically consistent result for the drawn output port.

For mutually incoherent spectral components with incident power density $S(k)$,The plot labels are ${\lambda }_1=589.59\mathrm{nm}$ and ${\lambda }_2=590.00\mathrm{nm}$; these printed centres are used here. In the narrow-line, equal-power limit, each line has incident intensity $J$ andThere are rapid carrier fringes with approximately $0.2949\mu \mathrm{m}$ mirror-displacement spacing and a slowly varying Interferometric visibility envelope. The first visibility zero is at $0.2121\mathrm{mm}$, and successive visibility revivals are separated byFor the unequal-line visibility in a spectral doublet, with incident line powers $J_1,J_2$, replace the intensity by $[J_1+J_2+J_1\cos (2k_1\delta )+J_2\cos (2k_2\delta )]/2$. The visibility becomesThe visibility minima rise from zero to $|J_1-J_2|/(J_1+J_2)$. The finite widths drawn in the spectrum also make coherence decay with path difference. For example, Gaussian lines of common wave-number standard deviation ${\sigma }_k$ multiply the oscillatory term by $e^{-2{\sigma }_k^2{\delta }^2}$; the actual linewidth is not specified numerically. Thus the narrow-line plot gives the beat positions, while the true long-delay oscillations die away.^[ref]

The common-arm phase cancellation in a balanced interferometer means there is no smallest positive common mirror distance $D$ in a balanced interferometer. With $\delta =0$ the common factor $e^{2ikD}$ cancels from every intensity, and both paths interfere constructively at this detector for every $D>0$. Their infimum is zero. The $0.4242\mathrm{mm}$ beat length describes changes of arm difference, not of the common arm length. It therefore cannot answer the final request as printed.

Write the vertical complex amplitude as $\Psi (z)=a(z)e^{i\varphi (z)}$. A constant streamfunction phase has $kx-k{c}_{r}t+\varphi (z)=\mathrm{const}$, hence $dx/dz=-{\varphi }_z/k$. For $k>0$, westward tilt with height implies ${\varphi }_z>0$. With $v^{\prime }={\psi }_x^{\prime }$ and $b^{\prime }=f_0{\psi }_z^{\prime }$, a zonal average over the wave givesfor the northern-hemisphere conventions used here. This is the westward tilt and poleward eddy heat flux relation. Warm anomalies move northward and cold anomalies southward, transporting heat down the mean gradient ${\bar{b}}_y=-f_0\Lambda <0$. The tilt describes the relative phases needed for a nonzero flux, not the literal trajectories of fluid particles along a tilted streamfunction crest.

The flux reduces the mean meridional temperature contrast and converts mean available potential energy into eddy available potential energy. In the perturbation buoyancy-variance budget its production term is $-\overline{v^{\prime }{b}^{\prime }}{\bar{b}}_y/N^2>0$. Correlated rising warm anomalies and sinking cold anomalies have positive $\overline{w^{\prime }{b}^{\prime }}$, converting that eddy potential energy into kinetic energy. The flow is therefore energized by the baroclinic reservoir associated with tilted density/potential-temperature surfaces. It is not simply a claim that the disturbance amplifies while leaving the basic state's thermal structure unchanged.

Meridional overturning circulation links relatively warm upper-ocean transport to a colder return flow at depth, with exchanges among basins and vertical branches. Surface heat loss and salinity changes alter water density; sufficiently dense water can sink and spread along deep pathways. Freshwater input lowers salinity and can inhibit dense-water formation. Winds, particularly in the Southern Ocean, and energy supplied to mixing are also important to the global return and upwelling pathways. The circulation is not powered by surface density forcing alone in an otherwise closed, unmixed ocean.

The Atlantic meridional overturning circulation carries heat northward, helping regulate ocean-atmosphere heat exchange and regional climate. Overturning also ventilates the deep ocean and redistributes carbon, oxygen and nutrients. Changes in its strength can alter regional sea-surface temperatures, air-sea heat release and atmospheric circulation. They need not eliminate every western-boundary current or every source of European warmth: wind-driven ocean circulation and atmospheric heat transport are distinct, interacting contributions. The two-box model below isolates a density/salinity feedback rather than representing all these mechanisms.

Illustration TODO. The source's two-volume circulation sketch on page 4 is unavailable. Its exact arrows and layout cannot be recovered. The surviving flow law, salinity/temperature definitions and printed steady polynomial allow the calculation below under an explicit freshwater sign convention; no unavailable diagram values are invented.

Use $\Delta T=T_2-T_1>0$, $\Delta S=S_2-S_1$, and let positive $q$ be the thermally driven state with cold box 1 denser than warm box 2. Take positive $F$ to freshen box 1, represented in a fixed-volume salt budget by an equivalent salt removal $F{S}_1$ from that box and addition to the other reservoir. This is the freshwater convention compatible with the printed $+BF$ and positive $B$. It can represent a compensating saline volume transfer or the usual virtual-salt-flux approximation. The exact pipe assignment is not claimed to have been read from the missing illustration.

Exchange in either circulation direction reduces a salinity difference at rate proportional to $|q|\Delta S$. Steady salt balance therefore givesThe supplied linear equation of state givesEliminating $\Delta S$ yields the two-box overturning model equationHere $A$ and $B$ both have transport units when $q,F,k$ do, while $\alpha \Delta T$ and $\beta S_1$ are dimensionless. Although the source calls $q$ a mass flow, its stated Sverdrup units are volume per time; the common reference density cancels in the Boussinesq salt/volume approximation. At $F>0$, $q=0$ cannot be a finite-salinity steady state in this model. At $F=0$, zero flow needs separate treatment instead of division by $|q|$.

The circulation weakens when freshwater freshens the cold reservoir, and the weaker circulation removes the resulting salinity contrast less effectively. This salt-advection feedback is the mechanism behind multiple steady states, not an additional externally imposed switch.

Use all coefficients exactly as printed and use Sv consistently for transports. The present salt contrast and density-driven friction coefficient areThusThe large numerical value of $k$ reflects division by a small dimensionless density contrast, about $2.07\times 10^{-6}$ with the supplied effective expansion coefficients. It is not the actual overturning transport, which remains $15$ Sv.

For $q>0$, the steady equation is a quadratic with rootsThe critical freshwater forcing in a two-box overturning model occurs where the two roots coalesce:There is also a negative-flow branch for every $F>0$:At the present forcing, the three roots are $15.000$, $6.774$ and $-3.950$ Sv. The prescribed present state is therefore the upper positive branch, not an arbitrary choice among the roots.

Stability needs a salinity evolution law, not just the algebraic steady equation. For fixed temperature contrast and reference salinity, the standard salt-contrast dynamics can be writtenwhere $V_{\mathrm{eff}}>0$ summarizes the reservoir volumes. The linear slope at a positive root is $(A-2q)/V_{\mathrm{eff}}$: the upper branch is stable and the lower branch unstable. At a negative root it is $(2q-A)/V_{\mathrm{eff}}<0$, so the reversed branch is stable. The two stable states below the fold have different attraction basins, demonstrating bistability. Their stability is conditional on this fixed-temperature, fixed-reference-salinity closure.

As forcing increases quasistatically from $0.5$ Sv, the normal branch weakens and reaches a saddle-node bifurcation at $F_c$. Beyond it no positive steady circulation remains, so the reduced dynamics move toward the stable reversed branch. The reversal is salinity dominated: the fresh cold reservoir is no longer the denser one. Recovery under a later forcing change can depend on the system's state and feedback history.

Signed radial balance. Let $r$ increase outward from the vortex centre and $U$ be signed azimuthal velocity, positive anticlockwise. At $45^{\circ }$N, $f=2{\Omega }_{\oplus }\sin 45^{\circ }>0$. The positive speed is $V=|U|$, $f$ is the Coriolis parameter, the density is $\rho$, and $p_r=\partial p/\partial r$. In circular motion the inward centripetal acceleration is $-U^2/r$, while the outward radial acceleration from the Coriolis force is $fU$. The pressure gradient acceleration is $-p_r/\rho$. HenceThe term $U^2/r$ is an acceleration, not an additional physical force in the rotating-frame momentum equation. It can be moved to the force-balance side and depicted as outward centrifugal force in a frame following the circular parcel.

For a cyclone, $U=+V$ and $p_r/\rho =V^2/r+fV>0$: pressure increases outward and the centre is a low. Its inward pressure force balances the outward Coriolis and centrifugal terms. For a conventional anticyclone, $U=-V$ and $G=-p_r/\rho >0$; the balance isThe inward Coriolis force balances the outward pressure and centrifugal terms, and pressure has a maximum at the centre.

Anticyclonic branches and near-centre regularity. The anticyclonic gradient-wind branches follow from the quadratic balance. Solving givesThe minus root is the normal branch connected continuously to geostrophic flow for small forcing. The plus root is an anomalous high-speed branch and must not be discarded algebraically, even though it often fails the stability conditions relevant to large-scale anticyclones. Reality requiresA finite smooth central pressure therefore has zero radial gradient at the origin, with $-p_r/\rho$ no larger than order $r$. For $p=p_c-\rho C{r}^2/2+\cdots$, the limiting curvature obeys $C\le f^2/4$. A pressure maximum with a finite nonzero limiting outward pressure force cannot support the required circular anticyclonic balance arbitrarily close to its centre.

For positive $\mathrm{Ro}=V/(fr)$ and geostrophic magnitudes,The cyclonic gradient wind is subgeostrophic and the normal anticyclonic wind supergeostrophic. On the latter branch $0\le \mathrm{Ro}\le 1/2$. At a point on the eastern side of a cyclone, geostrophic and actual winds point north but the actual arrow is shorter; the ageostrophic difference points south. On the eastern side of a normal anticyclone both geostrophic and actual winds point south, with the actual arrow longer and its ageostrophic addition also south.

At lower latitude, the same $V$ and radius give a larger Rossby number and larger differences from geostrophic wind. The contour topology remains a central low or a conventional central high, but for fixed $V,r$ the required radial pressure gradients change with $f$: the cyclonic gradient is $\rho (V^2/r+fV)$ and the anticyclonic outward pressure force is $fV-V^2/r$. The latter conventional branch cannot be maintained once the reduced Coriolis contribution is too weak. Actual spatial gradients, rather than an arbitrary spacing of contour values, determine how tightly drawn contours should be.

The high-centre conclusion is conditional on the conventional anticyclonic branch. Gradient balance alone also allows a sufficiently rapid clockwise circulation about a low: if $V>fr$, then $p_r/\rho =V^2/r-fV>0$. For example, $V=2fr$ gives $p_r/\rho =2f^{2}r>0$. Such a configuration is outside the usual geostrophic anticyclone regime. Tornado-scale vortices have very large Rossby number, so the dominant balance is nearly cyclostrophic, $p_r/\rho \simeq V^2/r>0$. Their centres are low pressure irrespective of the sign of rotation; low pressure supplies the inward acceleration needed to bend the fast flow.

Use $u$ as a spectroscopic wavenumber, in the same units throughout. Write the mass absorption coefficient of the absorbing species asHere $S$ is the spectral line strength, the integrated mass absorption coefficient; it is not the line-centre coefficient. If ${\kappa }_{\nu }$ has units ${\mathrm{m}}^2\mathrm{k}{\mathrm{g}}^{-1}$ and $\nu$ is in $\mathrm{c}{\mathrm{m}}^{-1}$, then $S$ has units ${\mathrm{m}}^2\mathrm{k}{\mathrm{g}}^{-1}\mathrm{c}{\mathrm{m}}^{-1}$. For absorber mass density ${\rho }_a$, the volume coefficient in Beer-Lambert law is ${\rho }_a{\kappa }_{\nu }$.

Lifetime broadening. Finite lifetimes of the upper and lower states destroy optical coherence, even for an isolated molecule. Exponential coherence decay gives a Lorentzian line shape. For population decay rates ${\Gamma }_u,{\Gamma }_l$ and no extra dephasing, its wavenumber half-width is $({\Gamma }_u+{\Gamma }_l)/(4\pi c)$, with $c$ expressed in the length unit reciprocal to $\nu$. This sets a natural linewidth and can matter at extremely low collision rates; it is generally much smaller than atmospheric collisional or thermal widths.

Thermal Doppler broadening. A Maxwell-Boltzmann distribution of molecular line-of-sight velocities givesThis Gaussian distribution has a $1/e$ half-width ${\Delta }_D$, determined by absorber molecular mass $m_a$ and temperature. It is especially relevant in the rarefied upper atmosphere, where collisions are infrequent.

Pressure broadening. Collisions perturb transition frequencies and phases. In the impact regime, exponential phase decorrelation givesHere ${\Delta }_L$ is the half-width at half maximum, proportional to the coherence-destroying collision rate. It depends on perturber density, molecular collision properties, and temperature; it is often parameterized by ${\Delta }_L\propto p(T_{\ast }/T{)}^{n_T}$ for a specified broadener. It dominates much of the dense lower atmosphere. The Lorentzian line shape has algebraic $\delta {\nu }^{-2}$ wings, whereas the Doppler profile has rapidly decreasing Gaussian wings. Where both mechanisms matter, their convolution is a Voigt profile. A pure pressure-broadened profile should not be extrapolated to arbitrarily low pressure while discarding the finite Doppler core.

Below cloud base, an initially unsaturated closed parcel expands as ambient pressure falls. In the absence of heat exchange its temperature falls approximately at the dry adiabatic lapse rate $g/c_{p,d}$. Its vapour mixing ratio is conserved, but its saturation mixing ratio decreases strongly as it cools. Hence relative humidity rises until the lifting condensation level is reached. The density normally decreases through expansion; the environmental profile determines the parcel's buoyancy and thus whether ascent can continue.

Near saturation, slight atmospheric supersaturation permits suitable cloud condensation nuclei to activate. Köhler theory describes the competition between curvature, which raises equilibrium vapour pressure, and dissolved solute, which lowers it. Cooling generates supersaturation, while condensation onto the new droplets consumes it. Their balance normally keeps in-cloud relative humidity close to, but not identically at, saturation. The activated droplet number depends on the aerosol population and the maximum supersaturation produced by the ascent.

Above cloud base, condensation lowers the vapour content and raises the liquid-water content. The latent heat released partly offsets expansion cooling, so the moist adiabat has a smaller lapse-rate magnitude than the dry one. In a dilute approximation,A closed adiabatic parcel conserves total water $q_t=q_v+q_l$ and approximately its moist entropy; retaining condensate makes this a reversible-cloud idealization, rather than a precipitating pseudoadiabat. In a shallow cloud, the saturation-mixing-ratio decrease and slow density variation yield nearly constant condensation per height, so liquid-water content is approximately linear above base. Once activation is complete, assuming nearly fixed droplet number per volume gives increasing droplet size as condensate accumulates. Entrainment, continued activation, collision-coalescence and atmospheric precipitation can all violate that simplified description; they are excluded from the stated retrieval model.

Put the rigid lid at $z=0$, the interface at $z=-h(x,y,t)$ and ${\rho }_1={\rho }_0-\Delta \rho$. Assume hydrostatic vertical balance and a resting deep layer. Its pressure is $p_2=C-{\rho }_{0}gz$, with horizontally uniform $C$. Pressure continuity at the interface then fixes the upper-layer pressure:Hence its horizontal pressure gradient at fixed $z$ is ${\nabla }_h{p}_1=\Delta \rho g{\nabla }_{h}h$, independent of depth. Taking the horizontal momentum equation and a vertically uniform upper-layer velocity givesThe exact density denominator for upper-layer acceleration is ${\rho }_1$; the familiar ${\rho }_0$ denominator is the small-density-contrast Boussinesq approximation. The restoring pressure is weaker than for a free surface because only the density difference contributes.

The horizontal shallow-layer equations use the vertical component of planetary rotation: $f=2\Omega \sin \phi$. Thus the source's vector notation is interpreted as the projected vector $f\mathbf{k}$ after the traditional approximation. Taking the full magnitude $2\Omega$ as a midlatitude Coriolis parameter would discard the latitude dependence. Together with layer-volume conservation, the model has ${\partial }_{t}h+{\nabla }_h\cdot (h\mathbf{u})=0$.

Both descriptions are depth-averaged hydrostatic shallow-layer models. In a one-layer free-surface shallow water equations, a homogeneous layer overlies a rigid bottom and its upper surface moves. Surface elevation changes give a restoring acceleration based on full gravity $g$, with external-wave speed $\sqrt{g{H}_0}$. Barotropic coastal storm surges, long surface gravity waves or a basin seiche are examples.

In the rigid-lid reduced gravity model, the movable boundary is a density interface beneath the active upper layer, while the much deeper lower layer remains at rest. The internal interface mode instead has speed $\sqrt{g^{\prime }{h}_0}$ and represents, for example, thermocline displacements and an upper-ocean baroclinic eddy. The mathematical shallow-water form is retained with $g$ replaced by $g^{\prime }$ and with the appropriate interface-depth convention. A freely moving two-layer surface would support both external and internal modes; the rigid lid filters out the fast external mode, and the resting abyss removes its independent dynamics.

A numerical scale estimate requires a reference ocean density and latitude not specified exactly. Taking ${\rho }_0\simeq 1025\mathrm{kg}{\mathrm{m}}^{-3}$, $g\simeq 9.81\mathrm{m}{\mathrm{s}}^{-2}$ and a representative midlatitude $|f|\simeq 10^{-4}{\mathrm{s}}^{-1}$ givesFor a small-amplitude interface displacement of scale $L=1\mathrm{km}\ll L_R$, gravity drives rapid horizontal acceleration and radiates most of the displacement in inertia-gravity waves. The gravity adjustment time is $L/c\sim 8$ minutes, much shorter than $1/|f|\sim 2.8$ hours. A weak balanced remnant remains in an unbounded inviscid adjustment problem; energy is carried away by waves, not dissipated merely by geostrophic adjustment.

For $L=100\mathrm{km}\gg L_R$, rotational adjustment produces a largely retained interface displacement with a geostrophic current. Wave emission and inertial oscillations accompany the adjustment, but the balanced structure is much more important. To quantify that distinction under an initially resting, linear single-Fourier-component assumption, potential-vorticity conservation givesThis is about $0.002$ for the kilometre scale and $0.95$ for the hundred-kilometre scale. These fractions use $L$ as an inverse wavenumber scale, not a full wavelength, and are illustrative rather than unique predictions for an unspecified disturbance or bounded basin.

Take an isothermal hydrostatic reference atmosphere, ${\rho }_0(z)={\rho }_s{e}^{-z/H}$, and linearize about rest. The horizontal acceleration and pressure-gradient force retain the equatorial Coriolis parameter $f\simeq \beta y$, with $\beta =2\Omega /a$. The traditional approximation retains only its vertical component. Quadratic perturbation advection is neglected. The vertical acceleration is balanced by perturbation buoyancy and pressure gradient after removing the reference hydrostatic balance.

The continuity constraint is the anelastic approximation, $\nabla \cdot ({\rho }_0\mathbf{u})=0$, giving $u_x+v_y+w_z-w/H=0$. It retains stratification-driven expansion while filtering acoustic density fluctuations. It is therefore not the full compressible continuity equation. Adiabatic displacement in a stably stratified background gives $b_t+N^{2}w=0$, with $N$ the Brunt–Väisälä frequency; upward displacement develops restoring negative buoyancy.

Pressure convention matters. With the entropy-buoyancy equation just described, a consistent reduced-pressure form writes the pressure acceleration as $-\nabla \Phi$, where $\Phi =p^{\prime }/{\rho }_0$. Its vertical term is $-{\partial }_z(p^{\prime }/{\rho }_0)$, not simply $-{\rho }_0^{-1}{p}_z^{\prime }$. The original source explicitly specifies variable ${\rho }_0(z)$ but prints the latter. Those differ by $p^{\prime }/({\rho }_{0}H)$, so the two cannot silently be equated. Parts b and c display the consequence and then solve the explicitly corrected reduced-pressure model that yields the requested wave equation. This is a mathematical inconsistency in the printed set, not missing source material.

First check the printed equations literally. Set $v=0$ and use ${\rho }_0^{\prime }/{\rho }_0=-1/H$. Differentiating continuity in time and using horizontal momentum givesVertical momentum and buoyancy give $w_{tt}+N^{2}w=-p_{zt}/{\rho }_0$. Combining them therefore yieldswhich differs from the requested equation. The factor two arises when differentiating both ${\rho }_0$ and $w_z-w/H$. Keeping variable density in continuity while treating it as constant during that differentiation would be inconsistent.

For the intended anelastic reduced-pressure Kelvin wave, make the explicit correction $w_t-b+{\Phi }_z=0$, while $u_t+{\Phi }_x=0$ and $\Phi =p^{\prime }/{\rho }_0$. Continuity then gives ${\Phi }_{xx}=w_{zt}-w_t/H$, and vertical momentum gives $w_{tt}+N^{2}w=-{\Phi }_{zt}$. Cross differentiation now producesThis proves the advertised equation for the corrected model, without claiming it follows from the incompatible printed set.

A concrete counterexample makes the distinction testable. In nondimensional units $H=N=k=1$, take $m=0$, $\omega =2/\sqrt{5}$ and $w=e^{z/2}{e}^{i(x-\omega t)}$, multiplied by its Gaussian meridional envelope. This satisfies the last boxed wave equation, but substitution into the literal printed-system wave operator gives $-2w/5\ne 0$. No choice of accompanying pressure and horizontal velocity can then make it a solution of all five printed equations. Conversely, the literal operator has vertical factors $e^{z/H+imz}$ and dispersion ${\omega }^2=N^2{k}^2/(k^2+m^2)$; its scaling is different from the requested half-scale-height factor.

Work with the explicitly corrected reduced-pressure system in b. Write $d=1/(2H)$ and use the complex wave amplitude, taking real parts at the end. The vertical factor satisfies $({\partial }_z^2-H^{-1}{\partial }_z)w=-(m^2+d^2)w$, so the wave equation givesFor a Kelvin wave, continuity, zonal momentum and the buoyancy equation then determine all polarizations:In particular $u,w,b,\Phi$ have amplitude proportional to $e^{z/(2H)}$, whereas pressure $p^{\prime }$ has amplitude proportional to $e^{-z/(2H)}$. Directly substituting these expressions into vertical momentum gives ${\Phi }_z=-i\omega (m^2+d^2)w/k^2=i(\omega -N^2/\omega )w$, verifying both the sign and the dispersion relation.

Meridional momentum is $\beta yu+{\Phi }_y=0$. Since $\Phi =cu$, it integrates to a Gaussian $\exp [-\beta y^2/(2c)]$. ThusThe eastward intrinsic branch $c>0$ is selected by decay away from the equator; choosing $c<0$ would make the meridional amplitude grow without bound. Although the squared dispersion relation has both frequency signs, a meridionally trapped Kelvin wave has this one phase-propagation direction.

For short vertical wavelengths, $m^2\gg k^2+d^2$, so ${\omega }^2\simeq N^2{k}^2/m^2$, $c\simeq N/|m|$ and $L_D^2\simeq N/(\beta |m|)$. The wave is a slowly oscillating, hydrostatic internal gravity wave with strong vertical variation and relatively narrow trapping. Here $u\simeq -(m/k)w$ and $p^{\prime }/{\rho }_0\simeq -(\omega m/k^2)w$; the buoyancy expression above remains valid.

For long vertical wavelengths, $|m|\ll \sqrt{k^2+d^2}$, the exact limiting speed is $c=N/\sqrt{k^2+d^2}$. To obtain the stated height-controlled limit also requires long horizontal wavelength, $|k|H\ll 1$. With both $|m|H\ll 1$ and $|k|H\ll 1$,The finite density scale height prevents the hydrostatic $1/|m|$ phase speed from diverging as $m\to 0$. The polarizations tend to $u\simeq -iw/(2Hk)$, $\Phi \simeq -iNw/k$ and $b\simeq -iNw/(2Hk)$. Long vertical wavelength alone at fixed large $kH$ would not justify $c=2NH$. In the opposite horizontally short limit the frequency approaches $N$ instead. The growth $e^{z/(2H)}$ compensates decreasing reference density in wave-energy density; a radiating wave train is a local modal description, not a finite-total-energy eigenfunction over an infinite vertical domain.

$S=e/e_s^l$ is the ambient saturation ratio over a flat liquid-water surface; $S-1$ is the liquid atmospheric supersaturation. Curvature increases the surface equilibrium vapour pressure, producing the Kelvin penalty $A/r$, while dissolved solute lowers it and gives the opposing Raoult contribution $B/r^3$. Thus the dilute Köhler equilibrium curve is $S_{\mathrm{eq}}-1=A/r-B/r^3$, and the numerator is the atmospheric supersaturation available to drive growth beyond surface equilibrium.

$F_d^l$ describes the resistance of vapour diffusion through the surrounding air. $F_k^l$ describes the thermal-conduction resistance to removing latent heat released by condensation: a warmed droplet surface has a higher equilibrium vapour pressure and hence a weaker mass-transfer drive. Common continuum approximations are $F_d^l\simeq R_{v}T/(D_v{e}_s^l)$ and $F_k^l\simeq L_v[k_{\mathrm{air}}T{]}^{-1}[L_v/(R_{v}T)-1]$, showing the roles of diffusivity, conductivity, saturation pressure and latent heat. Both resistances have the supplied units $\mathrm{s}\mathrm{m}\mathrm{k}{\mathrm{g}}^{-1}$ and add for this coupled transport problem.

For a spherical droplet of constant density, $m=4\pi \rho r^3/3$. Multiplying the radius equation by $4\pi \rho r$ givesAt sufficiently large radius, curvature and solute corrections are negligible, soThis limit also requires sustained positive atmospheric supersaturation. It need not describe the earliest submicrometre stage or an unactivated haze particle.

In the quasisteady diffusion approximation, the vapour-density disturbance outside a small stationary particle satisfies Laplace's equation, with a specified equilibrium value on its surface and a uniform far-field value. This has the same mathematical boundary-value problem as the electrostatic potential outside a conductor held at fixed potential. The integrated normal gradient is proportional to the conductor's geometrical capacitance.

Here $C$ is a length: a sphere of radius $r$ has $C=r$. The corresponding electrical capacitance in SI units would be $4\pi {ϵ}_{0}C$, not $C$ farads. Replacing the spherical diffusive conductance $4\pi r$ by $4\pi C$ accounts for the crystal's non-spherical geometry. Coupling the analogous heat-conduction problem to vapour diffusion supplies the ice-specific thermal and diffusional resistances.

The approximation assumes continuum diffusion, equilibration of the surface vapour concentration, an approximately isothermal particle, a slowly changing shape, weak perturbations and a far field sufficiently distant relative to the crystal size. It neglects ventilation by falling motion, interactions with neighbouring particles and facet-dependent attachment kinetics. Shape affects diffusive growth through capacitance, not merely through projected collecting area as in a collision calculation.

The target liquid-water mass isThe thin-plate ice-growth time uses the prescribed ice mass law, which therefore requires $r_{i\ast }=(m_{\ast }/\alpha {)}^{1/3}=45.31\mu \mathrm{m}$, not the same radius as the water droplet. With $C=2r_i/\pi$ and $F_i=F_k^i+F_d^i$,Hence, for negligible seed mass and constant conditions,This is about $8.24$ times faster than the liquid estimate to the same mass. Both laws have $dm/dt\propto m^{1/3}$ under the given shape assumptions. Ice benefits here from five times the atmospheric supersaturation and a larger geometrical capacitance for the same mass, partly offset by its slightly greater transport resistance. The mass coefficient $\alpha$ encodes the assumed plate shape; replacing it by bulk ice density would change the supplied model.

In a mixed-phase cloud, saturation vapour pressure over ice is lower than over supercooled liquid water. Air can therefore be subsaturated with respect to liquid yet supersaturated with respect to ice: droplets evaporate while ice crystals grow, transferring water through the vapour phase. This is the Wegener–Bergeron–Findeisen mechanism. Subsequent aggregation, riming and melting can produce atmospheric precipitation much faster than condensation alone.

The large difference is not a universal fixed ice/liquid time ratio. Atmospheric supersaturation evolves, latent heating increases thermal resistance, crystal habits and capacitance change, and surface attachment or ventilation may matter. Ice nucleation and the availability of freezing nuclei also precede the calculated growth. The boxed time describes the specified already-present crystals under maintained conditions, not a complete prediction of the time until rain begins.

Linearize about a uniform horizontal flow $U$ and a hydrostatic, vertically stratified basic state. The basic-flow material derivative is $D_0={\partial }_t+U{\partial }_x$. The first two equations are the horizontal and vertical perturbation momentum balances: pressure gradients accelerate the disturbance, while the vertical equation retains the buoyancy force $-g{\rho }^{\prime }/{\rho }_0$ after the basic hydrostatic balance is subtracted.

The Boussinesq approximation neglects fractional density variations in inertia and enforces incompressibility, $u_x+w_z=0$, but retains them in buoyancy. Advecting the stable background density gradient gives $D_0{\rho }^{\prime }+wd\bar{\rho }/dz=0$. With $N^2=-(g/{\rho }_0)d\bar{\rho }/dz>0$, this is the fourth supplied equation. A displaced parcel then has a restoring buoyancy tendency.

Small amplitude allows products of perturbations to be discarded; two-dimensional motion removes $y$ dependence and transverse flow. Rotation, viscosity, diffusion, diabatic heating and mean shear are neglected. ${\rho }_0$ is the reference inertial density, whereas the small background density gradient still sets $N$. In a compressible real atmosphere the appropriate stability measure involves potential temperature; this Boussinesq system is a local idealization that filters acoustic waves rather than a full compressible atmospheric model.

Take $S$ as the integrated line mass absorption coefficient, and $q$ as the absorber mass fraction, assumed constant; for a pure absorbing gas set $q=1$. All spectral intervals below use the paper's wavenumber coordinate. The normalized Lorentz line shape isIf instead an absorption coefficient per length is wanted, multiply by absorber density. Neglecting temperature dependence, pressure broadening gives $\Delta (p)={\Delta }_{s}p/p_s$. In the line wings, with $d=|\nu -{\nu }_0|\gg \Delta (p)$,Define the line-opacity parameter and characteristic wavenumber widthAt an emitting level ${\tau }_{\nu }\sim 1$, the wing estimate implies $d/\Delta (p)\sim \sqrt{C}$. Thus the strong-line limit relevant to outgoing radiation isIt ensures that the opaque-to-transparent transition lies in the Lorentz wings, rather than the narrow central core, and the characteristic absorption width is much larger than ${\Delta }_s$. At a grazing ray angle the criterion is even easier to satisfy because of the larger optical depth. This is a criterion for the radiatively important wings, not a claim that collision broadening remains physical at arbitrarily small pressure.