Earth Energy Imbalance & Radiative Transfer Simulator

Documentation Overview

This document was inspired by a denier citing an argument in a paper co-authored by a former MIT Professor of Physics. I'm no professor. The simulator is based on mathematics from the referenced IPCC report.

It is possible that Prof Lindzen is right, and this is wrong. The atmosphere's sensitivity to increases in carbon would then be much lower and the peril much less.

However, the simulator uses the IPCC assumptions from the start of the Keeling measurements of actual atmospheric carbon, dating back to the 1950s. It then arrives at roughly the measured atmospheric warming today.

It is possible that the warning will tail off, and the worst predictions start to match the Lindzen theory. But we won't know that until we're there, and if the model is right, we will then be in deep trouble, and beyond any possibility of returning our climate to the one in which human civilization has thrived.

The only option available is to stop polluting our atmosphere with waste carbon molecules.

For a detailed critique of Lindzen's paper start here:
https://www.carbonbrief.org/a-disservice-to-the-scientific-method-climate-scientists-take-on-richard-lindzen#

                Authorship Verification Stamp (SHA-256): 2308c58eff8d7bacfcbfcb553eb4985f25abf8dc1a55234c163a31a5b243d02f
              

Select a chapter below or use the navigation tabs above to review the physical foundations, closed-form mathematical proofs, and empirical calibration data powering this simulation.

📐 Page 1: Mathematical Proofs & Model Verification
Complete analytical proofs: Equilibrium Climate Sensitivity (ΔT_eq = ΔF / λ), Stefan-Boltzmann radiating altitude inversion in SI metres (z_eff), and eigenvalue decomposition of ocean thermal lag.
☀️ Page 2: Vertical Radiative Transfer & The "Saturation Myth" Refutation
Detailed physical breakdown of the five atmospheric pins: Solar shortwave transparency and albedo reflection, greenhouse 15 µm molecular resonance, barometric optical depth thinning (τ ≈ 1), and radiative choking at altitude (σT⁴).
🌊 Page 3: IPCC AR6 Calibrated Two-Layer Ocean EBM IPCC AR6 Calibrated 2-Box EBM
Governing differential equations (Held et al. 2010; Geoffroy et al. 2013), parameter calibration, historical hindcast integration from 1750 unforced equilibrium, and harmonization with the standard 1850–1900 observational baseline (+0.12°C offset).
🧊 Page 4: 1-Million-Year Paleoclimate Reconstructions & Sources EPICA / LR04 / Temp12k
Benchmark ice-core and marine sediment datasets synthesized in the illustrative 1,000,000-year / 500-year rate comparator (EPICA Dome C, Bereiter et al. 2015, Lisiecki & Raymo LR04, Shakun et al. 2012, Kaufman Temp12k), temporal resolution limits, and glacial termination rates.
🛢️ Page 5: Fossil Fuel Metrics — Barrels of Oil & Gas to Atmospheric CO₂ (ppm) & Temperature Rise Combustion Stoichiometry & ppm Mapping
The direct physical and mathematical chain linking physical fuel commodities (million/billion barrels of crude oil, trillion cubic feet of gas, tonnes of coal) to atmospheric concentration (ppm) and warming. Explains atmospheric mass (5.15 × 10⁶ Gt), the universal 7.82 Gt CO₂ / ppm factor, airborne fraction (AF ≈ 46%), why it takes ~39.5 billion barrels of oil to add +1.0 ppm, and how each +10 ppm commits Earth to +0.10 °C of long-term equilibrium warming.
📉 Page 6: Global Carbon Budgets & The "Unburnable Carbon" Reality IPCC AR6 / TCRE / Paris Limits
The Transient Climate Response to Cumulative Emissions (TCRE ≈ 0.45 °C per 1,000 Gt CO₂), remaining carbon budgets for +1.5 °C (~200 Gt CO₂) and +2.0 °C (~1,150 Gt CO₂), translating budgets into remaining years of current extraction, and comparing finite atmospheric absorptive capacity against the ~3,400 Gt CO₂ of proven fossil reserves (the "Unburnable Carbon" paradox where >90% of reserves must stay unextracted).
🏭 Page 7: Global Emissions Anatomy — Sectors, Fuels & Non-CO₂ Warming Gases (Our World in Data) OWID / Climate Watch / GCP
Complete quantitative breakdown of global annual emissions (54.6 Gt CO₂e total, ~40.9 Gt CO₂). Analyzes the 4 primary economic sectors: Energy (73.2% — Industry 24.2%, Transport 16.2%, Buildings 17.5%, Fugitive 5.8%), Agriculture & Land Use (AFOLU, 18.4% — livestock methane, fertilizer N₂O, deforestation), Direct Industrial Processes (5.2% — cement calcination, chemicals), and Waste (3.2%). Maps fossil fuels by source (Coal 41%, Oil 33%, Gas 21%) and pairs each emission category with its direct technological replacement.
💡 Page 8: The Clean Energy & Storage Transition — Replacing 41 Gt/yr of Fossil Energy Electrification Factor & Storage Physics
The physical mathematics of replacing 41 Gt/yr of fossil energy. Explains the "Electrification Efficiency Dividend": why 135,000 TWh of fossil primary energy (combustion loss >60%) only requires ~75,000 TWh/yr of clean useful electricity. Examines the generation mix (Solar, Wind, Nuclear baseload, Hydro, Geothermal) and evaluates six grid storage technologies across their governing physical laws: Chemical Batteries (Li/Na/Flow), Pumped Hydro (m·g·h), Compressed Air (P·dV), Gravity Storage, Flywheels (½·I·ω²), and Hydrogen seasonal storage.
🌡️ Addendum: Present-Day Forced Warming (+1.26°C) vs. Headlines Reporting +1.5°C Breached Climate Benchmark Clarification
Why the simulator models +1.26°C of forced warming for 2026 while news media reported 2023–2024 breached +1.5°C. Details the scientific distinction between single-year weather spikes (El Niño, solar maximum) and the multi-decadal forced climate response (IPCC AR6 standard), pre-industrial baseline differences (1750 vs. 1850–1900), and ocean heat uptake inertia.

Page 1: Mathematical Proofs & Model Verification

The mathematical validity of the simulation can be verified directly using closed-form analytical solutions of the governing differential equations.

Proof 1: Equilibrium Climate Sensitivity (ECS) Closed-Form Limit

At thermal equilibrium, time derivatives vanish (dΔT_s/dt = 0, dΔT_d/dt = 0). Consequently, the mixed layer and deep ocean equilibrate to the exact same temperature (ΔT_s = ΔT_d = ΔT_eq). Substituting this into the mixed-layer ordinary differential equation:

C_m · (dΔT_s / dt) = ΔF - λ·ΔT_s - γ·(ΔT_s - ΔT_d) = 0
ΔF - λ·ΔT_eq - γ·(ΔT_eq - ΔT_eq) = 0
ΔF - λ·ΔT_eq = 0 ⟹ ΔT_eq = ΔF / λ

For a doubling of atmospheric CO₂ (C / C₀ = 2), the effective radiative forcing (IPCC AR6 Working Group I, Chapter 7) is ΔF_2× = 3.93 W/m² (giving α = 3.93 / ln(2) ≈ 5.670 W/m²). With Equilibrium Climate Sensitivity ECS = 3.00 °C, the climate feedback parameter is derived from F_2× / ECS as λ = 3.93 / 3.00 = 1.310 W/m²/K:

ΔT_eq = 3.930 W/m² / 1.310 W/m²/K = 3.000 °C

Verification: Clicking "3. Instant 2× CO₂ Pulse" and then "⚡ Instant Ocean Equilibrium" sets ΔT to exactly +3.00 °C and Net Imbalance (EEI) = 0.00 W/m².

What "Instant Ocean Equilibrium" Means Physically:
In the real Earth system, surface air warming lags behind greenhouse gas radiative forcing by decades to centuries because the deep ocean acts as an immense planetary thermal buffer (absorbing over 90% of excess trapped heat, a delay termed ocean thermal lag or transient warming).

The "⚡ Instant Ocean Equilibrium" control mathematically evaluates the asymptotic steady-state limit (t → ∞):

Equalization of Thermal Layers: Forces both the upper mixed layer (C_m) and deep ocean abyss (C_d) to reach identical temperature anomalies (ΔT_s = ΔT_d = ΔT_eq = ΔF / λ).
Elimination of Downward Ocean Heat Flux: When surface and deep ocean temperatures equalize, downward heat uptake ceases entirely (γ·(ΔT_s − ΔT_d) = 0).
Planetary Radiative Balance Restored: Outgoing Longwave Radiation (OLR) increases by exactly λ·ΔT_s, bringing Top-of-Atmosphere net imbalance back to zero (EEI = 0.00 W/m²).
Committed Warming in the Pipeline: At current 2026 conditions (427 ppm CO₂, ΔF ≈ 2.40 W/m²), transient warming is ~+1.26 °C while ~0.78 W/m² continues flowing into the ocean sink. Clicking Instant Ocean Equilibrium reveals the +1.83 °C of warming that is already physically locked into the pipeline if atmospheric CO₂ were held constant today.

