Codex @codex  0

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]

TODO: user-deferred biological component. Biological carbon uptake, storage and release are outside the current scope. The physical and ocean-chemistry explanation is given here; the deferred component means this is not a complete all-process answer.

There is no single first-order atmospheric lifetime for an added carbon dioxide pulse. Air–sea gas exchange and adjustment of the surface mixed layer act on years-to-decades scales; redistribution into the ocean interior takes centuries to millennia. These transfers partition carbon between atmosphere and ocean rather than making the entire added carbon inventory disappear. Dissolution also changes the ocean's carbonate speciation: conversion among dissolved carbon dioxide, bicarbonate and carbonate buffers the uptake, but provides a finite capacity at a given alkalinity. A substantial perturbation therefore remains airborne after the faster exchanges have adjusted.

Further removal involves dissolution of carbonate sediments and carbonate compensation over several thousand years, followed by the much slower silicate-weathering and burial pathway, typically tens to hundreds of thousands of years. The exact fractions and times depend on pulse size and model conditions. A thousand years exceeds the fast exchange times but not the complete ocean–solid-Earth adjustment time, so a residual atmospheric perturbation is expected. An intercomparison found a substantial long tail after ocean equilibration, with further carbonate neutralization on roughly 3–7 kyr scales.^[ref]

Much stratospheric heating is supplied by absorption of solar ultraviolet radiation by oxygen and ozone, whereas carbon dioxide is an important emitter of thermal infrared radiation. In the usual stratospheric balance, increasing carbon dioxide enhances net infrared cooling at the original temperature without supplying a matching increase in the solar heating. Temperature must fall until infrared emission and absorption again balance the heating.

A useful optically thin schematic is $Q_{\mathrm{solar}}=c_{\mathrm{C}{\mathrm{O}}_2}K[B(T)-\overline{J}]$, where $B(T)$ represents the relevant thermal source function and $\overline{J}$ the incident infrared radiation field. With positive fixed solar heating and an approximately fixed radiation field, increasing $c_{\mathrm{C}{\mathrm{O}}_2}$ lowers the required $B(T)-\overline{J}$, hence lowers $T$. This explains why the surface greenhouse argument cannot simply be applied locally to the stratosphere. Spectral blocking of upwelling radiation contributes too; the real response depends on altitude, ozone heating and radiative transfer, rather than emissivity alone.^[ref]

The solar constant measures the incident power per area normal to the sunlight. A spherical Earth of radius $R$ intercepts sunlight over the disc $\pi R^2$, but radiates over its surface $4\pi R^2$. For uniform surface albedo $a$ in the global energy-budget approximation, absorbed power is $(1-a)S\pi R^2$. 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, $T_E^4=⟨T^4⟩$; 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 $R=6.371\times 10^6\mathrm{m}$, so $A_E=4\pi R^2=5.10\times 10^{14}{\mathrm{m}}^2$. The changed area is $A_i=10^{12}{\mathrm{m}}^2$, givingThe initial planetary equilibrium temperature is $T_0=254.8\mathrm{K}$. ThusThe differential estimate $\Delta T\simeq -T_0\Delta a/[4(1-a_0)]$ 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 ${\rho }_i=917\mathrm{kg}{\mathrm{m}}^{-3}$ and the given $L_f=336000\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}$. 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 $92$ days, so the mean ocean heat flux to the ice isUsing ${\rho }_i\simeq 900\mathrm{kg}{\mathrm{m}}^{-3}$ instead gives about $76\mathrm{W}{\mathrm{m}}^{-2}$, within the intended accuracy. Around the date of this examination, the estimated net anthropogenic radiative forcing since preindustrial times was about $1.6\mathrm{W}{\mathrm{m}}^{-2}$; doubling carbon dioxide gives about $3.7\mathrm{W}{\mathrm{m}}^{-2}$. 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 $100\mathrm{m}$ water column, take ${\rho }_w=1025\mathrm{kg}{\mathrm{m}}^{-3}$ and liquid-water specific heat capacity $c_w\simeq 4000\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$. Assuming all the same energy remains in that uniformly mixed column,If the water initially lies near the seawater freezing point, about $-1.8{}^{\circ }\mathrm{C}$, its final temperature is approximately $\boxed{-0.3{}^{\circ }\mathrm{C}}$. 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.

The case for anthropogenic global warming rests on a physically consistent combination of forcing, energy storage and observed response, rather than on the claim that every unusual event has the same cause. For this examination, the observational comparison is the evidence available around 2010, and the two forecast horizons are approximately 2060 and 2110. Separating the robust mechanisms from uncertain regional details is essential to a coherent attribution of recent climate change.

