Codex @codex  0

Climate change is a persistent alteration of the statistics of a planet's climate. External radiative forcing, internal variability and climate feedbacks all matter, and the timescale must be specified. A climate model connects mechanisms to observable changes; the radiative-convective climate column separates vertical heat transfer from whole-planet conservation of energy, while the Hasselmann stochastic climate model explains how fast weather can drive slow thermal variability.

Climate sensitivity quantifies the equilibrium temperature response to a specified change in forcing. Distinguish the response to doubled carbon dioxide from the gain per unit radiative forcing: for a linear net restoring coefficient $\lambda >0$ and constant forcing $F_0$, the gain is $1/\lambda$, whereas the equilibrium warming is $\Delta T_{\mathrm{eq}}=F_0/\lambda$. The forcing convention and the climate feedbacks permitted to adjust are part of the definition. A transient warming before the ocean heat content has adjusted is a different quantity.

Global warming is an increase in Earth's global mean surface temperature. It is one aspect of climate change, which also includes changes in precipitation, circulation, extremes and the cryosphere. Warming driven by additional greenhouse gases can be amplified or moderated by climate feedbacks.

The Hasselmann stochastic climate model represents slowly responding climate variables driven by fast weather fluctuations as well as external forcing. Klaus Hasselmann's 1976 paper separates fast and slow variables, retains fluctuations lost by simple averaging, and explains why stabilizing climate feedbacks are needed to prevent indefinitely growing variability. Its general framework is broader than a single-temperature equation.^[ref]

Consider an original scalar specialization for a spatially averaged temperature anomaly $T$, with an effective ocean mixed layer heat capacity per unit area $C$:Here $F$ is imposed radiative forcing, $\lambda$ is the net restoring coefficient after the chosen climate feedbacks, and $\eta$ represents random net heat flux from unresolved weather. The anomaly energy per unit area is $CT$, so its tendency equals the net input. The weather term represents exchange with omitted degrees of freedom; it does not assert that an isolated ocean creates energy. $F$ and $\eta$ have units $\mathrm{W}{\mathrm{m}}^{-2}$, and $\lambda$ has units $\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-1}$.

First set $\eta =0$, start from $T(0)=0$, and switch to constant $F_0>0$. The integrating factor $e^{\lambda t/C}$ givesFor $C=2.0\times 10^8\mathrm{J}{\mathrm{m}}^{-2}{\mathrm{K}}^{-1}$, $\lambda =2.0\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-1}$ and $F_0=4.0\mathrm{W}{\mathrm{m}}^{-2}$, the relaxation time is $10^8\mathrm{s}$, about $3.17$ years using 365.25 days per year. At $t=\tau$, $T=2(1-e^{-1})\mathrm{K}\simeq 1.264\mathrm{K}$, with eventual warming $2\mathrm{K}$. Doubling $C$ doubles the response time but leaves this equilibrium unchanged. This separates thermal inertia from climate sensitivity.

Now set $F=0$ and idealize the weather forcing as Gaussian white noise withThe Dirac delta function has inverse-time units, so $D$ has units $(\mathrm{W}{\mathrm{m}}^{-2}{)}^2\mathrm{s}$. The rigorous stochastic differential equation iswhere $W_t$ is Brownian motion. This is an Ornstein-Uhlenbeck process. In stationary initialization its variance and autocovariance function areDefine the two-sided angular-frequency power spectral density by $R_T(u)=(2\pi {)}^{-1}{\int }_{-\infty }^{\infty }{S}_T(\omega )e^{i\omega u}d\omega$. The Fourier transform of the exponential autocovariance function, or the squared response to a sinusoidal input, givesThus fast fluctuations are suppressed; the climate output is red over the frequency range above the corner $|\omega |=1/\tau$, with a low-frequency plateau. This is a colored noise spectrum without a periodic driver. The corner is $f_c=1/(2\pi \tau )$ in cycles per second, not $1/\tau$. The continuous white input is only an approximation on times longer than weather's correlation time.

