Гидрологические последствия изменения климата и антропогенного воздействия в криолитозоне (РНФ)НИР

Hydrological outcomes of climate change and anthropogenic impacts in the cryolithozone

Источник финансирования НИР

грант РНФ

Этапы НИР

# Сроки Название
1 7 мая 2024 г.-31 декабря 2024 г. Гидрологические последствия изменения климата и антропогенного воздействия в криолитозоне
Результаты этапа: В 2024 г. проведен комплекс исследований, направленных на обобщение ранее полученных результатов в области изучения динамики стока рек криолитозоны и механизмов связи ее с изменением метеорологических характеристик, морфологического строения водосборов и параметров вечной мерзлоты, создание новых и расширение существующих баз данных, установление основных особенностей и закономерностей изменений стока и водного режима, температурных условий, стока наносов рек под влиянием меняющего климата и физико-географических условий. Выполнен обзор многочисленной отечественной и зарубежной литературы по оценке влияния изменений климата на сток рек и водный режим рек криолитозоны, термический режим, сток наносов, опасность наводнений и заторов льда на реках. Актуализированы имеющиеся базы данных гидрологических и метеорологических характеристик суточного и месячного разрешения, данных GRACЕ. Существенно расширена имеющаяся база данных по средним месячным, годовым и максимальным расходам взвешенных наносов, мутности воды. Для бассейнов арктических рек была собрана база данных различных спутниковых продуктов для оценки снегозапасов, испарения, влажности почвы, общих влагозапасов для сравнения с наземными наблюдениями. Анализ данных наблюдений за стоком рек включал оценки периодов нарушения стационарности в рядах годового и максимального стока, изменений средних и дисперсий за выбранные периоды времени, анализ трендов и построение соответствующих карт. Нарушение стационарности в рядах годового стока для Западной Сибири и значительной части бассейна Енисея произошло примерно в 1985 г., для бассейна Байкала, Амура, Лены – в 1996 г. Для рядов максимальных расходов воды получена схожая картина - 1980 и 1988 гг. в пределах ЕТР сменяются на 1985 г. для Западной Сибири, среднего и нижнего течения Енисея и 1996 г. - для бассейна Амура, Лены и северо-востока России. В пределах АТР реки с существенно положительными аномалиями годового стока (10-50% и больше) занимают нижнюю часть бассейна Иртыша и Обь-Иртышское междуречье, огромную территорию от Хатанги до Колымы. Существенно меньшую площадь занимают два региона со снижением годового стока - верховья Лены (южнее 60° с.ш.), бассейн Селенги, верховья Амура, а также Чукотку, сев. побережье Охотского моря, большую часть Камчатки. Ситуация с максимальными за год расходами воды рек более разнообразная. Здесь явно преобладают реки с положительными аномалиями. Причем реки на юго-востоке АТР со значительной ролью летне-осенних дождевых паводков демонстрируют уменьшение максимальных расходов, тогда как у рек с основной фазой в весеннее половодье доминирует их увеличение. Результаты калибровки и верификации модели формирования стока ECOMAG для ключевых бассейнов криолитозоны показали, что для исследуемых створов получены хорошие результаты расчета суточных гидрографов как по данным метеостанций, так и на основе баз данных реанализов, Сравнение данных спутниковых продуктов о всех доступных компонентах водного баланса по 20 бассейнам в арктической зоне с данными сетевых наблюдений продемонстрировало высокую корреляционную связь месячных осадков IMERG с осадками на 43 метеостанциях, значимую корреляцию спутниковых данных о снегозапасах с данными натурных измерений на 31 снеговом профиле, значимую корреляционную связь между стоком и суммарным испарением для Печеры, Енисея и Ангары. Связь между годовым стоком и интегральными запасами воды на водосборе (TWC), включающими запасы грунтовых вод и водохранилищ, наблюдалась для 13 водосборов. Результаты сравнения годовых величин показали, что на большинстве водосборов наилучшее совпадение модельной (ECOMAG) и спутниковой межгодовой изменчивости наблюдается для снегозапасов (SWE) и интегральных запасов воды (TWC). В бассейне Печоры наблюдается высокая связь межгодового хода модельной и спутниковой эвапотранспирации. При анализе среднемесячных величин было выявлено, что модельные и спутниковые параметры SWE и TWC демонстрируют хороший согласованный сезонный ход и его межгодовую изменчивость на большинстве модельных водосборов, месячные значения эвапотранспирации модельно-спутникового продукта GLEAMS показывают высокую корреляцию с эвапотранспирацией, рассчитанной по модели ECOMAG. В то же время для влажности почвы модельная изменчивость межгодового и сезонного цикла слабо связана с изменчивостью, полученной по спутниковым измерениям. Спутниковые данные о суммарных влагозапасах (данные GRACE), как интегральные характеристики влагозапасов, могут применяться для итоговой верификации модели. Основные результаты изучения термического режима рек связаны с исследованием роли различных географических факторов. Показано, что роль температуры воздуха в формировании температуры воды (в среднемесячном масштабе времени) растёт с северо-востока на юго-запад и с юга на север в горных районах восточнее р. Лена. В то же время показана незначительная роль атмосферных осадков и стока воды в формировании термического режима рек, дана количественная оценка степени изменчивости связи температуры воды с различными её факторами. Получены количественные оценки температуры воды в реках в современный период и ее изменение. Кроме того, оценён тепловой сток, формируемый реками различного размера. На основе дополнения и актуализации баз данных по стоку взвешенных наносов и мутности воды рек региона проведен анализ материалов многолетних наблюдений, установлены факторы его временной изменчивости, в том числе посредством построения многомерных зависимостей, проведен расчет основных характеристик. В 2024 г. такая работа проводилась для рек бассейнов Индигирки и Алдана, имеющих общую границу, схожие антропогенные условия формирования и руслового транзита наносов, но различающихся по величине и естественному режиму стока наносов. Детально изучена хронология, структура и география горнодобычи в бассейнах рек, влияние климата и его изменения, стока и водного режима рек, температуры воды и теплового стока, состояние вечной мерзлоты и ландшафтного строения данных водосборов.
2 1 января 2025 г.-31 декабря 2025 г. Гидрологические последствия изменения климата и антропогенного воздействия в криолитозоне
Результаты этапа: В 2025 г. планировалось продолжение комплекса научных исследований, направленных на обобщение полученных результатов в области изучения динамики стока рек криолитозоны и механизмов связи ее с изменением метеорологических характеристик, морфологического строения водосборов и параметров вечной мерзлоты, создание новых и расширение существующих баз данных, установление основных особенностей и закономерностей изменений стока и водного режима, температурных условий, стока наносов рек под влиянием меняющего климата и физико-географических условий. Работа велась по нескольким направлениям: 1. Построены карты изменения годового и максимального стока рек криолитозоны, включая даты нарушения стационарности рассматриваемых рядов. Проанализированы величины изменений и их статистическая значимость с помощью различных критериев. Оценена синхронность колебаний годового и максимального стока. С помощью программы CALM проведен анализ данных стационарных наблюдений за состоянием многолетнемерзлых пород, изменений мощности сезонно-талого слоя (СТС). По данным наблюдений за период 1990–2024 гг. рассчитаны средние за декаду скорости изменения глубины СТС. Собраны полученные для регионов России результаты о динамике метеорологических параметров и глубине СТС. Анализ корреляционных зависимостей выявил существенные пространственные различия влияния климатических факторов на динамику СТС в различных регионах криолитозоны Северного полушария. Для некоторых рек российской криолитозоны рассмотрена зависимость между изменением гидрологических характеристик, метеорологических параметров и глубиной СТС. Для этого выделены две группы рек: азиатской части России, водосборы которых находятся в зоне распространения преимущественно сплошной многолетней мерзлоты, и реки Кольского полуострова, где мерзлота островная. Рассматривались изменения средних значений в 2000–2024 гг. относительно 1979–1999 гг. Был выполнен анализ связи между годовым слоем осадков, речного стока, потенциального испарения и ряда физико-географических характеристик водосборов рек криолитозоны. В качестве источника данных об осадках и потенциальном испарении использовались две крупномасштабные (30 секунд) базы данных – CHELSA 2.1 [Karger et al., 2017] и WorldClim [Fick and Hijmans, 2017]. В обеих базах данных для расчета потенциального испарения (PET) используется уравнение Пенмана-Монтиса. Осадки в CHELSA рассчитываются на основе реанализа ERA5 и базы GPCC. WorldClim использует лишь данные наземных наблюдений за осадками. Были выбраны 172 речных бассейна криолитозоны. Источником информации о физико-географических особенностях речных водосборов послужили базы данных HydroATLAS (https://www.nature.com/articles/s41597-019-0300-6), включающие в себя озерность, лесистость, глубину зеркала грунтовых вод, уклон речной сети, уклон водосбора, среднюю влажность почвы, долю водосбора, занятую многолетнемерзлыми породами, карстом; глобальная километровая база данных (https://essd.copernicus.org/articles/16/2007/2024/essd-16-2007-2024.html) – экспозиция склонов, высота водосбора, стандартное отклонение высот в пределах водосбора, содержание песка и глины в верхнем 10 см слое почвы, листовой индекс, индекс сезонности осадков, число дней со снежным покровом; а также данные архива CHELSA – минимальная температура за год, среднегодовая температура, число случаев перехода температуры через ноль градусов, сумма положительных температур, индекс сезонности климата (рассчитан по данным о сезонном ходе осадков и потенциального испарения). С использованием байесовского подхода была выполнена корректировка данных о многолетнем водном балансе речных водосборов. Был оценен вклад различных особенностей водосборов в изменчивость параметра w уравнения Фу [Zhang et al., 2004], отвечающего за величину испарения при заданных осадках, речном стоке и потенциальном испарении. Масштабы и степень современных изменений генетических составляющих стока рек севера и средней части Сибири и Дальнего Востока в разных физико-географических условиях оценены с помощью автоматического расчленения гидрографа по алгоритму GrWat (Киреева и др., 2019). С использованием этой программы проанализированы показатели, характеризующие все фазы водного режима рек. В рамках исследования была проведена типизация и анализ трансформации гидрографов стока для замыкающих створов 13 рек арктической зоны по годам, а также в среднем за современный и исторический период. Более подробно был изучен бассейн Колымы. Методика типизации заключалась в разделении гидрографов по двум основным признакам: а) превышение максимального расхода относительно среднегодового; б) доля объема стока половодья в общем годовом стоке. По данным признакам гидрографы делились на две группы: с ярко выраженным половодьем (максимальный расход более 50 % обеспеченности за весь период наблюдений, объем половодья более 50% от среднего годового объема стока) и с более сглаженным половодьем. Затем рассчитывалась доля тех и иных случаев за период с 1960 по 1988 и с 1988 по 2017 гг. Дополнительно был проведен анализ максимальных и минимальных расходов воды, дат их прохождения и продолжительности характерных фаз водности, а именно основной из них – весенне-летнего половодья. Был выполнен анализ факторов изменения минимального месячного стока рек криолитозоны. В качестве величины изменения выступала разница минимального месячного стока за период до и после 1980 г. (1940–1980 гг. и 1981–2023 гг.). Для периода с наличием данных наблюдений за речным стоком рассчитывались величины суммы положительных температур за год и слой осадков. Для анализа выбрано 122 гидрологических поста. Помимо основных физико-географических характеристик, аналогичных тем, что использовались для анализа годового стока, в годы с наличием данных наблюдений были рассчитаны величины годового слоя осадков и температуры воздуха по данным реанализа ERA5. Реанализ ERA5 имеет систематическую ошибку, однако в работе он использовался лишь для сравнения слоя осадков и температуры воздуха за два периода. Поскольку для многих рек криолитозоны наблюдается перемерзание, для характеристики изменения минимального стока было выбрано изменение слоя стока, в мм/мес., а не в процентах от базового периода. 2. Для изучения состояния и многолетней динамики растительного покрова на водосборах рек Арктики были использованы спутниковые продукты, полученные по измерениям спутника MODIS. Продукт MOD13C2 (Версия 061) содержит данные о двух растительных индексах NDVI (Normalized Diferential Vegetation Index) и EVI (Enhanced Vegetation Index). Считается, что последний более чувствителен к изменениям свойств растительного покрова в регионах его высокой плотности. Пространственное разрешение продукта 5.6 х 5.6 км, временное - 1 месяц. Период наблюдений составляет 24 года (2000 - 2024 гг.). Изменения растительного покрова в результате пожаров оценивались по данным продукта ESA CCI-L4_FIRE-BA (версии 5.1). В продукте представлены данные по классификации растительного покрова (18 классов) и оценки площадей гарей в каждом из растительных классов. Для исследований использовалась сеточные данные, агрегированные на регулярной основе 0.25х0.25°. Продукт также имеет месячное разрешение, период наблюдений 2001-2022 гг. Для оценки годовых значений гарей, месячные значения были просуммированы. Годовые же значения растительных индексов представляют собой среднемесячные оценки. Значения растительных индексов и оценки площадей гарей были экстрагированы и осреднены для 22 бассейнов крупных и средних рек, расположенных в криолитозоне. В качестве дополнительного источника информации об изменениях состояния наземных экосистем был привлечен спутниковый продукт, содержащий оценки температуры земной поверхности (Land Surface Temperature, LST). Используемый продукт (SSMI_SSMIS_L3C проекта CCI_LST) получен по микроволновым измерениям инструментов SSM/I и SSMIS (спутников серии DMSP), имеет дневное разрешение в интервале лет 1996–2020. Пространственное разрешение продукта 0.25°х0.25°. Для исследования взаимосвязей с водно-балансовыми характеристиками привлекались базы спутниковых данных, собранных на первом этапе выполнения проекта. Для водосборов со значительным количеством гарей, выявленных по спутниковым данным, были проведены численные эксперименты с моделью формирования стока ECOMAG. Для этого в среднем течении р. Лены в подбассейнах р. Вилюй и Алдан были выделены элементарные водосборы, для которых на месте гарей, изначально покрытых таежной растительностью, в модель формирования стока вводились другие типы ландшафтов, также характерные для территории бассейна. Были рассмотрены створы на притоках первого порядка р. Вилюй (р.Тонгуо, Илин-Джели, Чилли, Вилюйчан, Оччугуй-Боттобуя) и контрольный створ на основном русле р. Вилюй – Сунтар, а также притоки в верховьях р. Алдан (верховья рек Амга, Тимптон, Гыным). В расчетах учитывалась максимальная площадь гарей до 20% от площади элементарного водосбора. Средняя площадь каждого элементарного водосбора, где учитывалась смена ландшафтов, составила 5тыс. км2. Была рассмотрена замена таежного леса на открытую незаросшую территорию, на горные луга, на редколесье и на лиственный лес. Расчеты проводились за период 2010-2020 г., в контрольных створах ниже элементарных водосборов с гарями оценивалось изменение модельных гидрографов при исходном задании типов ландшафтов и при их смене на 20% площади элементарных водосборов (для этого редактировался файл stat.bas модели ECOMAG, учитывающий доли тип почв и ландшафтов в каждом элементарном водосборе). Для бассейнов рек Пур, Надым и Колыма в ИМК ECOMAG были подготовлены модельные выдачи о снегозапасах, испарении и суммарных влагозапасах за основной период, соответствующий наличию спутниковых данных - с 2001 по 2020 г. Для сравнения со спутниковыми продуктами об испарении (база данных GLEAM3.6b), снегозапасов (по базе данных Copernicus CCI_Snow SWE), и данных о суммарных влагозапасах спутниковой миссии GRACE проводилось осреднение снегозапасов, испарения и суммарных влагозапасов по бассейнам с месячным шагом. Качество моделирования оценивалось по совпадению смоделированных и полученных по спутниковым данным величин снегозапасов, испарения и суммарных влагозапасов, а также по данным о стоке в ключевых створах р. Надым – Надым, р. Пур – Уренгой, р. Колыма - Колымское с применением критерия качества NSE и систематической ошибки расчетов BIAS. Все модельные расчеты проводились на основе базы данных реанализа W5E5 v2.0, основанного на реанализах ERA-INTERIM и ERA5. По базе данных Copernicus CCI_Snow SWE и результатам моделирования дополнительно сравнивались максимальные снегозапасы. После калибровки и верификации модели для задач дальнейшего анализа компонентов водного баланса и суммарных влагозапасов было проведено моделирование стока исследуемых бассейнов на основе ИМК ECOMAG за весь современный период 1979-2019 г. Для проведения моделирования в ИМК ECOMAG на основе суточных спутниковых данных об осадках (продукт IMERG) были написаны дополнительные скрипты на языке R. Спутниковые данные имеют пространственное разрешение 0.1°х0.1°. Данные о температуре воздуха брались по базе данных реанализа W5E5 v2.0. Было проведено сопоставление сезонного хода осадков по данным реанализа и спутниковым данным IMERG, оценены систематические ошибки в спутниковых осадках при сравнении с реанализом. Моделирование стока проводилось как на основе исходных спутниковых данных, так и с учетом коррекции систематических ошибок в спутниковых осадках отдельно для холодного и теплого периода. Определены тенденции изменения гидрологических показателей (включая оценки компонентов водного баланса, влагозапасов в бассейнах) для ключевых бассейнов арктических рек, полученных на основе спутниковых измерений за общий период наблюдений. Для расчетов линейных трендов наблюденных расходов воды и показателей водного баланса были использованы базы данных натурных и спутниковых наблюдений, собранных на предыдущем этапе выполнения проекта. 3. В 2025 г. проведен детальный анализ ситуации с затоплениями поймы и речными наводнениями на Северо-Востоке азиатской части страны, охватывающем северо-восточные улусы республики Саха (Якутия), Чукотский АО и Магаданскую обл., бассейны рек Яны, Индигирки, Алазеи, Колымы, Анадыря, бассейнов Чукотского и Охотского морей. Здесь наводнения сопровождаются затоплениями в населенных пунктах, аэродромов и портовой инфраструктуры, приводят к размыву автомобильных дорог и обрушению мостов, прорыву дамб хвостохранилищ, нарушают работу горнодобывающих предприятий, увеличивая экологические риски. С 1990-х гг. высокие наводнения были в 1991, 1995-2000, 2004, 2006-2009, 2013, 2014, 2017, 2018, 2019, 2022, 2023, 2025 гг. Были подготовлены несколько информационных массивов: 1) по расходам и уровням воды на 71 постах региона (с 1936 по 2023 гг.), 2) критическим и неблагоприятным отметками, 3) многослойная ГИС региона, 4) База данных детализированных сведений о событиях, связанных с речными наводнениями (91 событие) и неблагоприятными маловодьями (19), 5) 50 спутниковых снимков в период с 2004 г. Проведен детальный статистический, графический, картографический анализ гидрометеорологических условий современного ухудшения ситуации с затоплениями и наводнениями в регионе. Далее на основе созданных баз гидрологических данных подготовлены многолетние ряды максимальных срочных уровней воды с датами (по месяцам и годовые), глубин и продолжительности затопления поймы, случаев превышения отметки ОЯ; для всех постов построены соответствующие хронологические графики, на весь регион – 132 карты соответствующих характеристик за годы с наводнениями (после 1960 г.). Проведен детальный анализ режима максимальных уровней воды и характеристик затоплений. Помимо затоплений рассмотрен вопрос опасных маловодий, связанных прежде всего с низкими расходами и уровнями воды; получены сведения о затруднении навигации в устье Яны, на рр. Колыма и Анадырь. 4. В 2025 г. продолжен анализ истории развития, масштабов и структуры горнодобывающей деятельности в ДФО, ее воздействия на сток наносов рек на основе формируемых баз данных стационарного мониторинга стока взвешенных наносов, ГИС разрабатываемых месторождений и других информационных источников. Новыми объектами исследования стали верхняя и средняя часть бассейна Лены, включая бассейн р.Олёкмы. Для нее подготовлена база данных расходов воды, взвешенных наносов и мутности воды для 57 постов, из которых 33 с наблюдениями за стоком наносов с середины XX в. по 2023 г. Второй информационный массив составила ГИС разрабатываемых месторождений полезных ископаемых, созданная на основе ГИС-Атлас «Недра России», анализа спутниковых снимков, сведений из литературных и интернет источников, связанных с горнодобычей. Всего было проанализировано более 1000 как действующих, так и выведенных из эксплуатации месторождений, прежде всего связанных с золотодобычей. Проанализировано развитие горнодобычи в регионе, выделены основные районы производства (подрайоны внутри них), их специализация и современное состояние. Проведен сравнительный анализ многолетних колебаний стока воды и наносов рек региона, характера связей между ними. В 2025 г. было проведено расширение существующей базы данных о термическом режиме и о стоке рек с привлечением данных по бассейну озера Байкал. Данный бассейн интересен распространением различных типов многолетнемёрзлых пород, большим разнообразием высотных условий и достаточно хорошей гидрологической изученностью. Были собраны данные о среднемесячной и среднедекадной температуре воды на 77 гидрологических постах, данные о месячном и декадном стоке воды на 34 постах с 1950-х гг. По аналогии с нашими более ранними работами был проведён подробный статистический анализ данных, вычислен тепловой сток в озеро Байкал. Расширена база данных о стоке воды в гидрографическом регионе Лена, что позволило проводить более подробный анализ гидрологического режима в целом и теплового стока в частности. На основании собранных данных было выполнено картирование температуры воды методом изолиний для рек Арктики и криолитозоны. Для этого была выработана новая методика построения изолиний на картах термического режима рек. Предложены рекомендации построения мелкомасштабных карт для различных регионов. С использованием реанализа ERA5 была оценена роль температуры воздуха за разные временные периоды (с использованием средних значений и сумм) в формировании средней месячной температуры воды. На основании собственных полевых съемок температуры воды на верхней Колыме уточнена роль температуры воздуха и других географических факторов, разработок россыпных месторождений в формировании температуры воды. На основе расширенной базы данных дана оценка роли новых, построенных в последнее десятилетие водохранилищ в термическом режиме рек, антропогенной трансформации теплового стока рек в вершинах устьевых областей. 5. Долгосрочное прогнозирование характеристик стока рек Арктической зоны на примере реки Надым (сверх плана) Основной целью работы в этом году являлась дальнейшая разработка методики долгосрочного прогнозирования среднемесячного уровня воды заблаговременностью один месяц. В качестве объекта исследования был выбрана река Надым, вдоль берегов которой находятся основные газовые месторождения севера Западной Сибири и являющейся значимой в структуре экономики Надыма. Для исследования были использованы гидрометеорологические данные с разрешением 1 месяц, осредненные по суточным данным и полученные из базы данных Гидрометцентра России и архивов ВНИИГМИ МЦД. В качестве предиктора был выбран среднемесячный уровень воды на гидрологическом посту г. Надым. Для исследования универсальности разработанной методики анализ качества работы модели проводился на том же самом периоде, что и для реки Колыма (с 2011 по 2019 гг.), который включал в себя многоводные и маловодные годы. В процессе разработки моделей нейронных сетей были выделены обучающая (1985-2003 гг.), валидационная (2004-2010 гг.) и тестовая (2011-2019 гг.) выборки. На основе данных 2 метеостанций в пределах бассейна были подготовлены временные ряды месячного разрешения, включающие информацию о температуре и влажности воздуха, скорости ветра, суммарном количестве осадков, а также высоте снежного покрова, а также ряды осредненной по водосбору метеорологической информации за 1960 – 2020 гг. Оценена эффективность пяти моделей - модели инерционного прогноза, модели SARIMAX, рекуррентных нейронных сетей LSTM и GRU, а также модели LS-SVM с различными наборами входных предикторов. Оценка качества моделей проводилась с использованием различных метрик, таких как RMSE и коэффициент детерминации R² и S⁄σΔ . 6. Оценка эмиссии метана с поверхности Колымского водохранилища (сверх плана) На основе результатов натурных измерений концентрации метана в воде и его удельных потоков с водной поверхности впервые проведена оценка эмиссии метана в теплый период года из Колымского водохранилища, водосбор которого полностью находится в зоне вечной мерзлоты. Измерения удельного потока проводились методом плавучих камер, определение содержания метана в пробах – методом headspace [Goldenfum et al., 2010]. В качестве основного инструмента для численного воспроизведения удельных потоков метана из Колымского водохранилища была выбрана модель LAKE 3.2 [Lomov, Stepanenko et al., 2024]. В отчете обобщены материалы сезонных наблюдений в 2021-2023 гг. Выявлена пространственная и сезонная изменчивость как содержания метана в воде, так и его эмиссии. Полный отчет со всеми иллюстрациями и полученными результатами находится в прикрепленном файле. 1. Построены карты изменения годового и максимального стока рек криолитозоны, включая даты нарушения стационарности рассматриваемых рядов. Проанализированы величины изменений и их статистическая значимость с помощью различных критериев. Оценена синхронность колебаний годового и максимального стока. С помощью программы CALM проведен анализ данных стационарных наблюдений за состоянием многолетнемерзлых пород, изменений мощности сезонно-талого слоя (СТС). По данным наблюдений за период 1990–2024 гг. рассчитаны средние за декаду скорости изменения глубины СТС. Собраны полученные для регионов России результаты о динамике метеорологических параметров и глубине СТС. Анализ корреляционных зависимостей выявил существенные пространственные различия влияния климатических факторов на динамику СТС в различных регионах криолитозоны Северного полушария. Для некоторых рек российской криолитозоны рассмотрена зависимость между изменением гидрологических характеристик, метеорологических параметров и глубиной СТС. Для этого выделены две группы рек: азиатской части России, водосборы которых находятся в зоне распространения преимущественно сплошной многолетней мерзлоты, и реки Кольского полуострова, где мерзлота островная. Рассматривались изменения средних значений в 2000–2024 гг. относительно 1979–1999 гг. Был выполнен анализ связи между годовым слоем осадков, речного стока, потенциального испарения и ряда физико-географических характеристик водосборов рек криолитозоны. В качестве источника данных об осадках и потенциальном испарении использовались две крупномасштабные (30 секунд) базы данных – CHELSA 2.1 [Karger et al., 2017] и WorldClim [Fick and Hijmans, 2017]. В обеих базах данных для расчета потенциального испарения (PET) используется уравнение Пенмана-Монтиса. Осадки в CHELSA рассчитываются на основе реанализа ERA5 и базы GPCC. WorldClim использует лишь данные наземных наблюдений за осадками. Были выбраны 172 речных бассейна криолитозоны. Источником информации о физико-географических особенностях речных водосборов послужили базы данных HydroATLAS (https://www.nature.com/articles/s41597-019-0300-6), включающие в себя озерность, лесистость, глубину зеркала грунтовых вод, уклон речной сети, уклон водосбора, среднюю влажность почвы, долю водосбора, занятую многолетнемерзлыми породами, карстом; глобальная километровая база данных (https://essd.copernicus.org/articles/16/2007/2024/essd-16-2007-2024.html) – экспозиция склонов, высота водосбора, стандартное отклонение высот в пределах водосбора, содержание песка и глины в верхнем 10 см слое почвы, листовой индекс, индекс сезонности осадков, число дней со снежным покровом; а также данные архива CHELSA – минимальная температура за год, среднегодовая температура, число случаев перехода температуры через ноль градусов, сумма положительных температур, индекс сезонности климата (рассчитан по данным о сезонном ходе осадков и потенциального испарения). С использованием байесовского подхода была выполнена корректировка данных о многолетнем водном балансе речных водосборов. Был оценен вклад различных особенностей водосборов в изменчивость параметра w уравнения Фу [Zhang et al., 2004], отвечающего за величину испарения при заданных осадках, речном стоке и потенциальном испарении. Масштабы и степень современных изменений генетических составляющих стока рек севера и средней части Сибири и Дальнего Востока в разных физико-географических условиях оценены с помощью автоматического расчленения гидрографа по алгоритму GrWat (Киреева и др., 2019). С использованием этой программы проанализированы показатели, характеризующие все фазы водного режима рек. В рамках исследования была проведена типизация и анализ трансформации гидрографов стока для замыкающих створов 13 рек арктической зоны по годам, а также в среднем за современный и исторический период. Более подробно был изучен бассейн Колымы. Методика типизации заключалась в разделении гидрографов по двум основным признакам: а) превышение максимального расхода относительно среднегодового; б) доля объема стока половодья в общем годовом стоке. По данным признакам гидрографы делились на две группы: с ярко выраженным половодьем (максимальный расход более 50 % обеспеченности за весь период наблюдений, объем половодья более 50% от среднего годового объема стока) и с более сглаженным половодьем. Затем рассчитывалась доля тех и иных случаев за период с 1960 по 1988 и с 1988 по 2017 гг. Дополнительно был проведен анализ максимальных и минимальных расходов воды, дат их прохождения и продолжительности характерных фаз водности, а именно основной из них – весенне-летнего половодья. Был выполнен анализ факторов изменения минимального месячного стока рек криолитозоны. В качестве величины изменения выступала разница минимального месячного стока за период до и после 1980 г. (1940–1980 гг. и 1981–2023 гг.). Для периода с наличием данных наблюдений за речным стоком рассчитывались величины суммы положительных температур за год и слой осадков. Для анализа выбрано 122 гидрологических поста. Помимо основных физико-географических характеристик, аналогичных тем, что использовались для анализа годового стока, в годы с наличием данных наблюдений были рассчитаны величины годового слоя осадков и температуры воздуха по данным реанализа ERA5. Реанализ ERA5 имеет систематическую ошибку, однако в работе он использовался лишь для сравнения слоя осадков и температуры воздуха за два периода. Поскольку для многих рек криолитозоны наблюдается перемерзание, для характеристики изменения минимального стока было выбрано изменение слоя стока, в мм/мес., а не в процентах от базового периода. 2. Для изучения состояния и многолетней динамики растительного покрова на водосборах рек Арктики были использованы спутниковые продукты, полученные по измерениям спутника MODIS. Продукт MOD13C2 (Версия 061) содержит данные о двух растительных индексах NDVI (Normalized Diferential Vegetation Index) и EVI (Enhanced Vegetation Index). Считается, что последний более чувствителен к изменениям свойств растительного покрова в регионах его высокой плотности. Пространственное разрешение продукта 5.6 х 5.6 км, временное - 1 месяц. Период наблюдений составляет 24 года (2000 - 2024 гг.). Изменения растительного покрова в результате пожаров оценивались по данным продукта ESA CCI-L4_FIRE-BA (версии 5.1). В продукте представлены данные по классификации растительного покрова (18 классов) и оценки площадей гарей в каждом из растительных классов. Для исследований использовалась сеточные данные, агрегированные на регулярной основе 0.25х0.25°. Продукт также имеет месячное разрешение, период наблюдений 2001-2022 гг. Для оценки годовых значений гарей, месячные значения были просуммированы. Годовые же значения растительных индексов представляют собой среднемесячные оценки. Значения растительных индексов и оценки площадей гарей были экстрагированы и осреднены для 22 бассейнов крупных и средних рек, расположенных в криолитозоне. В качестве дополнительного источника информации об изменениях состояния наземных экосистем был привлечен спутниковый продукт, содержащий оценки температуры земной поверхности (Land Surface Temperature, LST). Используемый продукт (SSMI_SSMIS_L3C проекта CCI_LST) получен по микроволновым измерениям инструментов SSM/I и SSMIS (спутников серии DMSP), имеет дневное разрешение в интервале лет 1996–2020. Пространственное разрешение продукта 0.25°х0.25°. Для исследования взаимосвязей с водно-балансовыми характеристиками привлекались базы спутниковых данных, собранных на первом этапе выполнения проекта. Для водосборов со значительным количеством гарей, выявленных по спутниковым данным, были проведены численные эксперименты с моделью формирования стока ECOMAG. Для этого в среднем течении р. Лены в подбассейнах р. Вилюй и Алдан были выделены элементарные водосборы, для которых на месте гарей, изначально покрытых таежной растительностью, в модель формирования стока вводились другие типы ландшафтов, также характерные для территории бассейна. Были рассмотрены створы на притоках первого порядка р. Вилюй (р.Тонгуо, Илин-Джели, Чилли, Вилюйчан, Оччугуй-Боттобуя) и контрольный створ на основном русле р. Вилюй – Сунтар, а также притоки в верховьях р. Алдан (верховья рек Амга, Тимптон, Гыным). В расчетах учитывалась максимальная площадь гарей до 20% от площади элементарного водосбора. Средняя площадь каждого элементарного водосбора, где учитывалась смена ландшафтов, составила 5тыс. км2. Была рассмотрена замена таежного леса на открытую незаросшую территорию, на горные луга, на редколесье и на лиственный лес. Расчеты проводились за период 2010-2020 г., в контрольных створах ниже элементарных водосборов с гарями оценивалось изменение модельных гидрографов при исходном задании типов ландшафтов и при их смене на 20% площади элементарных водосборов (для этого редактировался файл stat.bas модели ECOMAG, учитывающий доли тип почв и ландшафтов в каждом элементарном водосборе). Для бассейнов рек Пур, Надым и Колыма в ИМК ECOMAG были подготовлены модельные выдачи о снегозапасах, испарении и суммарных влагозапасах за основной период, соответствующий наличию спутниковых данных - с 2001 по 2020 г. Для сравнения со спутниковыми продуктами об испарении (база данных GLEAM3.6b), снегозапасов (по базе данных Copernicus CCI_Snow SWE), и данных о суммарных влагозапасах спутниковой миссии GRACE проводилось осреднение снегозапасов, испарения и суммарных влагозапасов по бассейнам с месячным шагом. Качество моделирования оценивалось по совпадению смоделированных и полученных по спутниковым данным величин снегозапасов, испарения и суммарных влагозапасов, а также по данным о стоке в ключевых створах р. Надым – Надым, р. Пур – Уренгой, р. Колыма - Колымское с применением критерия качества NSE и систематической ошибки расчетов BIAS. Все модельные расчеты проводились на основе базы данных реанализа W5E5 v2.0, основанного на реанализах ERA-INTERIM и ERA5. По базе данных Copernicus CCI_Snow SWE и результатам моделирования дополнительно сравнивались максимальные снегозапасы. После калибровки и верификации модели для задач дальнейшего анализа компонентов водного баланса и суммарных влагозапасов было проведено моделирование стока исследуемых бассейнов на основе ИМК ECOMAG за весь современный период 1979-2019 г. Для проведения моделирования в ИМК ECOMAG на основе суточных спутниковых данных об осадках (продукт IMERG) были написаны дополнительные скрипты на языке R. Спутниковые данные имеют пространственное разрешение 0.1°х0.1°. Данные о температуре воздуха брались по базе данных реанализа W5E5 v2.0. Было проведено сопоставление сезонного хода осадков по данным реанализа и спутниковым данным IMERG, оценены систематические ошибки в спутниковых осадках при сравнении с реанализом. Моделирование стока проводилось как на основе исходных спутниковых данных, так и с учетом коррекции систематических ошибок в спутниковых осадках отдельно для холодного и теплого периода. Определены тенденции изменения гидрологических показателей (включая оценки компонентов водного баланса, влагозапасов в бассейнах) для ключевых бассейнов арктических рек, полученных на основе спутниковых измерений за общий период наблюдений. Для расчетов линейных трендов наблюденных расходов воды и показателей водного баланса были использованы базы данных натурных и спутниковых наблюдений, собранных на предыдущем этапе выполнения проекта. 3. В 2025 г. проведен детальный анализ ситуации с затоплениями поймы и речными наводнениями на Северо-Востоке азиатской части страны, охватывающем северо-восточные улусы республики Саха (Якутия), Чукотский АО и Магаданскую обл., бассейны рек Яны, Индигирки, Алазеи, Колымы, Анадыря, бассейнов Чукотского и Охотского морей. Здесь наводнения сопровождаются затоплениями в населенных пунктах, аэродромов и портовой инфраструктуры, приводят к размыву автомобильных дорог и обрушению мостов, прорыву дамб хвостохранилищ, нарушают работу горнодобывающих предприятий, увеличивая экологические риски. С 1990-х гг. высокие наводнения были в 1991, 1995-2000, 2004, 2006-2009, 2013, 2014, 2017, 2018, 2019, 2022, 2023, 2025 гг. Были подготовлены несколько информационных массивов: 1) по расходам и уровням воды на 71 постах региона (с 1936 по 2023 гг.), 2) критическим и неблагоприятным отметками, 3) многослойная ГИС региона, 4) База данных детализированных сведений о событиях, связанных с речными наводнениями (91 событие) и неблагоприятными маловодьями (19), 5) 50 спутниковых снимков в период с 2004 г. Проведен детальный статистический, графический, картографический анализ гидрометеорологических условий современного ухудшения ситуации с затоплениями и наводнениями в регионе. Далее на основе созданных баз гидрологических данных подготовлены многолетние ряды максимальных срочных уровней воды с датами (по месяцам и годовые), глубин и продолжительности затопления поймы, случаев превышения отметки ОЯ; для всех постов построены соответствующие хронологические графики, на весь регион – 132 карты соответствующих характеристик за годы с наводнениями (после 1960 г.). Проведен детальный анализ режима максимальных уровней воды и характеристик затоплений. Помимо затоплений рассмотрен вопрос опасных маловодий, связанных прежде всего с низкими расходами и уровнями воды; получены сведения о затруднении навигации в устье Яны, на рр. Колыма и Анадырь. 4. В 2025 г. продолжен анализ истории развития, масштабов и структуры горнодобывающей деятельности в ДФО, ее воздействия на сток наносов рек на основе формируемых баз данных стационарного мониторинга стока взвешенных наносов, ГИС разрабатываемых месторождений и других информационных источников. Новыми объектами исследования стали верхняя и средняя часть бассейна Лены, включая бассейн р.Олёкмы. Для нее подготовлена база данных расходов воды, взвешенных наносов и мутности воды для 57 постов, из которых 33 с наблюдениями за стоком наносов с середины XX в. по 2023 г. Второй информационный массив составила ГИС разрабатываемых месторождений полезных ископаемых, созданная на основе ГИС-Атлас «Недра России», анализа спутниковых снимков, сведений из литературных и интернет источников, связанных с горнодобычей. Всего было проанализировано более 1000 как действующих, так и выведенных из эксплуатации месторождений, прежде всего связанных с золотодобычей. Проанализировано развитие горнодобычи в регионе, выделены основные районы производства (подрайоны внутри них), их специализация и современное состояние. Проведен сравнительный анализ многолетних колебаний стока воды и наносов рек региона, характера связей между ними. В 2025 г. было проведено расширение существующей базы данных о термическом режиме и о стоке рек с привлечением данных по бассейну озера Байкал. Данный бассейн интересен распространением различных типов многолетнемёрзлых пород, большим разнообразием высотных условий и достаточно хорошей гидрологической изученностью. Были собраны данные о среднемесячной и среднедекадной температуре воды на 77 гидрологических постах, данные о месячном и декадном стоке воды на 34 постах с 1950-х гг. По аналогии с нашими более ранними работами был проведён подробный статистический анализ данных, вычислен тепловой сток в озеро Байкал. Расширена база данных о стоке воды в гидрографическом регионе Лена, что позволило проводить более подробный анализ гидрологического режима в целом и теплового стока в частности. На основании собранных данных было выполнено картирование температуры воды методом изолиний для рек Арктики и криолитозоны. Для этого была выработана новая методика построения изолиний на картах термического режима рек. Предложены рекомендации построения мелкомасштабных карт для различных регионов. С использованием реанализа ERA5 была оценена роль температуры воздуха за разные временные периоды (с использованием средних значений и сумм) в формировании средней месячной температуры воды. На основании собственных полевых съемок температуры воды на верхней Колыме уточнена роль температуры воздуха и других географических факторов, разработок россыпных месторождений в формировании температуры воды. На основе расширенной базы данных дана оценка роли новых, построенных в последнее десятилетие водохранилищ в термическом режиме рек, антропогенной трансформации теплового стока рек в вершинах устьевых областей. 5. Долгосрочное прогнозирование характеристик стока рек Арктической зоны на примере реки Надым (сверх плана) Основной целью работы в этом году являлась дальнейшая разработка методики долгосрочного прогнозирования среднемесячного уровня воды заблаговременностью один месяц. В качестве объекта исследования был выбрана река Надым, вдоль берегов которой находятся основные газовые месторождения севера Западной Сибири и являющейся значимой в структуре экономики Надыма. Для исследования были использованы гидрометеорологические данные с разрешением 1 месяц, осредненные по суточным данным и полученные из базы данных Гидрометцентра России и архивов ВНИИГМИ МЦД. В качестве предиктора был выбран среднемесячный уровень воды на гидрологическом посту г. Надым. Для исследования универсальности разработанной методики анализ качества работы модели проводился на том же самом периоде, что и для реки Колыма (с 2011 по 2019 гг.), который включал в себя многоводные и маловодные годы. В процессе разработки моделей нейронных сетей были выделены обучающая (1985-2003 гг.), валидационная (2004-2010 гг.) и тестовая (2011-2019 гг.) выборки. На основе данных 2 метеостанций в пределах бассейна были подготовлены временные ряды месячного разрешения, включающие информацию о температуре и влажности воздуха, скорости ветра, суммарном количестве осадков, а также высоте снежного покрова, а также ряды осредненной по водосбору метеорологической информации за 1960 – 2020 гг. Оценена эффективность пяти моделей - модели инерционного прогноза, модели SARIMAX, рекуррентных нейронных сетей LSTM и GRU, а также модели LS-SVM с различными наборами входных предикторов. Оценка качества моделей проводилась с использованием различных метрик, таких как RMSE и коэффициент детерминации R² и S⁄σΔ . 6. Оценка эмиссии метана с поверхности Колымского водохранилища (сверх плана) На основе результатов натурных измерений концентрации метана в воде и его удельных потоков с водной поверхности впервые проведена оценка эмиссии метана в теплый период года из Колымского водохранилища, водосбор которого полностью находится в зоне вечной мерзлоты. Измерения удельного потока проводились методом плавучих камер, определение содержания метана в пробах – методом headspace [Goldenfum et al., 2010]. В качестве основного инструмента для численного воспроизведения удельных потоков метана из Колымского водохранилища была выбрана модель LAKE 3.2 [Lomov, Stepanenko et al., 2024]. В отчете обобщены материалы сезонных наблюдений в 2021-2023 гг. Выявлена пространственная и сезонная изменчивость как содержания метана в воде, так и его эмиссии.
3 1 января 2026 г.-31 декабря 2026 г. Гидрологические последствия изменения климата и антропогенного воздействия в криолитозоне
Результаты этапа: -

Прикрепленные к НИР результаты

Для прикрепления результата сначала выберете тип результата (статьи, книги, ...). После чего введите несколько символов в поле поиска прикрепляемого результата, затем выберете один из предложенных и нажмите кнопку "Добавить".