Past exam of the mathematics course of the University of Cambridge / 2013 / iii / Paper 72 (Polar Oceans and Climate Change)

Write $w(x,t)$ for the vertical displacement of the sea ice, and take the undisturbed water surface as $z=0$, with water occupying $z<0$. The elastic plate has areal mass $m={\rho }_{i}h$ and bending stiffnessHere $E$ is Young's modulus and $\eta$ is Poisson's ratio. We neglect in-plane prestress, viscosity and plate shear deformation, and linearize about hydrostatic equilibrium. These are important assumptions: perfect elasticity alone does not specify every term in a floating-plate model.

For a plane wave $w=\hat{w}{e}^{i(kx-\omega t)}$, $k>0$, potential flow in deep water has velocity potential $\varphi =\hat{\varphi }{e}^{kz}{e}^{i(kx-\omega t)}$. This solves Laplace's equation and decays downwards. The kinematic boundary condition gives $k\hat{\varphi }=-i\omega \hat{w}$. Linearizing the water pressure at the displaced interface gives an upward excess loadThe elastic plate equation is $m{w}_{tt}+D{w}_{xxxx}=p$. Substitution and multiplication by $k$ therefore give the flexural-gravity wave dispersion relationEquivalently, with $\alpha ={\rho }_{i}h/{\rho }_w$ and $\beta =D/{\rho }_w$,This determines the positive wavenumber implicitly for prescribed positive angular frequency. It is unambiguous: the derivative of ${\omega }^2$ is $(g+5\beta k^4+4\alpha \beta k^5)/(1+\alpha k{)}^2>0$, while ${\omega }^2$ runs from zero to infinity.

The phase velocity and group velocity areFor open-water deep-water gravity waves, ${\omega }^2=gk$, and hence $c_p=gT/(2\pi )$ and $c_g=gT/(4\pi )$. At large period the ice-covered curves approach these straight lines. At shorter period, plate bending raises the speeds, so both ice-covered curves turn upward as period decreases. In the bending regime with negligible plate inertia, $\omega \propto k^{5/2}$ and $c_g\simeq (5/2)c_p$; in the formal plate-inertia-dominated limit, $\omega \propto k^2$ and $c_g\simeq 2c_p$. The latter extrapolation eventually leaves thin-plate validity and should not be read as a prediction at arbitrarily small wavelength.

The plotted parameters are illustrative rather than measured at the observation site. They give a group-velocity minimum of a flexural-gravity wave near $15.1\mathrm{m}{\mathrm{s}}^{-1}$ at period $17.6\mathrm{s}$. There is no arbitrarily slow wave-energy branch under a continuous elastic sheet. Energy put into a localized disturbance travels away at at least this minimum group velocity; a slowly moving wind system cannot retain a wave packet indefinitely beneath itself. This reduces the opportunity for sustained local growth compared with slow, short open-water waves, and the continuous cover also prevents direct wind forcing of an exposed water surface. Incoming long swell can still propagate.

A group-velocity minimum of a flexural-gravity wave is not by itself a universal minimum wind speed for wave generation. A steadily translating forcing pattern requires a phase velocity matching its translation speed, so the minimum of $c_p$, not $c_g$, supplies the corresponding resonance threshold. Random wind forcing, dissipation and aerodynamic coupling must be specified before making an absolute generation claim.

The sea ice acts as a frequency-selective filter. Over a long path, wave attenuation in sea ice is generally much stronger for shorter surface gravity waves, through repeated scattering, internal ice losses and water-side dissipation. Their amplitude can fall below the tiltmeter's detection level even if the source initially generated them. Thus the long-period swell survives while the shorter-period tail does not.

The disappearance near 14 s is an attenuation and detectability limit, not a forbidden-frequency interval of the ideal elastic plate. The flexural-gravity wave dispersion relation admits a real positive wavenumber for every positive angular frequency. In an impulsive-source picture, shorter-period deep-water gravity waves also arrive later because their group velocity is smaller, but delayed arrival alone does not explain a persistent observed cutoff.