Proof 2: Analytical Derivation of Effective Radiating Altitude (z_eff in Metres)

Earth's Outgoing Longwave Radiation (OLR) escaping to space originates from the altitude where the atmosphere becomes optically thin. Inverting the Stefan-Boltzmann law:

OLR = σ · (T_eff)⁴ ⟹ T_eff = (OLR / σ)^1/4

Coupling this with the standard tropospheric lapse rate (Γ = 0.0065 K/m, or 6.5 °C/km), where T(z) = T_s - Γ·z:

T_eff = T_s - Γ · z_eff ⟹ z_eff (metres) = (T_s - T_eff) / 0.0065

At pre-industrial equilibrium (OLR = 239.0 W/m², T_s = 13.85 °C = 287.00 K, σ = 5.670374 × 10⁻⁸ W/m²/K⁴):

T_eff = (239.0 / 5.670374×10⁻⁸)^1/4 = 254.80 K (-18.35 °C)
ΔT = 287.00 K - 254.80 K = 32.20 K
z_eff = 32.20 K / 0.0065 K/m = 4,953.8 metres (displayed as 4,954 m)

Verification: In the pre-industrial baseline state, the app displays 4,954 m (-18.4 °C), matching the analytical calculation.

Proof 3: Eigenmode Decomposition of Ocean Thermal Lag

The two-layer ocean differential system exhibits two distinct physical eigenvalues corresponding to the fast upper-ocean response and slow deep-ocean heat uptake:

τ_fast ≈ 3.6 years (Upper ocean mixed layer warming)
τ_slow ≈ 230 years (Deep ocean thermal diffusion)

Because human emissions ramp up year after year (+2.5 ppm/yr), surface warming continuously lags radiative forcing, sustaining a modeled Earth Energy Imbalance of ≈ +0.78 W/m² (consistent with modern satellite- and Argo-observed estimates of ~0.8 to 1.0+ W/m²).

Page 2: Radiative Transfer & The "Saturation Myth" Refuted

The five interactive markers across the atmospheric column trace the physical journey of solar and terrestrial energy:

Step	Layer / Altitude	Physical Mechanism
1	Top of Atmosphere (13,000 m)	Solar Shortwave Transparency: The Sun's intense temperature (~5,778 K) concentrates solar radiation into the shortwave spectrum (spanning ultraviolet, visible light, and near-infrared). The atmosphere is broadly transparent to incoming solar radiation across this spectrum because greenhouse gases like CO₂ lack vibrational resonance modes at these high frequencies, allowing incoming solar energy to reach the planetary surface unimpeded (~239 W/m² post-albedo).
2	Earth Surface (0 m)	Surface Wavelength Conversion (Longwave IR): Dark ground and oceans absorb incoming shortwave light and warm to ~288 K (~15 °C). In accordance with Wien's law, a 288 K planetary surface re-emits in the longwave thermal infrared (peaking at 10–15 µm). Because this ~15 µm wavelength matches the exact molecular bending resonance of CO₂, greenhouse gases trap energy on the way out, not on the way in. (Bright ice reflects shortwave light back to space unchanged. Without greenhouse gases, Earth would be ~33 °C colder, plunging surface temperatures to ~18 °C below zero; Earth froze into a snowball twice in geological history).
3	Lower Troposphere (2,500 m)	Photons excite molecules; energy equilibrates into surrounding air via thermal collisions, percolating upward.
4	Upper Troposphere (5,000 m)	Exponential barometric density thinning allows optical depth to space to drop (τ ≈ 1), escaping freely.
5	Mid/Lower Column	Adding greenhouse gases pushes the escape layer higher into colder air. Radiating less energy (σT⁴) traps excess heat below.

☀️ Fundamental Physics: Why CO₂ Traps Outgoing Heat but Ignores Incoming Sunlight
A common point of confusion is why greenhouse gases do not simply absorb heat from the Sun as sunlight travels down through the atmosphere. The reason is the profound difference in wavelength between solar radiation and terrestrial heat (governed by Planck's Law and Wien's Displacement Law, λ_peak ∝ 1 / T):

Inbound Solar Radiation (Shortwave Spectrum, ~0.2 to 4 µm): The Sun's surface is blisteringly hot (~5,778 K / ~5,500 °C). High temperature compresses radiation into short, energetic wavelengths spanning ultraviolet, visible light (0.4 to 0.7 µm), and near-infrared. The atmosphere is transparent to incoming solar radiation across this entire spectrum because the linear carbon dioxide molecule (O=C=O) has electronic transitions only in the far ultraviolet and vibrational modes in the thermal infrared; it has zero absorption bands across the solar shortwave spectrum. Thus, incoming solar radiation passes through the atmospheric column unimpeded.
Surface Energy Transformation: Sunlight strikes Earth's oceans and continents. Absorbed shortwave energy warms the planet to an average surface temperature of ~288 K (~15 °C).
Outbound Terrestrial Radiation (Longwave Thermal Infrared, ~15 µm): Because Earth is ~20 times cooler than the Sun, it radiates at wavelengths roughly ~20 to ~30 times longer—in the thermal infrared band (4 to 50 µm).
The 15 µm Molecular Resonance: The fundamental bending vibration (ν₂ vibrational mode) of the CO₂ molecule resonates at an optical wavenumber of 667 cm⁻¹, which corresponds to a wavelength of exactly 14.99 µm (~15 µm). This aligns precisely with the peak of Earth's thermal emission curve. Therefore, while incoming solar radiation slipped through without interaction, outgoing thermal infrared is vigorously absorbed by CO₂ and thermalized into the atmosphere.
The Essential Greenhouse Blanket & Snowball Earth: Greenhouse warming is not intrinsically bad—it is essential to maintaining liquid water and complex life. Without carbon dioxide and other greenhouse gases, Earth's bare blackbody radiating equilibrium (with planetary albedo A ≈ 0.30) is T_eff = (S₀(1 - A) / 4σ)^1/4 ≈ 255 K (-18 °C / 0 °F). The natural greenhouse effect provides approximately +33 °C of warming, lifting average global surface temperatures to a comfortable +15 °C. When atmospheric greenhouse gas levels plummeted in geological deep time, Earth experienced runaway ice-albedo freezes into a global snowball twice: the Huronian Glaciation (~2.4–2.1 Ga) ↗ (when early photosynthesizers oxygenated the atmosphere and destroyed ancient methane) and the Cryogenian Snowball Earth (~720–635 Ma) ↗ (comprising the Sturtian and Marinoan glaciations, when the breakup of supercontinent Rodinia spurred rapid silicate weathering; see also Britannica: Snowball Earth Hypothesis ↗).

Page 3: Two-Layer Ocean EBM Formulation

Based on Geoffroy et al. (2013) and the IPCC AR6 Chapter 7 Report (Cross-Chapter Box 7.1) ↗:

C_m · (dΔT_s / dt) = ΔF - λ·ΔT_s - γ·(ΔT_s - ΔT_d)
C_d · (dΔT_d / dt) = γ·(ΔT_s - ΔT_d)

Parameter	Symbol	Model Value	Source & Status	Physical Meaning
Mixed Layer Heat Capacity	C_m	7.3 W·yr/m²/K	Published: Geoffroy et al. (2013) Table 4 CMIP5 multi-model mean (7.3 ± 1.4)	Thermal capacity of the atmosphere and upper ~75m ocean mixed layer
Deep Ocean Heat Capacity	C_d	106.0 W·yr/m²/K	Published: Geoffroy et al. (2013) Table 4 CMIP5 multi-model mean (106 ± 47)	Thermal capacity of the intermediate and deep ocean abyss
Heat Uptake Efficiency	γ	0.73 W/m²/K	Published: Geoffroy et al. (2013) Table 4 CMIP5 multi-model mean (0.73 ± 0.20)	Downward turbulent and thermohaline heat transfer
Net Climate Feedback	λ	1.310 W/m²/K	Derived: IPCC AR6 Working Group I, Chapter 7 (F_2× / ECS = 3.93 / 3.00)	Net radiative restoration response per degree of surface warming

Historical Hindcast & The Baseline Question: Why 1750 vs. 1850–1900?

