Спотлайт #31 – Наука про клімат: радіаційне форсування, петлі зворотного зв'язку, вуглецевий цикл та точки неповернення

Наука про клімат застосовує ті самі фізичні закони, що й решта цього блогу — перенесення випромінювання, гідродинаміку, термодинаміку, геохімію — до системи, яка найбезпосередніше впливає на людську цивілізацію. Атмосфера — це тонка оболонка (99% маси нижче 30 км на планеті радіусом 6371 км), що вловлює інфрачервоне випромінювання, перерозподіляє тепло через гігантські конвективні комірки й обмінюється вуглецем з океаном та наземною біосферою на часових масштабах від днів до мільйонів років.

1. Радіаційне форсування та парниковий ефект

Кліматична система зрештою рухається дисбалансом між вхідним сонячним короткохвильовим випромінюванням та вихідним довгохвильовим випромінюванням, яке випромінює Земля. Парниковий ефект — поглинання й повторне випромінювання інфрачервоного випромінювання атмосферними газами — підвищує температуру поверхні Землі з замерзлого −18 °C (рівновага чорного тіла голої планети) до спостережуваного +15 °C. Будь-яке зовнішнє збурення, що порушує цей баланс, є радіаційним форсуванням.

Нульвимірна модель енергетичного балансу:
  C dT/dt = S_0(1 − α)/4 − εσT&sup4;
  S_0 = 1361 Вт/м² (сонячна стала)
  α = 0,30 (планетарне альбедо)
  ε ≈ 0,78 (ефективна випромінювальна здатність; <1 через парниковий ефект)
  σ = 5,67 × 10^−8 Вт м^−2 К^−4 (стала Стефана-Больцмана)

Рівноважна температура (при dT/dt = 0):
  T_рівн = [S_0(1−α)/(4εσ)]^{1/4}
  Без парникового ефекту (ε=1): T_рівн = 255 К = −18 °C
  З реальною атмосферою (ε=0,78): T_рівн ≈ 288 К = +15 °C
  Внесок парникового потепління: ~33 °C

Радіаційне форсування ΔF [Вт/м²]:
  Зовнішньо накладений дисбаланс у Вт/м² до дії зворотних зв'язків.
  Подвоєння CO&sub2;:   ΔF = 5,35 ln(C/C_0) ≈ 3,7 Вт/м²  (найкраща оцінка IPCC)
  Доіндустріальний CO&sub2; = 280 ppm; значення 2024 року = 422 ppm → ΔF ≈ 2,2 Вт/м²
  Рівноважна кліматична чутливість (ECS): ΔT = λ ΔF
  λ ≈ 0,8 К/(Вт/м²) (найкраща оцінка IPCC AR6)
  ECS = λ · 3,7 ≈ 3 °C на подвоєння CO&sub2; (ймовірний діапазон 2,5–4 °C)

Інші агенти форсування (значення AR6, Вт/м²):
  CH&sub4;:              +0,54
  N&sub2;O:              +0,28
  Галогенвуглеці:      +0,36
  Аерозолі (сумарно): −1,06 (охолодження)
  Зміна землекористування:  −0,20
  Загальне антропогенне 1750–2019: +2,72 Вт/м²

2. Кліматичний зворотний зв'язок — підсилювачі та стабілізатори

Зворотні зв'язки підсилюють або гасять початкову реакцію на форсування. Кліматична система має як позитивні (підсилювальні), так і негативні (стабілізувальні) зворотні зв'язки. Сумарний зворотний зв'язок визначає рівноважну кліматичну чутливість. IPCC AR6 (2021) визначив зворотний зв'язок Планка як домінантний стабілізувальний механізм, зворотні зв'язки льодового альбедо й водяної пари як домінантні підсилювачі, а хмарні зворотні зв'язки як найбільше джерело невизначеності в кліматичних прогнозах.

Рамка зворотного зв'язку:
  ΔT = ΔF / (−λ^−1)  де λ = λ_Планка + Σ_i f_i
  Позитивне f_i → підсилення; негативне f_i → гасіння.