The greenhouse effect follows from radiative transfer. Sunlight supplies energy, while the Earth loses energy mainly through long-wave radiation. Carbon dioxide, methane and other greenhouse gases absorb and emit in parts of that outgoing spectrum. Increasing their concentrations alters the levels and temperatures from which radiation escapes, initially reducing the net loss of energy for a given surface state. The surface and atmosphere must then warm until outgoing radiation again balances the input. Strong absorption near the centres of some bands does not eliminate the influence of additional gas: absorption in the wings and the changing height of emission still matter. For scale, the approximate carbon dioxide forcing law $\Delta F\simeq 5.35\log (C/C_0)\mathrm{W}{\mathrm{m}}^{-2}$ gives about $3.7\mathrm{W}{\mathrm{m}}^{-2}$ for doubling. It describes a forcing, not the ultimate temperature response, which additionally involves climate feedbacks and ocean heat uptake.

Measured increases in atmospheric carbon dioxide are consistent with fossil-fuel combustion and changes in land use. The carbon isotope signature and atmospheric oxygen consumption provide independent evidence of the combustion source, while the ocean and land absorb part of the emitted carbon. Methane and other climatically active gases add to the forcing. Aerosols exert partially compensating effects, so net human forcing is smaller than the sum of the warming gases alone. Uncertainty in that compensation affects the quantitative rate of warming without removing the greenhouse mechanism.

A sustained rise in global temperature is accompanied by increased ocean heat content, glacier retreat and sea level rise. This combination is harder to explain as an artefact of a few surface thermometers than a temperature series in isolation: the ocean is a large energy reservoir, and thermal expansion and melting have independent observational consequences. Internal variations of ocean circulation can move heat between regions and alter warming over particular decades, but cannot indefinitely supply a positive global energy accumulation without a change in the energy budget. Solar changes and volcanic eruptions also influence climate, and must be included in attribution; they do not reproduce the full recent warming pattern and sustained heat gain. Volcanic aerosol episodes chiefly produce temporary cooling. Broad troposphere warming together with stratosphere cooling is another useful fingerprint: additional greenhouse gases and ozone changes have a different vertical response from simply increasing solar input. No one fingerprint is unique to one cause, but their joint agreement with greenhouse forcing strengthens the attribution.

Polar observations are especially informative because they connect atmosphere, ocean and ice. In the Arctic, reduced summer sea ice extent, thinning and loss of multi-year sea ice indicate a diminished ice reservoir. Warmer air, earlier melt and delayed freeze-up contribute, while wind and current changes redistribute ice and export thick floes. Pacific heat inflow and changes in Atlantic-origin waters supply oceanic routes to melting. The protection afforded by the cold halocline means that warmer Atlantic water does not automatically heat the ice: the intervening mixing and density stratification must be considered. The 2007 retreat illustrates how long-term thinning, anomalous winds and greater solar absorption can combine. It neither proves that weather variability has ceased nor provides a reason to ignore the longer-term energy imbalance.

Loss of grounded ice has a different global implication. Increased surface melt and runoff on Greenland, and changes in outlet-glacier discharge, transfer water from land to the ocean. In Antarctica, warming of the Peninsula and disintegration of some ice shelves demonstrate strong regional change. Warm ocean water reaching shelf cavities can thin floating shelves from below; loss of ice-shelf buttressing can then speed grounded-ice discharge. The floating shelf's disappearance is not itself equivalent to adding its whole volume to sea level. The important subsequent response is increased loss from grounded ice sheets. Increased snowfall can compensate some losses, so a mass budget must distinguish accumulation, melting and ice flow.

The two polar regions should not be assumed to behave identically. Around 2010, Antarctic sea-ice changes included strong regional contrasts and an overall satellite-era behaviour quite different from Arctic summer retreat. Winds, ozone-related atmospheric changes, freshwater effects and the circulation of the Southern Ocean influence the Antarctic pattern. A large cold continent surrounded by an ocean is physically different from an ocean surrounded by land, and vigorous Southern Ocean heat uptake can delay surface warming. A regional increase in Antarctic sea ice is therefore compatible with global ocean heat gain and losses from some grounded Antarctic ice; it still needs its own dynamical and thermodynamic explanation. Credible attribution accounts for these differences instead of requiring uniform warming everywhere.

Several mechanisms produce polar amplification. The ice-albedo feedback is the most direct: retreating snow and ice expose darker surfaces that absorb more sunlight, causing further warming and melt. Much of the extra energy over newly open Arctic water is stored during summer and released to the atmosphere in autumn and winter, when an ice-covered ocean would be better insulated. Consequently the largest near-surface atmospheric response need not occur in the season of greatest solar input. Sea-ice loss changes both reflection and the timing of ocean heat release.

