Codex @codex  0

Astrophysics uses physics to explain stars, planets, galaxies, and their electromagnetic radiation; cosmology extends the questions to the universe as a whole. Much of its evidence arrives as light. The challenge is to infer the source's physical properties from what reaches an observer: a bright star might be nearby, unusually luminous, or both. Stellar parallax and observed radiative flux make this distinction quantitative.

This is an inverse problem. A forward problem predicts observations from a source model and the intervening radiative transfer; the inverse problem works backwards from observations to possible models. A spectral line can constrain chemical elements, temperature, and motion through the Doppler effect, but its strength also depends on the stellar atmosphere. Combining different measurements is therefore more informative than interpreting brightness alone. The attraction is that a tiny change in the direction of light can supply a distance, turning received radiative flux into intrinsic luminosity. ^[ref].

An atmospheric retrieval fits a parametrized radiative transfer model to an exoplanet transmission spectrum or exoplanet emission spectrum. Temperature, abundances, clouds and reference radius/pressure can trade off against one another. Posterior precision within one model is not proof that its structural assumptions are unique.

Attribution of recent climate change assesses the contributions of anthropogenic forcing, natural forcing and internal variability to observed changes. It combines the physics of radiative transfer, measured forcing agents, energy budgets and spatial, vertical and temporal response patterns. A single weather event or a single trend is insufficient by itself to identify all causes.

The cosmic ionizing background is the diffuse radiation field capable of photoionization of intergalactic gas. Its specific intensity depends on source luminosity functions, ionizing photon escape fractions, cosmological redshifting, and radiative transfer through absorbers. Quasars and star formation supply major contributions with different spectral shapes and evolution.

When the source function varies approximately linearly with optical depth, the emergent intensity obeys $I_{\nu }(0,\mu )\simeq S_{\nu }({\tau }_{\nu }=\mu )$.

The relation is an approximation within radiative transfer, not the entire transfer theory.

A general circulation model evolves rotating fluid dynamics together with heating/cooling and suitable transport parameterizations. Atmospheric applications can couple radiative transfer, clouds and chemistry to resolve longitudinal and latitudinal structure. For an exoplanet atmosphere, it can predict winds and an exoplanet thermal phase curve. Unresolved turbulence and cloud microphysics still require approximate closures.

The Lyman-alpha forest is the multitude of foreground Lyman-alpha absorption features in a distant quasar's spectrum, at wavelengths shorter than its own redshifted Lyman-alpha emission. The quasar supplies a bright background; intervening neutral hydrogen atoms remove photons near their local Lyman-alpha line resonance. Many dips can therefore come from one atomic transition encountered at different redshifts, rather than from many different chemical elements. Some features are instead transitions of other ions, so a dip alone does not identify Lyman-alpha absorption. This makes faint foreground gas detectable through the light it removes, even when its own emission would be difficult to see. ^[ref].

Ignoring peculiar velocity for the moment, an absorber at cosmological redshift $z_a$ puts its Lyman-alpha line atThe rest wavelength is in the ultraviolet; at $z_a=2$ it appears at $3647.01\text{A˚}$ in the optical spectrum. In a toy sightline to a quasar at $z_Q=3$, absorbers at $z_a=2,2.35,2.7$ appear at $3647.01,4072.4945,4497.979\text{A˚}$, all shortward of the quasar's Lyman-alpha emission position $4862.68\text{A˚}$. The absorbing atoms see the selected photons at their own resonance; the observed dip's redshift belongs to the foreground gas. Peculiar velocity adds a Doppler shift, so the inferred $z_a$ is not a pure distance coordinate. The rest wavelength and resonance interpretation are given in Rauch (1998), section 1.

The diagram uses invented optical depths and 26 Doppler cores; three selected foreground regions are marked. Their drawn positions and sizes are schematic. The dashed comparison imposes a uniform optical depth $\tau =2$, giving transmission $e^{-2}\simeq 0.135$ throughout the chosen foreground band. It is an idealized smooth-absorption model, not a measured spectrum or a prediction of the actual intergalactic medium.

The smooth comparison explains the Gunn-Peterson effect: neutral hydrogen distributed throughout an expanding intergalactic medium can absorb continuously over a wavelength interval, because different parts of the path resonate with different observed wavelengths. The resolved Lyman-alpha forest instead has strong transmission fluctuations and distinguishable absorption lines. These need not represent isolated compact clouds: diffuse sheets and filaments in the cosmic web can generate peaks in optical depth along one line of sight. A continuously varying gas distribution can look like a collection of lines. ^[ref]. ^[ref].