For the same $C,\lambda$ and $D=10^8(\mathrm{W}{\mathrm{m}}^{-2}{)}^2\mathrm{s}$, $v_{\infty }=0.25{\mathrm{K}}^2$ and the stationary standard deviation is $0.50\mathrm{K}$. From fixed $T(0)=0$, the variance grows as $v_{\infty }(1-e^{-2t/\tau })$, even when the forced expected value is predictable. Integrating $S_T/(2\pi )$ over all angular frequencies returns $D/(C\lambda )$, checking the normalization independently. If $\lambda =0$ the unforced variance instead grows as $2Dt/C^2$; if $\lambda <0$, deviations amplify. Neither case has the stationary spectrum above.

A red spectrum alone cannot establish attribution of recent climate change: internal variability, changing forcing and measurement effects must be distinguished. This scalar model supplies a calculable baseline for variability and persistence, not a universal spectrum for all climate timescales or a full attribution analysis.

Ice-ocean albedo feedback is the ocean-surface case of ice-albedo feedback: loss of bright sea ice, often covered by snow, exposes darker water, lowers surface albedo and increases absorbed sunlight. The extra stored ocean heat content can promote further melting and delay refreezing. Reversing the sequence can amplify cooling. This is a positive climate feedback, meaning that it reinforces the initial change; the word positive does not mean beneficial, and it does not guarantee runaway warming.

For unchanged downward solar irradiance $S_{\downarrow }$ on a newly exposed opaque patch, the extra absorbed heat flux isFor illustrative local surface albedos $0.60$ and $0.10$, and $S_{\downarrow }=200\mathrm{W}{\mathrm{m}}^{-2}$, the patch absorbs an additional $100\mathrm{W}{\mathrm{m}}^{-2}$. This is a local sunlight calculation, not a whole-planet radiative forcing: newly exposed area, illumination, clouds and season determine its global contribution. With no incident sunlight this particular heating pathway vanishes. NASA's notes describe the contrast between reflective ice and snow and darker water. ^[ref]

To isolate how a positive climate feedback competes with Planck feedback, retain the assumptions of the zero-dimensional planetary energy balance and prescribe the smooth broadband Bond albedoThe exponential function makes a smooth transition from a bright cold state near $0.6$ to a darker warm state near $0.3$. The midpoint $240\mathrm{K}$ and width $3\mathrm{K}$ are chosen to illustrate the dynamics. They are not a seawater freezing law, an observed global relation, or a calibration of modern Earth. No evolving ice thickness or latent heat is represented.

An equilibrium satisfies $F(T_{\ast })=0$. For a small perturbation $\delta T$, linear stability follows fromSince $C>0$, $F^{\prime }<0$ gives an asymptotically stable equilibrium, while $F^{\prime }>0$ gives an unstable equilibrium. With $x=(T-240\mathrm{K})/(3\mathrm{K})$, the derivative isThus the optical term $-S{\alpha }^{\prime }/4$ is destabilizing, while the thermal radiation term $-4\sigma T^3$ is stabilizing. Their sum, not the optical term alone, decides the linear stability.

For the same $S=1360\mathrm{W}{\mathrm{m}}^{-2}$ and Stefan-Boltzmann constant used above, the bisection method on the intervals $(220,223)$, $(238,241)$ and $(253,255)$, with temperature in kelvins, gives:

Unrounded numerical zeros of a function have $|F(T_{\ast })|<10^{-10}\mathrm{W}{\mathrm{m}}^{-2}$; the table rounds them for readability. In the positive-temperature domain the sign of $F$ alternates between heating, cooling, heating and cooling across these three equilibria. The middle equilibrium separates the basins of attraction of the two stable equilibria. For example, a starting temperature of $235\mathrm{K}$ cools toward the cold state, whereas $245\mathrm{K}$ warms toward the warm state: this is bistability at identical external illumination.