The atmospheric vertical structure also matters. In the tropics, convective adjustment often produces appreciable warming aloft, which increases radiative cooling. Stable polar boundary layers can concentrate warming near the surface, producing a different lapse-rate feedback and less compensating cooling from upper layers. Furthermore, the idealized Planck feedback scales as $4\sigma T^3$: colder emitting regions have a smaller increase of blackbody radiation per kelvin, although a quantitative atmospheric comparison needs full radiative transfer. Water-vapour and cloud changes modify both short-wave and long-wave fluxes, with cloud effects depending on season and surface properties. Changes in poleward heat and moisture transport provide further coupling. Thus polar amplification has several interacting causes; it is not simply the same local ice-albedo process repeated throughout the year.

Polar change affects the rest of the planet. Grounded-ice loss contributes to eustatic sea level change, while ocean warming contributes to thermosteric sea level change; together they increase coastal flooding and erosion risks. Redistribution of ice and freshwater changes density stratification and dense-water formation, potentially altering the Atlantic meridional overturning circulation and other oceanic heat and carbon transports. A weakening is physically plausible under warming and freshening, but an immediate complete shutdown is not a necessary prediction. Permafrost carbon feedback can release additional carbon dioxide and methane, increasing warming outside the Arctic. Habitat loss, changing primary production and altered food webs affect marine ecosystems and people, while a longer open-water season changes navigation, coastal exposure and resource access. These impacts are not exhausted by a single global-mean temperature number.

Forecasts require an emissions scenario as well as physical understanding. On the roughly fifty-year horizon, continued emissions imply further ocean and atmospheric warming, a longer Arctic melt season, less late-summer ice and a greater fraction of young seasonal ice. Mountain glaciers generally continue to lose mass, and sea level continues rising through thermal expansion and land-ice loss. Variability can interrupt a regional trend for several years, and the exact first nearly ice-free Arctic summer cannot be fixed from the 2007 minimum. Ocean heat storage and slow ice-sheet adjustment also mean that some further change is committed even if emissions are reduced.

By approximately 2110, sustained high emissions would imply a substantially warmer world, potentially several degrees above early twenty-first-century global temperatures, with stronger Arctic warming and likely recurrent nearly ice-free late summers. This is a seasonal prediction, not a claim that the entire Arctic Ocean must be ice-free in winter. Losses of grounded ice and thermal expansion would produce larger cumulative sea level rise, plausibly measured in several decimetres and potentially more if ice dynamics accelerate. The largest uncertainty is not whether adding heat expands water, but the future emissions trajectory and the response of major ice sheets. Complete loss of Greenland on this century-long horizon is not required by that prediction. With strong mitigation, substantially smaller late-century warming and ice loss are possible, although sea level can continue to rise beyond stabilization of surface temperature.

Anthropogenic forcing explains the broad warming signal; feedbacks amplify polar change, while emissions and ice dynamics govern future outcomes. This distinction permits both a firm physical conclusion and honest uncertainty about dates, regional details and the magnitude of rapid ice-sheet change.

Take $x$ increasing upslope and let $u$ be the superficial Darcy velocity, so the mobile discharge is $uh$. Its pore transport speed is $u/\varphi$. This interpretation is required by the factors of porosity in the printed equations; a literal pore-speed convention would instead replace $u$ there by $\varphi u$.

For capillary residual trapping, let $M(x,t)={\max }_{0\le \tau \le t}h(x,\tau )$ be the maximum invaded thickness, including the initial state. The stored carbon-dioxide volume per plan area isDuring advance, $M=h$ and new pores fill with mobile carbon dioxide; during recession, $M$ is fixed and the newly vacated pores retain saturation $s$. Conservation of mobile plus trapped fluid isConsequentlyAssume $a,L,u,\varphi >0$ and $0<s<1$. The two speeds differ because the advancing front fills a whole pore volume while the receding tail removes only its mobile fraction.

The method of characteristics keeps height constant on $x=x_0+v_{A}t$ on the leading face and $x=x_0+v_{R}t$ on the trailing face. Matching the two linear profiles gives the triangular current with capillary retention:Before extinction, the rear, front, crest position and crest height areEvery positive height has one trailing and one leading characteristic; their intersections trace the crest. The mobile current disappears at

The final trapped region is the maximum-invasion envelope of a retained current. For $0<x<L$, the initial height is already the maximum, $M_{\infty }=ax$. For $L<x<x_{\ast }$, the maximum occurs when the crest passes $x$, at $t_m=2(x-L)/(v_R+v_A)$. Substituting into the trailing profile givesIf $z$ is distance into the aquifer measured normally from its upper boundary, trapped carbon dioxide occupies $0<z<M_{\infty }(x)$ at residual pore saturation $s$. The occupied pore volume per unit transverse width isexactly the initial mobile volume. This is an independent mass-conservation check of both the envelope and extinction distance. In the limit $s=0$, the current translates without shrinking and nothing is trapped; the finite-extinction formula is not used in that limit.

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.