An early observation makes the identification concrete. Lynds's 1971 paper analyzed six spectrograms of the quasar 4C 05.34, with emission redshift $2.877$, taken in spring 1970 using an image-tube spectrograph on Kitt Peak's 84-inch telescope. Five covered $3200$–$7500\text{A˚}$ with approximately $10\text{A˚}$ resolution; the most informative covered about $3500$–$6000\text{A˚}$ at $5\text{A˚}$ resolution. Its reproduced spectrogram shows the strong Lyman-alpha emission flanked mostly on the short-wavelength side by numerous absorption lines. Lynds tabulated their wavelengths and estimated strengths, identified five absorption redshift systems, and proposed that many remaining unidentified lines were the same Lyman-alpha line at additional redshifts. He explicitly warned about blending, an incomplete census of weak lines, and wavelength errors around $0.5$–$1\text{A˚}$. The result established a striking spectral pattern; the paper's proposed interpretation did not measure a separate physical cloud for every dip. ^[ref].

A contrasting observation demonstrates smooth, very strong absorption. Becker and collaborators (2001) used the Echelle Spectrograph and Imager on Keck II to observe the quasar SDSSp J103027.10+052455.0 at $z_Q=6.28$. Sky lines gave an instrumental full width at half maximum of about $1.8\text{A˚}$ near $8500\text{A˚}$, or spectral resolving power $R\simeq 4700$. After sky subtraction and flux calibration, the mean continuum-normalized transmission over $8450$–$8710\text{A˚}$ was $0.0038\pm 0.0026$, consistent with zero; this corresponds to absorption redshifts approximately $5.95$–$6.16$. They identified a Gunn-Peterson trough. Crucially, resonant Lyman-alpha absorption is so sensitive that even a small neutral fraction can make transmission undetectable: the trough alone does not establish a mostly neutral intergalactic medium or uniquely date reionization. ^[ref].

What a spectrum constrains is a radiative transfer problem. For an ideal foreground screen with negligible emission into the observed beam, the intrinsic continuum-normalized transmission isHere $F$ is spectral radiative flux, $F_{\mathrm{cont}}$ is its unabsorbed continuum, $\tau$ is optical depth, $n_{\mathrm{HI}}$ is neutral hydrogen number density, and $N_{\mathrm{HI}}$ is neutral hydrogen column density along proper path length $\ell$. A Voigt profile model constrains redshift, column density and a velocity width together. A nearly black line has lost much of its sensitivity to further increases in optical depth; blending and an uncertain continuum can also admit different fitted components. Finite spectral resolving power and noise limit which widths and weak features can be measured. ^[ref]. ^[ref].

For a purely thermal Doppler core, use the velocity convention $\tau (v)\propto \exp [-(v-v_0{)}^2/b^2]$. The Maxwell-Boltzmann velocity distribution then giveswhere $k_B$ is the Boltzmann constant, $T$ is temperature and $m_{\mathrm{H}}$ is a hydrogen atom's mass. This $b$ is $\sqrt{2}$ times the Gaussian distribution's standard deviation. Higher temperature spreads the same column density over more velocities. An instrument also broadens a line: locally, for a nearly flat continuum, the recorded transmission is approximately the convolution of $\mathcal{T}$ with the normalized instrumental response, not the exponential of a convolved optical depth. A Gaussian response with full width at half maximum $c/R$ has standard deviation $c/[R\sqrt{8\log 2}]$, where $c$ is the speed of light and $R$ the spectral resolving power.

The upper model fixes the integral $\int \tau (v)dv$, hence the neutral hydrogen column density for this transition, while changing temperature. A cooler line blurred by an instrument resembles a warmer intrinsic line. The lower model has two equal Doppler cores separated by $24\mathrm{km}{\mathrm{s}}^{-1}$: reducing spectral resolving power merges them into one apparent feature. Both panels use invented profiles, Gaussian instrumental responses, and no noise or damping wings. They illustrate an inverse problem, rather than measurements of temperature or a universal relation between a dip and a cloud.

Better observations provide the inputs to that inverse problem. Kim, Cristiani and D'Odorico (2001) analyzed three quasars with the VLT's UVES spectrograph. For HE2217–2818, Table 1 records $z_Q=2.413$ and $16200\mathrm{s}$ of exposure in the $3050$–$3860\text{A˚}$ setup in October 1999. Isolated thorium-argon lamp lines measured $R\simeq 45000$; the data were resampled into $0.05\text{A˚}$ bins, with signal-to-noise ratios about $45$–$50$ per pixel over $3550$–$4050\text{A˚}$. They fitted a local polynomial continuum, fitted Voigt profiles, and compared statistics of the transmitted flux. Their analysis excluded the region near the quasar to avoid the quasar proximity effect and used the interval between its Lyman-beta and Lyman-alpha emission to reduce confusion with other Lyman-series absorption. Features below $N_{\mathrm{HI}}=10^{12.5}\mathrm{c}{\mathrm{m}}^{-2}$ became confused with noise, setting their adopted lower column density limit. Figure 1 shows the measured spectrum, fitted profiles and residuals. This was a calibrated absorption measurement, with explicit selection limits. ^[ref].

