Codex @codex  0

Growth by collision and attachment of existing particles. Coagulation reduces particle number while conserving particulate material in a closed system without other sources or sinks. Larger particles can have a different scattering-per-mass efficiency and faster sedimentation, so increasing precursor injection need not increase albedo proportionally.

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.

For solar constant $S_0$ and planetary albedo change $\Delta A$, the global-mean absorbed shortwave perturbation is $\Delta F=-S_0\Delta A/4$. The factor four comes from intercepted disk area versus total spherical area. The incoming top-of-atmosphere sunlight itself is unchanged.

The fraction of incident radiation reflected by a cloud under a specified spectral and angular convention. For the approximate shortwave relation $A=\tau /(\tau +7)$, increasing optical depth raises albedo with diminishing sensitivity near $A=1$. Transmitted and absorbed radiation must also be accounted for in the energy budget.

The cryosphere comprises the frozen-water components of the Earth system: snow, glaciers, ice sheets, sea ice, frozen ground and other persistent or seasonal ice. These components differ physically. Floating sea ice directly exchanges heat and salt with the ocean, whereas grounded ice can store water above sea level.

Melting consumes latent heat, while snow and ice can strongly affect the reflected fraction of incoming sunlight through albedo. Measurements therefore combine area, thickness, mass and temperature rather than using ice-covered area as a complete measure of stored ice. A change in one component need not imply an identical response in every other component.

Geophysics applies physics and mathematics to the structure, composition, and evolution of planets, including Earth. It connects measurements with mechanisms: plate tectonics explains moving crust, the geodynamo explains a planetary magnetic field, and albedo connects reflected sunlight to a planetary surface energy balance. Even a small model can reveal a useful organizing principle, such as the competition between ice-ocean albedo feedback and stabilizing thermal radiation.

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.

A melt pond is liquid water collected in a depression on melting ice. Its lower albedo than snow-covered ice can strengthen ice-albedo feedback. Pond drainage and connected pores can redistribute an oil spill under sea ice during warming.

For the albedo-to-forcing conversion, Earth intercepts sunlight over a disk but distributes the global mean over its full sphere, giving incident mean flux $S_0/4$. The absorbed shortwave flux is $F_{\mathrm{abs}}=(1-A)S_0/4$, and henceThe cooling magnitude is $1.37\times 10^{-5}\mathrm{W}{\mathrm{m}}^{-2}$ for the assumed one-year effect of one mass unit of injection. Albedo does not change the sunlight incident at the top of the atmosphere; it changes the absorbed/net incoming flux. The minus sign therefore accompanies the radiative cooling.

Injected sulfur dioxide is oxidized, principally through hydroxyl radical chemistry, to sulfuric acid, which nucleates or condenses into aqueous sulfate particles. A suitably sized, weakly absorbing particle population scatters some incoming sunlight back to space, increasing planetary albedo and decreasing absorbed solar energy. The sulfur dioxide molecules are precursors, not the persistent reflective population.

The albedo lifetime reflects aerosol formation, growth/aerosol coagulation, sedimentation, stratospheric transport and eventual return to efficiently scavenged tropospheric air. The stratosphere has little atmospheric precipitation, so the effect can outlast analogous tropospheric aerosol. The assumed year is an effective persistence scale, not a universal lifetime of sulfur dioxide oxidation or an exact clock for every particle.

Hold liquid-water content and cloud thickness fixed, as required for a partial sensitivity to droplet number. Then $\tau =K{N}_d^{1/3}$, so $d\tau /d{N}_d=\tau /(3N_d)$. Differentiating the albedo relation gives $dA/d\tau =7/(\tau +7{)}^2$, henceSince $A=\tau /(\tau +7)$ and $1-A=7/(\tau +7)$,The factor $A(1-A)$ peaks at $A=1/2$, so the albedo response to a fixed fractional droplet-number change is greatest at intermediate albedo and weaker near either saturation limit.

The left side is atmospheric heat storage per unit surface area. Incoming solar flux $Q$ is averaged over the spherical surface by the factor $1/4$; albedo $a$ reflects a fraction of it. The outgoing term is a reduced Stefan-Boltzmann flux, with effective emissivity $\gamma$ representing atmospheric greenhouse trapping. Ice-albedo feedback usually makes $a$ decrease as temperature rises through an ice-melting range; greater carbon-dioxide pressure increases infrared opacity and therefore decreases this effective $\gamma$. The simple estimates below hold albedo fixed.

