Стаття
Клімат і Земля · ⏱ ≈ 13 хв читання · Останнє оновлення: 23 червня 2026 р.

Енергетичний баланс Землі — фізика клімату

Клімат — це не метеорологія. Якщо погода — це хаотичний щоденний стан атмосфери, то клімат — це довгостроковий енергетичний бюджет планети, баланс між поглинутим сонячним світлом та інфрачервоним теплом, випроміненим у космос. Починаючи з випромінювання Стефана-Больцмана та нульвимірної моделі енергетичного балансу, ми виводимо, чому Земля на 33 K тепліша, ніж «мала б» бути, як парникові гази зміщують баланс і звідки беруться точки незворотності.

Автор: Команда MySimulator · Редакційна перевірка: Редакція MySimulator

1. Сонячне випромінювання та альбедо

Світність Сонця становить L☉ = 3,846 × 10²⁶ Вт. На середньому орбітальному радіусі Землі d = 1,496 × 10¹¹ м сонячна стала (повний потік на верхній межі атмосфери) дорівнює:

S₀ = L☉ / (4πd²) ≈ 1361 Вт м⁻² Поглинуте сонячне випромінювання (ASR): ASR = S₀ · πR² · (1 − α) / (4πR²) = S₀(1−α)/4 α = планетарне альбедо ≈ 0.30 (30% сонячного світла відбивається) ASR ≈ 342 · 0.70 ≈ 239 Вт м⁻²

Множник 4 виникає тому, що Сонце освітлює диск площею πR², але планета випромінює зі своєї повної сфери 4πR². Внески в альбедо: свіжий сніг ≈ 0,85, океан ≈ 0,06, хмари ≈ 0,60 — у планетарному середньому домінують хмари.

2. Стефан-Больцман та ефективна температура

Абсолютно чорне тіло випромінює потужність, пропорційну T⁴. Прирівнявши випромінену потужність до поглинутої сонячної, отримуємо ефективну температуру випромінювання Tₑ:

OLR = σ T⁴ (вихідне довгохвильове випромінювання) σ = 5.67×10⁻⁸ Вт м⁻² K⁻⁴ (стала Стефана-Больцмана) Енергетичний баланс: σ Tₑ⁴ = S₀(1−α)/4 Tₑ = [S₀(1−α)/(4σ)]^¼ = [(1361 · 0.70)/(4 · 5.67×10⁻⁸)]^¼ ≈ 255 K (≈ −18°C)

Спостережувана середня приземна температура Землі — ~288 K (15°C). Різниця в 33 K — це природний парниковий ефект: атмосферні гази поглинають і перевипромінюють вихідне інфрачервоне випромінювання, нагріваючи поверхню вище за прогноз для голого чорного тіла.

3. Нульвимірна модель енергетичного балансу

0D-EBM розглядає всю планету як єдиний добре перемішаний резервуар з теплоємністю C:

C · dT/dt = ASR − OLR = S₀(1−α)/4 − σεT⁴ ε = ефективна випромінювальна здатність атмосфери (0 < ε ≤ 1) ε = 1 → чорне тіло (без парникового ефекту) ε < 1 → часткове парникове потепління Рівновага (dT/dt = 0): T* = [S₀(1−α)/(4σε)]^¼

Підлаштувавши ε ≈ 0,61, отримуємо T* = 288 K. Зменшення ε (сильніший парниковий ефект) збільшує T*. Стала часу EBM дорівнює τ = C/λ, де λ = dOLR/dT = 4σεT³ ≈ 3,3 Вт м⁻² K⁻¹ — це зворотний зв'язок Планка. За теплоємності перемішаного шару океану C ≈ 10⁸ Дж м⁻² K⁻¹, τ ≈ 30 років.

4. Парниковий ефект та радіаційне форсування

Парникові гази (CO₂, H₂O, CH₄, N₂O, O₃) поглинають вихідне довгохвильове випромінювання у певних спектральних смугах, зменшуючи OLR за заданої T. Утворений енергетичний дисбаланс — це радіаційне форсування ΔF (Вт м⁻²):