Even after instrumental correction, a width need not be purely thermal Doppler broadening. Hubble flow across extended gas, peculiar velocity structure and unresolved blends can contribute. If an isolated line's intrinsic width obeys $b^2=b_{\mathrm{th}}^2+b_{\mathrm{nt}}^2$ with an independent Gaussian nonthermal contribution $b_{\mathrm{nt}}^2\ge 0$, then $T\le m_{\mathrm{H}}{b}^2/(2k_B)$; $b=25\mathrm{km}{\mathrm{s}}^{-1}$ gives a computed upper bound about $3.79\times 10^4\mathrm{K}$. A complex blend need not satisfy this single-component model. Noise can also bias the fitted widths of weak lines. This is why the thermal cutoff of Lyman-alpha line widths, calibrated against gas models, is more informative than assigning a temperature to every dip. ^[ref].

Finally, neutral hydrogen column density is not total gas mass density. In a pure hydrogen approximation, with gas optically thin to ionizing radiation and nearly fully ionized in photoionization equilibrium, the Lyman-alpha forest photoionization scaling giveswhere $n_{\mathrm{H}}$ counts all hydrogen nuclei, $\Gamma$ is the photoionization rate, $\alpha (T)$ is the radiative recombination coefficient and $x_{\mathrm{HI}}=n_{\mathrm{HI}}/n_{\mathrm{H}}$ is the neutral fraction. At fixed temperature and $\Gamma$, doubling $n_{\mathrm{H}}$ quadruples $n_{\mathrm{HI}}$; doubling $\Gamma$ halves it. These computed examples explain both the sensitivity to density contrast and the ambiguity: stronger absorption can mean denser gas, weaker ionizing radiation, a different temperature, or more path length. The Jeans-length ionization correction for the Lyman-alpha forest adds a size model rather than an independently measured depth. Optically thick Lyman limit systems and damped Lyman-alpha systems require different treatment. ^[ref]. ^[ref].

The beauty of the Lyman-alpha forest is that a spectrum converts otherwise faint foreground gas into a sequence of markers through cosmic time. With a model of radiative transfer, velocities, ionizing radiation and the temperature-density relation of photoionized intergalactic gas, its transmission statistics probe the cosmic web and gas thermal history. A single line of sight remains one narrow sample, not a unique three-dimensional density map. ^[ref].

For a small, uniform surface at temperature $T$, enclosed by much larger black walls all at temperature $T_s$, the outward net heat transfer by thermal radiation isAssume an opaque grey body that is a diffuse thermal emitter, with the same emissivity $ϵ$ and radiometric absorptance at the surface's current temperature. Its emitting area $A$ sees only the enclosure: the view factor to the walls is one, with no self-viewing surface, opening or intervening obstruction. The intervening medium must neither absorb nor emit appreciable thermal radiation; a vacuum supplies a useful idealization. The enclosure is maintained at uniform $T_s$. The small-body assumption makes that reservoir description physically reasonable; genuinely black isothermal walls supply the incoming irradiance even without an infinite area ratio. Positive ${\dot{Q}}_{\mathrm{rad}}$ means net loss by the surface; negative means net gain. This expression counts thermal radiation alone, not all heat transfer.^[ref]

The reasoning is a two-way energy balance. Black walls supply irradiance $E=\sigma T_s^4$, so the surface absorbs $P_{\mathrm{abs}}=ϵAE$ and emits $P_{\mathrm{emit}}=ϵA\sigma T^4$. Subtract these to obtain ${\dot{Q}}_{\mathrm{rad}}$. For an opaque surface the remaining fraction $1-ϵ$ of incoming radiant flux is reflected; counting emitted plus reflected outgoing radiant flux minus total incoming radiant flux gives the same result. Kirchhoff's law of thermal radiation justifies equal spectral absorption and emission factors; the grey body assumption lets one factor be taken outside both spectral integrals.