Let $K$ be intrinsic permeability of a porous medium, $\varphi$ the mobile pore fraction, and let $h$ measure thickness normal to the sloping cap. Define cold and reservoir-temperature properties byThe imposed inequality makes ${\delta }_c>0$, so the cold fluid is still buoyant. Assume $K>0$, $\sin \theta >0$, positive densities and viscosities, and the usual positive thermal-expansion and viscosity coefficients. Away from localized fronts, neglect the along-current thickness gradient. Hydrostatic water pressure and Darcy's law then give the upslope Darcy velocitiesTypically $u_h>u_c$: warming increases buoyancy and decreases dynamic viscosity.

Current depths and the moving moving thermal front in a porous current. For this slope-driven porous gravity current, the injected mass flux per unit well length fixes the cold plateau thickness:The moving thermal front in a porous current moves at $V_T=\Gamma u_c$. Since $u_c$ is explicitly a Darcy velocity, the fluid mass stored per unit horizontal area is $\varphi \rho h$. Apply mass conservation across the moving moving thermal front in a porous current:Writing $\gamma =\varphi \Gamma$ and $r=u_c/u_h$, this gives the distal warm depthFor a retarded moving thermal front in a porous current take $0\le \gamma <1$; otherwise this two-plateau geometry needs reconsideration. The factor $(1-\gamma )/(1-\gamma r)$ is important: the distal mass flux is not generally $Q$, because mass is being stored as the cold region replaces the warm region. Setting the two fluxes equal would silently assume a stationary moving thermal front in a porous current.

If the model absorbs porosity into storage and uses pore velocity as its transport speed, the same formulas use $\gamma =\Gamma$. The explicit Darcy convention in the PDF instead gives $\gamma =\varphi \Gamma$; porosity must be specified or absorbed consistently. In the stationary-front limit $\gamma =0$, the two steady depths reduce to $Q/({\rho }_c{u}_c)$ and $Q/({\rho }_h{u}_h)$. In the Boussinesq approximation, use a common reference density in the mass factors while retaining ${\delta }_c,{\delta }_h$ in the driving force.

Leakage thresholds. A normal thickness $h_i$ generates cap overpressure ${\Pi }_i={\delta }_{i}g{h}_i\cos \theta$. Before leakage modifies the current, the two plateau values areFor the caprock leakage threshold $\Delta p>0$, the candidate critical mass-injection rates are thereforeCold-only leakage occurs for $Q_c^{\ast }\le Q<Q_h^{\ast }$, provided $Q_c^{\ast }<Q_h^{\ast }$. Both plateaus can leak once $Q\ge \max (Q_c^{\ast },Q_h^{\ast })$. In the usual common-density approximation, ${\mu }_c>{\mu }_h$ and $r<1$ ensure $Q_c^{\ast }<Q_h^{\ast }$, giving the expected sequence: no leakage, cold-only leakage, then leakage from both cold and warm regions. Without that approximation the ordering must be checked; the given inequality $\Delta \rho >\alpha \Delta T$ alone does not establish it. These are onset criteria computed on the nonleaking current, not a post-leakage mass budget.

For small slopes $\cos \theta \simeq 1$ and $\tan \theta \simeq \sin \theta \simeq \theta$. If another thickness convention is used, its hydrostatic column and projected flux must be changed consistently; one should not mix a normal thickness with a vertical-pressure formula lacking the cosine.

Cross-current heat conduction. Heat transfer from the warm formation makes temperature vary across the carbon dioxide depth and introduces a warming time controlled by thickness, thermal diffusivity and the surrounding rock's heat capacity. It smooths the sharp thermal adjustment and causes gradual changes in buoyancy, viscosity and the velocity profile. The cold fluid generally warms, becomes more mobile and requires less depth to carry a prescribed mass flux; the enhanced cold-region overpressure and its distinct leakage zone tend to shrink. Warming may occur before a parcel reaches the idealized advective moving thermal front in a porous current, especially for a thin current. A single supplied $\Gamma$ no longer describes all heat transport. Quantitative depths and thresholds then require a coupled temperature equation and thermal boundary data; their exact changes cannot be inferred from the linear property laws alone.

The usual major carbon/oxygen reservoirs in a hydrogen-rich hot Jupiter are water, carbon monoxide and methane, with their relative importance set by thermochemical equilibrium. This is not a universal ranking for every temperature and composition: nitrogen molecules or carbon dioxide can exceed a strongly depleted member of this trio.

