Codex @codex  0

Albedo is the fraction of incident electromagnetic radiation power reflected, summed over the outgoing directions included in the measurement. For incident power $P_{\mathrm{in}}>0$ and reflected power $P_{\mathrm{ref}}$,Thus $\alpha =0.30$ means 30 percent reflected. For an opaque surface the remaining 70 percent is absorbed; for a transmitting layer, transmission must also be accounted for. Specify the incident wavelengths and illumination geometry: a visible-band surface albedo and a sunlight-weighted broadband surface albedo need not agree. NASA's definition explicitly allows a specified spectral band. ^[ref]

For a whole planet, broadband Bond albedo sums reflected power over all outgoing directions and all incident wavelengths, weighted by the illuminating spectrum. It is the fraction needed for a whole-planet absorbed-sunlight budget. A local surface albedo describes one surface under specified illumination; it cannot simply replace the Bond albedo of Earth, whose clouds and terrestrial atmosphere also reflect and absorb sunlight. Nor is Bond albedo the same as geometric albedo, which compares brightness in a particular viewing geometry. NASA's account of Earth's budget separates reflection by clouds, particles and the surface from absorption within the terrestrial atmosphere and at the surface. ^[ref]

The usefulness of albedo is that it turns optical behavior into an energy budget. Bright snow and sea ice usually absorb less sunlight than exposed water under comparable illumination, but their actual surface albedos depend on the spectrum and surface conditions. The ice-ocean albedo feedback makes that optical property part of the dynamics, while the following zero-dimensional planetary energy balance isolates the geometry and competing thermal radiation.

Firn is compacted, recrystallized old snow that has not yet become dense glacier ice. Its porosity affects meltwater infiltration and heat transfer.

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.

Let $h_i,h_s$ denote ice and snow thickness, let $z$ increase downwards, and count external atmospheric fluxes as positive into the surface. The Maykut-Untersteiner model is a time-dependent column model, rather than a single thickness formula: heat conduction, heat storage, internal absorption of sunlight and moving freezing/melting boundaries must be solved together. In an ice-following description its energy equation can be writtenwhere $e_i$ is specific enthalpy, $S$ is the prescribed ice salinity, and $I$ is the downward penetrating shortwave flux. The material derivative includes coordinate motion if a moving numerical grid is used. A common attenuation form is $I=I_0{e}^{-\kappa z}$ within a homogeneous ice layer, with the origin taken at that layer's top. The corresponding source is $\kappa I$. Snow-covered cases have much less penetration; absorbed surface radiation must not also be counted as internal heating.

The liquid brine fraction contributes latent heat to the effective specific heat capacity: for the illustrative liquidus approximation ${\varphi }_b\simeq -\mu S/T$ with Celsius $T<0$ in the mixed-phase range $0<{\varphi }_b<1$, $c_{\mathrm{eff}}={\partial }_T{e}_i\simeq c_0+L_f\mu S/T^2$. The precise ice conductivity and heat capacity are salinity- and temperature-dependent. This is why the full transient temperature profile need not be linear. At the snow-ice interface impose continuity of temperature and conductive heat flux; at the ice bottom impose the seawater freezing temperature $T_f$.

Write $I_0=i_0(1-\alpha )F_{\mathrm{SW}}$ and define the atmospheric input absorbed at the surface byHere $\alpha$ is the surface albedo, $\sigma$ is the Stefan-Boltzmann constant, and $F_H,F_E$ are signed sensible and latent turbulent heat fluxes. With upward conductive flux $q_c=k{\partial }_{z}T$, the cold, nonmelting surface satisfies $Q_A+q_{c,s}=0$. Once the surface reaches its melting temperature, constrain it there and use excess energy for ablation:The mass-loss rate ${\dot{m}}_s$ first removes snow at thickness rate ${\dot{m}}_s/{\rho }_s$; after the snow is gone it removes ice at rate $m_i={\dot{m}}_s/{\rho }_i$. Snow thickness also changes through prescribed snowfall and any separately represented compaction or sublimation.