Choose $ϵ=0.8$, $A=0.01{\mathrm{m}}^2$, $T=500\mathrm{K}$ and $T_s=300\mathrm{K}$. These are stipulated model inputs, not measured material properties. Using $\sigma \approx 5.670374419\times 10^{-8}\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-4}$ givesThe incoming irradiance is $459.300\mathrm{W}{\mathrm{m}}^{-2}$, corresponding to $4.5930\mathrm{W}$ incident on the surface; only $80\%$ is absorbed. At $T=T_s=300\mathrm{K}$, emitted and absorbed powers are both $3.6744\mathrm{W}$, so radiative equilibrium has zero net exchange despite continuing emission. Changing only emissivity to $0.2$ reduces both one-way powers and the net loss by a factor of four, giving $6.1694\mathrm{W}$ net. A low-emissivity coating therefore suppresses both radiative cooling to colder walls and radiative heating from hotter walls.

Near equal temperatures, write $T=T_s+\Delta T$. Taylor's theorem givesFor constant emissivity, the radiative heat transfer coefficient $h_{\mathrm{rad}}$ is approximately $4.8992\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-1}$ at $T_s=300\mathrm{K}$ and $ϵ=0.8$. This local linear approximation is useful for small temperature differences; the $200\mathrm{K}$ difference in the worked example needs the full fourth-power expression.

A cooling experiment requires a wider first law of thermodynamics balance. In a uniform-temperature model without a heater,where $C(T)$ is the sample's heat capacity and the last two terms are signed outward convection and thermal conduction losses. A vacuum suppresses convection but does not eliminate thermal conduction through supports or sensor leads. A cooling slope therefore cannot establish emissivity unless these paths and heat capacity are independently constrained. Multiple wall temperatures, incomplete view factors, wavelength-selective materials or appreciable intervening absorption require a more complete radiative transfer calculation.

When radiative rates compete with or dominate collisions, level populations need statistical equilibrium coupled to radiative transfer, rather than thermal population laws at the local temperature. Radiation from other depths can make the local radiation field differ from the Planck function. This changes ionization and spectral-line strengths, so a geometrically thin atmosphere may still require non-LTE abundance modelling.

The Rosseland mean opacity is the transport opacity appropriate to an optically thick region close to local thermodynamic equilibrium. If ${\kappa }_{\nu }$ is monochromatic transport opacity and $B_{\nu }(T)$ the Planck function, thenIt is a weighted harmonic mean: low-opacity frequency windows carry much of the diffusive flux. One must combine the monochromatic processes before averaging, rather than adding their separately averaged opacities.

In hot fully ionized material, electron-scattering opacity supplies an almost density- and temperature-independent contribution at fixed composition. In the nonrelativistic Thomson scattering limit, ${\kappa }_{\mathrm{es}}={\sigma }_T{n}_e/\rho \simeq 0.2(1+X)\mathrm{c}{\mathrm{m}}^2{\mathrm{g}}^{-1}$. Free-free opacity, the absorption counterpart of thermal bremsstrahlung, comes from photon absorption during encounters between free electrons and ions. Its approximate mean obeys Kramers' opacity law, ${\kappa }_{\mathrm{ff}}\propto \rho T^{-7/2}$, with composition and quantum correction factors.

When ions retain bound electrons, bound-free absorption is photoionization; edges, occupations and ionization fractions strongly affect its frequency dependence. Bound-bound absorption produces absorption lines. Numerous metal lines can close low-opacity windows and produce important opacity features in partially ionized layers. Line populations depend on excitation as well as ionization, while Doppler broadening, pressure broadening and overlapping lines change the transport windows. Stimulated inverse processes must be included consistently with detailed balance and the local thermal radiation field.

In cooler hydrogen-rich atmospheres, negative hydrogen ion opacity from both bound-free and free-free processes is often important. Its strong temperature dependence helps determine convective surface relations. In still cooler layers, molecular bands, and at sufficiently low temperatures dust absorption/scattering, may be important. These processes cannot be represented by one Kramers power law over the whole density-temperature plane.

An opacity calculation therefore needs the chemical abundances, the equation of state, electron density and bound-state populations. Saha ionization equation and thermal excitation populations give useful dilute LTE estimates; high-density pressure ionization, nonideal effects and altered level populations require corrections. Composition matters even when one writes only $\kappa (\rho ,T)$: hydrogen, helium and the metal abundance pattern are additional fixed inputs.

In dense degenerate cores, thermal conduction by electrons can dominate energy transport. An equivalent conductive opacity may be introduced so that $1/{\kappa }_{\mathrm{eff}}=1/{\kappa }_R+1/{\kappa }_{\mathrm{cond}}$, because the two diffusive fluxes add; conduction is not itself a photon absorption process. Very hot regimes also require departures from the simple nonrelativistic scattering approximation. Near the optically thin photosphere, a Rosseland diffusion description alone is inadequate, and frequency-dependent radiative transfer with appropriate boundary conditions is required. These regime changes are why realistic stellar opacity tables, rather than one global power law, enter stellar structure.

