Climate change Created 2026-10-07 Updated 2026-10-18
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 model Created 2026-10-18 Updated 2026-10-20
A climate model is a mathematical description of the mechanisms and statistics of a planet's climate change and variability. A useful hierarchy starts with conservation of energy in a zero-dimensional planetary energy balance, adds vertical radiative transfer and convection in a radiative-convective climate column, and can resolve three-dimensional fluid mechanics and ocean circulation. Each level answers different questions. A column can explain why changing an absorber changes surface temperature; it cannot locate tomorrow's storm.
Three questions should be kept distinct: what sets an average temperature, how fast does it respond to radiative forcing, and how large are its unforced fluctuations? Climate sensitivity concerns the first, heat capacity controls response times, and the Hasselmann stochastic climate model illustrates the third. An observational test must compare a predicted observable with an instrument's actual readout, as in the satellite test of greenhouse spectral change.
Climate sensitivity Created 2026-10-18 Updated 2026-10-25
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 and constant forcing , the gain is , whereas the equilibrium warming is . 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.
The incremental radiative forcing per increment of a specified atmospheric gas burden. State whether the abundance unit is emitted mass, mole fraction or ppb, and specify the radiation-adjustment convention. For equal emitted masses in a fixed atmospheric air inventory, a per-ppb coefficient must be divided by the gas molar mass before comparing efficiencies. Spectral overlap with existing absorption makes efficiency depend on the background state.
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 , with an effective ocean mixed layer heat capacity per unit area :
Here is imposed radiative forcing, is the net restoring coefficient after the chosen climate feedbacks, and represents random net heat flux from unresolved weather. The anomaly energy per unit area is , 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. and have units , and has units .
First set , start from , and switch to constant . The integrating factor gives
For , and , the relaxation time is , about years using 365.25 days per year. At , , with eventual warming . Doubling doubles the response time but leaves this equilibrium unchanged. This separates thermal inertia from climate sensitivity.
Now set and idealize the weather forcing as Gaussian white noise with
The Dirac delta function has inverse-time units, so has units . The rigorous stochastic differential equation is
where is Brownian motion. This is an Ornstein-Uhlenbeck process. In stationary initialization its variance and autocovariance function are
Define the two-sided angular-frequency power spectral density by . The Fourier transform of the exponential autocovariance function, or the squared response to a sinusoidal input, gives
Thus fast fluctuations are suppressed; the climate output is red over the frequency range above the corner , with a low-frequency plateau. This is a colored noise spectrum without a periodic driver. The corner is in cycles per second, not . The continuous white input is only an approximation on times longer than weather's correlation time.
Figure 1.
Thermal inertia delays forced warming and filters random weather forcing
. Original one-box model: the step-forced mean approaches two kelvins, with one-standard-deviation spread from a fixed initial temperature and chosen white-noise strength. The right panel shows the analytic normalized two-sided spectrum, not measured climate data. Details.
For the same and , and the stationary standard deviation is . From fixed , the variance grows as , even when the forced expected value is predictable. Integrating over all angular frequencies returns , checking the normalization independently. If the unforced variance instead grows as ; if , 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 2026-10-13
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 on a newly exposed opaque patch, the extra absorbed heat flux is
For illustrative local surface albedos and , and , the patch absorbs an additional . 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 albedo
The exponential function makes a smooth transition from a bright cold state near to a darker warm state near . The midpoint and width 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 . For a small perturbation , linear stability follows from
Since , gives an asymptotically stable equilibrium, while gives an unstable equilibrium. With , the derivative is
Thus the optical term is destabilizing, while the thermal radiation term is stabilizing. Their sum, not the optical term alone, decides the linear stability.
For the same and Stefan-Boltzmann constant used above, the bisection method on the intervals , and , with temperature in kelvins, gives:
Equilibrium temperature ()Bond albedo ()Behavior
Asymptotically stable cold equilibrium
Unstable equilibrium
Asymptotically stable warm equilibrium
Unrounded numerical zeros of a function have ; the table rounds them for readability. In the positive-temperature domain the sign of 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 cools toward the cold state, whereas warms toward the warm state: this is bistability at identical external illumination.
Figure 1.
Two stable equilibria and one unstable equilibrium in an illustrative smooth ice-albedo model
. Left: net heating versus temperature at fixed illumination, with arrows pointing toward the stable zeros. Right: the equilibrium curve as illumination varies; solid branches are stable and the dashed branch is unstable. The dotted vertical line selects the same three equilibria shown on the left. Details.
The same model illustrates hysteresis. Its equilibrium curve can be written as . The two turning points, where , are saddle-node bifurcations: approximately and . For illumination strictly between these thresholds there are two stable equilibria and one unstable equilibrium. Under sufficiently slow increase of , 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 bounds the incoming heating, whereas outgoing thermal radiation grows as ; 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.
A blackbody absorbs all incident radiation and, in thermal equilibrium, emits the Planck law spectrum with integrated flux . Earth reflects some incoming sunlight, while its surface and atmosphere absorb and emit selectively with wavelength. With solar irradiance and Bond albedo A, the equilibrium planetary effective radiating temperature is
Reflection lowers this temperature relative to the zero-albedo value. The greenhouse effect makes the surface warmer than the effective emitting temperature: in a simple diagnostic , . This effective emissivity summarizes atmospheric transfer, not a universal material constant. At fixed absorbed sunlight, greenhouse warming does not require a different equilibrium outgoing flux. [ref]
The stated is positive net energy uptake, not necessarily the imposed radiative forcing itself. A linear climate forcing-response energy balance is
Here F is the perturbing forcing, is the net stabilizing radiative response and C is an effective areal heat capacity. The residual heats the Earth system, mainly through ocean heat uptake, with additional storage and melting. Finite heat capacity and ocean exchange delay equilibration. The imbalance alone determines neither F, the warming rate without C, nor the eventual equilibrium warming . Feedbacks change the response coefficient. The given number is treated as the examination's premise, not a present-day observational estimate. [ref]
The logarithmic forcing relation implies decreasing incremental forcing per additional ppm, explaining the two annual-increment results. It does not imply that a higher total carbon dioxide burden requires less compensation relative to a fixed original climate. The chosen reference and whether an end-of-year, mean-year or steady maintained aerosol burden is intended must remain explicit.
The numerical mass estimates use an unusually simple, stipulated constant albedo-per-mass calibration. Real particle size and lifetime change with loading, aerosol coagulation, injection altitude/latitude and circulation. Sulfate absorbs radiation as well as scattering it and changes stratospheric chemistry. A global-mean shortwave/longwave radiative forcing balance also need not reproduce regional temperature, rainfall or transient ocean heat uptake; it does not undo CO2-driven acidification. Chemistry–climate simulations explicitly resolve some of these competing aerosol-size, temperature and ozone effects.[ref]
For a liquid cloud at fixed water content and geometrical thickness, more cloud condensation nuclei can produce more, smaller droplets. The increased total extinction area raises optical depth and reflected sunlight: the Twomey effect gives negative shortwave radiative forcing. It is an indirect aerosol effect mediated by cloud microphysics, rather than direct aerosol scattering.[ref]
Smaller droplets can also delay collision-coalescence and drizzle, potentially increasing liquid-water retention, cloud lifetime or coverage; that adjustment would usually strengthen cooling for low reflective clouds. It is not universal: evaporation and mixing with dry air can reduce water content or cloud cover, and dynamical changes can oppose the simple lifetime argument. Absorbing particles can heat an atmospheric layer, sometimes evaporating cloud and producing a positive semi-direct forcing, but their vertical position can instead stabilize a cloud layer and reverse the response. Aerosols can additionally alter ice nucleation and mixed-phase precipitation, with no single forcing sign independent of cloud type. The net response combines shortwave reflection and longwave thermal effects, and depends on cloud height, optical depth and adjustments. Thus the fixed-water Twomey sign is robust within its assumptions; adjustment and ice-cloud signs require context.[ref]
The solar constant measures the incident power per area normal to the sunlight. A spherical Earth of radius intercepts sunlight over the disc , but radiates over its surface . For uniform surface albedo in the global energy-budget approximation, absorbed power is . 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 give
For a nonuniform surface the same power balance instead defines an effective planetary equilibrium temperature, ; 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 , so . The changed area is , giving
The initial planetary equilibrium temperature is . Thus
The differential estimate 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 and the given . 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 is
July through September has days, so the mean ocean heat flux to the ice is
Using instead gives about , within the intended accuracy. Around the date of this examination, the estimated net anthropogenic radiative forcing since preindustrial times was about ; doubling carbon dioxide gives about . 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 water column, take and liquid-water specific heat capacity . Assuming all the same energy remains in that uniformly mixed column,
If the water initially lies near the seawater freezing point, about , its final temperature is approximately . 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.
Denote aerosol optical depth by to distinguish it from cloud optical depth. Retrieve as in d, collocate it with , and estimate an aerosol-to-droplet susceptibility
A regression should control for water content, updraught and cloud regime; aerosol optical depth is not itself a direct measurement of activated nuclei. For an approximately constant over a factor two in aerosol, .
Integrate the supplied cloud albedo susceptibility rather than treating a finite doubling as infinitesimal:
For a small response, . If is the appropriate mean incoming shortwave flux and the fraction of the scene occupied by these clouds, the leading top-of-atmosphere perturbation, positive downward, is
The two transmissions account approximately for the atmosphere above the cloud and may be absorbed into an effective illuminating flux. In the simplest global transparent-overburden estimate . Positive gives increased reflection and negative shortwave forcing, the Twomey effect. Numerical forcing cannot be determined from the given susceptibility alone: it needs baseline albedo, cloud fraction, illumination and an empirically constrained . Humidity swelling, cloud contamination of aerosol retrievals and common meteorological influences can create spurious correlations. Changes in cloud amount, water path or lifetime are additional adjustments, not the fixed-water Twomey term derived here.
is the imposed radiative forcing; is the restoring change in outgoing radiation; and transfers heat from the surface layer to the deep ocean, whose equal and opposite gain appears in its own energy balance. The additional is a radiative adjustment correlated with the surface–deep-ocean temperature contrast, representing a pattern-dependent effective feedback rather than another conserved interlayer heat flux.
For constant forcing, final equilibrium has . During transient warming , so the radiative response inferred from a recent global energy budget is . If that response were mistaken for the final feedback, the inferred equilibrium climate sensitivity would be too small. Summing the two equations makes the distinction explicit: .
A linearly increasing radiative forcing produces a fast adjustment over about four years superposed on a slow deep-ocean response over about four hundred years. Once forcing stops increasing at year 70, surface warming continues and approaches the final equilibrium. The right panel of solution gives a model response normalized to the supplied endpoint and equilibrium.
Write after stabilization and , with . For a 70-year forcing ramp, the normalized step response is , . The ramp endpoint implies
This gives , and . Thus
at years respectively. The usual order estimate after the fast tail has largely disappeared is , with a still visible fast contribution near year 10.
To hold fixed, the deep layer obeys and the surface equation requires . Hence . Using gives and . With and initial forcing , the initial fractional decline is
The rate declines slowly as the ocean contrast disappears. A coarser estimate treating the post-fast warming as an effective equilibrium change gives per year; distinguishing the instantaneous surface feedback from the equilibrium feedback explains that difference.
Wording qualification: the printed text says a “1% increase … over 70 years”. The conventional transient climate response experiment uses per year, which nearly doubles concentration in 70 years. A literal total increase would not define the usual equilibrium climate sensitivity or transient climate response. The supplied and are used here as the endpoint normalizations, and a logarithmically uniform concentration growth gives the assumed linear forcing ramp.
The transient climate response to cumulative carbon emissions is , where is cumulative emitted carbon dioxide in a stated mass unit. Here the unit is tonnes of , rather than tonnes of carbon. Increasing concentration weakens the incremental radiative forcing per added atmospheric tonne because forcing is approximately logarithmic. Over a substantial range this is partly offset by a weakening of carbon sinks and by the changing ratio between transient warming and forcing as ocean heat uptake evolves. These approximate compensations make warming nearly proportional to cumulative emissions; no exact cancellation follows from the logarithmic forcing law alone. The approximation is range- and pathway-dependent, and is not a theorem valid for arbitrary emissions histories.
Sulfate aerosol 2026-10-20
Atmospheric particles containing sulfate, commonly in sulfuric-acid/water or neutralized-salt mixtures. Particle size, water uptake, acidity and accompanying material determine scattering, absorption, cloud activation and heterogeneous chemistry; sulfate mass alone does not specify radiative forcing.
Twomey effect 2026-10-20
Increasing cloud droplet number at fixed liquid-water content and cloud thickness reduces mean droplet size and increases total extinction area, raising shortwave optical depth and cloud albedo. This idealized response produces negative shortwave radiative forcing. Changes in cloud amount, water content and dynamics are separate adjustments.[ref]