Use a common emission time $t_0$ and a path length $L$, and treat each observed period as a narrow wave packet. For deep-water gravity waves, its group velocity is $gT/(4\pi )$, not the phase velocity $gT/(2\pi )$. Therefore the dispersive swell source inversion givesThus period decreases hyperbolically, while frequency $1/T$ increases linearly. The elapsed time between the two observations is $40\mathrm{h}40\min =146400\mathrm{s}$, soThe frequency slope is $2.6021\times 10^{-7}\mathrm{Hz}{\mathrm{s}}^{-1}$. Initially $dT/dt=-0.843\mathrm{s}{\mathrm{h}}^{-1}$; at 14 s it has slowed to $-0.184\mathrm{s}{\mathrm{h}}^{-1}$. The inferred emission was $C/30=35\mathrm{h}35\min$ before the first arrival, namely 19 March at 1710 UT. Both observed endpoints give this same time.

A northward-propagating swell with this distance scale points to a remote energetic storm south of the high Arctic rather than local wind acting on continuous sea ice. Along a meridional route, 3,000 km is about $27^{\circ }$ of latitude, putting the source on the scale of the northern North Atlantic and adjacent open seas. A route through Fram Strait is plausible, but longitude and refraction are not supplied, so no particular storm or unique source position follows. Long swell arriving first and steadily increasing frequency are the expected signatures of remote wave dispersion. The distance and time are conditional on an approximately impulsive source and open-water propagation speeds; passage through ice, currents and finite storm duration produce corrections.

It is important to distinguish loss of forward-going surface-gravity-wave energy from conversion of mechanical energy into heat. Scattering attenuation by ice floes redirects energy; it can attenuate a coherent transmitted wave without dissipating the total energy.

For fixed floe geometry, increasing frequency usually increases attenuation over the relevant swell range: shorter wavelengths respond more strongly to the contrast between water and the elastic plate, and to repeated floe edges. Long surface gravity waves have weak curvature and often penetrate much farther. This is a trend over a specified frequency range, not a theorem excluding resonances.

The diameter dependence is governed by $d/\lambda$. An ice floe much smaller than the wavelength moves nearly with the water and scatters weakly. Scattering becomes appreciable when floe size is comparable with the wavelength, and interference between its two edges can give maxima and minima. At fixed ice concentration, larger ice floes also mean fewer edges per unit propagation distance, roughly proportional to $p/d$. Consequently the attenuation coefficient need not increase monotonically with diameter: the single-floe reflection and the number of encounters must both be considered. Thickness increases areal inertia as ${\rho }_{i}h$ and bending stiffness as $h^3$, generally increasing wave mismatch and reflection, although detailed frequency-dependent resonances again prevent a universal monotonic law.

When $d\ll \lambda$, particularly for frazil ice and pancake ice, weak individual scattering leaves other processes dominant. Relative crystal and water motion causes viscous dissipation; an aggregate layer can behave as a viscous or viscoelastic material, and pancake ice collisions, rubbing and overwash remove energy. Their importance depends on concentration and wave amplitude.

For a uniform continuous sheet with horizontal dimensions much greater than the wavelength, there are no repeated floe edges in its interior. A perfectly elastic sheet over inviscid water supports undamped flexural-gravity waves, so internal scattering is not an explanation of decay there. Real attenuation can instead arise from internal ice anelasticity or viscoelasticity, a dissipative sub-ice viscous boundary layer, turbulence, cracks and brine-related processes. Small-floe mixtures and continuous sheets require dissipation models beyond the isolated-floe scattering picture.

Take $x$ positive in the off-ice wind direction, from the compact pack towards the open sea. The wind initially separates the outer ice floes, creating irregular polynyas. A larger opening gives more open-water fetch, so stronger short wind waves develop before reaching its downwind edge. Their reflection supplies a positive force on the downwind ice floes; these catch their neighbours and compact into an ice-edge band. Incoming longer swell exerts force in the opposite direction. The short waves can exert substantial force despite their smaller amplitude because a small ice floe reflects them much more effectively than it reflects the long swell.

