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.
Greenhouse gas radiative efficiency 2026-10-20
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.
Hasselmann stochastic climate model 2026-10-18
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]
Stochastic climate models Part I. Theory by Klaus Hasselmann. Hasselmann, Stochastic climate models. Part I. Theory, 1976, Section 2; Section 7, summary points 1–4
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 givesFor , 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 withThe Dirac delta function has inverse-time units, so has units . The rigorous stochastic differential equation iswhere 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 . 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 , 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.
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 isFor 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]
https://science.nasa.gov/learn/heat/resource/explore-albedo-activity/. NASA, Explore Albedo Activity
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 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 fromSince , gives an asymptotically stable equilibrium, while gives an unstable equilibrium. With , the derivative isThus 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.
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]
https://cdn.repub.ch/s3/republik-assets/assets/pdf/budyko_1969.pdf. M. I. Budyko, The effect of solar radiation variations on the climate of the Earth, Tellus 21(5), 611–619, 1969, pp. 614–616
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]https://acoustique.ec-lyon.fr/chaos/Sellers69.pdf. William D. Sellers, A Global Climatic Model Based on the Energy Balance of the Earth-Atmosphere System, Journal of Applied Meteorology 8, 392–400, 1969
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. Past exam of the atmospheric physics course of the University of Cambridge 2012 Paper 1 1 a Solution 2026-10-20
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 isReflection 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]
Climate and Earth's Energy Budget by NASA Earth Observatory. Absorbed/reflected sunlight and effective radiating versus surface temperature
The stated is positive net energy uptake, not necessarily the imposed radiative forcing itself. A linear climate forcing-response energy balance isHere 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 Earth's Energy Budget, Climate Feedbacks, and Climate Sensitivity by Piers Forster and others. §7.4.1 andBox7.1, forcing–response budget; sign convention here uses positive damping lambda
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]
The impact of geoengineering aerosols on stratospheric temperature and ozone by P. Heckendorn and others. Environmental Research Letters 4 (2009), 045108, Abstract: coupled aerosol-size and chemistry–climate simulations
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]
The Influence of Pollution on the Shortwave Albedo of Clouds by S. Twomey. Journal of the Atmospheric Sciences 34 (1977), pp. 1149–1152, fixed-water-content droplet-number and reflectance effect
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]
Past exam of the mathematics course of the University of Cambridge 2010 iii Paper 69 1 Solution Created 2026-10-03 Updated 2026-10-14
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 giveFor 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 , givingThe initial planetary equilibrium temperature is . ThusThe 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 isJuly through September has days, so the mean ocean heat flux to the ice isUsing 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 susceptibilityA 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, isThe 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.
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.
The Earth’s Energy Budget, Climate Feedbacks, and Climate Sensitivity by Piers Forster and others 2026-10-20
Chapter 7 of Climate change 2021: The Physical Science Basis, pp. 923–1054. §§7.3.2–7.3.3 contrast greenhouse-gas and aerosol effective 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]
The Influence of Pollution on the Shortwave Albedo of Clouds by S. Twomey. Fixed-water-content droplet-number and albedo response