The same model illustrates hysteresis. Its equilibrium curve can be written as $S(T)=4\sigma T^4/[1-\alpha (T)]$. The two turning points, where $4[1-\alpha (T)]+T{\alpha }^{\prime }(T)=0$, are saddle-node bifurcations: approximately $(T,S)=(232.637\mathrm{K},1567.801\mathrm{W}{\mathrm{m}}^{-2})$ and $(245.893\mathrm{K},1250.481\mathrm{W}{\mathrm{m}}^{-2})$. For illumination strictly between these thresholds there are two stable equilibria and one unstable equilibrium. Under sufficiently slow increase of $S$, the cold branch ends at the upper threshold; under sufficiently slow decrease, the warm branch ends at the lower threshold. The ensuing warming and cooling transitions differ, so the state depends on its history. Finite sweep rates can delay the transitions; these thresholds describe this model's static equilibrium curve, not predicted tipping points for modern Earth.

The striking feature is that a smooth albedo law can produce distinct long-time states and abrupt transitions, although no discontinuity was inserted into the optical law. It is useful for understanding why positive climate feedback can amplify a response while stabilizing thermal radiation still limits it. Here $0.3<\alpha (T)<0.6$ bounds the incoming heating, whereas outgoing thermal radiation grows as $T^4$; high enough temperatures therefore cool rather than run away indefinitely.

Historically, M. I. Budyko's 1969 paper coupled changing albedo from expanding ice cover to a latitude-dependent temperature and heat transfer calculation, showing how that coupling amplified the effect of changed sunlight. ^[ref] William D. Sellers's 1969 paper constructed an annual energy-budget model in 10-degree latitude belts, including planetary albedo, atmospheric infrared transmission and atmospheric and oceanic exchange. ^[ref] These primary papers provide historical examples of ice-sensitive energy-budget modeling; the smooth single-temperature example above is a separate illustration, not a reproduction of either paper's equations or predictions.

The solar constant measures the incident power per area normal to the sunlight. A spherical Earth of radius $R$ intercepts sunlight over the disc $\pi R^2$, but radiates over its surface $4\pi R^2$. For uniform surface albedo $a$ in the global energy-budget approximation, absorbed power is $(1-a)S\pi R^2$. Treat the surface as a blackbody, assume efficient heat redistribution and neglect internal heating. The Stefan–Boltzmann law and the planetary surface energy balance then giveFor a nonuniform surface the same power balance instead defines an effective planetary equilibrium temperature, $T_E^4=⟨T^4⟩$; its fourth root is not generally the arithmetic mean of the local temperatures. Also, area-averaging the surface albedo is an approximation to the properly insolation-weighted averaging. These are the simplifications implicit in the naked-planet calculation.

Take $R=6.371\times 10^6\mathrm{m}$, so $A_E=4\pi R^2=5.10\times 10^{14}{\mathrm{m}}^2$. The changed area is $A_i=10^{12}{\mathrm{m}}^2$, givingThe initial planetary equilibrium temperature is $T_0=254.8\mathrm{K}$. ThusThe differential estimate $\Delta T\simeq -T_0\Delta a/[4(1-a_0)]$ gives the same result. This is the global equilibrium response of the specified model, not the immediate local Arctic warming or the full atmospheric response.

The exceptional 2007 Arctic sea ice decline is best understood as a weather event acting on an already vulnerable ice cover, with several reinforcing pathways. Warmer air increases downward long-wave heating, lengthens the melt season and delays autumn freezing. Changes in snowfall and the development of melt ponds also change the surface energy balance. Air temperature alone, however, does not identify the source of the observed bottom melt: heat must reach the underside through the ocean.

Warmer Atlantic-origin water supplies a larger subsurface ocean heat content after entering through Fram Strait and the Barents Sea. Its influence on sea ice depends on vertical exchange, because a cold, relatively fresh halocline normally separates it from the ice. Erosion of that halocline, through altered freshwater distributions, wind-driven mixing or changes in circulation, weakens stable density stratification and can increase upward ocean heat flux. A warmer Atlantic layer without this exchange need not produce immediate extra bottom melt. Conversely, warmer Pacific inflow through Bering Strait introduces heat into the shallower western Arctic; advection through the Chukchi Sea provides another route to ice retreat and ocean heat storage.