Зворотний зв'язок Планка (стабілізатор чорного тіла):
  f_Планк = −d(4εσT³)/dT ≈ −3,3 Вт м^−2 К^−1  (завжди негативний — стабілізуючий)
  Фізичний механізм: тепліша планета випромінює більше, відновлюючи баланс.

Зворотний зв'язок водяної пари:
  f_ВП ≈ +1,8 Вт м^−2 К^−1   (підсилювальний)
  Механізм: тепліше повітря утримує більше водяної пари (Клаузіус-Клапейрон);
  H&sub2;O — парниковий газ → більше потепління → більше H&sub2;O (підсилювальна петля).
  Приблизно подвоює лише чутливість зворотного зв'язку Планка (ECS зростає з ~1 °C до ~2 °C).

Зворотний зв'язок льодового альбедо:
  f_ЛІД ≈ +0,3 Вт м^−2 К^−1   (підсилювальний)
  Механізм: потепління плавить лід → оголюється темніший океан/суша (α з 0,9 до 0,06)
  → більше поглинання → більше потепління. Особливо сильний в Арктиці.

Зворотний зв'язок градієнта температури (лапс-рейт):
  f_ГТ ≈ −0,4 Вт м^−2 К^−1   (слабко стабілізує в тропіках, підсилює на полюсах)
  Верхня тропосфера тропіків нагрівається швидше за поверхню → збільшує вихідне ДХ-випромінювання.
  Полярні регіони: слабкий градієнтний зворотний зв'язок → полярне підсилення (Арктика теплішає в 3 рази швидше за глобальне середнє).

Хмарні зворотні зв'язки (найбільша невизначеність):
  Ймовірний діапазон IPCC AR6: −0,2 до +0,8 Вт м^−2 К^−1
  Низькі морські хмари (стратокумулюс): зменшуються з потеплінням → позитивний зворотний зв'язок
  Високі перисті хмари: складно; чистий знак невизначений
  Невизначеність хмарного зворотного зв'язку охоплює ~1 °C невизначеності ECS.

Сумарний зворотний зв'язок та ECS:
  Сумарне λ = f_Планк + f_ВП + f_ГТ + f_ЛІД + f_хмари ≈ −1,3 Вт м^−2 К^−1
  ECS = 3,7 / 1,3 ≈ 2,85 °C на подвоєння CO&sub2; (IPCC: ймовірно 2,5–4,0 °C)

3. Глобальний вуглецевий цикл

Вуглець циркулює між атмосферою, океаном, наземною біосферою та геологічними резервуарами на часових масштабах від днів (фотосинтез/дихання) до мільйонів років (вивітрювання порід та вулканізм). Антропогенні викиди (~11 ГтC/рік у 2023 році) збурюють швидкий цикл із швидкістю приблизно у 100 разів вищою за будь-яку геологічну подію вкидання вуглецю за останні 66 мільйонів років.

Основні резервуари (ГтC = гігатонни вуглецю):
  Атмосфера:                    880 ГтC   (2023; доіндустріальна 600 ГтC)
  Наземна біосфера (жива):  600 ГтC
  Органічний вуглець ґрунту:         1 500 ГтC
  Вічна мерзлота (CH&sub4;+CO&sub2;):       1 700 ГтC   (вразлива до потепління)
  Поверхневий океан:                 900 ГтC
  Глибокий океан:                 37 000 ГтC
  Літосфера (осадові породи):   5 000 000 ГтC  (надзвичайно повільний обмін)

Поточні річні потоки (ГтC/рік, приблизно 2020-ті):
  Викиди від викопного палива + цементу:    +10,5
  Зміна землекористування:                   +1,2
  Наземне поглинання (чисте):          −3,1
  Океанічне поглинання (чисте):                −3,0
  Атмосферне накопичення:          +5,6  (близько 51% викидів залишаються в повітрі)