Radiative equilibrium givesDifferentiating $\gamma T^4=\text{constant}$ gives $\Delta \gamma \simeq -4{\gamma }_0\Delta T/T_0\simeq -0.02474$. A doubling from $36$ Pa changes pressure by $36$ Pa, henceReference weathering balance requires $A_L{W}_0=v_0$, so $\boxed{W_0=2.0\times 10^{-3}\mathrm{kg}{\mathrm{m}}^{-2}\mathrm{y}{\mathrm{r}}^{-1}}$.

For the requested approximate system, take $(T_0,p_0)$ as the reference radiative/weathering equilibrium, set $\theta =(T-T_0)/\Delta T_c$, $P=p/p_0$, and $\hat{t}=t/t_{\mathrm{C}{\mathrm{O}}_2}$, whereLinearize the radiative loss, retaining terms linear in $T-T_0$ and $p-p_0$. Its two terms are $4\sigma {\gamma }_0{T}_0^3(T-T_0)$ and $-\sigma {\gamma }_1{T}_0^4(p-p_0)$. The carbon balance scales exactly with $v=v_{\mathrm{C}{\mathrm{O}}_2}/v_0$. ThereforeNumerically $\boxed{t_{\mathrm{C}{\mathrm{O}}_2}\simeq 7.41\times 10^3\mathrm{yr},t_T\simeq 0.0958\mathrm{yr},ϵ\simeq 1.29\times 10^{-5},\Gamma \simeq 0.1795}$. The time conversion uses seconds for the radiative constants and years for the supplied carbon production.

Because $ϵ\ll 1$, temperature rapidly approaches $\theta =\Gamma (P-1)$, while carbon changes slowly. A sustained $v\simeq 40$ has equilibrium determined byIts left side is strictly increasing on positive pressure, giving $\boxed{P_{\ast }\simeq 16.83,\Delta T_{\mathrm{eventual}}=\Delta T_c\Gamma (P_{\ast }-1)\simeq 36.9\mathrm{K}}$ within this extrapolated toy model. Counting volcanic production as well gives $v=41$ and $37.2$ K instead. An inventory estimate for the rise is $t_{\mathrm{C}{\mathrm{O}}_2}(P_{\ast }-1)/(v-1)\simeq 3.0\times 10^3$ years. Linearization near the final equilibrium gives the e-folding timeThus the carbon/weathering response is millennial, despite the short thermal adjustment time.

There are two numerical-reference inconsistencies in the printed setup. The requested expansion treats $(T_0,p_0)$ as its reference equilibrium although earlier labels pair $T_0$ with present-day $36$ Pa rather than $p_0=28$ Pa. Also, the calibrated $\Gamma =0.1795$ gives $13^{0.3}{e}^{13\Gamma }\simeq 22.3$, not the stated hint. That hint instead corresponds to $\Gamma \simeq 3/13$, obtained by calibrating a doubling of $28$ Pa. Such a convention gives a pressure around $13$–$14$ and a similar rough warming of several tens of kelvin. The displayed results above consistently use the stated $36$ Pa doubling calibration and explicitly declared reference-model approximation. A large extrapolation of a linearized radiative law is a property of this exercise, not a quantitative climate projection.

This is a weathering feedback with fast thermal adjustment.

Use dimensional quantities $T_d,t_d$ and set $T_d=T_{0}T$, $t_d=t_{0}t$. Dividing the energy balance by $\sigma \gamma T_0^4$ and the albedo equation by its time scale givesThus $\epsilon$ is the ratio of the atmospheric thermal response time to the much longer albedo response time. Converting kilometres to metres and years to seconds gives $C=10^7\mathrm{J}{\mathrm{m}}^{-2}{\mathrm{K}}^{-1}$, $T_0\simeq 266.3\mathrm{K}$, $t_0\simeq 3.156\times 10^{10}\mathrm{s}$, andThe thermal time is about half a year, compared with a millennium for albedo adjustment. This is a singular perturbation: setting $\epsilon =0$ removes the independent initial condition for temperature.

Set $x=T^4/T_0^4-1$. Since $\sigma \gamma T_0^4=Q/8$, the steady state energy balance is $2(1-a_0)=1+x$. Substitution givesThe two sides have slopes $\mu$ and $2$ at the origin. For $x>0$, $\tanh (\mu x)$ is strictly concave, its initial slope is $\mu$, and it remains below $1$. Consequently no positive root exists for $\mu \le 2$, while for $\mu >2$ there is exactly one positive root $x_{+}\in (0,1/2)$. Oddness supplies the corresponding negative root. At $\mu =2$ the only root is still zero. Hence multiple steady temperatures occur precisely when $\mu >{\mu }_c=2$.