A common question when interpreting climate models is why the differential equations integrate from 1750 CE when virtually all public targets (the UNFCCC Paris Agreement 1.5 °C / 2.0 °C goals, IPCC policy summaries, and tipping-point literature) are referenced to the 1850–1900 baseline:

🌍 1750: The Physical Equilibrium Baseline

Why climate physics starts in 1750: The fundamental greenhouse forcing equation ΔF = α·ln(C / C₀) and the two-layer differential equations require an unperturbed thermodynamic rest state. In 1750, atmospheric CO₂ was stable at ~278–280 ppm (from Antarctic ice cores), ocean heat uptake was in balance, and the planetary energy imbalance was zero (EEI = 0, ΔT_s = 0, ΔT_d = 0). Starting here avoids having to guess arbitrary non-zero initial ocean heat fluxes.

🌡️ 1850–1900: The Observational Thermometer Baseline

Why public policy starts in 1850–1900: Reliable worldwide instrumental thermometer networks (HadCRUT, Berkeley Earth, NOAA) only began around 1850. Because global surface thermometers did not exist in 1750, the IPCC and Paris Agreement operationally defined "pre-industrial" as the 1850–1900 multi-decadal observational mean.

Bridging the ~0.12 °C Early-Industrial Offset: Between 1750 and the late 19th century, early coal burning and agricultural land clearance nudged atmospheric CO₂ from 280 ppm to ~291 ppm, warming the planet slightly before global thermometer networks existed. In IPCC AR6 Working Group I (Chapter 2), paleoclimate reconstructions estimate this early warming at ~0.1 °C (likely range -0.05 °C to +0.23 °C).

When our two-box EBM hindcast runs forward from 1750 through historical Law Dome and Mauna Loa Keeling data, it calculates an average 1850–1900 surface warming of exactly +0.121 °C above 1750 (in near-perfect agreement with the IPCC central estimate). The simulator harmonizes both frameworks:

Physical Radiative & Ocean Integration: Solved forward from the unforced 1750 equilibrium (ΔT = 0, EEI = 0) with zero artificial initialization fudge factors.
Tipping Points & Chart Reference Lines: Evaluated relative to the standard 1850–1900 Paris baseline by accounting for the model's +0.12 °C offset, ensuring all threshold crossing dates (e.g. Greenland Ice Sheet, Coral Reefs) align with published empirical science (Armstrong McKay et al. 2022).

Transient Warming vs. Instant Ocean Equilibrium (Committed Warming)

The two-box ocean energy balance model distinguishes between two fundamental timescales of climate change:

Transient Climate Response (In-Flight Simulation): While greenhouse gas concentrations are actively changing, the upper ocean mixed layer warms relatively quickly (timescale τ_fast ≈ 3.5 years), but heat is constantly drawn downward into the massive intermediate and deep ocean (γ·(ΔT_s − ΔT_d)). This deep ocean heat uptake creates a persistent planetary energy imbalance (EEI > 0) that acts as a brake, delaying surface warming by decades to centuries.
Instant Ocean Equilibrium (t → ∞): The "⚡ Instant Ocean Equilibrium" feature calculates what happens when the entire ocean column has finally caught up to atmospheric CO₂ (ΔT_s = ΔT_d = ΔF / λ). Downward heat transport drops to zero, Top-of-Atmosphere net imbalance ceases (EEI = 0.00 W/m²), and the full "committed warming in the pipeline" (Equilibrium Climate Sensitivity, ECS = 3.00 °C per doubling) is revealed.

Climate State	Simulator Action	Physical Mechanism	Earth Energy Imbalance (EEI)	Response to 2× CO₂ (3.93 W/m²)
Transient Response (TCR)	Active dynamic stepping (`▶ Start Simulation`)	Ocean mixed layer warms while deep abyss absorbs ~90% of heat uptake via γ·(ΔT_s − ΔT_d).	EEI > 0 (~0.7 to 1.0+ W/m² today)	+1.96 °C (at year 70 under 1%/yr CO₂ increase)
Instant Ocean Equilibrium (ECS)	Clicking `⚡ Instant Ocean Equilibrium`	Deep ocean abyss catches up completely (ΔT_d = ΔT_s); downward ocean heat uptake drops to 0.	EEI = 0.00 W/m² (exact TOA radiative balance)	+3.00 °C (ΔT_eq = ΔF / λ = 3.93 / 1.310)

Page 4: 1-Million-Year Paleoclimate Sources & Resolution

The pop-out 1-Million-Year / 500-Year Rate-of-Change Comparator provides an illustrative reconstruction parameterized from published benchmark glacial–interglacial stages, terminations, and Holocene syntheses, rather than loading continuous raw observational datasets. Because no single physical archive spans 1,000,000 years at decadal resolution, paleoclimatology synthesizes polar ice cores with global marine sediment stacks.

1. Primary Scientific Syntheses & Benchmark Proxy Sources

Time Span (Years BP)	Archive & Dataset Type	Primary Reference	Temporal Resolution & Nature of Proxy
0 – 12,000 BP (Holocene)	Multi-proxy global surface temperature reconstruction (Temp12k / Marcott) & Law Dome ice cores	Kaufman et al. (2020), Scientific Data; Marcott et al. (2013), Science; Rubino et al. (2019)	10 – 50 years (multi-proxy surface temperature stack & high-resolution firn/ice CO₂)
11,000 – 22,000 BP (Termination I / Last Deglaciation)	Global multi-proxy deglaciation stack (80 globally distributed records) + WAIS Divide & NGRIP ice cores	Shakun et al. (2012), Nature; Tierney et al. (2020), Nature (LGM ΔT = -6.0 °C)	20 – 100 years (globally distributed marine/terrestrial surface proxies & polar ice)
0 – 800,000 BP (8 Glacial Cycles)	EPICA Dome C Antarctic Ice Core deuterium (δD) local temperature & composite trapped-gas atmospheric CO₂ record	Jouzel et al. (2007), Science; Bereiter et al. (2015), Geophysical Research Letters (NOAA/NCEI)	50 – 500 yrs. Records local East Antarctic site temperature (scaled to GMST) and direct continuous trapped-gas CO₂ back to 800,000 BP (the oldest continuous ice core CO₂).
800,000 – 1,000,000+ BP (Mid-Pleistocene Transition)	LR04 Global Benthic Stack (57 globally distributed deep-sea cores) & derived surface temperature reconstructions	Lisiecki & Raymo (2005), Paleoceanography; Snyder (2016), Nature; Hansen et al. (2013)	500 – 1,000 years. LR04 records benthic foraminiferal δ¹⁸O (reflecting deep-ocean temperature and continental ice volume, not direct surface temperature or CO₂); surface ΔT is indirectly derived, and pre-800 kyr CO₂ values are proxy/orbital estimates.

2. Polar Amplification & Global Mean Surface Temperature (GMST) Scaling

Antarctic ice-core deuterium isotopes (δD in Jouzel et al. 2007) record local East Antarctic site temperature swings of ~10 °C to 12 °C between deep glacial minimums and interglacial peaks. Due to polar amplification, polar temperature changes are roughly ~2× larger than the global average. Following Tierney et al. (2020), Shakun et al. (2012), and IPCC AR6 Chapter 2, Antarctic isotope anomalies are scaled by the empirical global-to-polar ratio (~0.52) so that the Last Glacial Maximum (21,000 BP) corresponds to a Global Mean Surface Temperature anomaly of -6.0 °C relative to 1750 CE, and warm interglacials (such as the Eemian / MIS 5e at 125,000 BP) peak at +1.0 °C to +1.4 °C.

3. Resolution Limits & Why the 500-Year Window Comparison Is Robust

• Fastest Prehistoric Global Warming (Termination I, ~15,000–14,500 BP): +0.35 °C per 500 yrs (+0.070 °C / century)
• Typical Holocene Natural Variability (6,000 BP): ±0.04 °C per 500 yrs (±0.008 °C / century)
• Observed & Projected Modern Trajectory (1958–2058 CE at +2.5 ppm/yr): +2.15 °C per 100 yrs (~31× faster than Termination I)

In deep ice cores (400,000 to 800,000 BP) and marine sediment cores (800,000 to 1,000,000 BP), molecular diffusion and core compression smooth out sub-centennial fluctuations. However, in the younger record (0 to 50,000 BP)—where WAIS Divide, NGRIP, and global deglacial stacks preserve 20-to-100-year resolution—we can directly resolve the fastest natural global warming events in Earth's recent history (Orbital Glacial Terminations). Even during the steepest 500-year phase of Termination I (15,000 to 14,500 BP), global temperature rose by ~0.35 °C over five centuries, driven by orbital Milankovitch insolation shifts and an ~8 ppm natural CO₂ outgassing from warming oceans. By contrast, human industrial emissions have added +149 ppm of CO₂ since 1750 (and over +110 ppm since 1958 alone), driving a warming slope approximately 31 times steeper than peak natural deglaciation.