The surface-gravity-wave energy of the short waves grows with fetch in the windward polynya, then decays rapidly across the band. Swell enters from the seaward side and usually attenuates more slowly. This gives inward forcing from the two sides: short-wave forcing is largest at the windward face and swell forcing is largest at the seaward face. In the next open polynya, short wind waves regrow from the weak transmitted component. Reflected waves also enhance energy locally on the incident side, with interference that is omitted from a smooth, phase-averaged sketch.

Each plotted component is normalized by its own incident energy; the curves do not assert equal absolute wind-wave and swell energies. The band is partially transmitting. For a perfectly opaque reflector, the transmitted component would instead vanish.

The initial bands have unequal floe inventories, widths and forcing. Differential drift and collisions merge some into composite bands. Larger bands tend to shield smaller downstream accumulations from the short-wave forcing needed to keep them separate, while sufficiently wide intervening polynyas can generate fresh wind waves and maintain separation. Finite available ice, available wave fetch, the opposing swell, wind strength and duration, floe size and thickness, and subsequent mergers limit the number of persistent bands. The wave-force formula alone does not select a universal count. This mechanism and its merger interpretation are supported by the original ice-band study.

There is a normalization issue in the requested stress calculation. Let $\mathcal{E}={\rho }_{w}g{a}^2/2$ be standard linear surface-gravity-wave energy for crest amplitude $a$. The deep-water wave radiation stress is $S_{xx}=\mathcal{E}/2$. Hence standard momentum balance givesFor lossless reflection with amplitude reflection coefficient $R$, $r^2=R^2{a}^2$ and $t^2=(1-R^2)a^2$, so $F_{\mathrm{phys}}={\rho }_{w}g{R}^2{a}^2/2$.

In the usual independent-floe, weak-reflection approximation, there are about $p/d$ effective layers per unit distance. Each removes a fraction $R^2$ of the forward energy, givingThis yields the intended wave-driven ice-band compaction:The attenuation approximation retains the leading term in $R^2$; an independent discrete-layer model instead gives an energy coefficient $-(p/d)\log (1-R^2)$. Neither coefficient follows from the single-object force formula without this additional scattering closure.

For ordinary crest amplitudes, the force printed in the paper is $4F_{\mathrm{phys}}$. With that printed normalization and the same attenuation law, its derivative is four times the requested stress. More generally, if $a^2=a_0^2{e}^{-\kappa x}$, the printed lossless force gives $-F^{\prime }=2{\rho }_{w}g{R}^2\kappa a^2$. Recovering the stated stress from it requires $\kappa =p{R}^2/(4d)$ instead, a different unspecified attenuation convention. The supplied force and stress cannot both be derived from standard amplitudes and the usual floe-layer attenuation law. Restoring the factor $1/4$ in the force gives a consistent intended model.

Using the printed force for the numerical question, perfect reflection gives $r=a$, $t=0$, and thereforeUsing the requested stress formula for partial reflection givesFor comparison, the consistently normalized perfect-reflection force is $50.3\mathrm{N}{\mathrm{m}}^{-1}$, and the printed force with standard attenuation would give $0.628\mathrm{N}{\mathrm{m}}^{-2}$ for the partial-reflection derivative. The approximate stress is a force per horizontal area, not a direct measure of the three-dimensional ice-skeleton stress.

The inward force gradient helps maintain a coherent, close-packed band, especially near the incident-wave faces. It is modest enough that wind, current, swell changes or mergers can disrupt the arrangement; it does not guarantee permanent mechanical stability. Perfect reflection estimates a bounding force, whereas the much smaller partial-reflection stress varies as $R^4$. No mechanical-strength law is supplied, so stability can be assessed only qualitatively.