The largest root corresponds monotonically to the largest positive absolute temperature:For large $\mu$, the positive root tends to $1/2$, so $a_{+}\simeq 1/4$ and $T_{+}/T_0\simeq (3/2{)}^{1/4}$. The negative branch has larger albedo and lower temperature.

The heat capacity per unit horizontal area is $C=\rho c_{p}d$: mass density times specific heat capacity times atmospheric depth. Thus $C\dot{T}$ is the rate of stored thermal energy per unit area. The solar constant $Q$ is the incident radiative flux on a surface normal to sunlight; the ratio of the Earth's intercepting disk area to its spherical area is $1/4$. The absorbed fraction is $1-a$, where albedo $a$ is the reflected fraction. The outgoing radiative flux is $\sigma \gamma T^4$, a Stefan–Boltzmann law with effective radiation factor $\gamma$ incorporating the simplified greenhouse effect. Every term therefore has units of power per area.

The albedo curve decreases smoothly, passes through $(T_0,1/2)$, and approaches $1/4$ at high temperature. At $T=0$ it is $(2+\tanh \mu )/4$, close to $3/4$ when $\mu$ is large. Its derivative isAn ice-covered planet reflects more sunlight than exposed land or ocean. Warming removes reflective ice, lowers albedo, and increases absorbed radiative flux: this is ice-albedo feedback. The parameter $\mu$ controls the sharpness of this transition; the formula is an idealized global response, rather than a local melting law.

Initially $a=1/2$ and $T=T_{+}/T_0>1$, so $\epsilon T_t=1-T^4<0$. On the fast time $s=t/\epsilon$, albedo changes by only $O(\epsilon )$ andThus temperature rapidly falls towards $1$. At the same time $a_t=a^{\ast }(T)-1/2<0$ while $T>1$, which seeds selection of the warm branch. Afterwards the approximate slow manifold is $T^4=2(1-a)$: decreasing albedo makes temperature rise slowly back towards $T_{+}/T_0$. The preceding figure separates these two time scales; a single time axis would compress the initial fall into a nearly vertical segment.

In the sharp-transition regime, after selection of the warm branch $a^{\ast }(T)\simeq 1/4$. Solving the linear ordinary differential equation with the leading matched value $a(0^{+})=1/2$ givesThese are outer solutions: their initial temperature is $1$, so the fast cooling layer is required to satisfy the actual initial temperature.

There is a useful qualification to the two-limit argument. For fixed finite $\mu >2$, put $x=1-2a$ on the slow manifold. Then $x_t=-x+\frac{1}{2}\tanh (\mu x)$, and $x=0$ is an unstable equilibrium with growth rate $\mu /2-1$. The small positive $O(\epsilon )$ seed from the fast layer can therefore produce a delay of order $\log (1/\epsilon )/(\mu /2-1)$. The simple exponential approximation assumes that the sharp albedo transition is sufficiently fast to make this delay negligible on the displayed slow scale. It is not a uniform approximation for every fixed $\mu >2$ as $\epsilon \to 0$.

Balances and meanings. The first equation is conservation of energy per unit area. Here $c$ is effective areal heat capacity, $Q$ is the solar constant, and $Q/4$ accounts for the ratio of the Earth's intercepting disc area to its total surface area. The albedo $a$ reflects a fraction of that input. The outgoing thermal radiation is the Stefan–Boltzmann law modified by $\gamma (p)$, an effective transmission factor for the greenhouse effect; it normally decreases as atmospheric carbon dioxide increases.

The albedo relaxes towards its temperature-dependent equilibrium on the ice-adjustment time $t_i$. The atmospheric mass inventory per unit pressure is $M_{{\mathrm{CO}}_2}A/(M_{a}g)$, using the molar mass ratio of carbon dioxide to air and hydrostatic atmospheric mass. Its sources and sinks are volcanic input $v$, removal by chemical weathering over land area $A_L$ at rate $W$, and the atmosphere-to-ocean exchange $h(p-p_s)$. The ocean equation is a molar mass conservation balance: ${\rho }_{{\mathrm{H}}_2\mathrm{O}}{V}_{\mathrm{oc}}$ is ocean water mass, division by $M_{{\mathrm{CO}}_2}$ converts the incoming carbon-dioxide mass to moles, and $bC$ removes dissolved carbon by burial/export. With $C$ in moles per kilogram, $b$ has units of kilograms of water-equivalent per year.

Times and reduction. Linearize the temperature equation at fixed $a,p$ about radiative equilibrium. Its exponential relaxation time isThe equality uses $\sigma \gamma T^4=Q(1-a)/4$. A stock-over-flux thermal time, sometimes used as a scale, is four times this, $0.381\mathrm{y}$; both describe a rapid thermal response relative to the other processes.