4. Direct Links to Public Archival Repositories

NOAA National Centers for Environmental Information (NCEI) Paleoclimatology Ice Core Portal:
https://www.ncei.noaa.gov/products/paleoclimatology/ice-core ↗
Bereiter et al. (2015) / EPICA Dome C 800,000-Year Composite CO₂ Data:
https://www.ncei.noaa.gov/access/paleo-search/study/17975 ↗
Lisiecki & Raymo (2005) LR04 57-Core Marine Benthic Stack:
https://lorraine-lisiecki.com/stack.html ↗

Page 5: Fossil Fuel Metrics — Barrels of Oil & Gas to Atmospheric CO₂ (ppm) & Temperature Rise

While international climate negotiations and scientific syntheses frame planetary heating in abstract metrics—such as gigatonnes of carbon dioxide (Gt CO₂) or radiative forcing in watts per square metre (W/m²)—the industrial economy operates in physical commodities: barrels of crude oil, cubic feet of natural gas, and tonnes of coal.

This chapter details the exact physical, chemical, and atmospheric formulas that map physical units of extracted fossil fuels directly to atmospheric concentration in parts per million (ppm), top-of-atmosphere radiative forcing (W/m²), and eventual global surface warming (°C).

1. Combustion Stoichiometry: From Hydrocarbons to Mass of CO₂

When fossil fuels are burned for power, heat, or transport, atmospheric oxygen combines with the fuel's carbon atoms in exothermic oxidation reactions. Because every molecule of CO₂ contains one carbon atom (atomic mass 12.011 u) and two oxygen atoms (15.999 u each), the molecular mass of CO₂ is 44.010 u. Consequently, burning 1.0 kg of elemental carbon produces:

44.010 g/mol CO₂ / 12.011 g/mol C = 3.664 kg of CO₂ per 1.0 kg of Carbon burned

Applying this stoichiometric principle across the primary fossil commodities yields precise emission factors (calibrated against US EPA, US Energy Information Administration, and IPCC Guidelines for National Greenhouse Gas Inventories):

Fossil Commodity	Standard Physical Unit	Carbon Content by Mass	Direct CO₂ Emitted per Unit	Emissions per 1 Billion Units
Crude Oil (Petroleum)	1 barrel (bbl = 42 US gal = 158.99 L)	~84% – 86% Carbon (~117 kg C / bbl)	0.430 metric tonnes CO₂ (430 kg)	0.430 Gt CO₂ per billion barrels (Gbbl)
Fossil Natural Gas	1,000 cubic feet (1 Mcf at STP)	~75% Carbon (predominantly CH₄)	0.0545 metric tonnes CO₂ (54.5 kg)	0.0545 Gt CO₂ per Trillion Cubic Feet (Tcf)
Natural Gas (Energy Equiv.)	1 barrel of oil equivalent (boe ≈ 5,800 cu ft)	Equivalent thermal energy to 1 bbl crude	≈0.316 metric tonnes CO₂ (5,800 cu ft × 0.0545 t / 1,000 cu ft)	≈0.316 Gt CO₂ per billion boe
Thermal Coal (Bituminous)	1 metric tonne of hard coal (1,000 kg)	~65% – 75% Carbon (~660 kg C / t)	2.42 metric tonnes CO₂ for hard bituminous coal; ≈1.80 t CO₂ per tonne for the global production mix (incl. lower-carbon lignite)	≈1.8 – 2.4 Gt CO₂ per billion tonnes (this guide uses 1.80 for annual global totals)

2. The Mass of Earth's Atmosphere & The Universal ppm Factor

Atmospheric concentration of greenhouse gases is measured in parts per million by volume or mole fraction (ppmv). To calculate how much physical mass corresponds to 1.0 ppm, we compute the total mass of Earth's atmosphere:

Mean global surface atmospheric pressure: P₀ = 985.5 hPa (accounting for land topography above sea level)
Earth surface area: A = 4πR² = 5.101 × 10¹⁴ m²
Mean acceleration of gravity: g = 9.80665 m/s²
Total atmospheric mass: M_atm = (P₀ · A) / g = 5.148 × 10¹⁸ kg = 5.148 × 10⁶ Gt

With the mean molar mass of dry air at M_air = 28.964 g/mol and carbon dioxide at M_CO₂ = 44.010 g/mol, the mass ratio is 44.010 / 28.964 ≈ 1.5194. Evaluating 1 part in 10⁶ of total atmospheric moles:

1.0 ppm atmospheric CO₂ ≡ 7.824 Gt CO₂ (or 2.134 Gt of pure Carbon)

In a sealed, non-absorbing container the size of Earth's atmosphere, injecting 7.824 gigatonnes of CO₂ would raise concentration by exactly 1.0 ppm.

3. The Airborne Fraction (AF ≈ 46%) & Earth's Carbon Sinks

In the real Earth system, the atmosphere is not an isolated glass flask. When human industry emits CO₂, Earth's dynamic biological and geochemical sinks immediately absorb a significant portion of the gas:

The Ocean Carbon Sink (~26% absorbed): Atmospheric CO₂ dissolves across the air-sea boundary, reacting with water molecules to form carbonic acid, bicarbonate, and carbonate ions (H₂O + CO₂ ⇌ H₂CO₃ ⇌ HCO₃⁻ + H⁺). This geochemical dissolution drives progressive ocean acidification while buffering atmospheric accumulation.
The Terrestrial Biosphere Sink (~28% absorbed): Elevated ambient CO₂ concentrations stimulate photosynthetic carboxylation in land plants (CO₂ fertilization effect) and enhance soil organic matter storage, sequestering roughly a quarter of annual anthropogenic emissions in forests and soils.
The Atmospheric Retention (~46% remains airborne): Over the multi-decadal instrumental era (1958–2026), the empirical Airborne Fraction (AF) has remained remarkably steady at AF ≈ 0.46 ± 0.04.

Because only ~46% of every emitted tonne accumulates in the atmosphere, the effective emission mass required to produce a net 1.0 ppm rise is scaled by 1 / AF:

Effective Mass per ppm = 7.824 Gt CO₂ / 0.46 ≈ 17.01 Gt of gross CO₂ emitted per +1.0 ppm atmospheric increase

The 17-Gigaton Principle:
As a robust rule of thumb across contemporary climate science: Every ~17 gigatonnes of gross CO₂ humanity emits from fossil fuel combustion and industrial processes produces a net +1.0 ppm permanent rise on the atmospheric Keeling Curve.

4. Direct Equivalence: How Many Barrels of Oil & Gas Make 1.0 ppm?

With combustion stoichiometry and the airborne fraction established, we can directly solve for the physical commodity quantities that drive atmospheric concentration changes:

A. Crude Oil (Petroleum):

Each 1 billion barrels (1 Gbbl) of crude oil produces 0.430 Gt CO₂ upon combustion. Factoring in the 46% airborne fraction:

ΔC per billion barrels = (0.430 Gt × 0.46) / 7.824 Gt/ppm = +0.02528 ppm per billion barrels of oil
Barrels required for +1.0 ppm = 1.0 / 0.02528 ≈ 39.55 billion barrels of crude oil

Global Reality: Global oil consumption currently stands at ~102 million barrels per day (mbpd), or ~37.2 billion barrels per year. Thus, burning oil alone adds:

Annual oil contribution = 37.2 Gbbl/yr × 0.02528 ppm/Gbbl ≈ +0.94 ppm / year from crude oil alone

B. Fossil Natural Gas:

Each 1 trillion cubic feet (1 Tcf) of natural gas emits 0.0545 Gt CO₂. With the 46% airborne fraction:

ΔC per Tcf = (0.0545 Gt × 0.46) / 7.824 Gt/ppm = +0.003204 ppm per Trillion Cubic Feet
Gas volume required for +1.0 ppm = 1.0 / 0.003204 ≈ 312.1 Tcf of natural gas (or 72.4 billion boe)

Global Reality: World natural gas consumption is currently ~140 Tcf per year (~4,000 billion m³), emitting ~7.63 Gt CO₂/year. Natural gas combustion adds:

Annual gas contribution = 140 Tcf/yr × 0.003204 ppm/Tcf ≈ +0.45 ppm / year from natural gas alone

C. Synthesis of Global Atmospheric Growth (≈ +2.6 ppm/yr):

Summing across all human fossil extraction and land-use activities gives a growth rate in the same range as the observed Keeling Curve rise at NOAA's Mauna Loa Observatory (≈2.4 ppm/yr decadal mean; individual years vary with El Niño and sink strength):