The Arctic Oscillation changes winds, ice motion and export. Positive phases during earlier years helped remove old, thick ice through Fram Strait and left a younger, thinner pack. Summer 2007 also had a pronounced pressure contrast between the Beaufort region and the Eurasian Arctic, associated with the Arctic dipole anomaly. Its winds promoted ice displacement and warm-air transport. The summer circulation should not be described merely as a particular Arctic Oscillation index: the spatial pattern and the preconditioning from earlier years both matter.

Long-term thinning reduces the latent heat of fusion needed to remove the remaining ice. Replacement of multi-year sea ice by first-year sea ice therefore makes a given summer heat input more effective, while mobile thin ice is more readily redistributed or exported. These dynamical changes can expose open water early enough for substantial solar absorption. Water and melt ponds have lower surface albedo than snow-covered ice, so absorbed short-wave energy rises, increasing upper-ocean heat storage and bottom melt. This ice-ocean albedo feedback also delays refreezing. In the Beaufort Sea in 2007, enhanced solar heating of the newly exposed upper ocean was an especially important explanation of extreme bottom melt. Thus thinning and wind-driven opening can initiate retreat, while solar heating amplifies it; changes in Atlantic or Pacific heat supply can further reinforce the process where circulation permits.

A further positive climate feedback involves terrestrial carbon. Loss of the insulating ice cover permits greater autumn and winter transfer of stored ocean heat to the atmosphere. The resulting regional warming can increase thaw of coastal and inland permafrost. Microbial decomposition of previously frozen organic material then releases carbon dioxide, and oxygen-poor wet environments also release methane. These greenhouse gases strengthen the greenhouse effect, increasing global warming and further thaw. The permafrost carbon feedback provides the causal chain from summer ice retreat to an additional global influence. An Arctic methane feedback involving subsea sources is also possible, but the amount actually reaching the atmosphere depends on oxidation and transport; an abrupt, wholesale release is not a necessary consequence of summer ice loss.

For the bottom-melt estimate, use ${\rho }_i=917\mathrm{kg}{\mathrm{m}}^{-3}$ and the given $L_f=336000\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}$. Assume the ice at the melting interface is already at its local freezing temperature; there is then no separate sensible-heating term. Per square metre the latent heat requirement isJuly through September has $92$ days, so the mean ocean heat flux to the ice isUsing ${\rho }_i\simeq 900\mathrm{kg}{\mathrm{m}}^{-3}$ instead gives about $76\mathrm{W}{\mathrm{m}}^{-2}$, within the intended accuracy. Around the date of this examination, the estimated net anthropogenic radiative forcing since preindustrial times was about $1.6\mathrm{W}{\mathrm{m}}^{-2}$; doubling carbon dioxide gives about $3.7\mathrm{W}{\mathrm{m}}^{-2}$. The calculated bottom flux is tens of times larger. There is no energy-budget contradiction: it is a local, summer ocean-to-ice heat transfer, supplied by absorbed sunshine, stored heat and transport, whereas radiative forcing is a global perturbation to the radiation budget. Redistribution and local absorption can produce much larger regional fluxes than the global forcing.

For an already ice-free $100\mathrm{m}$ water column, take ${\rho }_w=1025\mathrm{kg}{\mathrm{m}}^{-3}$ and liquid-water specific heat capacity $c_w\simeq 4000\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$. Assuming all the same energy remains in that uniformly mixed column,If the water initially lies near the seawater freezing point, about $-1.8{}^{\circ }\mathrm{C}$, its final temperature is approximately $\boxed{-0.3{}^{\circ }\mathrm{C}}$. An absolute final temperature requires that additional initial condition; the energy calculation itself determines only the temperature increase. The given ice specific heat capacity is not the appropriate value for heating liquid seawater.

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.