Механізми поглинання вуглецю океаном:
  1. Розчинювальний насос: CO&sub2; розчиняється в холодній поверхневій воді (закон Генрі: s ∝ 1/T)
  2. Біологічний насос: фітопланктон фіксує CO&sub2; → органічний вуглець осідає з мертвими клітинами
     Ефективність < 25% (більшість органічної речовини мінералізується до досягнення дна)
  Закислення океану: pH знизився з 8,25 (доіндустріальний) до 8,06 (2023)
  ΔpH = −0,19 → [H&sup+;] зросла на 55%  (логарифмічна шкала pH)

Тривалість перебування CO&sub2; в атмосфері:
  Частка імпульсу, що залишається в повітрі після:
    10 років:     ~64%
   100 років:     ~40%
  1 000 років:    ~15–25%
  10 000+ років:  ~5–10%
  Немає єдиної часової шкали: рівновага з океаном (століття), вивітрювання (десятки тисячоліть).

4. Елементи неповернення — нелінійні переходи

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

Оцінка порогів елементів неповернення McKay та ін. (2022) (глобальне потепління вище доіндустріального):

Кріосфера:
  Західноантарктичний льодовий щит (WAIS): 1,5–3 °C  (колапс упродовж століть; підйом рівня моря 3–5 м)
  Гренландський льодовий щит:             0,8–3 °C  (майже незворотний вище 1,5 °C)
  Втрата багаторічного арктичного морського льоду:   ~1,5–2 °C  (спочатку сезонна, потім багаторічна літня втрата)
  ТН гірських льодовиків:             ~2 °C        (майже повна втрата; безпека водних ресурсів)

Океанічна циркуляція:
  Колапс AMOC (Атлантична термохалінна циркуляція):  3–5 °C  (можливо, вже дестабілізується)
  Різке ослаблення AMOC приносить похолодання Північній Європі, посилює посуху в Сахелі.

Екосистеми:
  Всихання Амазонії:                  3,5 °C  (взаємодія посухи й вирубки; досягнуто 17%)
  Зсув бореального лісу на північ:   4 °C
  Знебарвлення коралових рифів (>90% втрати): ~1,5–2 °C  (уже 67% пережили масове знебарвлення у 2023)

Метан вічної мерзлоти:
  Межа погано визначена; поступове танення вічної мерзлоти може не бути жорсткою точкою неповернення,
  але додає значну (0,1–0,4 ГтC/рік додатково до 2100 року за сценаріями з високими викидами)

Ризик каскаду неповернення:
  Lenton та ін. (2019); Armstrong McKay та ін. (2022): перетин однієї ТН при 2 °C
  приблизно подвоює ймовірність запуску іншої впродовж десятиліть.
  При 2 °C: ~9 ТН із порогами нижче 2 °C можуть наближатись або бути перевищеними.

5. Атмосферна циркуляція — комірки Гадлі, Феррела та полярна

Диференційне сонячне нагрівання між екватором і полюсами приводить у рух глобальну атмосферну циркуляцію, що перерозподіляє тепло до полюсів. Ефект Коріоліса від обертання Землі руйнує однокоміркову конвекцію, очікувану від чистого теплового приводу, і створює три меридіональні циркуляційні комірки в кожній півкулі: комірки Гадлі, Феррела та полярну. Ці комірки визначають великомасштабні моделі опадів — тропічні дощові пояси, субтропічні пустелі, помірні шторм-траки та полярні зони високого тиску.

Меридіональні циркуляційні комірки:
  Комірка Гадлі (0–30°):
    Висхідна гілка: ВТЗК (внутрішньотропічна зона конвергенції) поблизу екватора
    Низхідна гілка: субтропічні пояси високого тиску (пустелі Сахара, Аравійська, Атакама)
    Пасати: східні вітри на поверхні, зумовлені відхиленим до полюса повітрям (Коріоліс)
  Комірка Феррела (30–60°):  непряма термічно-керована комірка; західні вітри на поверхні
  Полярна комірка (60–90°):  холодне низхідне повітря на полюсах; східні вітри на поверхні