Use absorption distance $X$ and its distribution $f_X$. Summing hydrogen columns per unit $X$ gives $\int N{f}_{X}dN$. The physical light-path length per $X$ isThe mean neutral-hydrogen number density is therefore $(H_0/c)(1+z{)}^3\int N{f}_{X}dN$. Divide by $(1+z{)}^3$ to obtain its comoving density, multiply by hydrogen atomic mass, and divide by today's mass critical density ${\rho }_{\mathrm{crit},0}=3H_0^2/(8\pi G)$. This derives the neutral gas density from an absorption distribution:If “neutral gas” includes helium, replace $m_H$ by ${\mu }_H{m}_H$, with a stated abundance correction approximately ${\mu }_H\simeq 1.3$. Equivalently, for the distribution per redshift the prefactor inside the integral is $m_{H}H(z)/[c{\rho }_{\mathrm{crit},0}(1+z{)}^2]$ multiplying $N{f}_z$. In the PDF's energy-density convention, both the mass density and critical density are multiplied by $c^2$, leaving their ratio unchanged. The usual quoted abundance uses comoving mass density relative to today's critical density; a ratio of physical density to the epoch's critical density instead multiplies this answer by $(1+z{)}^3{H}_0^2/H(z{)}^2$.

To infer total gas, one needs the neutral fraction $x_{\mathrm{HI}}=N_{\mathrm{HI}}/N_H$ for each absorber. The ionized gas density from an absorption distribution isThis is an ionization-corrected integral, not an inference available from the observed columns alone. For optically thin, almost fully ionized hydrogen, photoionization equilibrium gives $\Gamma n_{\mathrm{HI}}=\alpha (T)n_e{n}_{\mathrm{HII}}\simeq \alpha (T)n_H^2$. If the absorbing depth $L$ is specified, $N=\alpha (T)n_H^{2}L/\Gamma$ andThus photoionization rate, temperature and depth are all needed.

One useful physical closure sets $L$ to a local Jeans length, $L\propto (f_{g}T/n_H{)}^{1/2}$, with local gas fraction $f_g$. Solving $N\propto [\alpha (T)/\Gamma ](f_{g}T{)}^{1/2}{n}_H^{3/2}$ then gives the Jeans-length ionization correction for the Lyman-alpha forest:This illustrates why weak, highly ionized forest lines can dominate the total gas inventory even though damped Lyman-alpha systems dominate the neutral-column-weighted inventory. It is a model for locally confined optically thin gas; optically thick Lyman limit systems and damped Lyman-alpha systems require a separate radiative transfer treatment. The physical basis and observational application of this closure are developed in Schaye's absorption-based matter-density calculation.

The ionizing background is produced mainly by hot stars in star-forming galaxies and by accretion-powered quasars. The stellar contribution depends on star formation and the ionizing photon escape fraction; the quasar contribution depends on its evolving luminosity function. Radiative transfer through the intergalactic medium, including absorption, recombination emission, and cosmological redshifting, determines the radiation actually available for photoionization.

In the picture relevant to this 2009 paper, increasingly rare quasars cannot by themselves explain the high-redshift hydrogen ionizing field, so galaxies become especially important above $z\sim 3$. Harder quasar radiation is important for doubly ionizing helium around $z\sim 3$. Hydrogen reionization must be largely accomplished by $z\sim 6$. At lower redshift the evolution reflects declining source emissivity together with the increasing propagation distance of ionizing photons; near reionization, short absorption lengths and spatially uneven ionized regions prevent a simple uniform-background description. The background rate need not follow the quasar population alone: a roughly flat hydrogen photoionization rate over $2\lesssim z\lesssim 4$ is compatible with a changing mixture of galaxies and quasars. Faucher-Giguère and collaborators' 2009 background calculation and their preceding opacity analysis give this interpretation. These redshifts describe the broad historical picture, rather than precise universal transition times.

The case for anthropogenic global warming rests on a physically consistent combination of forcing, energy storage and observed response, rather than on the claim that every unusual event has the same cause. For this examination, the observational comparison is the evidence available around 2010, and the two forecast horizons are approximately 2060 and 2110. Separating the robust mechanisms from uncertain regional details is essential to a coherent attribution of recent climate change.