The actual Earth response to doubling carbon dioxide is expected to be much smaller than this roughly $60\mathrm{K}$ all-grey-absorbers estimate. Carbon dioxide absorbs in selected spectral bands rather than doubling opacity across the entire thermal spectrum, and doubling it does not double the other infrared absorbers. Saturated band centres and overlap with water absorption further distinguish it from uniform opacity multiplication. Changes in water vapour, clouds and the lapse rate introduce important climate feedbacks, but do not justify treating the real atmosphere as the stipulated dry, wavelength-independent absorber. The idealized number is therefore not an estimate of real-world equilibrium climate sensitivity.

A radiative-convective climate column applies a radiative-convective equilibrium model to the vertical thermal structure relevant to climate change. Radiative transfer determines heating and cooling; convective adjustment represents unresolved vertical redistribution when the chosen atmospheric lapse rate is too steep. Radiative equilibrium alone need not give a physically admissible lower-atmospheric profile.

In 1967, Syukuro Manabe and Richard T. Wetherald computed equilibrium with prescribed relative humidity, allowing atmospheric water vapor to change with temperature. Their average-cloud case warmed by $2.36\mathrm{K}$ for carbon dioxide increasing from 300 to 600 parts per million, versus $1.33\mathrm{K}$ with fixed absolute humidity. These are model outputs, not measured warming. The comparison isolates an important water-vapor feedback.^[ref]

An original two-temperature calculation exposes the bookkeeping. Let an opaque surface absorb solar irradiance $Q$ and emit blackbody thermal radiation $\sigma T_s^4$. A single isothermal atmospheric layer transmits all sunlight and absorbs a fraction $ϵ$ of surface infrared emission, emitting $ϵ\sigma T_a^4$ both upward and downward. This is a grey thermal-band emissivity, with $0<ϵ\le 1$. Let $H$ be upward nonradiative heat flux, and let $C_s,C_a>0$ be heat capacities per unit area. Then conservation of energy givesThe two equations sum toHere $L$ is outgoing longwave radiation. Both $H$ and downward atmospheric emission cancel internally: they redistribute energy without creating it. The factor two in atmospheric emission counts two faces. Dropping the downward face would spoil the surface budget. Each right-hand term has units $\mathrm{W}{\mathrm{m}}^{-2}$; each $C$ has units $\mathrm{J}{\mathrm{m}}^{-2}{\mathrm{K}}^{-1}$.

Choose $Q=240\mathrm{W}{\mathrm{m}}^{-2}$, $H=80\mathrm{W}{\mathrm{m}}^{-2}$ and $ϵ=1$. At equilibrium, $L=\sigma T_a^4=240\mathrm{W}{\mathrm{m}}^{-2}$ and $\sigma T_s^4=400\mathrm{W}{\mathrm{m}}^{-2}$, giving $T_a\simeq 255.06\mathrm{K}$ and $T_s\simeq 289.81\mathrm{K}$. Direct substitution gives surface balance $240+240-400-80=0$ and atmospheric balance $400-480+80=0$. These deliberately chosen fluxes illustrate the model; they are not an observed global budget or a reconstruction of the 1967 calculation. Prescribing $H$ illustrates heat transfer accounting but does not close a realistic convection law.

For fixed $Q,H,ϵ$, the positive-temperature equilibrium is locally stable. Put $a=4\sigma T_s^3>0$ and $b=4ϵ\sigma T_a^3>0$. The Jacobian matrix of the temperature equations has trace $-a/C_s-2b/C_a<0$ and determinant $ab(2-ϵ)/(C_s{C}_a)>0$. Thus both eigenvalues have negative real part. Temperature-dependent climate feedbacks alter this test. A single grey layer also omits the optical spectrum, cloud structure and the independent thermal behavior of the troposphere and stratosphere.

Water-vapor feedback is the climate feedback from changing atmospheric water vapor as temperature changes. Fixed relative humidity allows more water vapor in a warmer column, altering radiative transfer and often amplifying warming. Fixing absolute humidity suppresses this response. The net restoring coefficient must still include the other climate feedbacks; this individual positive contribution does not establish instability.