Форсування CO₂ (IPCC AR6): ΔF = α · ln(C/C₀) α ≈ 5.35 Вт м⁻² (коефіцієнт радіаційного форсування) C₀ = 280 ppm (доіндустріальний 1750) C = 421 ppm (2023) ΔF = 5.35 · ln(421/280) ≈ 2.17 Вт м⁻² Рівноважна кліматична чутливість (ECS): ΔT = ΔF / λ_eff λ_eff ≈ 1.0–1.2 Вт м⁻² K⁻¹ (Планк + зворотні зв'язки) ECS для подвоєння CO₂ ≈ 3 K (імовірний діапазон 2.5–4 K)

Водяна пара +

Потепління збільшує атмосферну водяну пару, яка й сама є парниковим газом — приблизно подвоюючи форсування CO₂. Найбільший позитивний зворотний зв'язок (~1,8 Вт м⁻² K⁻¹).

Градієнт температури −

Тропіки: вологий адіабатичний градієнт спричиняє більше нагрівання верхньої тропосфери, ніж поверхні — верхні шари випромінюють більше, частково компенсуючи водяну пару.

Крига-альбедо +

Потепління топить кригу з високим альбедо → темніші океан/суша поглинають більше сонячного → більше потепління. Особливо сильний в Арктиці.

Хмари −/+

Низькі хмари відбивають сонячне (охолодження), але також затримують довгохвильове (нагрівання). Чистий хмарний зворотний зв'язок ≈ +0,42 Вт м⁻² K⁻¹ — найбільш невизначений член.

5. Зворотні зв'язки альбедо-температура

Щоб урахувати зворотний зв'язок крига-альбедо, α стає функцією від T. Проста параметризація:

α(T) = α_ice якщо T < T_ice (≈ 263 K — повністю вкрита кригою) = α_ocean якщо T > T_melt (≈ 293 K — без криги) = лінійна інтерполяція посередині dOLR/dT − dASR/dT = λ_Planck + λ_WV + λ_LR + λ_alb + λ_cld Сумарний зворотний зв'язок λ_eff = Σ λᵢ λ_eff > 0 → стабільний (відновлювальний) λ_eff < 0 → некерований (нестабільний)

На Марсі зворотний зв'язок крига-альбедо настільки сильний, що усунення парникового форсування могло б перекинути його у стан «снігової кулі». На Венері некерований парниковий ефект зафіксував планету на 737 K.

6. Точки незворотності та бістабільність

За нелінійної α(T) рівняння енергетичного балансу C dT/dt = F(T) може мати кілька рівноваг. Якщо зобразити енергетичний бюджет як функцію від T, виявляється до трьох нерухомих точок: дві стабільні (тепла Земля, Земля-сніжка) та одна нестабільна (між ними).

F(T) = S₀[1−α(T)]/4 − σεT⁴ Стабільні рівноваги: F(T*) = 0 та dF/dT|_{T*} < 0 Нестабільна: F(T*) = 0 та dF/dT|_{T*} > 0 Незворотність: якщо форсування перевищує нестабільну точку, система незворотно «падає» у протилежний стабільний стан.

Реальні елементи незворотності Землі включають Західноантарктичний льодовиковий щит, Гренландський льодовиковий щит, Атлантичну меридіональну циркуляцію (AMOC) та відмирання Амазонії — кожен зі своїм пороговим форсуванням та часовим масштабом незворотності.

7. Симуляція EBM на JavaScript

// Нульвимірна модель енергетичного балансу зі зворотним зв'язком крига-альбедо
const sigma = 5.67e-8;  // Стефан-Больцман
const S0    = 1361;     // Вт/м², сонячна стала
const C     = 4e8;      // Дж/(м²·K) теплоємність перемішаного шару

function albedo(T) {
  if (T < 263) return 0.62; // снігова куля
  if (T > 293) return 0.28; // без криги
  return 0.62 + (0.28 - 0.62) * (T - 263) / (293 - 263);
}

function dTdt(T, eps, dForcing = 0) {
  const ASR = S0 * (1 - albedo(T)) / 4;
  const OLR = sigma * eps * T ** 4;
  return (ASR - OLR + dForcing) / C;
}

function runEBM(T0, eps, dForcing, years = 500, dt = 3.15e7) {
  // dt = 1 рік у секундах
  let T = T0;
  const result = [{year: 0, T}];
  const steps = years * 12; // місячні часові кроки
  const dtM = dt / 12;
  for (let i = 1; i <= steps; i++) {
    T += dTdt(T, eps, dForcing) * dtM;
    if (i % 12 === 0) result.push({year: i / 12, T});
  }
  return result;
}

// Стабільна сучасна Земля
const modern = runEBM(288, 0.61, 0);
console.log('Рівноважна T:', modern.at(-1).T.toFixed(2), 'K');

// Подвоєння CO₂: форсування +3.7 Вт/м²
const co2x2 = runEBM(288, 0.61, 3.7);
console.log('ΔT для CO₂×2:', (co2x2.at(-1).T - 288).toFixed(2), 'K');

// Земля-сніжка із сучасного стану (збільшення альбедо)
const snowball = runEBM(288, 0.61, -40); // −40 Вт/м² → запускає крижану лавину
console.log('T сніжки:', snowball.at(-1).T.toFixed(2), 'K');

8. Ієрархія кліматичних моделей

Симуляція: EBM вище демонструє ту саму емерджентну бістабільність, що й SIR-модель епідемії — нагадування, що мультистабільність є універсальною рисою нелінійних динамічних систем.
🌍 Відкрити «Атмосферу» →