A sea-ice pressure ridge forms where horizontal convergence compresses sea ice. The initially thinner sheet or colliding ice floes fracture and ride over or under one another. Continued convergence piles broken blocks into an emergent sail and a submerged keel. The submerged volume is normally larger because buoyancy supports the pile. Pores and brine-filled gaps initially make the rubble unlike a solid intact sheet; refreezing can consolidate it.

A sea-ice shear ridge develops along a fracture where neighbouring ice moves tangentially in opposite directions or at different speeds. Rough edges interlock, crush and locally converge, producing chains of piled blocks along the shear boundary. Thus the large-scale strain is mainly shear, but the actual production of ridge rubble involves local compression. Pure sliding of perfectly smooth parallel surfaces need not create a ridge. Pressure ridging is driven by convergence; shear ridging is driven by relative tangential motion with local crushing and convergence.

Let $H$ denote a ridge's peak sea-ice draft, reserving $h$ for the draft at a randomly sampled position. In the exponential ridge-draft model, normalization by the line density $\mu$ givesThe mean peak sea-ice draft isConsequentlywith $h_m>h_0$, $b$ having dimensions inverse length and $B$ inverse length squared. The normalized peak probability density function is a shifted exponential distribution.

For the triangular argument, interpret the common ridge shape as geometrically similar triangles with common along-track slope $\tan \delta$ and variable peak height. Literal congruence would require identical sizes and could not coexist with an exponential peak-draft distribution. Each side of a triangle has $dx=|dh|/\tan \delta$. A ridge reaching draft $H\ge h$ therefore contributes $2\cot \delta dh$ of horizontal track in the interval $[h,h+dh]$. Summing this occupation length over all qualifying peaks proves the triangular ridge occupation identity:ThusThis is a tail relation for sampled draft occupation, not an instruction to normalize $n$ and $g$ identically. Below $h_0$, the ideal triangles contribute $2\mu \cot \delta$ rather than the same exponential; level ice and gaps contribute their own draft distributions. If triangular keels are referenced to a level-ice base, the vertical coordinate must be shifted consistently. We also require nonoverlapping occupation: arbitrary choices of $\mu$, mean draft and slope can otherwise demand more than the available track length.

Observed mean keel slopes are typically of order $20^{\circ }$–$30^{\circ }$, with broad individual variation rather than a single universal angle. Orientation matters: if a track crosses a straight ridge at angle $\psi$ to the crest, $\tan {\delta }_{\mathrm{track}}=\tan {\delta }_{\perp }|\sin \psi |$. The track slope can therefore approach zero at a grazing crossing. A sonar morphology study found location-dependent mean slopes about $22^{\circ }$–$27^{\circ }$ after correcting for ridge orientation.

Young sea-ice pressure ridges often have recognizably triangular sections with angular, porous rubble and comparatively continuous crests. Melting, refreezing and repeated cracking modify older ridges: their blocks can become rounded and consolidated, and their keel or crest can fragment into separated hummocks rather than retain one triangular shape. A pre-exam multibeam study found first-year sea ice ridge slopes averaging roughly $27^{\circ }$, while multi-year sea ice ridges often consisted of irregular separated smooth blocks. Multi-year sections can be broader or locally shallower, but age alone does not determine one slope angle. The constant-slope triangle is a useful statistical idealization, not a faithful shape for every old ridge.

Use the exponential ridge-draft model along the direction of motion relative to the seabed. During time $dt$, the stationary point samples track length $|V|dt$. Peaks deeper than the seabed have line densityThe candidate seabed gouging by ice encounter rate is thereforeFor the usual case $D\ge h_0$, replace $\max (D,h_0)$ by $D$. For $D<h_0$, every counted ridge is deep enough and ${\nu }_D=|V|\mu$.

The units are inverse time. This is the number of potentially scouring keel passages per unit time, not sediment volume removed per unit time. The latter also requires keel width, sediment properties and an erosion law. The expression assumes the ridge statistics remain applicable during drift and grounding does not arrest or reshape the deep keels before they reach the point. If speed fluctuates, use mean relative speed under independence from ridge occurrence, not the magnitude of a vector-mean drift that might cancel reversing motion.