Потенційна температура та CAPE:
  Висхідна гілка Гадлі: CAPE ~ 1 000–3 000 Дж/кг (енергія тропічної конвекції)
  Волога комірка Гадлі: масштабування Голда-Суареса → межа комірки до полюса ∝ T_s / Ω  (Ω = швидкість обертання Землі)

Струменеві течії:
  Субтропічна струменева течія: ~12 км висоти, ~30–35° широти, v ~ 60 м/с
  Полярна фронтальна струменева течія: ~9 км, ~55–65°, v ~ 30–100 м/с (західні вітри, шторм-траки)
  Термічний вітровий баланс: ∂u/∂z = −(g/fT)(∂T/∂y)   (f = коефіцієнт Коріоліса)
  Потепління (арктичне підсилення) → зменшення ∂T/∂y → слабша струменева течія → більш хвиляста, повільніше меандрує

Розширення комірки Гадлі (зсув до полюса при потеплінні):
  Спостережено: ~1–2° широти зсуву на °C глобального потепління
  Наслідок: субтропічні сухі зони розширюються до полюсів (висихання Середземномор'я, ПЗ США, ПЗ Австралії)

ENSO (Ель-Ніньйо — Південна осциляція):
  Циркуляція Уокера в тропіках: зі сходу на захід уздовж екваторіальної частини Тихого океану
  Ель-Ніньйо: ослаблені пасати → тепла водна маса зміщується на схід → глобальні телезв'язки
  ΔT (регіон Ніньйо-3.4) > +0,5 °C впродовж 5 послідовних місяців = подія Ель-Ніньйо
  PDO та AMO (декадні осциляції) модулюють амплітуду ENSO та тренд глобальної температури.

6. Кліматичні моделі — від нульвимірної до глобальних зв'язаних GCM

Кліматичні моделі охоплюють ієрархію від простої нульвимірної моделі енергетичного балансу (одне рівняння, одна змінна) до повністю зв'язаних моделей загальної циркуляції (GCM) із сотнями мільйонів комірок сітки в атмосфері, океані, суші та морському льоді. Кожен щабель ієрархії додає фізичні процеси ціною обчислювальних витрат та структурної невизначеності.

Ієрархія моделей:
  0-В МЕБ:       C dT/dt = F_вх − F_вих   (навчання концепцій; швидко)
  1-В МЕБ:       додає широту; меридіональне перенесення тепла Q = −D d(T)/dφ
  2-В РКМ:       радіаційно-конвективна модель; вертикальна структура; атмосферні стовпці
  AOGCM:         зв'язана GCM атмосфера-океан (моделі CMIP5/6)
  ESM:           модель земної системи = AOGCM + інтерактивний вуглецевий цикл, хімія тощо

Горизонтальна роздільна здатність (типова CMIP6):
  Атмосфера: 100–200 км (реконструкція), 25–50 км (HighResMIP)
  Океан:      25–100 км (з параметризацією вихорів), 10 км (з розв'язанням вихорів: дуже витратно)
  Субсіткові процеси (мікрофізика хмар, конвекція, турбулентність) параметризуються.

Ансамбль моделей CMIP6 (Zelinka та ін., 2020):
  ~ 40 моделей із 20+ інститутів по всьому світу
  Історичне потепління 1850–2014: усі моделі відтворюють спостережуваний запис з точністю ~0,2 °C
  Діапазон ECS у CMIP6: 1,8–5,6 °C (ширший за CMIP5: 2,1–4,7 °C)
  Моделі з високою ECS (вище 5 °C) тепер оцінюються як малоймовірні (обмеження IPCC AR6).

Прогнози IPCC AR6 (2021) для 2100 року порівняно з базовим рівнем 1850–1900:
  SSP1-1.9 (агресивне пом'якшення): +1,0 – +1,8 °C
  SSP2-4.5 (проміжний):          +2,1 – +3,5 °C
  SSP5-8.5 (інтенсивне використання викопного палива):  +3,3 – +5,7 °C
  Підйом рівня моря (ймовірний діапазон):     0,32–1,01 м (SSP1–SSP5)
  Примітка: спостережуване потепління 2024 року ~ +1,45 °C вище доіндустріального середнього.