Сезонна інсоляція · Вплив CO₂ · Зворотний зв'язок альбедо · Площа льоду
Ця симуляція моделює динаміку арктичного морського льоду на спрощеній полярній сітці. Кожна комірка може бути відкритим океаном або льодом, а її стан залежить від локального енергетичного балансу: вхідна сонячна радіація (яка змінюється з сезоном і широтою), відбита радіація (визначається альбедо — лід відбиває ~60 % сонячного світла, тоді як відкрита вода поглинає ~94 %) та вихідне довгохвильове випромінювання, модифіковане парниковим ефектом. Ключова ідея — позитивний зворотний зв'язок лід–альбедо: коли лід тане, темніший океан поглинає більше тепла, прискорюючи подальше танення, що потенційно може спричинити різкий зламний момент.
Площа арктичного морського льоду зменшується приблизно на 13 % за десятиліття (вересневий мінімум) з початку супутникових спостережень у 1979 році. Зворотний зв'язок лід–альбедо є одним з найсильніших підсилюючих механізмів кліматичної системи. Під час Еоценового термального максимуму (~50 млн років тому) CO₂ перевищував 1 000 ppm, й Арктика була вільна від льоду — крокодили мешкали за Полярним колом.
Ця симуляція моделює сітку 40×60 арктичних комірок, кожна з яких зберігає частку льоду та температурну аномалію, що змінюються місяць за місяцем за спрощеним локальним енергетичним балансом. Поглинута сонячна енергія залежить від сезонної інсоляції та альбедо кожної комірки (лід ≈0,62, відкритий океан ≈0,06), тоді як вихідне довгохвильове випромінювання послаблюється логарифмічним членом парникового впливу CO₂, ln(CO₂/280)/ln(2). Отримана зміна температури спричиняє замерзання нижче −2 °C та танення вище 0 °C, а проста дифузія згладжує кожну комірку відносно сусідніх. Оскільки танення оголює темніший, більш теплопоглинальний океан, модель відтворює петлю зворотного зв'язку лід–альбедо, здатну спричинити неконтрольоване відступання арктичного льоду.
Полярна сітка, де кожна з 2400 комірок відстежує власну частку льоду (0–1) та температуру. Вкриті льодом комірки відображаються біло-блакитними, відкритий океан — темно-синім з температурним відтінком. Додатковий лінійний графік показує загальну площу льоду (частку сітки, що досі замерзла) протягом симульованих десятиліть, роблячи зламні точки зворотного зв'язку видимими.
Перетягніть повзунок Рівня CO₂ (280–1200 ppm), щоб змінити парниковий вплив, та Множник Сонця (80–120%), щоб імітувати зміни інсоляції на кшталт циклів Міланковича. Швидкість (1×–10×) визначає, скільки симульованих років минає за секунду. Увімкніть Сітку, щоб побачити окремі комірки, та Графік площі, щоб показати/приховати історію площі льоду. Пауза/Скинути перезапускають симуляцію із новим випадковим розподілом льоду.
Реальні супутникові дані з 1979 року показують, що вересневий мінімум площі арктичного морського льоду зменшується приблизно на 13% за десятиліття — тенденція, яку якісно відтворює петля зворотного зв'язку цієї симуляції. Під час Еоценового термального максимуму атмосферний CO₂ перевищував 1000 ppm, а Арктика була настільки вільною від льоду, що крокодили мешкали біля полюса.
Лід відбиває більшість вхідного сонячного світла (високе альбедо, близько 0,6), тоді як відкритий океан поглинає більшість цього світла (низьке альбедо, близько 0,06). Коли лід тане, темніша поверхня океану поглинає більше сонячної енергії, ще більше нагріваючи воду й розтоплюючи ще більше льоду. Ця самопідсилювана петля — один із найсильніших підсилювачів у кліматичній системі Арктики, і саме її відтворює розрахунок альбедо для кожної комірки в симуляції.
Повзунок CO₂ впливає на коефіцієнт парникового ефекту, який обчислюється як 1 + 0,7·log₂(CO₂/280). Вищий CO₂ збільшує цей коефіцієнт, що зменшує змодельоване вихідне довгохвильове випромінювання кожної комірки. Оскільки зміна температури кожної комірки залежить від поглинутої сонячної енергії мінус вихідне випромінювання, підвищення CO₂ схиляє більше комірок до чистого потепління, сприяючи таненню замість замерзання.
Кожен рядок представляє смугу широти, де рядок 0 моделює полюс, а останній рядок — край області поблизу екватора. Початкова ймовірність льоду масштабується широтою у степені (lat^1,2), тож комірки ближче до полюса мають вищий шанс бути льодом від початку. Сезонна інсоляція також зважена за широтою, тому полярні комірки отримують менше середнього сонячного світла й мають тенденцію залишатися холоднішими.
Кожна комірка має власну температурну аномалію, що оновлюється за локальним енергетичним балансом щомісяця симуляції. Якщо температура комірки опускається нижче −2 °C, її частка льоду зростає; якщо піднімається вище 0 °C, частка льоду зменшується. Між цими порогами крижаний покрив комірки залишається приблизно стабільним, тому вузька смуга температур відповідає прикордонній зоні льоду, де площа найчутливіша до впливу.
Це спрощена концептуальна модель, а не модель загальної циркуляції. Вона відтворює основні механізми сезонної інсоляції, зворотного зв'язку альбедо та парникового впливу на грубій сітці 40×60 з базовою дифузією між сусідами, чого достатньо для відтворення якісної поведінки, як-от підсилення зворотного зв'язку та зламних точок. Вона не враховує перенесення тепла океаном, зворотні зв'язки хмар, товщину льоду та динаміку, тож абсолютні цифри слід сприймати як ілюстративні, а не прогностичні.