The independent, spatially homogeneous random-placement idealization gives a Poisson process of intersections along the sampling line. With line intensity $\mu$, an interval of length $s$ contains no ridge with probability $e^{-\mu s}$. The spacing $X$ thus has exponential distributionRandom orientation changes the intersection intensity, which has already been represented by the measured $\mu$. Merely saying that positions are random does not prove a Poisson process: independence and homogeneity are additional assumptions. Finite segments, clustering or excluded widths can invalidate them.

For a shifted lognormal distribution, write $Z=\log (X-\theta )\sim N(m,{\sigma }^2)$ with $\sigma >0$. Applying the change of variables formula, $dZ/dX=1/(X-\theta )$, givesThe threshold $\theta$ is a minimum observable separation. Finite keel widths impose geometric exclusion, and the sonar's footprint and ridge-identification criterion can suppress a shallow peak near a deeper peak. This sonar ridge shadowing or resolution effect means that $\theta$ need not be a universal physical distance between every pair of ridges; it can partly reflect how the profiles were processed.

A lognormal distribution is compatible with products of many positive factors: taking logarithms turns multiplicative changes into sums, which can approach a normal distribution. Repeated deformation, breakup, convergence and merging across scales are plausible contributors. The observed form therefore suggests correlated or multistage ridging rather than the simplest independent-intersection model. It does not identify a unique ridging mechanism. For example, directly generating $X=\theta +\exp Z$ with a normal distribution for $Z$ produces the same spacing law without specifying any particular mechanics. Distributional agreement must be supplemented by dynamical and spatial evidence.

Use historical changes up to 2012 rather than the present-day Arctic state, and separate extent, actual ice-covered area, thickness and age composition. They do not have interchangeable rates.

For summer horizontal coverage, September minimum extent fell from roughly $7$ million square kilometres around 1980 to $3.4$ million in 2012: about half the earlier extent. The fitted September monthly-mean trend through 2012 was about $-91,600\mathrm{k}{\mathrm{m}}^2\mathrm{y}{\mathrm{r}}^{-1}$, or $-13\%$ per decade relative to the 1979–2000 mean. Extent includes the whole area of grid cells above the specified ice-concentration threshold; actual covered area additionally weights fractional cover. A separate concentration-weighted September ice-area analysis for 1979–2012 gives a decline of roughly $0.95\times 10^6\mathrm{k}{\mathrm{m}}^2$ per decade, or $95,000\mathrm{k}{\mathrm{m}}^2\mathrm{y}{\mathrm{r}}^{-1}$; this retrospective historical-period estimate is reported in an observational-area comparison. It should not be confused with the extent trend, even though these two absolute slopes are similar. These figures and definitions are documented in the 2012 Arctic sea-ice observations.

For thickness, the submarine and satellite record in the declassified submarine-data region, covering about $38\%$ of the Arctic Ocean, gives a winter mean of $3.64\mathrm{m}$ in 1980 versus $1.89\mathrm{m}$ in 2008: a $48\%$ reduction, averaging about $0.063\mathrm{m}\mathrm{y}{\mathrm{r}}^{-1}$. It is a regional winter comparison, not a basin-wide summer measurement. The same combined analysis reported recent 2003–2008 declines around $0.10\mathrm{m}\mathrm{y}{\mathrm{r}}^{-1}$ in winter and $0.20\mathrm{m}\mathrm{y}{\mathrm{r}}^{-1}$ in summer. The summer record is shorter and cannot justify extrapolating one constant summer-thickness slope back to 1980. See the original thickness analysis.

For composition, repeated summer loss and export depleted the thick multi-year sea ice reservoir and increased the relative importance of young and first-year sea ice. A directly comparable age indicator is the March fraction aged at least four years: about $26\%$ in 1988, $19\%$ in 2005 and only $7\%$ in 2012. This winter age measure records the loss of ice that had survived earlier summers; it is not the fraction of surviving September ice that is first-year ice. The youngest ice disproportionately melts in summer, so the age mix of survivors differs from that of the preceding winter cover. Overall, summer cover became smaller, thinner and supported by a much depleted reservoir of older ice.