Emissions Source	Annual Extraction / Activity	Gross Emissions (Gt CO₂/yr)	Net Airborne CO₂ (Gt/yr)	Annual Growth Added	Share
Crude Oil Combustion	~37.2 billion barrels / yr (~102 mbpd)	~16.00 Gt	~7.36 Gt	+0.94 ppm / yr	35.9%
Coal Combustion	~8.5 billion metric tonnes / yr	~15.30 Gt	~7.04 Gt	+0.90 ppm / yr	34.4%
Fossil Natural Gas	~140 trillion cu ft / yr (~4,000 bcm)	~7.63 Gt	~3.51 Gt	+0.45 ppm / yr	17.2%
Cement & Industrial Flaring	Clinker production & upstream flaring	~1.80 Gt	~0.83 Gt	+0.11 ppm / yr	4.2%
Land-Use Change & Deforestation	Tropical clearing & peat degradation	~3.80 Gt	~1.75 Gt	+0.22 ppm / yr	8.4%
TOTAL ANTHROPOGENIC FORCING	All Global Industrial Activity	~44.5 Gt (gross, summed from rows above)	~20.5 Gt	≈ +2.6 ppm / year (decadal observed mean ≈ 2.4; 2023–24 ≈ 3.4–3.8 as sinks weakened)	100.0%

5. The Final Temperature Link: From Barrels of Fuel to °C of Warming

To connect physical fuel extraction to actual surface warming, we evaluate the simulator's core radiative forcing formula (Myhre et al. 1998; IPCC AR6 Chapter 7):

ΔF = α · ln(C / C₀) = 5.67 · ln(C / 278) W/m²

Taking the marginal derivative with respect to CO₂ concentration (dF / dC) evaluated at present-day background levels (C ≈ 427 ppm):

dF / dC = 5.67 / 427 ≈ 0.01328 W/m² of radiative forcing per 1.0 ppm

Under the IPCC AR6 calibrated net climate feedback parameter λ = 1.31 W/m²/K (corresponding to Equilibrium Climate Sensitivity ECS = 3.00 °C per CO₂ doubling), the eventual equilibrium warming committed by a 1.0 ppm increase is:

d(ΔT_eq) = (dF / λ) = 0.01328 W/m² / 1.310 W/m²/K ≈ 0.01014 °C of equilibrium warming per 1.0 ppm

The "One-Hundredth of a Degree" Golden Rule:
At current atmospheric concentrations, every +1.0 ppm increase in atmospheric CO₂ commits the planet to about 0.010 °C (1/100th of a degree Celsius) of long-term equilibrium warming.
• An increase of +10 ppm commits Earth to +0.10 °C of eventual heating.
• An increase of +100 ppm commits Earth to +1.00 °C of eventual heating.

6. Master Fuel-to-Warming Conversion Matrix

The table below synthesizes the complete chain: physical fuel volume → gross emissions → airborne ppm added → top-of-atmosphere radiative forcing → transient surface warming (realized within decades) → long-term committed equilibrium warming:

Fuel Volume / Quantity	Gross CO₂ Emitted	Atmospheric ppm Added	Incremental Forcing (ΔF)	Transient Warming (Decadal)	Committed Equilibrium Warming (ΔT_eq)
1 Million Barrels of Oil (1 mmbbl)	0.00043 Gt	+0.0000253 ppm	+0.00000034 W/m²	+0.00000019 °C	+0.00000026 °C
100 Million Barrels of Oil (~1 Day Global Oil Burn)	0.0430 Gt	+0.00253 ppm	+0.0000336 W/m²	+0.0000195 °C	+0.0000256 °C
1 Billion Barrels of Oil (1 Gbbl)	0.430 Gt	+0.0253 ppm	+0.000336 W/m²	+0.000195 °C	+0.000256 °C
3.95 Billion Barrels of Oil	1.70 Gt	+0.100 ppm	+0.00133 W/m²	+0.00077 °C	+0.00101 °C (1/1000th °C)
37.2 Billion Barrels of Oil (~1 Year Global Oil Burn)	16.00 Gt	+0.941 ppm	+0.0125 W/m²	+0.0073 °C	+0.00954 °C (~0.01 °C / yr)
1 Trillion Cubic Feet of Gas (1 Tcf)	0.0545 Gt	+0.00320 ppm	+0.0000425 W/m²	+0.0000247 °C	+0.0000324 °C
140 Trillion Cubic Feet of Gas (~1 Year Global Gas Burn)	7.63 Gt	+0.449 ppm	+0.00596 W/m²	+0.0035 °C	+0.00455 °C
1 Billion Tonnes of Hard Coal (1 Gt coal)	2.42 Gt	+0.142 ppm	+0.00189 W/m²	+0.00110 °C	+0.00144 °C
8.5 Billion Tonnes of Coal (~1 Year Global Coal Burn)	15.30 Gt	+0.899 ppm	+0.0119 W/m²	+0.0069 °C	+0.00911 °C

Page 6: Global Carbon Budgets & The "Unburnable Carbon" Reality

Because carbon dioxide persists in the planetary atmosphere and surface ocean for millennia, global warming is not governed by the rate of annual emissions in any single year, but by the cumulative total mass of carbon ever emitted since the pre-industrial era. This immutable physical reality gives rise to the scientific framework of the Planetary Carbon Budget.

1. The Physical Foundation: Why Temperature Scales Linearly with Cumulative Emissions (TCRE)

In complex Earth System Models and observational records alike, global mean surface temperature rise exhibits a near-perfect linear proportionality to cumulative CO₂ emissions. This linearity is the result of a remarkable physical compensation between two opposing non-linear geophysical processes:

Logarithmic Radiative Saturation (Decelerating Heating): As CO₂ concentration rises, greenhouse absorption bands become increasingly opaque, so each additional ton of CO₂ exerts slightly less marginal radiative forcing (ΔF ∝ ln(C/C₀)).
Ocean Heat & Carbon Uptake Saturation (Accelerating Heating): As the upper ocean warms and acidifies, the seawater Revelle chemical buffer factor increases, reducing the ocean's ability to absorb further carbon. Simultaneously, upper-ocean stratification reduces downward heat diffusion into the deep abyss, leaving a larger fraction of trapped energy at the surface.

Over century timescales, the decelerating radiative effect is almost exactly counterbalanced by the accelerating ocean saturation effect. This relationship is quantified by the Transient Climate Response to Cumulative Carbon Emissions (TCRE), established in IPCC AR6 Working Group I (Chapter 5):

TCRE ≈ 0.45 °C of surface warming per 1,000 Gt CO₂ (likely range: 0.27 °C to 0.63 °C per 1,000 Gt CO₂)
Equivalently: 1.0 °C of global warming requires ~2,220 Gt of cumulative CO₂ emissions

2. Remaining Global Carbon Budgets for Paris Agreement Limits

From the onset of the Industrial Revolution (1750) through the end of 2024, human activity emitted approximately 2,600 Gt of CO₂ (roughly 1,800 Gt from fossil fuel combustion and cement manufacturing, plus ~800 Gt from deforestation and land-use change). This historical burn has pushed forced surface temperatures to +1.26 °C above 1750 (+1.14 °C above 1850–1900).

With global emissions continuing at approximately 41.0 Gt CO₂ per year (37.4 Gt fossil + 3.6 Gt land use), the remaining allowable emissions quotas (carbon budgets) from January 2025/2026 forward are vanishingly small (synthesized from IPCC AR6 WGI Chapter 5, Forster et al. 2024 Earth System Science Data, and the Global Carbon Project 2024):

Warming Target	Probability of Staying Below	Remaining Carbon Budget (from 2025/2026)	Time Remaining at Current Burn Rate (~41 Gt/yr)	Equivalent Oil Combustion
+1.5 °C (Paris Target)	50% (Even Odds)	~200 Gt CO₂	~4.9 years (Depleted ~2030)	~465 billion barrels
+1.5 °C (Paris Target)	67% (Two-Thirds Chance)	~130 Gt CO₂	~3.2 years (Depleted ~2028)	~302 billion barrels
+1.7 °C (Intermediate)	50% (Even Odds)	~650 Gt CO₂	~15.9 years (Depleted ~2041)	~1,510 billion barrels
+2.0 °C (Upper Guardrail)	50% (Even Odds)	~1,150 Gt CO₂	~28.0 years (Depleted ~2053)	~2,674 billion barrels
+2.0 °C (Upper Guardrail)	67% (Two-Thirds Chance)	~950 Gt CO₂	~23.2 years (Depleted ~2048)	~2,209 billion barrels
+2.0 °C (Upper Guardrail)	83% (High Confidence)	~750 Gt CO₂	~18.3 years (Depleted ~2043)	~1,744 billion barrels