At about one bar, the useful net reaction isThe rightward reaction is exothermic. Cooler gas favors methane and water; warming favors carbon monoxide and suppresses methane. The CO/CH4 crossover is of order $10^3\mathrm{K}$ and shifts with pressure, elemental inventory and metallicity; it is not a universal temperature. Water remains an important oxygen reservoir for oxygen-rich compositions but can dissociate at sufficiently high temperature.

Increasing the atmospheric metallicity of a giant planet raises the available carbon and oxygen. In a dilute, H2-dominated regime, major CO and H2O abundances roughly increase with the enrichment factor, and cool-regime CH4 does likewise. CO2 can rise faster, approximately quadratically in suitable warm regimes. At very high enrichment, the H2 fraction and mean molecular weight also change, invalidating simple linear scalings.

At high temperature and atmospheric carbon-to-oxygen ratio below unity, CO binds much of the carbon, leaving excess oxygen for H2O. As C/O approaches or exceeds one, CO consumes nearly all available oxygen and H2O is strongly depleted; excess carbon enhances CH4, hydrogen cyanide and acetylene. Cooler CH4-dominated chemistry uses less oxygen in CO and can leave more water even at the same elemental ratio. These trends are carbon partition and atmospheric water abundance, and assume equilibrium rather than vertical or horizontal quenching.

The three parameter groups are local temperature, local pressure, and elemental composition. The elemental composition includes both atmospheric metallicity of a giant planet and atmospheric carbon-to-oxygen ratio; these are independent compositional controls, so a specification by $T$ and $P$ alone is incomplete. Given these inputs, thermochemical equilibrium minimizes Gibbs free energy with element conservation.

Three useful chemical observables in a hydrogen-rich hot Jupiter are:

CO/CH4, H2O, and CO2 bands provide complementary temperature and elemental-abundance diagnostics, but their observed strengths also depend on the atmospheric pressure-temperature profile, exoplanet clouds and observing geometry. A spectral feature is not by itself a direct equilibrium abundance.

Take $z$ upwards and write $\zeta =z+Wt$, so that the unperturbed interface is $\zeta =0$, with carbon dioxide above and water below. Use the sharp-interface, incompressible flow version of Darcy's law, with constant porosity and permeability of a porous medium, no capillary pressure, and effectively unbounded layers. The distinction between Darcy flux and interface speed matters: the base flux is $q_z=-\varphi W$, since the printed $W$ is the speed of the interface. Thus the base pressure gradients satisfyThe mass density ordering is stabilizing, whereas injection of the less viscous fluid into the more viscous fluid produces the Saffman–Taylor instability.

For linear stability analysis, displace the interface by $\eta =\hat{\eta }{e}^{\sigma t+i\alpha x}$. Taking the divergence of Darcy's law shows that the perturbation pressure satisfies the Laplace equation. Decay away from the interface gives, for $a_{\alpha }=|\alpha |>0$,The kinematic boundary condition and continuity of normal Darcy flux giveMeanwhile pressure continuity, evaluated on the displaced interface, givesEliminating $A,B$ yields the buoyancy-modified Darcy fingering dispersion relation:Every nonzero wavenumber decays when $W<W_c$, grows when $W>W_c$, and is neutral at equality, whereThe zero wavenumber is a uniform displacement and is neutral. Without a short-scale regularization such as surface tension, this model has no fastest-growing finite wavenumber: on the unstable side its growth rate increases without bound as $|\alpha |$ increases. Finite layer depths would change the dispersion relation; none are specified here.

Taking $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$ and a year of $365.25$ days, the data giveThe numerical paragraph supplies no porosity, so it cannot determine a unique numerical interface speed. If $W$ were instead intended to mean Darcy velocity, the numerical critical value would be $103.2\mathrm{m}\mathrm{y}{\mathrm{r}}^{-1}$, and the interface would move at $W/\varphi$. These two velocity conventions must not be conflated.

For solar elemental abundances, take a molecular hydrogen- and helium-dominated atmosphere with atmospheric carbon-to-oxygen ratio near $0.55$. The principal oxygen-bearing molecule is generally water; the carbon and nitrogen carriers depend on both temperature and pressure.

At the cool $500\mathrm{K}$ observable upper layers, methane is the expected main carbon reservoir under chemical equilibrium, with abundant water and much less carbon monoxide. At the hotter layers approaching $1400\mathrm{K}$, carbon monoxide becomes the main carbon carrier and water contains much of the oxygen not bound in it. The relevant law of mass action follows fromLower temperature favors the exothermic methane-forming direction, while increasing pressure favors the side with fewer molecules. This explains why the change of dominant carrier cannot be specified by temperature alone.

