virial pressure-radius curve with uniform-sphere binding energy. Put . Isothermality makes , so the virial theorem gives
The curve crosses zero at , rises to its maximum at
then decreases toward zero at large radius. Since , a cloud displaced outward at fixed external pressure has a restoring pressure deficit only on the descending branch. Thus the homologous, isothermal radial mode is stable for , unstable for , and marginal at . The negative-pressure branch is not a physical equilibrium with a positive confining pressure. This conclusion uses the stipulated gravitational-energy approximation, not the exact isothermal density profile.
Figure 1.
Original pressure-radius curve and stable versus unstable branches of the stipulated virial cloud model
. Details.
critical mass from a virial pressure turning point. At the turning point, ; hence
Using the cloud's mean density , the condition becomes
Using the boundary density instead changes the numerical coefficient, not this scaling. The coefficient should not be identified with the exact Bonnor-Ebert sphere threshold because the assumed potential energy is different.
mass and luminosity weighting of a stellar initial mass function. Impose finite birth-mass cutoffs . The stellar birth mass and main-sequence luminosity scale as
For a broad mass interval, the low-mass end controls the first integral and the high-mass end the second. Per logarithmic mass interval the weights are respectively and . Finite lifetimes remove massive stars from an aging population, so the luminosity conclusion applies while those stars remain on the main sequence.
exponential closed-box gas consumption. The gross formation rate is , but the permanent locking rate is . The instantaneous recycling approximation and mass conservation give , . Therefore
Here means surviving locked stellar mass in the galactic gas fraction, not the cumulative gross mass that has passed through star formation.
locked-mass and returned-mass yield conventions. Pre-existing metals are locked at rate ; the stipulated newly produced metal source is . Thus
For and initially pristine gas,
The effective yield is per unit permanently locked mass; it is not the source parameter . For , and the direct time equation gives . For , the explicitly stipulated metal source vanishes, so this particular enrichment model gives . Constant yield is a trace-enrichment approximation; the formal logarithmic divergence as gas is exhausted must not be interpreted as a physical metal fraction exceeding one.
backward extrapolation of exponential closed-box star formation. Denote the present age of the universe by . With the last paragraph's ,
Adopting a present age of the universe for the estimate gives about , requiring of gas at that epoch in the same model. The boxed expression makes the dependence on the adopted present age explicit. This is an extreme history for a present Milky-Way-like system, although intense high-redshift star formation by itself is not forbidden.
The extrapolation forces a fixed gas reservoir to decline exponentially over many gas-depletion times. Gas accretion, outflows, evolving gas density and star-formation efficiency, delayed recycling, and galaxy assembly break that assumption. The empirical disk-averaged relation is not a universal constant-efficiency law for all gas: the original measurements combined Balmer series H-alpha, atomic gas, carbon monoxide and infrared data, with tracer calibration and dust corrections entering the inference.[ref] Cosmological calculations explicitly support continuing cold-gas supply to rapidly forming high-redshift galaxies, providing a physical alternative to this closed-box backward extrapolation.[ref]
For a monatomic ideal gas, write the thermal energy density as . At low optical depth, radiative loss uses an astrophysical cooling function, often written , while contains heating. At fixed density the net radiative gas cooling time is approximately when the denominator is positive; a cooling flow closer to constant pressure must include enthalpy and compressional work. Comparing this time with the free-fall time of a uniform sphere determines whether pressure can be lost during gravitational contraction.
Hot ionized gas cools by bremsstrahlung and recombination; collisional excitation followed by photon escape supplies atomic and metal-line cooling at lower temperatures. molecular line cooling through rotational or vibrational lines and dust thermal emission become important in cold dense gas. Abundances, ionization state, density-dependent collisional de-excitation, and photon trapping all change the cooling function. Optically thick gas requires radiative transfer, not a simple local density-squared loss law.
Heating comes from photoionization, grain photoelectric heating driven by ultraviolet light, cosmic rays, shocks, turbulent dissipation and compression. photoionization heating deposits the photon energy above the ionization threshold as electron kinetic energy. A thermal equilibrium satisfies ; its stability depends on how the net loss changes under the relevant pressure or density constraint, not merely on equality. Different stable balances permit a multiphase interstellar medium. Spectral-line ratios constrain temperature and ionization, but interpreting them requires density diagnostics, interstellar extinction corrections and a radiation-field model.
Stellar radiation, stellar winds and supernovae return energy, momentum, mass and metals to surrounding gas. Heating, ionization, turbulent support and gas removal can reduce subsequent star formation; swept-up shells can also compress gas locally. Whether thermal energy drives an outflow depends on its radiative gas cooling time relative to expansion: radiative losses can remove energy while leaving substantial momentum transfer. Thus deposited energy alone does not specify the dynamical effect.
The galactic gas regulator has mass budget , where is the returned fraction and is the mass-loading factor. In a quasi-steady state, . This connects feedback to the available gas supply. A wind's kinetic power and momentum flux must fit the source budget and gravitational binding cost. Gas that fails to escape may form a galactic fountain and later return; preferential metal loss changes chemical evolution.
Absorption against a galaxy's stellar continuum provides a concrete outflow test: compare interstellar sodium absorption with the systemic velocity and stellar absorption spectrum. A primary starburst survey found blueshifted interstellar absorption and an inclination dependence consistent with minor-axis outflows. Its saturated lines also illustrate why velocity detection is more direct than a unique column density or mass flux.[ref] Deprojection, covering fraction, ionization, radius and the unseen hot phase limit total outflow-rate estimates.
An accreting supermassive black hole adds radiation, winds and astrophysical jets. Coupling these to the interstellar or circumgalactic medium can suppress cooling or expel gas, especially in deep potentials where stellar driving alone may be inadequate. Jet cavities, line-emitting gas and multiphase outflows test different channels; coexisting activity and low star formation alone do not prove a causal feedback history.
A spiral density wave is an organized density and gravitational-potential perturbation through which stars and gas can move; it is not generally a permanently labelled set of stars. Differential rotation winds a material arm rapidly. In a wave description, a perturbation has phase , with pattern speed distinct from the local stellar angular speed.
For a thin fluid disk the local tightly wound dispersion relation is
Rotation and pressure resist self-gravity. A collisionless stellar disk requires its stellar response function instead of simply replacing by a dispersion; the familiar local axisymmetric stability parameter is . This is distinct from the gas .
Waves exchange angular momentum at Lindblad resonances and near corotation radius. Swing amplification can strengthen a leading disturbance as shear turns it into a trailing one. Bars or companions can force a coherent spiral; recurrent self-excited disturbances also occur. Thus a single indefinitely steady wave is one model, not a description established for every spiral galaxy. Gas compression and shocks along arms can organize star formation, while gravitational torques redistribute angular momentum. Arm brightness, recent star formation and stellar mass density need not trace precisely the same pattern.[ref]
Figure 1.
Differentially orbiting stars can cross a spiral density pattern with its own angular speed
. Source. Details.
In a matter-dominated Einstein-de Sitter universe, the growing linear density contrast is proportional to the scale factor. When , linear evolution no longer describes collapse. An overdense region initially expands more slowly than the background, reaches cosmological turnaround, and then contracts under self-gravity.
The spherical top-hat collapse model supplies a useful idealization: pressureless shells obey . At formal collapse, the linearly extrapolated growing mode is
This is a linear extrapolation, not the true finite nonlinear density at that instant. In the ideal energy-conserving spherical model, virial equilibrium gives a virial radius about half the turnaround radius; its mean density at the collapse time is times the cosmic matter density.
Real perturbations are nonspherical: tides and anisotropic collapse produce sheets, filaments and haloes. Collisionless dark matter develops shell crossing and phase mixing rather than radiatively losing its orbital energy. Baryonic gas can shock, cool and collect within the resulting potential wells. Hierarchical merging then assembles larger systems where small-scale fluctuations have sufficient initial power. Press-Schechter formalism connects a collapse threshold with the statistics of the smoothed initial field, but halo abundances and internal structure need more than the single spherical model. These ideas apply after perturbations become nonlinear; one must not extend pressureless subhorizon equations through all radiation-era or horizon-scale conditions.[ref]
With constant isothermal sound speed , the equation of state is . Spherical hydrostatic equilibrium and mass continuity require
Try . With no central point mass, , and hydrostatic balance becomes . The positive-density branch is therefore
The constant external pressure fixes the truncation radius through , giving
This is the pressure-truncated singular isothermal gas sphere, one hydrostatic solution, not a claim that every isothermal cloud has a singular centre. Its central density diverges, but its enclosed mass tends to zero as .
First specify the perturbation: take fixed gas density, fixed ionic abundance and the same temperature-independent heating, so thermal stability requires at an equilibrium. Let , be the collider density and the atomic density. A two-level optically thin balance is
The critical density of an atomic transition is approximately . At low density, each rare collisional excitation is followed by an escaping photon, so
The isochoric thermal stability of a collisionally excited line therefore gives
At the linear restoring coefficient vanishes; the cooling maximum is marginal to first order and is not robust two-sided stability.
At high density, collisional de-excitation competes successfully with radiative decay, and collisions establish local thermodynamic equilibrium level populations. Detailed balance gives , where . The cooling is then
Thus the high-density fixed-density branch is stable for finite , becoming weakly restoring as the upper fraction saturates at high temperature. It is incorrect to retain the low-density factor in the radiative cooling law once populations are thermalized.[ref]
The isobaric thermal stability of a two-level coolant is a distinct question. If the intended heating is a fixed rate per particle, then and stability tests the cooling per particle. The low-density boundary becomes ; the high-density cooling per particle remains proportional to , so remains stable. If instead heating is constant per volume, the isobaric low-density boundary is , and the high-density stability condition is . The source gives no pressure constraint or explicit heating normalization: the boxed result is the fixed-density answer, with these conditional alternatives stated to avoid asserting a universal boundary.
Figure 1.
Low-density collisional line cooling has a temperature maximum, whereas high-density LTE emission rises toward saturation
. Source. Details.
Let be the separation of galaxy centres and a small displacement from the satellite centre along the line of centres. Expanding the central point-mass field gives the radial differential acceleration
Inside a uniform-density satellite, the restoring self-gravity is
Comparing the two at the same displacement yields the differential-gravity disruption limit of a uniform satellite:
The satellite radius cancels because both accelerations are linear in . This is the instantaneous tidal-gravity estimate requested by the given criterion, not a full rotating Roche calculation. Including centrifugal terms for a circular corotating orbit replaces the radial coefficient two by three, giving while preserving the required scaling. An eccentric orbit requires checking the stronger tide near periapsis, and extended density profiles require a more detailed restoring-force calculation.
Put and let , be constant background matter fractions. Define the total matter perturbation and the relative mode . Adding with the background weights and subtracting the two given equations respectively gives
The total and relative density modes of two pressureless matter fluids are consequently
These are the general consequences of the supplied equations; a single ratio additionally needs initial-mode assumptions.
For the usual dark-matter-dominated growing-mode approximation, take , neglect the decaying modes, and set the relative-velocity mode . Then . With at recombination, , so the baryon catch-up approximation after recombination is
At redshifts , replacing by gives the printed . The printed relation is a high-redshift approximation with a specified relative mode, not an exact identity for arbitrary solutions of the two-fluid equations.
For example, retaining a finite baryon fraction with a growing total mode and constant relative mode gives, with ,
Alternatively, if initially both and while the dark-matter growing mode is already present, then in the dark-matter-dominated limit the relative transient is not zero and
These different solutions show why the initial-velocity convention matters. The quoted catch-up law is the dark-matter-dominated branch with the relative decaying mode neglected, together with the high-redshift approximation stated above. The general solution retains the transient and a finite baryon fraction, so the additional assumptions needed for that special case remain explicit.
The physical conclusion is robust: after decoupling, baryonic gas can fall into pre-existing dark-matter potential wells and its fractional density perturbation catches up. The absolute constant offset need not vanish; it becomes small relative to the growing common mode. The pressureless approximation applies on scales above the baryonic Jeans scale and while perturbations remain linear; subsequent collapse requires nonlinear treatment.[ref]
Figure 1.
Baryons approach the dark-matter density contrast, with transient behaviour set by the initial relative velocity
. Source. Details.
An observed catalogue is the result of a population, an instrument and a selection rule. Its astronomical selection function can depend on flux, color, size, surface brightness, sky location, cadence and follow-up success. The relevant cut is often on noisy measured quantities. A model that selects on noiseless true flux instead can predict a different sample even when it uses the same nominal threshold.
Malmquist bias occurs because brighter members of a population are visible through a larger volume in a flux-limited survey. For a uniform Euclidean population, and . The selected luminosity distribution is therefore weighted toward large . For an intrinsic Gaussian absolute-magnitude distribution with mean and variance , multiplying by and completing the square gives the Gaussian absolute-magnitude Malmquist shift
This assumes the stated simple population, a complete flux cut and no additional boundaries. Exactly identical noiseless standard candles have a distance truncation, not a luminosity spread to reweight. Scatter, extinction, cosmological volume and source evolution need their own model.
Lutz-Kelker bias concerns parallax errors, spatial population weighting and selected calibration samples. Uniform spatial density gives locally, hence on the physically valid range. More distant objects can scatter into a measured-parallax cut. In addition, inversion is nonlinear and gives misleading estimates at low signal-to-noise; a negative measured parallax is possible even for a positive true distance. A posterior instead combines the parallax likelihood with a justified population prior. An unbounded uniform-volume prior is improper here; a finite population model or a proper distance prior for parallax inference is necessary.[ref]
Bayes' theorem makes the assumptions explicit. If denotes population parameters, latent properties and the catalogue inclusion event, a selected-sample likelihood has the form
This selection-normalized astronomical likelihood accounts for the model-dependent probability of entering the catalogue. For individual inference, selection information already determined by the fully conditioned measured data must not be multiplied in a second time. The sampling design, inferential target and actual conditioned information determine the appropriate expression.
Other effects include Eddington bias, where steep number counts and measurement scatter move more faint objects upward than bright objects downward, and surface-brightness selection that misses diffuse galaxies despite substantial total flux. Dust extinction affects both detectability and inferred luminosity. Variable-source detection depends on cadence and duration; spectroscopic follow-up can favor easy-to-measure redshifts or particular colors. Blends, calibration drift and redshift errors can mimic population changes.
Examples include magnitude-limited supernova samples used in distance-redshift fits, where luminosity scatter and follow-up must be modeled, and radial-velocity planet surveys, which favor large-amplitude, short-period systems and leave inclination degeneracies. Forward modeling through the survey pipeline, injection-and-recovery estimates of completeness, calibrated noise and properly modeled censored or nondetected observations address these effects. A universal correction or an unexamined inverse-volume weight does not replace an explicit generative and selection model.
Thermal bremsstrahlung arises from electron acceleration in ion Coulomb fields. For a thermal plasma its optically thin spectral emissivity has the dependence
The exponential high-energy slope helps constrain temperature in X-ray spectra, with composition, absorption and other emission components included. At radio frequencies with , the thin spectrum is much flatter and does not alone determine a unique temperature; the Gaunt factor supplies a weak frequency dependence.
For an approximately isothermal layer, radiative transfer gives . Radio free-free opacity scales approximately as , where the hydrogen-plasma emission measure is . A resolved optically thick brightness temperature approaches ; an unresolved filling factor lowers the observed value. The turnover between thin and thick regimes helps constrain emission measure.
An angular extent and distance give a physical size, and an assumed geometry supplies a line-of-sight depth. Combining that depth with emission measure estimates an rms density, . Clumping, filling factor, composition and unknown depth preserve a density-depth degeneracy of free-free radiation. These measurements constrain different combinations of temperature, density and size rather than independently fixing all three without assumptions.
The observed population is the parent population filtered by an astronomical selection function and measurement errors. Malmquist bias arises because a flux-limited galaxy selection admits luminous objects from larger volumes: in Euclidean space and . Thus observed luminosity counts are weighted by . Even a nominal standard candle with intrinsic scatter is preferentially observed on its bright side near the limit; ignoring this can underestimate distances. A volume-limited sample or a justified weighting helps, subject to completeness, geometry and evolution assumptions. Separately, Eddington bias occurs when noise scatters a steeply rising population across a detection threshold: more faint objects enter than bright objects leave.
Lutz-Kelker bias concerns noisy stellar parallaxes and population/selection weighting. Locally uniform spatial density gives and, by changing variables, over the model's valid finite range. More distant stars are available to scatter into a parallax-limited sample than nearby stars to scatter out. Inverting an uncertain parallax also introduces a distinct nonlinear-estimation problem. There is no universal magnitude correction for every star: spatial density, luminosity, extinction and how the sample was selected all matter. The classical treatment motivates this issue; it should not replace an explicit generative model.[ref]
Bayes' theorem combines the measurement likelihood function with a justified prior distribution. For example, with a Gaussian parallax error and consistent distance/parallax units,
A proper distance prior for parallax inference might have with a justified scale ; extending to infinity gives an improper posterior here. Negative measured parallaxes remain valid noisy data. For selection conditioning in Bayesian inference, distinguish individual inference from population inference: if selection is determined entirely by the exact data already conditioned on, it adds no further information about that individual's parameters. For population parameters , the selected-data likelihood must instead include the detection normalization . Applying a cut twice is also an error.
Selection can depend on surface brightness, angular resolution, dust, sky coverage, colour and observing cadence as well as total flux. Crowding and unresolved binaries bias distance calibrators. Transit surveys favour large planets and short periods, with an additional geometric transit probability approximately for small planets on circular orbits; injection/recovery tests quantify the remaining transit survey completeness. Cluster catalogues select noisy mass proxies from a steep abundance distribution. Planck's 2013 SZ catalogue explicitly characterized detection completeness and reliability rather than identifying every detection with an unbiased mass census.[ref] Calibrators, artificial-source tests and a forward model of detection, measurement and follow-up are the practical checks on any correction.
The thermal Sunyaev-Zeldovich effect is inverse Compton scattering of cosmic microwave background photons by hot cluster electrons. For an optically thin, nonrelativistic thermal plasma,
The Compton y parameter measures the electron pressure column. The distortion is a decrement below roughly 217 GHz, an increment above it, and has in the Rayleigh-Jeans law limit. Relativistic electron corrections modify the spectrum and null.[ref]
Multifrequency maps separate this spectrum from primary Cosmic microwave background anisotropy, radio sources and dust; atmospheric subtraction, beams and calibration require control. The kinetic Sunyaev-Zeldovich effect has and a CMB-temperature-like spectrum, so it is both a contaminant and a peculiar-velocity probe. Planck demonstrated blind multifrequency cluster detection with external confirmation and a measured selection function.[ref]
For a fixed physical pressure column, has no ordinary cosmological surface-brightness dimming: the distorted spectrum and CMB redshift together. The integrated Compton parameter still scales as for fixed integrated physical pressure. SZ observations therefore find distant clusters and measure thermal pressure; integrated is a useful mass proxy, requiring calibration and a selection model for cosmological counts. X-ray emission scales approximately as , while SZ scales as . Combining them and a measured temperature constrains density, depth and, with angular geometry, distance. gas clumping factor, asphericity, temperature structure and nonthermal support limit that inference.
thermal bremsstrahlung is continuum radiation from electrons accelerated in ion electric fields. For Maxwellian electrons,
At high frequency the exponential cutoff constrains ; X-ray fitting must also model foreground absorption, line emission, detector response and possible multitemperature or nonthermal components. At radio frequencies, , an optically thin free-free spectrum is almost flat, with the weak frequency dependence set mainly by the Gaunt factor.
For pure ionized hydrogen, radio free-free opacity gives
The turnover constrains emission measure; a resolved, optically thick source has brightness temperature approaching . An unresolved source has an unknown beam filling factor, so its observed brightness need not equal its electron temperature. Angular size and an independent distance give a geometric depth ; a uniform model then yields . Clumping instead measures a density-squared average, not the mean density.[ref]
A free-free spectrum alone does not uniquely determine density and size. The density-depth degeneracy of free-free radiation is explicit: for a resolved uniform column, , leaves both emissivity-integrated brightness and optical depth unchanged at fixed temperature. With independent SZ information and known , the electron column is ; then
These equalities require consistent composition and geometry, and uniform density and temperature. Relating to angular extent can give distance; asphericity and the gas clumping factor are important limitations of this combined cluster method.
An electron accelerated in an ion's Coulomb field emits bremsstrahlung. Both electron states are unbound, giving free-free radiation rather than a discrete bound-state line. Averaging over a Maxwellian electron distribution gives an optically thin spectral emissivity with dependence
Here depends on the emissivity convention, and the Gaunt factor supplies the quantum correction. No numerical emissivity normalization is needed for these diagnostics. For a mixture, sum over ion species. The bolometric emissivity is proportional to with a frequency-averaged Gaunt factor. Thermal electron speeds control the spectrum while the encounter rate supplies its density-squared dependence.
temperature can be measured from the high-energy exponential cutoff at photon energy of order , provided the spectrum spans that turnover and foreground absorption and nonthermal contamination are modelled. At radio frequencies , the exponential is nearly one, so the optically thin continuum is a poor thermometer by slope alone. Kirchhoff's law, , and the Rayleigh-Jeans law give for an isothermal thermal plasma. The transfer solution is
in the radio limit. A resolved optically thick region therefore measures its electron temperature through brightness temperature, if it fills the beam. The turnover between thick and thin emission constrains the line-of-sight emission measure for a hydrogen-dominated fully ionized plasma. In the thin regime , rather than . The slowly varying Gaunt factor gives a nearly flat optically thin radio flux-density spectrum, in contrast to the thick limit at fixed projected area.[ref]
To obtain density, one also needs geometry. An angular extent and independent distance give a transverse size; an assumed shape provides a depth or volume. For approximately uniform emitting gas with astrophysical volume filling factor , , so . Alternatively a measured column emission measure and assumed path length give an rms density. Clumping raises the mean squared density relative to the square of the mean: the inferred value is not automatically the mean gas density. Free-free brightness alone does not uniquely determine temperature, density and physical size. Resolved imaging, a distance, optical-depth information or an independent column-density diagnostic is needed to remove those degeneracies. Cluster SZ data are complementary because their leading signal depends on , whereas free-free emission depends on density squared. Metal lines, unresolved temperature structure, foreground absorption, synchrotron emission and background subtraction must be controlled before attributing the continuum to a single thermal component.
Type Ia supernovae are standardizable candles. Repeated multiband photometry measures peak apparent magnitude, light-curve width and colour; spectra identify the supernova class and determine redshift, often from the host galaxy. Nearby objects calibrate the relation between light-curve shape, colour and absolute luminosity. Applying those relations produces a corrected apparent magnitude , not an assumption that every event has identical luminosity. The distance modulus then gives
A common uncertainty in is degenerate with the absolute distance scale, but does not erase the redshift-dependent shape of the Hubble diagram.
luminosity distance probes the expansion history. In a spatially flat model,
For fixed low-redshift calibration, a negative increases the second-order distance relative to a decelerating model. Objects at a given redshift are therefore fainter. A Hubble diagram spanning enough redshift can distinguish this curvature from the common magnitude intercept. In matter-plus-cosmological-constant models, ; fitting the distances can test whether . The two original supernova teams found distances supporting an accelerating expansion in that framework.[ref][ref]
The comparison requires matched rest-frame passbands and K corrections, correction for redshift time dilation in the light curves, a common photometric calibration, and explicit uncertainty propagation. Dust extinction can imitate dimming; multicolour and infrared observations test reddening and the adopted extinction law. Evolution in the progenitor population or host-galaxy dependence can alter standardization, so one compares matched populations and fits host-dependent terms where warranted. Flux-limited discovery introduces Malmquist bias and other selection effects, which must be modelled through the survey detection efficiency. peculiar velocities matter at low redshift, while gravitational lensing and intrinsic scatter broaden the distant sample. Classification and outlier tests reduce contamination by other transient types.
Independent data reductions, nearby anchors, spectroscopic checks and a fit over the full redshift range distinguish an expansion-history signal from one instrumental offset. Complementary cosmological observations help separate matter density, curvature and dark-energy parameters. The acceleration inference comes from the relative distance-redshift relation, conditional on standardization and an expansion model. Distant supernovae integrate the past expansion history; they do not directly measure the instantaneous second time derivative of today's scale factor without such modelling.
For relativistic dust, , the relativistic Euler equation gives . The nowhere-vanishing density therefore makes the four-acceleration zero: is a timelike geodesic vector field. Let be the differential one-form obtained by the musical isomorphism. Its exterior derivative satisfies
Here is the interior product. By hypersurface orthogonality and the Frobenius theorem, . Taking the interior product and using gives
Thus . By the Poincaré lemma, locally , and then . Choose spatial coordinate functions on one level set of and extend them constantly along the integral curves of a vector field , using its smooth flow. This constructs . Since ,
The metric on the orthogonal complement of a timelike vector is positive-definite, so the synchronous coordinates for irrotational dust have
This construction is local; crossing geodesics can prevent these synchronous coordinates from covering a larger region.
Let denote the actual horizontal velocity field. With constant nonzero Coriolis parameter , geostrophic balance and hydrostatic balance, using pressure per reference density , give
Differentiate the first equation in to obtain thermal-wind balance:
Taking the dot product with proves
Thus the horizontal velocity is tangent to the horizontal buoyancy contours at every height. For a fluid whose thermal equation of state makes horizontal buoyancy contours coincide with horizontal temperature contours, there is no horizontal advection across those contours. In a dry terrestrial atmosphere, use potential temperature as the materially conserved thermal variable. A salinity contribution to buoyancy would require separate treatment before interpreting the result as a statement about temperature.
The usual statement about net meridional heat transport also follows explicitly. Introduce a horizontal streamfunction such that . Integration of thermal-wind balance gives
and hence
For a zonally periodic function , or vanishing endpoint contributions, the zonal integral is zero. Likewise, the anomalous eddy buoyancy flux is
The required correlation between meridional velocity and thermal anomaly is absent. A nonzero vertical phase tilt, as in the next part, permits that correlation.
Qualification of the wording. The result does not say that the local heat flux vector vanishes, or that an arbitrary open-domain integral vanishes without boundary conditions. For a concrete counterexample to that literal reading, take and . These satisfy the displayed thermal-wind balance, but is generally nonzero. What the hypotheses alone guarantee is zero transport across horizontal buoyancy contours; the vanishing zonally integrated heat transport uses the stated endpoint condition.
The equations describe the linear reduced gravity shallow water equations for a uniform upper layer over a deep, resting lower layer, or equivalently one separated internal vertical mode. The layer is shallow compared with its horizontal scales, permitting hydrostatic balance, and its horizontal velocity is independent of depth. The Boussinesq approximation retains the small density contrast through . We neglect friction, external forcing, background currents, and quadratic perturbation terms. On an equatorial beta plane, with , and the equatorial Kelvin wave branch imposes . The meridional Cauchy momentum equation then reduces exactly within this linear model to a balance between the Coriolis force and pressure gradient force.
Strictly, the printed must denote the thickness perturbation , with total thickness and . A signed oscillation around zero cannot represent the total positive layer thickness. Define . For a nonzero zonal wavenumber , the zonal Cauchy momentum equation and conservation of mass give
On the eastward branch , the wave polarization is . Meridional geostrophic balance gives
The branch instead grows as and is excluded by decay away from the equator. Matching the paper's Gaussian function therefore gives
The usual equatorial deformation radius is , with the envelope written . Thus the printed exponent uses . Calling that width a deformation radius is a convention; using the standard radius without the factor would fail the printed meridional equation.
Figure 1.
Equatorial Kelvin-wave thickness anomaly and zonal velocity at one instant, using the paper's Gaussian-width convention
. Details.
At an instant, alternating positive and negative thickness anomalies have their largest magnitudes on the equator and decay symmetrically away from it. Eastward arrows occupy the positive anomalies, westward arrows the negative anomalies; the velocity vanishes at the nodes and weakens away from the equator. The whole pattern travels eastward at speed . Meridionally, is geostrophic balance. For a positive anomaly north of the equator, the northward pressure gradient force balances the southward Coriolis force; both directions reverse south of the equator. Zonally, and : the zonal pressure gradient force accelerates the flow, so it is not in zonal geostrophic balance. On the equator the meridional terms both vanish smoothly.
Use eastward , northward , constant layer depth , constant reference density , and a depth-independent horizontal velocity field . Assume a flat bottom and a rigid upper boundary, giving the two-dimensional incompressibility condition . The traditional approximation retains only the vertical component of planetary rotation in the horizontal equations. Take the Coriolis parameter , a constant linear drag rate , and no other torques. The component Cauchy momentum equations are
Define relative vorticity . Apply to the second equation and subtract of the first. The pressure gradient force has zero curl because is constant and mixed partial derivatives commute. Direct expansion of the material-acceleration terms gives
The curl of the rotation terms is
The imposed zonal wind stress contributes , and constant linear drag contributes . Hence
Since , defining absolute vorticity gives
Without the incompressibility condition, the left-hand side would contain the stretching term . Variable layer depth, variable density or spatially varying linear drag would introduce further terms. No small-Rossby number assumption was needed to obtain this constant-depth vorticity equation.