For pressure in atmospheres, write $I=M_{{\mathrm{CO}}_2}A10^5/(M_{a}g)$ for the atmospheric carbon-dioxide mass per atmosphere. The factor $10^5$ is essential because $h$ is given per atmosphere. The pressure and slow ocean exponential relaxation times areThus $t_T\ll t_p\ll t_i=10^4\mathrm{y}<t_C$. The leading quasi-steady approximation givesSubstitution in the two remaining equations yieldsThe coefficient of $\dot{a}$ remains $t_i$; the printed $t_T$ in the reduced equation is a typographical error. Eliminating a fast temperature variable does not replace the ice-response time. The ocean time here is the slow time after atmospheric balance cancels the exchange flux; holding pressure fixed after imposing a carbon-dependent $p_s$ would describe a different conditional relaxation.

The supplied volcanic input givesabout $41.1\%$ of the stated present value $2.0\times 10^{-3}$; that value is about $2.43C_0$.

Carbon-dioxide levels. Since $A_{L}W=3.0\times 10^{11}\mathrm{kg}{\mathrm{y}}^{-1}=v$, the atmospheric balance with $p_s=C/K$ givesThis uses the specified present ocean concentration, as requested, rather than replacing it by the model's smaller volcanic-only equilibrium.

After sustained additional input $v_i$, $C$ barely changes during the atmospheric adjustment. At this stageThis is an intermediate balance, approached over several centuries. On the slower ocean time,with pressure following $p\simeq C/K+v_i/h$. Since $C(0)<C_{\infty }$, pressure continues to rise slowly. Its eventual level in this stipulated constant-parameter model isThe two-reservoir carbon adjustment therefore has a rapid rise towards the frozen-ocean balance followed by a much longer increase. The figure solves the coupled linear pressure–ocean equations, retaining their small inventory correction to the reduced time scales.

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.

Open water absorbs more sunlight than bright sea ice, and melt ponds also lower the surface albedo. More absorbed heat promotes further melt and can be stored in the upper ocean, delaying autumn freezing. This ice-albedo feedback amplifies an initially forced loss during the sunlit season. Perovich and collaborators' 2008 measurements identify solar heating of the upper ocean as the main heat source for enhanced bottom melting in their observed Beaufort Sea case. Feedback is especially important locally after open water develops; it does not supply solar heat during polar winter.

The Maykut-Untersteiner model solves a seasonally forced, one-dimensional sea ice thermodynamics problem with a snow layer above saline ice. Let $z$ increase downward from the upper surface, with snow–ice interface $z=h_s$ and ice bottom $z=h_s+h_i$. Write $q=-k{T}_z$ for conductive heat flux, positive downward, and $I(z)$ for downward shortwave radiation penetrating the material. Local energy conservation givesFor $I=I_0{e}^{-\kappa z}$, the radiation term is $+\kappa I$: absorption heats the interior. Here $c_{\mathrm{eff}}$ is the derivative of specific enthalpy with temperature at prescribed ice salinity, and includes the latent heat of changing brine fraction. Using only pure-ice specific heat capacity would omit this storage. A usual approximation has ${\rho }_i{c}_{\mathrm{eff}}={\rho }_i{c}_{i,0}+a{S}_i/T_C^2$ and $k_i=k_{i,0}+b{S}_i/T_C$, with $T_C<0$ in degrees Celsius and positive, unit-dependent constants $a,b$. Thus warmer saline ice has larger effective heat capacity and lower thermal conductivity. Layerwise constant $k_i$ gives the usual finite-layer heat equation; full differentiation retains its spatial variation. Snow uses its own density and thermal conductivity, and its stronger extinction can make shortwave penetration negligible beneath it.