Nitrogen is distributed between molecular nitrogen and ammonia. Cooling favors ammonia, but low pressure favors molecular nitrogen; it is therefore unsafe to call ammonia dominant throughout the $500\mathrm{K}$ low-pressure region. The exact partition requires the equilibrium constant for ${\mathrm{N}}_2+3{\mathrm{H}}_2\rightleftharpoons 2{\mathrm{NH}}_3$.

Carbon dioxide is usually a minor constituent at solar composition; hydrogen cyanide and acetylene are much less abundant than the principal carbon carriers in this oxygen-rich equilibrium case. Condensation can remove refractory species where a condensation curve is crossed. The robust cool-atmosphere expectation is a molecular hydrogen–helium background with methane and water, changing toward carbon monoxide in hotter layers. Precise mixing ratios require thermodynamic data and an explicitly specified elemental inventory after rainout.

At fixed temperature and pressure, the atmospheric carbon-to-oxygen ratio controls the division of the elemental inventory among carbon monoxide, methane, water and other molecules. For an oxygen-rich mixture, water can remain after carbon monoxide consumes its share of oxygen. Increasing the atmospheric carbon-to-oxygen ratio reduces the oxygen available for water, particularly in the hotter layers where carbon monoxide is stable. Near or above unity there, water can become strongly depleted while excess carbon feeds methane, hydrogen cyanide and acetylene. This familiar hot carbon-rich behavior should not be applied unchanged to every cool layer: at $500\mathrm{K}$, methane formation can leave substantial water even at high atmospheric carbon-to-oxygen ratio.

Increasing the atmospheric metallicity of a giant planet at fixed elemental ratios increases the heavy-element inventory relative to molecular hydrogen and helium. Abundances of water and the major carbon-bearing molecules generally rise while molecular hydrogen remains dominant. In that regime, the law of mass action for $\mathrm{CO}+{\mathrm{H}}_2\mathrm{O}\rightleftharpoons {\mathrm{CO}}_2+{\mathrm{H}}_2$ givesWhere carbon monoxide and water each scale roughly linearly with enrichment $Z$, carbon dioxide consequently scales approximately as $Z^2$. This scaling changes when molecular hydrogen ceases to dominate or chemical partitioning changes. Higher enrichment can also favor carbon monoxide over methane at fixed conditions, increase mean molecular weight, and reduce the atmospheric scale height and transmission-feature amplitudes. Cloud formation and atmospheric condensate rainout alter the observable elemental ratios relative to the bulk ones.

Use displacements $x=(x_1,x_2,x_3{)}^T$ of the left oxygen, carbon and right oxygen along a common positive direction. In this harmonic model the kinetic energy and elastic potential energy areThus the generalized eigenvalue problem for small oscillations is $Kv={\omega }^2\mathsf{M}v$, whereThe normal modes of three masses connected by two equal springs give the three independent solutionsThe second row follows by holding the centre fixed and moving the end masses oppositely. In the remaining internal mode the end displacements are equal, and zero displacement of the centre of mass gives $2m{v}_1+M{v}_2=0$, yielding the third row; substitution verifies its frequency. Both nonzero modes have fixed centre of mass. In the symmetric stretch the two bond-length changes are equal; in the antisymmetric stretch they are opposite.

The zero mode occurs because a uniform translation changes neither spring length, so $K{v}_0=0$. It represents free center-of-mass motion, not an additional restoring vibration. Its general time amplitude is $A+Bt$. Only one translational zero mode exists in the restricted one-dimensional model; transverse and rotational motion are outside its specified coordinates.

For the initial impulse, all displacement amplitudes at zero time vanish. Resolve the initial velocity into these three mode shapes:Writing ${\omega }_s=\sqrt{k/m}$ and ${\omega }_a=\sqrt{k/m+2k/M}$ gives the impulse on an end mass of a free symmetric three-mass chain:In particularwhere $V_{\mathrm{cm}}=mu/(M+2m)$. The initial derivatives are $(u,0,0)$, as required. The carbon dioxide model therefore drifts uniformly while undergoing both internal stretches. Its total momentum is the initial $mu$; omitting the zero mode would incorrectly hold the centre of mass fixed.

The useful mechanism is radioactive decay during retention, not chemical destruction of carbon dioxide. With the time in seconds,Here $a$ is a radioactive activity concentration; the supplied units are per volume rather than per mass. The original plot follows this concentration over four hours.

A practical illustrative choice is a three-hour guaranteed residence time, at which $a/a_0=0.00214175$, giving about a factor $467$ reduction. On a plug-flow basis, the carbon-dioxide-equivalent hold-up volume isUse a contained delay line or a validated retention arrangement that prevents early breakthrough, followed by monitored discharge. A retention bed could instead hold the radioactive carbon dioxide, but its capacity and breakthrough performance would have to be established. Gas retention and monitored discharge are waste-treatment functions, independent of the eventual use of the isotope.^[ref]