3. The "Unburnable Carbon" Paradox: Known Reserves vs. Planetary Sink

The central tension of modern climate economics lies in the severe arithmetic mismatch between the physical carbon budget of the atmosphere and the proven fossil fuel reserves currently claimed as commercial assets on the balance sheets of energy corporations and sovereign nation-states.

According to the BP Statistical Review of World Energy, the International Energy Agency (IEA World Energy Outlook), and IPCC AR6 Working Group III (Chapter 2), total global proven extractable fossil fuel reserves consist of:

Fossil Resource	Global Proven Reserves	Emission Factor	Potential CO₂ Embedded in Reserves
Proven Crude Oil Reserves	~1,750 billion barrels (Gbbl)	0.430 Gt CO₂ / Gbbl	~752 Gt CO₂
Proven Natural Gas Reserves	~7,300 trillion cubic feet (Tcf) (~207 trillion m³)	0.0545 Gt CO₂ / Tcf	~398 Gt CO₂
Proven Thermal Coal Reserves	~1,074 billion metric tonnes (Gt coal)	2.00 – 2.42 Gt CO₂ / Gt	~2,250 Gt CO₂
TOTAL PROVEN RESERVES	Economically extractable with existing tech	Combined fossil fuels	~3,400 Gt CO₂

The Mathematical Paradox of Unburnable Carbon:
• The remaining budget to hold warming to +1.5 °C (50% probability) is ~200 Gt CO₂.
• Total proven fossil reserves hold ~3,400 Gt CO₂.
• Conclusion: To achieve the +1.5 °C Paris target, over 94% of proven coal, oil, and gas reserves must remain permanently unextracted in the ground.
• Even for the upper +2.0 °C guardrail (budget ~1,150 Gt CO₂), roughly two-thirds (~66%) of all existing proven fossil fuel reserves are completely "unburnable".

4. What Happens if All Proven Reserves Are Combusted?

If human industry extracts and burns the ~3,400 Gt CO₂ of known fossil reserves, we can compute the resulting climate state directly using the physical equations in this simulator:

Emissions pulse: ΔE = 3,400 Gt CO₂
Net airborne retention (AF = 0.46): ΔC = (3,400 × 0.46) / 7.824 ≈ +200 ppm of CO₂
New atmospheric concentration: C = 427 ppm + 200 ppm = 627 ppm
Top-of-Atmosphere Radiative Forcing: ΔF = 5.67 · ln(627 / 278) = +4.61 W/m²
Committed Equilibrium Warming: ΔT_eq = 4.61 W/m² / 1.31 W/m²/K = +3.52 °C above 1750 (+3.40 °C above 1850–1900)

Warming of +3.5 °C to +4.0 °C would cross multiple catastrophic planetary tipping points identified by Armstrong McKay et al. (2022) in Science: the total collapse of the Greenland and West Antarctic Ice Sheets (committing 10+ metres of sea-level rise), the complete dieback of the Amazon rainforest, collapse of low-latitude coral reef ecosystems, and possible cessation of the Atlantic Meridional Overturning Circulation (AMOC).

5. Net Zero, The Zero Emissions Commitment (ZEC), & Drawdown Rebound

A common misconception in climate discussions is: "If emissions drop to zero, will temperatures continue to rise for centuries?"

The physical answer, verified by the IPCC AR6 and the international Zero Emissions Commitment Model Intercomparison Project (ZECMIP, MacDougall et al. 2020), is no:

The Zero Emissions Commitment (ZEC ≈ 0.0 °C): When anthropogenic CO₂ emissions cease entirely (Net Zero), ocean heat uptake continues absorbing excess energy from the lower atmosphere. At the same time, the ocean and terrestrial biosphere continue drawing down atmospheric CO₂, causing atmospheric radiative forcing to decline. In physical climate models, these two opposing curves almost exactly offset: global surface warming essentially halts within a few years of reaching true zero emissions.
The Carbon Drawdown "Sink Rebound" Effect: While stopping emissions halts further warming, restoring atmospheric CO₂ back down to pre-industrial (278 ppm) or safe historic levels (350 ppm) requires active Carbon Dioxide Removal (CDR).
However, as atmospheric CO₂ partial pressure drops, the oceans and land biosphere begin outgassing a fraction of the carbon they previously absorbed. Consequently, removing 1.0 ppm from the atmosphere requires scrubbing not just 7.82 Gt CO₂, but approximately ~17 Gt of CO₂ from the active Earth system to overcome the outgassing rebound.

6. Key References & Authoritative Scientific Sources

IPCC AR6 Working Group I, Chapter 5 (Global Carbon and Other Biogeochemical Cycles and Feedbacks):
https://www.ipcc.ch/report/ar6/wg1/chapter/chapter-5/ ↗
Forster et al. (2024), Earth System Science Data — Indicators of Global Climate Change 2023: Annual update of key indicators of the state of the climate system and human influence:
https://essd.copernicus.org/articles/16/2625/2024/ ↗
Global Carbon Project (2024) — Global Carbon Budget:
https://globalcarbonbudget.org/ ↗
MacDougall et al. (2020), Biogeosciences — Is there warming in the pipeline? A multi-model analysis of the Zero Emissions Commitment from CO₂:
https://bg.copernicus.org/articles/17/2987/2020/ ↗
McGlade & Ekins (2015), Nature — The geographical distribution of fossil fuels unused when limiting global warming to 2 °C:
https://www.nature.com/articles/nature14016 ↗

Page 7: Global Emissions Anatomy (Our World in Data & Climate Watch)

To eliminate the ~41 Gt/yr of global CO₂ emissions driving planetary energy imbalance, we must understand precisely where greenhouse gas emissions originate across the global economy. This synthesis integrates data from the Global Carbon Project (Friedlingstein et al. 2023) and Climate Watch / World Resources Institute (2023), as standardized by Our World in Data (Ritchie, Roser & Rosado 2020).

1. The Distinction Between Total Greenhouse Gases and Carbon Dioxide

Climate discussions frequently conflate total greenhouse gas emissions (measured in CO₂-equivalents, Gt CO₂e) with pure carbon dioxide emissions (Gt CO₂):

• Total Anthropogenic Greenhouse Gases: ~54.6 Gt CO₂e / yr (100%)
• Carbon Dioxide (CO₂): ~40.9 Gt CO₂ / yr (74.9% of total GWP₁₀₀)
— Fossil Fuel Combustion & Industry: 37.4 Gt CO₂ / yr
— Land Use Change & Deforestation: ~3.5 Gt CO₂ / yr
• Methane (CH₄, GWP₁₀₀ = 28): ~8.8 Gt CO₂e / yr (16.1% — fugitive gas, cattle, rice)
• Nitrous Oxide (N₂O, GWP₁₀₀ = 265): ~3.3 Gt CO₂e / yr (6.0% — agricultural soils, fertilizers)
• Fluorinated Gases (HFCs, PFCs, SF₆): ~1.6 Gt CO₂e / yr (2.9% — refrigeration, switchgear)

2. The Four Macro-Economic Sectors (Our World in Data)

Sector	Share of GHG	Annual Volume	Key Activities & Emission Sources	Primary Clean Abatement Pathways
1. Energy	73.2%	40.0 Gt CO₂e	Industry Energy (24.2%): Steel blast furnaces, chemical heat, paper/pulp. Energy in Buildings (17.5%): Residential & commercial space heating, boilers, cooking. Transport (16.2%): Road cars & trucks (11.9%), aviation (1.9%), shipping (1.7%). Fugitive Energy (5.8%): Oil/gas venting & flaring (3.9%), coal mine methane (1.9%). Other Combustion (9.5%): Tractors, fishing, unallocated.	Clean electricity generation (Solar, Wind, Nuclear, Hydro); grid storage; heat pumps for buildings; battery electric vehicles (BEVs) for road; green hydrogen for steelmaking (H₂-DRI).
2. Agriculture, Forestry & Land Use (AFOLU)	18.4%	10.0 Gt CO₂e	Livestock & Manure (5.8%): Methane from ruminant enteric fermentation. Agricultural Soils (4.1%): Nitrous oxide from synthetic nitrogen fertilizers. Deforestation (2.2%): Tropical forest clearing for pastures & soy. Crop Residues & Savanna (3.5%): Biomass burning. Rice Paddies (1.3%): Anaerobic soil methanogenesis.	Halting tropical deforestation; precision fertilizer application; dietary shifts toward plant-rich protein; alternate wetting & drying (AWD) in rice paddies; methane-inhibiting cattle feed additives.
3. Direct Industrial Processes	5.2%	2.8 Gt CO₂e	Cement Calcination (3.0%): Chemical limestone decarbonation (CaCO₃ → CaO + CO₂), separate from kiln fuel. Chemical Processes (2.2%): Non-combustion emissions from ammonia synthesis, ethylene steam cracking, solvents.	Limestone clinker substitution with calcined clay and pozzolans (LC3 cement); alternative geopolymer binders; point-source Carbon Capture & Storage (CCUS).
4. Waste	3.2%	1.7 Gt CO₂e	Landfills (1.9%): Anaerobic digestion of organic municipal waste. Wastewater (1.3%): Industrial & domestic sewage treatment methane.	Landfill methane gas capture for power generation; municipal organic composting; anaerobic digestion with biomethane upgrade; modern aerobic sewage aeration.