The greenhouse effect follows from radiative transfer. Sunlight supplies energy, while the Earth loses energy mainly through long-wave radiation. Carbon dioxide, methane and other greenhouse gases absorb and emit in parts of that outgoing spectrum. Increasing their concentrations alters the levels and temperatures from which radiation escapes, initially reducing the net loss of energy for a given surface state. The surface and atmosphere must then warm until outgoing radiation again balances the input. Strong absorption near the centres of some bands does not eliminate the influence of additional gas: absorption in the wings and the changing height of emission still matter. For scale, the approximate carbon dioxide forcing law $\Delta F\simeq 5.35\log (C/C_0)\mathrm{W}{\mathrm{m}}^{-2}$ gives about $3.7\mathrm{W}{\mathrm{m}}^{-2}$ for doubling. It describes a forcing, not the ultimate temperature response, which additionally involves climate feedbacks and ocean heat uptake.

Measured increases in atmospheric carbon dioxide are consistent with fossil-fuel combustion and changes in land use. The carbon isotope signature and atmospheric oxygen consumption provide independent evidence of the combustion source, while the ocean and land absorb part of the emitted carbon. Methane and other climatically active gases add to the forcing. Aerosols exert partially compensating effects, so net human forcing is smaller than the sum of the warming gases alone. Uncertainty in that compensation affects the quantitative rate of warming without removing the greenhouse mechanism.

A sustained rise in global temperature is accompanied by increased ocean heat content, glacier retreat and sea level rise. This combination is harder to explain as an artefact of a few surface thermometers than a temperature series in isolation: the ocean is a large energy reservoir, and thermal expansion and melting have independent observational consequences. Internal variations of ocean circulation can move heat between regions and alter warming over particular decades, but cannot indefinitely supply a positive global energy accumulation without a change in the energy budget. Solar changes and volcanic eruptions also influence climate, and must be included in attribution; they do not reproduce the full recent warming pattern and sustained heat gain. Volcanic aerosol episodes chiefly produce temporary cooling. Broad troposphere warming together with stratosphere cooling is another useful fingerprint: additional greenhouse gases and ozone changes have a different vertical response from simply increasing solar input. No one fingerprint is unique to one cause, but their joint agreement with greenhouse forcing strengthens the attribution.

Polar observations are especially informative because they connect atmosphere, ocean and ice. In the Arctic, reduced summer sea ice extent, thinning and loss of multi-year sea ice indicate a diminished ice reservoir. Warmer air, earlier melt and delayed freeze-up contribute, while wind and current changes redistribute ice and export thick floes. Pacific heat inflow and changes in Atlantic-origin waters supply oceanic routes to melting. The protection afforded by the cold halocline means that warmer Atlantic water does not automatically heat the ice: the intervening mixing and density stratification must be considered. The 2007 retreat illustrates how long-term thinning, anomalous winds and greater solar absorption can combine. It neither proves that weather variability has ceased nor provides a reason to ignore the longer-term energy imbalance.

Loss of grounded ice has a different global implication. Increased surface melt and runoff on Greenland, and changes in outlet-glacier discharge, transfer water from land to the ocean. In Antarctica, warming of the Peninsula and disintegration of some ice shelves demonstrate strong regional change. Warm ocean water reaching shelf cavities can thin floating shelves from below; loss of ice-shelf buttressing can then speed grounded-ice discharge. The floating shelf's disappearance is not itself equivalent to adding its whole volume to sea level. The important subsequent response is increased loss from grounded ice sheets. Increased snowfall can compensate some losses, so a mass budget must distinguish accumulation, melting and ice flow.

The two polar regions should not be assumed to behave identically. Around 2010, Antarctic sea-ice changes included strong regional contrasts and an overall satellite-era behaviour quite different from Arctic summer retreat. Winds, ozone-related atmospheric changes, freshwater effects and the circulation of the Southern Ocean influence the Antarctic pattern. A large cold continent surrounded by an ocean is physically different from an ocean surrounded by land, and vigorous Southern Ocean heat uptake can delay surface warming. A regional increase in Antarctic sea ice is therefore compatible with global ocean heat gain and losses from some grounded Antarctic ice; it still needs its own dynamical and thermodynamic explanation. Credible attribution accounts for these differences instead of requiring uniform warming everywhere.

Several mechanisms produce polar amplification. The ice-albedo feedback is the most direct: retreating snow and ice expose darker surfaces that absorb more sunlight, causing further warming and melt. Much of the extra energy over newly open Arctic water is stored during summer and released to the atmosphere in autumn and winter, when an ice-covered ocean would be better insulated. Consequently the largest near-surface atmospheric response need not occur in the season of greatest solar input. Sea-ice loss changes both reflection and the timing of ocean heat release.

