Арктичний сніг зовсім не схожий на м’яку пухку ковдру, до якої звикли мешканці помірних широт. У тундрі це тонкий, вітром вичищений і різко шаруватий шар кристалів льоду, від якого залежить, чи промерзне ґрунт, чи розтане мерзлота, а також чи переживуть зиму північні олені та лемінги. Саме тому команда дослідників суттєво оновила сніговий блок у моделі CLASSIC, що є складовою Канадської системної кліматичної моделі. Про роботу розповідає видання Scienmag.

Модель CLASSIC використовується не лише для наукових експериментів: це стандартний інструмент Канади для оцінки наслідків зміни клімату, у тому числі в межах доповідей МГЕЗК. Якщо снігові процеси в Арктиці в ній описані неправильно, помилки поширюються на моделювання температури ґрунту, відтавання мерзлоти та поведінки величезних запасів вуглецю в арктичних ґрунтах, яких майже вдвічі більше, ніж у всій атмосфері.
У чому помилялася стара модель снігу
Команда під керівництвом Мікаеля Лаланда (Mickaël Lalande) з Університету Квебека в Труа-Рів’єр проаналізувала роботу моделі на трьох арктичних майданчиках: острові Білот, на ділянці Umiujaq TUNDRA та в районі Trail Valley Creek. Базова версія CLASSIC систематично завищувала глибину снігу приблизно на 10–20 сантиметрів, що є дуже великою похибкою для снігового покриву, який рідко перевищує пів метра. Водночас модель давала занадто високі зимові температури ґрунту — на 5–10 °C тепліші за виміряні.
Головною причиною виявився надто «пухкий» і товстий сніг у розрахунках. Схема ущільнення обмежувала максимальну густину сухого снігу значеннями, значно нижчими за реальні в умовах вітряної тундри, де об’ємна густина 250–400 кілограмів на кубічний метр є звичайним явищем. Такий надмірно пухкий сніг у моделі працював як сильний утеплювач і штучно підігрівав ґрунт унизу.
Нова схема ущільнення: вітер переформатовує сніг
Ключова інновація дослідження — нова схема ущільнення снігу, що залежить від вітру. Коли швидкість вітру перевищує приблизно 2,5 метра за секунду, модель тепер дозволяє тонким сніговим покривам ущільнюватися до вищих граничних густин. Це імітує формування щільної вітрової кірки, яка в природі накриває арктичний сніг зверху.
Спеціальна гаусова функція ваг розташовує максимальний ефект вітру навколо глибини снігу близько 80 сантиметрів. Завдяки цьому глибокі снігові покриви в середніх широтах, де домінують гравітаційне осідання та метаморфізм кристалів, практично не зачіпаються. Перевага підходу в тому, що одношарова модель може неявно відтворити типову для тундри двошарову будову: низькогустинний нижній шар (depth hoar) та щільну вітрову кірку зверху — без явного розділення їх у коді.
Результат виявився разючим: середнє зміщення в оцінках глибини снігу на арктичних майданчиках скоротилося з 12,0 сантиметра до лише 0,1 сантиметра.
Виправлення «теплової сантехніки» снігу й ґрунту
Втім, однієї нової схеми ущільнення було замало. Дослідники помітили приховану проблему в розрахунку температури на межі «сніг–ґрунт». Стара формула надто сильно орієнтувалася на усереднену температуру снігу, недооцінюючи температуру його нижнього шару. Це спричиняло надмірний теплопотік від ґрунту до снігу, штучно підігріваючи сніг і водночас маскуючи теплу похибку ґрунту.
Новий підхід до обчислення температури нижнього шару снігу враховує температуру снігу і ґрунту з вагами, оберненими до товщини відповідних шарів. Фізика стала коректнішою, але на перших порах це навіть погіршило похибки ґрунтової температури, оскільки стара параметризація теплопровідності снігу «тихо» компенсувала помилки.
Тоді команда замінила широко вживану залежність теплопровідності за Sturm та співавторами, що базувалася на вимірюваннях нагрітим голчастим зондом (нині відомо, що вони недооцінювали теплопровідність крихкого шаруватого depth hoar), на новішу параметризацію Calonne та колег. Остання була відкалібрована за допомогою високороздільних тривимірних зображень мікроструктури снігу і дає вищі значення теплопровідності. У моделі це дозволяє надлишковому теплу з ґрунту виходити крізь сніг в атмосферу, а не накопичуватися в землі.
У сукупності з іншими змінами це скоротило середнє кореневе середньоквадратичне відхилення змодельованих температур ґрунту на арктичних майданчиках з 5,3 до 3,0 °C. На острові Білот зимове зміщення ґрунтової температури зменшилося майже з 7 до менш ніж 1 °C.
Як модель навчили краще відтворювати поверхню снігу
Окрема проблема стосувалася температури снігової поверхні. В умовах стійкої зимової стратифікації атмосфери стандартна теорія турбулентності Моніна–Обухова фактично «відключає» теплообмін, через що моделі часто дають надто холодну поверхню снігу — похибка могла сягати близько 10 °C у попередніх версіях.
Дослідники активували так званий коефіцієнт безвітрового обміну, запропонований Россом Брауном (Ross Brown) та колегами ще два десятиліття тому. Він додає невеликий базовий турбулентний теплообмін у 2 ватти на квадратний метр на кельвін, коли поверхня холодніша за повітря. На семи середньоширотних та альпійських майданчиках SnowMIP цей простий крок знизив середню похибку температури снігової поверхні з 3,50 до 1,93 °C, практично не змінюючи глибину снігу й температуру ґрунту.
Роль знесення снігу та перевірка на різних ділянках
Команда також випробувала арктично специфічний процес — сублімацію знесеного вітром снігу, коли частинки випаровуються просто в повітрі. Цей механізм реалізували за емпіричною формулою Gordon та колег. Виявилося, що він помітно впливає лише на найбільш вітряних станціях — Umiujaq, Trail Valley Creek та альпійському басейні Senator Beck, де зменшував змодельовану глибину снігу на 2–4 сантиметри. На більшості ділянок ефект був мінімальним, що нагадує: не кожен арктичний процес потребує складної параметризації, аби бути відтвореним справедливо.
Важливо, що покращення арктичного блоку не погіршили роботу моделі в інших умовах. Оцінки на майданчиках SnowMIP, серед яких класичні альпійські та середньоширотні локації на кшталт Col de Porte у Франції та Weissfluhjoch у Швейцарії, показали: фінальна конфігурація зберегла або навіть покращила точність, а помилки в глибині снігу зменшилися. Інтегральний показник якості моделі за всіма ділянками й змінними зріс з 0,284 до 0,380.
Станцію Trail Valley Creek спеціально не використовували для налаштування параметрів, залишивши її незалежним майданчиком перевірки. Дослідники також опублікували змінену версію параметрів ущільнення, що враховує й цю ділянку, підкреслюючи великі невизначеності в вимірюванні снігопадів на такій вітряній станції.
Чому сніг критичний для арктичного вуглецевого балансу
Наслідки оновлення моделі виходять далеко за межі «правильної» глибини снігу. Уточнені зимові температури ґрунту змінили змодельований чистий обмін вуглекислого газу між екосистемою та атмосферою на арктичних ділянках. У деяких випадках річний баланс вуглецю зсувався від джерела в бік поглинача, адже холодніший ґрунт означає менше зимове дихання ґрунтових організмів.
Ця чутливість демонструє, що снігова фізика в Арктиці — не вузька технічна деталь, а своєрідний «вартовий» вуглецевого зворотного зв’язку, який може посилити глобальне потепління.
Автори відкрито говорять і про обмеження. У моделі все ще немає явних шарів моху та лишайників, не враховано теплопровідні «містки» через чагарники, а арктичні польові вимірювання ускладнюють іній на сенсорах, полярна ніч і знесений сніг, що додає невизначеностей і в дані, і в спостереження.
Наступний етап роботи — перенести точкові покращення в масштаб панарктичних просторових симуляцій, де параметризації частки снігового покриву мають впоратися з повною мозаїкою умов у межах кожної клітинки сітки. Наразі ж дослідження показує, що навіть відносно проста одношарова модель снігу, якщо надати їй правильну фізику, може впевнено працювати в суворих умовах Арктики.
Оновлена модель снігу точніше відтворює арктичну мерзлоту з’явилася спочатку на Цікавості.

1140