3. The 10:1 Trophic Calorie Inefficiency & Carbon Opportunity Cost of Livestock

While biogenic methane (CH₄) from enteric fermentation is an acute short-term warming driver, an even larger systemic climate impact of livestock agriculture is its thermodynamic calorie loss and disproportionate land-use footprint:

• Thermodynamic Calorie Loss (Lindeman's Trophic Rule): Animals lose ~90% of ingested plant calories to basal metabolism, body heat, and inedible tissue. Converting plant feed into animal meat requires 10 to 25 kg of plant feed per 1 kg of beef (~3 to 4 kg for pork and poultry).
• The 80% / 18% Land Disparity: According to Poore & Nemecek (2018, Science), livestock rearing (grazing land plus animal feed croplands) monopolizes ~80% of all agricultural land on Earth (~40 million km²), yet delivers only ~18% of global human calories and 37% of global protein. The remaining 82% of human calories are grown on just ~20% of farmland.
• 75% Land Reduction Potential: If human calories were consumed directly from plants, global agricultural land demand would shrink by up to 75%, freeing approximately 30 million km² (an area the size of the entire African continent).
• Carbon Opportunity Cost: Monopolizing 40 million km² in pasture and feed crops actively prevents natural forest and grassland regeneration. Allowing this land to naturally rewild would sequester an estimated 300 to 500+ Gt CO₂ in vegetative and soil biomass over several decades (Searchinger et al. 2018, Nature; Hayek et al. 2021, Nature Sustainability)—equivalent to nearly a decade of all global fossil fuel emissions combined.

4. Fossil Fuel CO₂ by Fuel Type

Looking exclusively at the combustion of fossil fuels (~37.4 Gt CO₂):

Coal (41.0% · 15.3 Gt CO₂/yr): Emits the highest carbon per unit of energy (~94 kg CO₂/GJ). Predominantly consumed in power plants (~65%) and metallurgical coke for ironmaking (~30%).
Crude Oil (32.6% · 12.2 Gt CO₂/yr): Emits ~73 kg CO₂/GJ. Dominates transport fuels (gasoline, diesel, kerosene jet A-1, heavy bunker fuel) with ~15% for petrochemical feedstocks.
Natural Gas (21.1% · 7.9 Gt CO₂/yr): Emits ~56 kg CO₂/GJ. Powers combined-cycle gas turbines (CCGT), residential/commercial building heating, and industrial steam.
Cement Process & Flaring (5.3% · 2.0 Gt CO₂/yr): Direct limestone chemical calcination and upstream oilfield safety flaring.

Page 8: The Clean Energy & Storage Transition: Replacing 41 Gt/yr Fossil Energy

Eliminating ~41 Gt/yr of global fossil CO₂ emissions requires deploying a clean energy system capable of delivering reliable power, heat, transport, and industrial feedstocks 24/7/365.

1. The Electrification Efficiency Dividend

A common misconception is that the world must replace all 135,000 TWh of fossil primary energy with 135,000 TWh of clean power. This overlooks the massive thermodynamic waste inherent to combustion engines and thermal power generation:

• Internal Combustion Engine (ICE) Vehicles: Convert only 20% to 30% of fuel energy into forward motion; ~75% is rejected as heat through radiators and exhaust pipes.
• Battery Electric Vehicles (BEVs): Convert 85% to 90% of electrical energy into wheel motion (~3× to 4× less energy needed per kilometre).
• Gas & Oil Boilers: Deliver 80% to 90% thermal efficiency.
• Electric Heat Pumps: Draw ambient environmental heat with a Coefficient of Performance (COP) of 3.0 to 4.0 (delivering 300% to 400% heat per unit of power consumed, cutting energy demand by ~65% to 75%).
• Thermal Power Plants (Coal/Gas): Reject 50% to 65% of chemical energy into cooling towers and rivers.

The Bottom Line: Because electrified devices waste far less energy, fully replacing 135,000 TWh of fossil primary energy requires only ~70,000 to ~80,000 TWh/yr of clean useful electricity (we adopt 75,000 TWh/yr as the central planetary benchmark).

2. Clean Generation Technologies, Costs & Capacities

Solar Photovoltaics (PV): Current global capacity ~1,600 GW (output ~1,600 TWh/yr). Levelized Cost of Electricity (LCOE) is $0.03 to $0.05 / kWh (down >85% since 2010; the cheapest source of new bulk energy in history). Capacity factor averages 18% to 22%, necessitating daytime solar storage firming.
Wind Power (Onshore & Offshore): Current global capacity ~1,020 GW (output ~2,300 TWh/yr). LCOE is $0.035 to $0.05 / kWh onshore and $0.07 to $0.09 / kWh offshore. Capacity factors reach 35% onshore and 45%–55% for offshore turbines. Produces higher output in winter and night, offering strong seasonal complementarity to solar PV.
Nuclear Fission (Baseload & SMRs): Current global capacity ~375 GW across 413 reactors (output ~2,650 TWh/yr). LCOE is $0.07 to $0.15 / kWh for new builds (amortized existing fleet $0.03 to $0.04 / kWh). Operates at a ~90% capacity factor, providing 24/7 firm, weather-independent clean power with the lowest land-use footprint of any clean energy source (~0.1 m²/MWh).
Hydroelectric (Conventional & Run-of-River): Current global capacity ~1,250 GW (output ~4,300 TWh/yr). LCOE is $0.04 to $0.08 / kWh. Operates at ~42% capacity factor, providing essential dispatchable clean inertia, synchronous grid stabilization, and seasonal reservoir storage.
Geothermal, Biomass & Marine: Current global capacity ~175 GW (Biomass ~145 GW, Geothermal ~16 GW, Marine ~0.5 GW; output ~800 TWh/yr). LCOE is $0.05 to $0.08 / kWh for geothermal and $0.07 to $0.12 / kWh for biomass. Operates at ~70% to 85% capacity factor, delivering continuous high-temperature industrial steam and dispatchable firm power.

3. Physics, Costs & Mechanics of the Six Grid Storage Technologies

Because solar and wind generation vary with weather and daylight, energy storage and grid firming are required to prevent outages (such as during multi-day winter wind/solar lulls, or Dunkelflaute):

Technology	Current Capacity	Levelized Cost (LCOS)	Round-Trip Efficiency (RTE)	Discharge Duration	Optimal Grid Application & Physics
1. Chemical Batteries (Li-ion, Na-ion, LFP, Flow)	~0.8 TWh (~800 GWh / ~150 GW)	$0.10 – $0.18 / kWh-cycle	85% – 92%	2 – 8 hours	Sub-second frequency regulation; evening solar peak shaving; distribution substation firming. Moderate degradation over 3,000–8,000 cycles. (ΔG = -nFE)
2. Pumped Hydro Storage (PHS)	~9.0 TWh (>90% of global storage)	$0.05 – $0.12 / kWh-cycle	75% – 82%	8 – 24+ hours	Bulk daily and multi-day dispatch. Lowest levelized storage cost for massive volume. Zero capacity degradation over 50–100 years. (E = ρVgh)
3. Compressed Air (CAES)	~0.2 TWh (~200 GWh / ~2 GW)	$0.07 – $0.14 / kWh-cycle	55% – 70% (adiabatic thermal recovery)	6 – 24 hours	Large-scale intraday bulk storage in salt caverns without rare battery metals. High cycling endurance with deep depth of discharge. (W = ∫P\,dV)
4. Gravity Energy Storage	~0.05 TWh (~50 GWh commercial)	$0.08 – $0.15 / kWh-cycle	75% – 85%	4 – 12 hours	Medium-duration storage utilizing abandoned mine shafts or vertical towers. Zero standby self-discharge; 40+ year mechanical design life. (E = mgh)
5. Flywheels	~0.02 TWh (~20 GWh / ~1.5 GW)	$0.15 – $0.25 / kWh-cycle	85% – 90%	Seconds to 1 hour	Instantaneous sub-cycle grid frequency regulation (<10 ms); synchronous inertia to stabilize transmission grids lacking spinning steam turbines; >100,000 cycle lifespan. (E = ½Iω²)
6. Hydrogen & Power-to-Gas	~0.5 TWh (~500 GWh in salt caverns)	$0.18 – $0.35 / kWh-cycle	35% – 48% (power-to-gas-to-power)	Weeks to Months	Multi-week seasonal energy storage; storing summer solar surplus in salt caverns for winter heating; zero-carbon fuel for green steel and long-haul transport. (2H₂O → 2H₂ + O₂)

4. Transition Pathways & Grid Balancing Trade-Offs

The Clean Energy Solutions view models five characteristic global transition pathways illustrating the trade-offs between generation diversity, land footprint, and required energy storage:

Current 2026 Fleet: Reflects today's installed operational baseline (~11,650 TWh/yr clean power: hydro 4,300, nuclear 2,650, wind 2,300, solar 1,600, bio/geo 800 TWh; backed by ~10.6 TWh of storage, primarily pumped hydro).
IEA Net Zero 2050: The international consensus pathway aligning with Paris Agreement targets: aggressive wind (18,000 TWh) and solar (24,000 TWh) deployment supported by expanded nuclear (5,500 TWh), hydro (7,000 TWh), geothermal/biomass (3,500 TWh), and ~40 TWh of diversified storage.
100% Renewables: Relying purely on wind (26,000 TWh) and solar (38,000 TWh) alongside existing hydro (6,500 TWh). Because wind and solar output fluctuate with daylight and multi-week weather systems (Dunkelflaute), this scenario requires massive overbuilding of generation and ~91 TWh of long-duration and seasonal grid storage (hydrogen, pumped hydro, CAES) to eliminate outage risks.
High Nuclear Baseload: Heavy buildout of firm, weather-independent nuclear fission (20,000 TWh / ~2,500 GW) paired with moderate renewables. Because continuous 24/7 firm baseload matches the bulk of baseline industrial and residential demand, required grid storage drops to ~30 TWh and land surface footprint is minimized.
Balanced (Solar/Wind + Baseload): Combines ultra-low-cost variable wind (16,000 TWh) and solar (22,000 TWh) for bulk energy with dispatchable 24/7 firm clean baseload (10,000 TWh nuclear, 6,000 TWh hydro, 3,000 TWh geothermal/biomass). By supplying the "last 15% to 20%" of inflexible demand from firm baseload, the system cuts required long-duration storage by ~50% compared to a pure variable renewable architecture while avoiding over-reliance on any single technology.

Addendum: Present-Day Forced Warming (+1.26 °C) vs. Headlines Reporting +1.5 °C Breached

A common point of inquiry when reviewing climate models is: "Why does this physical simulator model a warming anomaly of approximately +1.26 °C for 2026, when international headlines in 2023–2024 widely reported that global temperatures had already breached +1.5 °C?"

This apparent divergence is grounded in the critical distinction between single-year weather variability (driven by chaotic internal ocean oscillations like El Niño) and the underlying multi-decadal forced climate response formally defined by the IPCC and the United Nations Paris Agreement.

1. Single-Year Anomaly Spikes vs. The Multi-Decadal Forced Trend

During calendar years 2023 and 2024, observational monitoring agencies (Copernicus / ECMWF ERA5, Berkeley Earth, NOAA GlobalTemp, NASA GISTEMP) reported that global surface air temperatures temporarily reached +1.48 °C to +1.55 °C above the 1850–1900 reference level.

However, these single-year records were amplified by a temporary confluence of natural internal variability stacked on top of human global warming:

A Powerful El Niño (ENSO): El Niño events temporarily vent vast quantities of ocean heat from the tropical Pacific into the troposphere, spiking global mean surface temperatures by +0.10 °C to +0.25 °C above the forced baseline. As El Niño decays into neutral or La Niña conditions, global temperatures naturally drop back toward the forced trend line.
Solar Cycle 25 Maximum: The 11-year solar cycle reached its peak in 2024–2025, temporarily adding ~0.05 W/m² of top-of-atmosphere solar irradiance.
Marine Aerosol Reductions & Volcanic Forcing: New international shipping fuel rules (IMO 2020) reduced sulfur emissions, temporarily thinning reflective marine stratocumulus clouds over shipping corridors, while the massive January 2022 eruption of the Hunga Tonga–Hunga Haʻapai submarine volcano injected unprecedented water vapor directly into the stratosphere.

2. The Official UNFCCC & IPCC Definition of "Breaching 1.5 °C"

Under the Paris Agreement, the World Meteorological Organization (WMO), and IPCC AR6, a long-term climate threshold is not considered "breached" by a single unusually warm year or month. Climate is scientifically defined as the 20-year to 30-year running average of global temperatures, which filters out short-term ENSO noise:

• Single-Year Observational Spikes (2023–2024 El Niño Peak): +1.48 °C to +1.55 °C (relative to 1850–1900; temporary weather record)
• Human-Induced Forced Warming (Forster et al. 2024, ESSD): +1.31 °C [1.1 to 1.7 °C] in 2023, ~1.36 °C in 2024 (relative to 1850–1900)
• Simulator Model State (2026 Forced Anomaly): +1.26 °C rel. to 1750 (+1.14 °C rel. to 1850–1900)
• Projected Permanent Crossing of the 20-Year Mean (+1.5 °C): Early-to-mid 2030s (under continued emissions)

The IPCC Sixth Assessment Report (AR6 Working Group I) and the authoritative Indicators of Global Climate Change 2023 (Forster et al. 2024, ESSD) define human-induced warming as the multi-decadal forced response to greenhouse gases and aerosols relative to the 1850–1900 reference period, estimated at +1.31 °C [1.1 to 1.7 °C] in 2023 and reaching approximately +1.36 °C in 2024. Relative to the unforced pre-industrial 1750 baseline (which includes ~0.12 °C of early warming between 1750 and 1900), the simulator's deterministic forced response reaches +1.26 °C in 2026 (equivalent to +1.14 °C above the 1850–1900 mean).

3. Pre-Industrial Baseline Differences: 1750 vs. 1850–1900

Another source of numerical variation in the literature is the choice of pre-industrial baseline:

The 1850–1900 Meteorological Baseline: Observational weather datasets (HadCRUT5, Berkeley Earth) use 1850–1900 because that was the earliest era with standardized, globally distributed thermometer networks.
The 1750 Pre-Industrial Baseline: In physical radiative transfer models (and Antarctic ice cores), the true pre-industrial baseline is 1750 CE, when atmospheric CO₂ was in steady-state equilibrium at 278 ppm prior to industrial coal combustion.

Between 1750 and 1880, early industrial coal burning and land-clearing had already elevated CO₂ from 278 ppm to ~290 ppm, producing approximately +0.12 °C of early warming before global thermometer networks were established. Consequently, tipping-point thresholds (Armstrong McKay et al. 2022) and Paris Agreement targets are referenced to 1850–1900, while physical radiative balance simulations are referenced to the 1750 unforced baseline. The simulator accounts for this by computing the model's 1850–1900 hindcast offset (mean +0.12 °C) and adjusting threshold crossings accordingly.

4. Why the Two-Layer Ocean EBM Models the Forced Signal

This simulator models the foundational physics of planetary energy balance and the Two-Layer Ocean Energy Balance Model (Held et al. 2010; Geoffroy et al. 2013; IPCC AR6 Chapter 7). Because it computes the deterministic thermal response to greenhouse gas forcing without synthesizing stochastic weather noise (like individual El Niño cycles or volcanic eruptions), it calculates the pure forced thermodynamic climate response.

At 427 ppm CO₂ in 2026, the model computes an Earth Energy Imbalance of ≈ +0.78 W/m² being absorbed by the ocean thermal sink (observational satellite & Argo ocean heat content estimates typically report ~0.8 to 1.0+ W/m² due to internal weather variability and aerosol changes), and the model's transient surface temperature anomaly lands precisely at +1.26 °C above the 1750 baseline (+1.14 °C above 1850–1900).

5. Key References & Authoritative Reading

Copernicus Climate Change Service (C3S) — How close are we to reaching a global warming of 1.5°C? (Global Temperature Trend Monitor):
https://climate.copernicus.eu/how-close-are-we-reaching-global-warming-15degc ↗
IPCC AR6 Working Group I, Chapter 7 (The Earth's Energy Budget & Climate Sensitivity):
https://www.ipcc.ch/report/ar6/wg1/chapter/chapter-7/ ↗
Forster et al. (2024), Earth System Science Data (Indicators of Global Climate Change 2023):
https://essd.copernicus.org/articles/16/2625/2024/ ↗