The atmospheric vertical structure also matters. In the tropics, convective adjustment often produces appreciable warming aloft, which increases radiative cooling. Stable polar boundary layers can concentrate warming near the surface, producing a different lapse-rate feedback and less compensating cooling from upper layers. Furthermore, the idealized Planck feedback scales as $4\sigma T^3$: colder emitting regions have a smaller increase of blackbody radiation per kelvin, although a quantitative atmospheric comparison needs full radiative transfer. Water-vapour and cloud changes modify both short-wave and long-wave fluxes, with cloud effects depending on season and surface properties. Changes in poleward heat and moisture transport provide further coupling. Thus polar amplification has several interacting causes; it is not simply the same local ice-albedo process repeated throughout the year.

Polar change affects the rest of the planet. Grounded-ice loss contributes to eustatic sea level change, while ocean warming contributes to thermosteric sea level change; together they increase coastal flooding and erosion risks. Redistribution of ice and freshwater changes density stratification and dense-water formation, potentially altering the Atlantic meridional overturning circulation and other oceanic heat and carbon transports. A weakening is physically plausible under warming and freshening, but an immediate complete shutdown is not a necessary prediction. Permafrost carbon feedback can release additional carbon dioxide and methane, increasing warming outside the Arctic. Habitat loss, changing primary production and altered food webs affect marine ecosystems and people, while a longer open-water season changes navigation, coastal exposure and resource access. These impacts are not exhausted by a single global-mean temperature number.

Forecasts require an emissions scenario as well as physical understanding. On the roughly fifty-year horizon, continued emissions imply further ocean and atmospheric warming, a longer Arctic melt season, less late-summer ice and a greater fraction of young seasonal ice. Mountain glaciers generally continue to lose mass, and sea level continues rising through thermal expansion and land-ice loss. Variability can interrupt a regional trend for several years, and the exact first nearly ice-free Arctic summer cannot be fixed from the 2007 minimum. Ocean heat storage and slow ice-sheet adjustment also mean that some further change is committed even if emissions are reduced.

By approximately 2110, sustained high emissions would imply a substantially warmer world, potentially several degrees above early twenty-first-century global temperatures, with stronger Arctic warming and likely recurrent nearly ice-free late summers. This is a seasonal prediction, not a claim that the entire Arctic Ocean must be ice-free in winter. Losses of grounded ice and thermal expansion would produce larger cumulative sea level rise, plausibly measured in several decimetres and potentially more if ice dynamics accelerate. The largest uncertainty is not whether adding heat expands water, but the future emissions trajectory and the response of major ice sheets. Complete loss of Greenland on this century-long horizon is not required by that prediction. With strong mitigation, substantially smaller late-century warming and ice loss are possible, although sea level can continue to rise beyond stabilization of surface temperature.

Anthropogenic forcing explains the broad warming signal; feedbacks amplify polar change, while emissions and ice dynamics govern future outcomes. This distinction permits both a firm physical conclusion and honest uncertainty about dates, regional details and the magnitude of rapid ice-sheet change.

In local thermodynamic equilibrium, thermal intensity samples the Planck function near an optical depth of order unity. The Eddington-Barbier relation makes this explicit: $I_{\lambda }(0,\mu )\simeq B_{\lambda }[T({\tau }_{\lambda }=\mu )]$. A molecular band has greater opacity than its adjacent continuum and therefore samples a higher layer. A band in emission relative to the continuum implies that this higher layer is hotter: the line-forming region has an atmospheric thermal inversion under the assumed LTE, thermal interpretation.

The continuum is thermal radiation from an optically thick, deeper photosphere, with comparatively smooth opacity. In an H/He hot Jupiter, collision-induced absorption by H2-H2 and H2-He collisions supplies an important continuum; weak overlapping molecular lines and opaque exoplanet clouds can contribute too. It is not a separate blackbody emitter floating above the gas. A strongly isothermal layer would erase LTE molecular contrast rather than generate emission peaks.

In the emitting inversion, $dT/dz>0$, whereas the dry adiabatic lapse rate has $dT/dz=-g/C_p<0$. Thereforewhich lies on the stable side of the Schwarzschild criterion. An upward-displaced parcel cools and becomes denser than the ambient hot upper gas. The region can thus carry and redistribute thermal energy by radiative transfer, not by unstable thermal convection. Winds may transport energy horizontally; stability rules out the specified buoyant vertical convection, not every possible motion.

For typical well-mixed structures, the dominant processes are:

Representative temperatures must specify the level: Earth has about $288\mathrm{K}$ at the surface (about $255\mathrm{K}$ effective emission temperature); Jupiter has about $165\mathrm{K}$ near one bar (about $125\mathrm{K}$ effective temperature); hot Jupiters commonly have photospheric temperatures of order $1000$–$2500\mathrm{K}$; and the Sun's photosphere is about $5800\mathrm{K}$. Upper layers, nightsides, deep interiors and the solar corona have different temperatures. These are characteristic values, not constant temperatures throughout each atmosphere.