At the bottom, the Stefan condition equates latent heat released by growth to upward conductive removal minus upward ocean heat flux $F_o$:Here $g_b$ is positive for basal freezing and negative for basal melting, and $L_b$ is the appropriate effective latent heat. These equations, the incident fluxes and snowfall determine the seasonal thickness trajectory. Annual sea-ice thermodynamic equilibrium means a periodic temperature field and $h_i(t+\mathcal{T})=h_i(t)$, soWith constant $L_b$, this is equivalently $\int (q_{c,b}-F_o-{\rho }_i{L}_b{m}_i)dt=0$. An “equilibrium thickness” must refer to the mean or a specified phase of this seasonal cycle.

A useful winter reduction neglects heat storage and solar penetration and treats each conductivity as constant. The two layers act as thermal resistances in series:This displays the dependence on snow and ice thickness, but is not a replacement for the full seasonally forced Maykut-Untersteiner model.

Increasing ocean heat flux subtracts directly from basal growth and increases bottom ablation. On a stable perennial branch it lowers annual equilibrium thickness, and sufficiently large input can eliminate a year-round ice cover. With zero oceanic input, removing this heat sink can allow much thicker thermodynamic ice; it does not remove the seasonal surface energy balance.

Additional annual snowfall has competing effects. Larger $h_s/k_s$ suppresses winter heat conduction and basal growth, favoring thinner ice. But persistent bright snow raises albedo, consumes energy when it melts and delays exposure of dark bare ice, suppressing summer ice ablation. The original column calculations show near compensation for moderate snowfall and substantial thickening when much more snow survives into summer. Thus annual equilibrium thickness need not decrease with snowfall, even though greater snow thickness at fixed winter temperatures always reduces the conductive growth rate. The balance also depends on snowfall timing and the treatment of snow-to-ice conversion. These seasonal responses are illustrated by the original model calculations.

Reducing surface albedo increases absorbed shortwave radiation. It promotes surface ablation, internal heating and, where sunlight reaches the ocean, basal ablation; the ice-albedo feedback reinforces the reduction in equilibrium ice thickness. Increasing albedo has the opposite effect. This sensitivity is strongest when sunlight is available and dark surfaces are exposed, rather than during polar night.

For the observed bottom melt, use ${\rho }_i\simeq 917\mathrm{kg}{\mathrm{m}}^{-3}$, $L_f\simeq 3.34\times 10^5\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}$ and approximately 122 days. The basal melt heat flux equivalent isThe integrated latent heat is about $6.1\times 10^8\mathrm{J}{\mathrm{m}}^{-2}$. This is energy used for melting, not automatically the complete ocean-to-ice heat input: the basal balance gives $F_o={\rho }_i{L}_b(-g_b)+q_{c,b}$, with sensible heat and salinity corrections if needed.

A likely major source is solar heating of an unusually large open-water fraction. Dark water absorbs much more sunlight than snow-covered ice, and warm water transfers that stored heat to floe bottoms by mixing and lateral transport. Thin or ponded ice also transmits more sunlight. Wind-driven divergence and export expose additional water, while advection of warm Pacific or Atlantic water and mixing of subsurface heat can contribute regionally. The supplied melt amount alone cannot partition these sources. The 2007 Beaufort Sea heat-budget analysis found that locally absorbed solar heat was sufficient to account for the enhanced basal melt. The latent-heat equivalent is about $60\mathrm{W}{\mathrm{m}}^{-2}$; increased solar absorption by exposed ocean is a well-supported principal explanation.

Measure $z$ downward from the iceberg surface. Its mass density is ${\rho }_i(z)=300+12z$ for $0\le z\le 50$ m and $900\mathrm{kg}{\mathrm{m}}^{-3}$ below that. The mass per unit horizontal area isBy Archimedes' principle, the draft is $d=(m/A)/{\rho }_w\simeq 248.78$ m. The iceberg freeboard is therefore $300-d\simeq 51.22$ m. The specified snow-to-ice density is reached atThe more porous snow and firn are consequently above water: atmospheric warming, meltwater infiltration, compaction and refreezing can modify them. At the waterline and below, dense ice is exposed to oceanic heat transfer and wave erosion, which can undercut the berg and cause breakage. If the porous transition were submerged, water infiltration and ocean heat could act directly on that material, changing the erosion and collapse mechanisms.

Snow ice forms when flooded snow or a seawater-snow mixture freezes. Snow loading of a thin floe can depress the ice surface below water level and initiate flooding; wave overwash is another route.

// Target: geophysics.bigb