The $4.00$ L result is a flow-delay volume, not a fill-and-immediately-dump tank specification. It credits radioactive decay during the guaranteed transit time and no subsequent storage after discharge from the delay system. Gas entering a filling tank near the end of filling has not aged for the entire filling time. For a literal tank filled for $T$ and immediately emptied, the radioactive activity at emptying is insteadThat alternative must not be analysed by multiplying the whole collected inventory by $e^{-\lambda T}$. If “no post filling radioactive decay” is intended to prescribe this literal fill-and-dump arrangement, the three-hour calculation does not establish compliance for it.

The nitrogen fraction is also unspecified. If $y$ is the carbon dioxide volume fraction at the stated volume reference conditions, a gas delay line requires $V_{\mathrm{mixture}}=V_{\mathrm{C}{\mathrm{O}}_2}/y$, not necessarily $4.00$ L of physical vessel volume. The graph alone does not define a unique optimum: a cost objective, operating calendar and residence-time model are needed. Parts (c) and (d) quantify the discharge constraint for the stated guaranteed-delay design.

Dating principle. Atmospheric cosmic rays produce secondary neutrons, which can generate carbon-14 through ${}^{14}\mathrm{N}(n,p){}^{14}\mathrm{C}$. Carbon then enters atmospheric carbon dioxide and living organisms. While an organism exchanges carbon with its environment, its carbon-14 to carbon-12 ratio approximately tracks that reservoir. After exchange stops, surviving carbon-14 undergoes beta-minus decay, whereas the stable carbon remains. On the closed-sample, constant-initial-ratio model, radioactive decay givesThus radiocarbon dating infers elapsed time from the loss of an isotope ratio or of radioactive activity per fixed carbon mass. It dates cessation of carbon exchange, which need not equal the date an object was manufactured or buried. Actual calendar dating also requires radiocarbon calibration and control of contamination and reservoir differences.

Mean lifetime. For exponential decay, the half-life obeys $e^{-\lambda T_{1/2}}=1/2$, while the mean lifetime is $\tau =1/\lambda$. HenceUsing $1\mathrm{yr}=365.25\times 86400\mathrm{s}$ gives $\lambda =3.83325\times 10^{-12}{\mathrm{s}}^{-1}$ and $\tau =2.60875\times 10^{11}\mathrm{s}$.

The mass minimum and beta stability. Treat the supplied masses as nuclear rest energies in $\mathrm{MeV}$, and omit nuclear pairing as instructed. At fixed mass number $A$, the volume and surface terms of the semi-empirical mass formula are constant. Expanding the remaining terms gives a convex nuclear mass parabola,The continuous minimum in atomic number is thereforeAt $A=14$, the quadratic coefficient is $7.01302\mathrm{MeV}$ and $Z_{\ast }=6.79450$. The nearest integer is $Z=7$; the nuclear isobar of minimum mass is nitrogen-14. Maximizing nuclear binding energy alone would omit the neutron-proton rest-mass difference; its continuous optimum is slightly different, but also selects integer $Z=7$ here.

The values most useful for checking neighboring decay thresholds areConsequently carbon-14 can lower its nuclear mass by changing a neutron into a proton:For bare nuclear masses, the Q value must include the emitted electron mass. Taking $m_e{c}^2=0.511\mathrm{MeV}$,The opposite proton-to-neutron conversion would lead to $Z=5$, whose nuclear mass is $18.1567\mathrm{MeV}$ higher than that of carbon-14. Neither beta-plus decay nor electron capture can bridge that gap by supplying an electron-mass correction. For the nitrogen-14 daughter, both neighboring nuclear masses are higher: $Z=6$ by $4.13070\mathrm{MeV}$ and $Z=8$ by $9.89535\mathrm{MeV}$. Even electron capture to $Z=6$ has $Q=0.511-4.13070<0$. Thus the daughter has no further energetically allowed single beta decay in this model; physical nitrogen-14 is stable.

These threshold numbers belong to the supplied smooth semi-empirical mass formula, not to an accurate prediction of the carbon-14 electron spectrum or its long half-life. Semi-empirical mass errors can change beta-decay thresholds. Also, the instruction to omit nuclear pairing is a model simplification: even $A=14$ does not make physical pairing vanish, and even-even carbon and odd-odd nitrogen need not have the same pairing correction. No omitted pairing coefficient is guessed here.

Activity of modern carbon. Neglect the small contribution of other stable carbon isotopes when counting atoms in a gram. The Avogadro constant givesThe expected radioactive activity of that one gram isor approximately $15.0$ decays per minute per gram of carbon. A becquerel is one decay per second; the normalization is per gram of carbon, rather than per gram of the entire organic sample.