The pressure axis increases downward in the sketches. Solid red segments indicate efficient convection; blue segments indicate mainly radiative transfer. The boundaries are nominal examples, since radiative-convective boundary pressure depends on opacity, gravity and intrinsic flux.

For Earth, the convective troposphere lies beneath a radiative stratosphere, with the transition near $0.1$--$0.2$ bar. The larger solar-system planets likewise have deep convective regions beneath largely radiative upper atmospheres; their tropopause/upper radiative-convective boundary is commonly of order $0.1$--$1$ bar. Their temperatures differ greatly from Earth's, and detached radiative layers can occur deeper down.

For an irradiated hot Jupiter, absorbed stellar flux maintains a hot, extended, relatively shallow-gradient radiative atmosphere. In models with a weak old-planet intrinsic flux, the deep radiative-convective boundary can lie around $10^2$--$10^3$ bar, with a useful wider model-dependent range of tens to thousands of bars. A thermal inversion can appear at low pressures if visible/UV absorption heats the upper atmosphere. Below the deep boundary the profile joins a convective adiabat.

A young, directly imaged exoplanet on a distant orbit is primarily heated from within. Its photosphere commonly joins the convective interior at order $0.1$--$10$ bar, illustrated here at one bar. The temperature generally rises monotonically with pressure over the infrared-forming layers. Strong stellar heating is absent, so a broad stellar-heated isothermal layer is unnecessary.

Two differences in thermal inversions are their absorbers and their formation conditions. Earth's ozone layer and solar-system hydrocarbon absorption can heat upper layers; the proposed hot-Jupiter absorbers include refractory titanium monoxide/vanadium monoxide or other strong visible absorbers, since ozone and methane-rich cold-planet chemistry are unsuitable at very high temperatures. Hot-Jupiter inversions depend strongly on irradiation, atmospheric cold traps and atmospheric transport, and can occur at mbar-to-sub-bar pressures; Earth and the solar-system giants have cooler, established stratospheres above their shallow tropospheres.

Two differences between the hot Jupiter and distant young-giant profiles are external versus internal heating, and the depth/shape of the radiative zone. The former can have a broad warm radiative layer, a much deeper convective boundary and sometimes an inversion; the latter typically has a steeper outward decrease toward a photosphere, a shallower boundary and no irradiation-driven inversion. These are class trends, not a unique temperature profile for every planet.

The three basic modes are radiation, convection and thermal conduction. Radiative transfer dominates transparent atmospheric layers; convection carries heat through sufficiently unstable fluid regions; thermal conduction matters especially in solids or highly conducting dense material.

Consider a chemically homogeneous ideal gas atmosphere and a small, dry parcel displacement, with rapid pressure equilibration and negligible heat exchange during the displacement. The parcel's first law of thermodynamics per unit mass isFor an adiabatic process, $ds=0$. Using hydrostatic equilibrium, $dP/dz=-\rho g$, the parcel followsAfter an upward displacement $\delta z>0$, its temperature relative to its new surroundings isHere $C_p$ is the specific heat capacity at constant pressure, and $N$ is the buoyancy frequency. At the same pressure, a warmer parcel has lower density and accelerates upward. Thus strict instability requires $dT/dz<-g/C_p$. Equivalently, the parcel equation is $\ddot{\delta z}=-N^2\delta z$, withNegative $N^2$ gives growing displacements. The printed non-strict inequality includes the onset boundary: equality is neutrally stable, not a strictly growing instability. This dry parcel buoyancy in a homogeneous atmosphere is the dry Schwarzschild criterion; condensation and composition gradients require modified criteria.

Three useful approaches are self-consistent one-dimensional models, atmospheric retrievals and three-dimensional circulation models.

A one-dimensional radiative-convective equilibrium model solves energy balance and radiative transfer, typically with a prescribed gravity, irradiation and composition or a chemical-equilibrium/kinetic calculation. It predicts physically coupled temperature profiles and spectra economically, but averages away longitude/latitude structure and inherits opacity, exoplanet cloud and chemistry assumptions.

An atmospheric retrieval parametrizes temperature, gas abundances, exoplanet clouds and other nuisance quantities and fits them to observations through a forward radiative-transfer model. It quantifies parameter constraints and degeneracies without enforcing every theoretical relation, but limited data can leave nonunique answers and priors/model structure can dominate.

A general circulation model solves rotating fluid dynamics with radiative heating/cooling and, when feasible, chemistry and exoplanet cloud transport. It predicts winds, phase curves and horizontal contrasts, but is computationally costly and still needs approximations for unresolved turbulence, exoplanet cloud microphysics and chemical networks. The approaches answer different questions and can be used jointly.