Several mechanisms can accelerate Arctic sea ice decline. The ice-albedo feedback increases solar absorption as dark water replaces bright ice. Additional ocean heat content delays autumn freeze-up, leaving less time for winter growth. Thinner ice needs less latent heat to disappear, and fractured mobile ice is more easily exported or redistributed by wind stress and currents. Melt ponds lower surface albedo; increased open-water fetch permits waves that break the ice further. Persistent atmospheric warming and warmer incoming water act on this weakened cover.

However, strict irreversibility is not implied by these positive feedbacks. Winter open water loses heat, and thin ice grows rapidly because its conductive resistance is low. As a concrete counterexample to an unavoidable one-way transition, a 2011 coupled-model experiment imposed an ice-free summer and found recovery of ice extent typically within two years. This establishes a physically consistent recovery mechanism, not a guarantee that every real loss reverses on that timescale.

Under continued warming, rebuilding the former multi-year sea ice cover is unlikely; loss of one summer's cover is nevertheless not intrinsically irreversible. Sustained greenhouse gas forcing changes the climatic state towards which ice recovers, and rebuilding several age classes takes multiple summers of survival. The qualified conclusion is persistence or worsening under the continuing forcing, not a proved thermodynamic prohibition of recovery at fixed or reduced forcing.

Use seawater specific heat capacity $c_w\simeq 4.0\times 10^3\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$, which is an additional standard material approximation: the paper supplies the ice heat capacity, not the seawater heat capacity. Relative to freezing, the assumed uniform column has $\Delta T=7-(-1.8)=8.8\mathrm{K}$. Its ocean heat content per horizontal area isThis is heat above the reference freezing state, not the absolute internal energy of seawater.

If all this heat reaches ice already at its fusion temperature, divide by the supplied latent heat:For a chosen ice mass density ${\rho }_i=917\mathrm{kg}{\mathrm{m}}^{-3}$, this corresponds to about $5.86\mathrm{m}$ of ice. Density was not specified for this conversion, so the mass per area is the result that needs no further ice-density assumption.

Cold ice must first be warmed. If its initial temperature is $T_i<T_m$, the corresponding ideal melt mass is $Q/[L_f+c_i(T_m-T_i)]$, using the supplied ice specific heat capacity $c_i=2100\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$. No initial ice temperature or salinity-dependent fusion temperature is supplied. The latent-only result is therefore an upper bound; actual melting is smaller when sensible warming, ocean or atmospheric losses and incomplete heat transfer matter.

Solar geometry must be included before multiplying by the summer duration. Let $\varphi =71^{\circ }$, solar declination $\delta$, and hour angle $u$, measured from local noon. The cosine of solar zenith angle isOnly positive values receive sunlight. Integrating through a day gives the daily mean solar irradiance at the top of the atmosphere:with $H_0$ clipped to $\pi$ for polar day and to zero for polar night. At the solstice, for example, polar-day averaging gives $\overline{S}=S_0\sin 71^{\circ }\sin 23.44^{\circ }\simeq 514\mathrm{W}{\mathrm{m}}^{-2}$. It would be wrong to apply $S_0$ continuously to a horizontal surface.

A simple seasonal approximation $\delta (n)=23.44^{\circ }\sin [2\pi (n-80)/365]$, with calendar day $n$, gives a June–August daily-mean average of about $427\mathrm{W}{\mathrm{m}}^{-2}$. There are 92 days. With surface albedo $\alpha =0.1$, the no-atmosphere absorbed-solar ceiling isThis calculation neglects the small seasonal change in Earth-Sun distance. The ceiling is larger than the $1.80\mathrm{GJ}{\mathrm{m}}^{-2}$ required in part (i), so geometry alone does not make a uniform $7^{\circ }\mathrm{C}$ column impossible.