Artefact age. Its remaining activity fraction is $0.13/0.25008=0.51983$. Inverting the exponential decay law,This is the radiocarbon age on the supplied constant-initial-activity, $5730$-year half-life model. It is not a calibrated calendar date; the sample's historical initial ratio and any later carbon exchange are not supplied.

The constants $R_a,R_b,R_d$ are effective reservoir capacities: stored carbon dioxide mass anomaly per unit effective concentration anomaly, in GtCO${}_2$ per ppmv. Only $C_a$ is a literal atmospheric mole fraction; $C_b,C_d$ are scaled reservoir variables chosen to share its equilibrium value. The exchange coefficients $r,s$ are conductances, in GtCO${}_2$ per ppmv per year: $r(C_a-C_b)$ is transfer from the atmosphere to the fast reservoir and $s(C_b-C_d)$ transfer from that reservoir to the deep ocean. They are not by themselves inverse times; division by a reservoir capacity gives an inverse-time rate.

Hydrostatic balance estimates total atmospheric mass as $M_{\mathrm{air}}=4\pi R_E^2\bar{p}/g$. One ppmv of a dilute ideal gas is a mole fraction $10^{-6}$, so, with dry-air molar mass $M_{\mathrm{mol},\mathrm{air}}\simeq 28.97$ g/mol and carbon dioxide molar mass $44.01$ g/mol,The final factor converts kilograms to gigatonnes. This is carbon dioxide mass, not the smaller carbon mass.

For the linear three-reservoir carbon model, with concentrations as a column vector, the linear ordinary differential equation isA homogeneous normal mode has $\mathbf{x}=\mathbf{v}{e}^{\lambda t}$, where $J\mathbf{v}=\lambda \mathbf{v}$. For negative real eigenvalues the e-folding relaxation time is $\tau =-1/\lambda$. The row sums vanish, so ${\mathbf{v}}_3=(1,1,1{)}^T$ has ${\lambda }_3=0$ andEqual effective concentrations generate no exchange. Summing the reservoir budgets givesThus injected carbon dioxide is conserved after emissions cease: the model redistributes it but contains no permanent removal. The zero normal mode records this conserved excess mass, rather than a long but finite decay.

To check the other normal modes are damped, set $D=\mathrm{diag}(R_a,R_b,R_d)$ and $L=-DJ$. ThenThe similarity transformation $D^{1/2}J{D}^{-1/2}=-D^{-1/2}L{D}^{-1/2}$ is a symmetric negative-semidefinite matrix with a one-dimensional kernel, so the remaining two eigenvalues are strictly negative.

Here $R=7.8$, $R_d=78$, $r=0.4$, and $s=0.04$, in the stated reservoir units. The exact characteristic coefficients are $P=0.10820513{\mathrm{yr}}^{-1}$ and $Q=0.000315582{\mathrm{yr}}^{-2}$, givingThe leading separation estimates give $9.75$ yr and $325$ yr if finite $R_d$ is retained, or $390$ yr if $R_d/R$ is also sent to infinity. Thus ${\tau }_1\ll 100\mathrm{yr}<{\tau }_2$; a ratio of ten warrants using the exact values for the requested sketch.

Assume initially zero anomalies in all reservoirs. The atmospheric impulse response, obtained by decomposing $e^{Jt}$ into the three normal modes, iswith $w_0=1/93.6=0.01068376$, $w_1=0.06087812$, and $w_2=0.05664324$, in ppmv/GtCO${}_2$. The checks $w_0+w_1+w_2=1/R_a$ and $w_1/{\tau }_1+w_2/{\tau }_2=r/R_a^2$ fix the instantaneous injection response and its initial decay.

The concentration is the convolution $C_a(t)={\int }_0^{t}G(t-t^{\prime })E(t^{\prime })d{t}^{\prime }$. To evaluate it transparently, define $H_0(q)=q^2/2$ for $q>0$, andsetting every $H_i(q)=0$ when $q\le 0$. Expressing the triangular emission history as $E(t)=0.5[t_{+}-2(t-100{)}_{+}+(t-200{)}_{+}]$ givesHere times are in years and $0.5$ has the corresponding emission-ramp units. Differentiating or evaluating the convolution givesThe maximum concentration occurs earlier than the end of emissions, at approximately year $172$, with $C_a\simeq 280.7$ ppmv. Initially there is a short adjustment over ${\tau }_1$; the curve then rises, turns over as emissions fall, and continues declining after year 200 as the deep ocean takes up carbon dioxide.

The total emitted mass is the area of the emission triangle, $\frac{1}{2}(200)(50)=5000$ GtCO${}_2$. All reservoirs eventually share one concentration, so conservation of mass gives