Define the net atmospheric input to the upper surface bywhere $\alpha$ is albedo, $H_{\uparrow }$ is turbulent sensible heat loss and $E_{\uparrow }$ is turbulent latent heat loss. Their signs can reverse. The outgoing longwave term uses the Stefan–Boltzmann law. Subtract the penetrating shortwave $I_0$, because it heats the interior rather than the infinitesimal surface. If $M_0\ge 0$ is surface ablation speed and ${\rho }_0$ the density of the material currently exposed, the upper boundary obeysBelow the melting point, $M_0=0$ and the flux balance determines $T_0$. Once it reaches $T_m$, the positive excess produces surface ablation; snow melts first, and only afterwards does that ablation remove ice. At the snow–ice interface, absent interfacial phase change,At the base, $T_i=T_f(S_w)$ and the Stefan condition isHere $G_b$ is signed bottom growth speed, and $F_o>0$ is ocean heat flux toward the ice. An upward conductive loss $-q_b>0$ freezes ice; an excess $F_o$ melts it. For a prescribed upward turbulent ocean heat flux, the model needs no separate ocean temperature solution. A water-side closure could write $F_o={\rho }_w{c}_{p,w}{K}_w{T}_{w,z}$ at the boundary in this downward coordinate. Snowfall is an independent mass input: with snowfall mass flux $P_s$, ${\dot{h}}_s=P_s/{\rho }_s-M_s$ in a constant-density approximation. The ice thickness obeys ${\dot{h}}_i=G_b-M_i$. Compaction or refreezing, if added, requires corresponding mass and energy terms rather than counting snowfall directly as ice growth. These boundary conditions and interior equations determine temperature, surface melt and basal growth together.

Under periodically repeated annual forcing, annual sea-ice thermodynamic equilibrium is a periodic annual state, not constant thickness:Winter accretion then balances net annual ablation. “Equilibrium thickness” normally denotes a mean or a specified point in this annual cycle. Different forcing can instead eliminate a perennial-ice solution.

Increasing ocean heat flux reduces bottom growth directly in the Stefan condition, lowering equilibrium thickness and eventually removing perennial ice. The stabilizing winter response is visible in the approximate conductive fluxwhen ice thins, conductive loss and winter regrowth increase. More snow lowers this flux and inhibits winter growth. However, snow also raises albedo and delays summer exposure of dark ice. The competing snow insulation and summer ice protection effects mean an annual snow sensitivity cannot be deduced from winter insulation alone. Under the original forcing, the Maykut-Untersteiner model gives an almost flat annual-thickness response below about $0.7\mathrm{m}$ annual snow depth, followed by strong thickening as snow persists into summer. This particular transition is not a universal threshold for other forcing or snowfall timing. Increasing summer albedo reduces absorbed solar energy, reduces surface and internal melting, and increases equilibrium thickness; sufficiently low albedo removes perennial ice. The original thermodynamic sensitivity calculations establish these three responses.

The three curves are original qualitative sketches, not a new numerical integration or digitization. The snow marker refers to that historical forcing; the unnumbered loss boundaries and ordinate indicate trends rather than calibrated universal thresholds.

For the summer retreat discussed around the time of this examination, the proposed mechanisms act as follows.

These processes overlap: circulation causes export, export changes thickness and open-water fraction, and those changes alter solar absorption. Without specified regional budgets there is no defensible unique percentage partition. For enhanced Beaufort summer basal melt, upper-ocean solar heating and western-Arctic heat pathways are more directly relevant than invoking unmodified heat from a deep, isolated Atlantic layer. ice-mass-balance observations of extreme Beaufort melt identify increased solar absorption in open water as an important source.

Take ${\rho }_i=900\mathrm{kg}{\mathrm{m}}^{-3}$ as an explicit ice-density approximation and use the supplied $L_f=336000\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}$. The latent heat for $1.5\mathrm{m}$ basal ablation isJune–September contains about $122$ days, soUsing $917\mathrm{kg}{\mathrm{m}}^{-3}$ would give $43.8\mathrm{W}{\mathrm{m}}^{-2}$. This is the net flux used for melting: the gross ocean heat flux can also supply conductive losses and changes in cold content. It greatly exceeds an annual central-pack background of a few $\mathrm{W}{\mathrm{m}}^{-2}$, but is compatible with a summer upper ocean absorbing solar radiation at tens to hundreds of $\mathrm{W}{\mathrm{m}}^{-2}$. For scale, an added open-water fraction $0.3$, a shortwave input $200\mathrm{W}{\mathrm{m}}^{-2}$ and an ice–water albedo contrast $0.6$ produce $0.3\times 200\times 0.6=36\mathrm{W}{\mathrm{m}}^{-2}$ additional area-mean absorption before losses. This is an illustrative budget, not a measurement. Wind mixing and circulation can deliver stored solar or advected heat to the underside.

The ice-albedo feedback closes the warming loop: retreat exposes darker water; extra solar absorption warms the upper ocean; basal and lateral melt increase retreat. The stored heat also delays autumn freeze-up, exposing water to atmospheric heat transfer for longer. Its later release warms the atmosphere, but must not be counted as a second source in addition to the solar energy that originally entered the ocean. In winter, greater heat loss from open water and rapid growth of thin ice partly oppose the summer positive feedback.