A reasonable conditional estimate includes an effective atmospheric short-wave transmission $\tau$ and net non-solar loss $\overline{L}$:before adding advection or subtracting ice melting. The parameters are scenario assumptions, not measurements supplied by the question. For example, $\tau =0.6$ gives $1.83\mathrm{GJ}{\mathrm{m}}^{-2}$ before other losses, barely enough; with a modest mean loss of $50\mathrm{W}{\mathrm{m}}^{-2}$, the retained amount is only $1.44\mathrm{GJ}{\mathrm{m}}^{-2}$. That would raise a uniform 50 m column from freezing by about $7.0\mathrm{K}$, reaching roughly $5.2^{\circ }\mathrm{C}$. With no other losses the transmission required for $7^{\circ }\mathrm{C}$ is $1.804/3.058\simeq 0.590$; with that illustrative loss it rises to about $0.720$. Clouds, emitted thermal radiation, evaporation, transfer to colder water and melting all affect the balance.

The satellite surface temperature is insufficient evidence for a $7^{\circ }\mathrm{C}$ seabed. A warm, shallow ocean mixed layer can overlie colder water because meltwater and salinity maintain stable density stratification. Heating only the upper 10 m through $8.8\mathrm{K}$ costs $0.361\mathrm{GJ}{\mathrm{m}}^{-2}$, much less than heating all 50 m. Warm Pacific-water advection can also raise surface temperature or supply additional heat. If measured net solar input were too small for full-depth warming, shallow surface heating would be the natural alternative; if mixing and additional heat supply were strong enough, full-depth warming remains possible. The missing transmission, loss, mixing and inflow information prevents a unique yes-or-no conclusion from the supplied surface observation.

Identified Arctic methane sources include anaerobic decomposition in wetlands and saturated tundra, bubbles from lakes and thermokarst lakes formed by thawing permafrost, geological gas seeps, and leakage from fossil-fuel extraction or transport. Warming can expose frozen organic carbon to microbes, promoting methanogenesis where oxygen is scarce; it does not convert all thawed carbon into methane. Field measurements of Siberian thaw lakes directly documented substantial methane bubbling associated with lake expansion into permafrost.

Shelf sediments can also generate or release methane, including gas stored beneath frozen sediment and, where pressure and temperature permit, gas in hydrates. The proposed Arctic methane feedback needs several links: seabed warming must reach the relevant material, release or production must occur, and enough methane must reach the atmosphere rather than dissolve or be oxidized in sediment and water. Visible plumes alone do not establish an imminent large atmospheric pulse.

Methane is a powerful greenhouse gas because it absorbs planetary thermal infrared radiation and alters the atmospheric emission balance. It also affects atmospheric chemistry, including ozone and stratospheric water vapour. Its relatively short atmospheric lifetime makes its warming contribution particularly strong over short horizons. The global warming potential compares time-integrated forcing from equal-mass pulses with that of carbon dioxide; the horizon must always be stated. Values available before this examination were about 72 over 20 years and 25 over 100 years, including the specified methane indirect effects, as recorded in the 2007 IPCC assessment. These are historical assessment values, not timeless material constants or instantaneous temperature multipliers.

A separate direct feedback is the ice-albedo feedback. When sea ice disappears, darker open water reflects less sunlight, stores more ocean heat content, and favours further melting and delayed freeze-up. Locally,For illustration, replacing ice of surface albedo $0.6$ by water of surface albedo $0.1$ under $200\mathrm{W}{\mathrm{m}}^{-2}$ incoming sunlight adds $100\mathrm{W}{\mathrm{m}}^{-2}$ of absorption on the newly exposed area. Global forcing requires weighting by affected area, season and clouds; this is a local example. Sea-ice loss can amplify warming through reduced reflection without requiring methane release. Increased evaporation and altered clouds supply further feedbacks, but their short-wave and long-wave effects can compete.