Моделирование процессов формирования стока рек горной криолитозоны северо-востока России с наледным питанием на примере р. Анмангында

Землянскова Анастасия Александровна ORCID: 0000-0002-2249-5708 лаборант, Северо-Восточный государственный университет; младший научный сотрудник, Санкт-Петербургский государственный университет 685017, Россия, Магаданская область, г. Магадан, ул. Портовая, 13 Zemlianskova Anastasiia Aleksandrovna Laboratory assistant, Laboratory of Monitoring and Forecasting Climate Change and the Environment, North-Eastern State University, Magadan, Russia; Junior Researcher and Postgraduate Student of the Institute of Earth Sciences, St. Petersburg University 685017, Russia, Magadan region, Magadan, Portovaya str., 13 anastasiazemlanskova@gmail.com Другие публикации этого автора

Макарьева Ольга Михайловна ORCID: 0000-0002-2532-4306 кандидат технических наук заведующий лабораторией, Северо-Восточный государственный университет; ведущий научный сотрудник, Санкт-Петербургский государственный университет 685000, Россия, Магаданская область, г. Магадан, ул. Портовая, 13 Makarieva Olga Mikhailovna PhD in Technical Science Head of Laboratory of Monitoring and Forecasting Climate Change and the Environment, North-Eastern State University, Magadan, Russia; Leading Researcher of Earth Sciences, St. Petersburg University 685000, Russia, Magadan region, Magadan, Portovaya str., 13 omakarieva@yandex.ru Другие публикации этого автора

Полевые исследования на Анмангындинской наледи в 2020-2023 гг. проведены при поддержке Санкт-Петербургского государственного университета в рамках проекта «Комплексная оценка естественных и антропогенных факторов интенсификации водообменных процессов криолитозоны в условиях изменения климата», обработка и анализ полевых данных, а также разработка блока «Наледь» выполнены в рамках проекта Российского научного фонда и Правительства Магаданской области No 23-17-20011 «Оценка рисков опасных явлений для инфраструктуры Магаданской области, вызванных климатически- и антропогенно-обусловленной трансформацией многолетнемерзлых пород».

Введение

Наледь — специфическая форма сезонного оледенения, характерная для горных условий многолетней мерзлоты ^[1-3]. Площадь одного объекта может измеряться десятками км², а толщина льда достигать 10—12 м ^[4]. В условиях глобального потепления климата наблюдается сокращение размеров наледей ^{[5, 6]} и их миграция, что оказывает влияние на гидрологический режим рек криолитозоны и геокриологические условия ^{[3, 7]}. Крупные размеры наледей обеспечивают высокую идентифицируемость этих объектов на космических снимках, что позволяет использовать их в виде качественных индикаторов трансформации водообмена на обширных, труднодоступных и недостаточно изученных территориях Северо-Востока России, где по последним оценкам насчитывается не менее 6683 наледей с общей площадью около 4500 км^{2 [8]}.

Наледи являются ключевым элементом в цепи водообменных процессов зоны распространения многолетнемерзлых пород. Согласно наблюдениям Б. Л. Соколова, в мае-июне до 50—90% наледного стока попадает в речной сток, в июле эта доля снижается до 30%, а в августе — до 7%. Гидрологическая роль наледей в формировании речного стока может быть сопоставима с такими снежно-ледовыми явлениями, как ледники, что подтверждается результатами многих исследований ^{[6, 9, 10]}. Наблюдения за наледями при строительстве Байкало-Амурской магистрали показали, что в годовом стоке рек доля наледного стока может достигать 35% для водосборов площадью до 500 км². В период весеннего половодья эта доля составляет от 5 до 40%, а в летнем стоке не превышает 10% ^[11].

Наледи могут играть важную роль в водообеспечении засушливых регионов ^{[2, 12]}. Несмотря на историю изучения этого природного явления (более 100 лет) и достигнутые результаты, отсутствуют методы прогноза развития наледных процессов, а также гидрологические модели, учитывающие долю наледного питания в речном стоке. На данный момент реализуются исследования по созданию концептуальных моделей образования наледей ^{[13, 14]}, что позволит лучше понять их природу. В известных гидрологических моделях, таких как HBV-Light model ^{[15, 16]}; PREVAH ^[17] и др. с помощью системы уравнений удовлетворительно описываются разрушение ледника и ледниковый сток, их можно также использовать в отношении наледей с помощью регулирования снежного индекса, но это теряет всякий физический смысл ^[18].

Проблема разработки математического аппарата, позволяющего прогнозировать размеры и характер наледных процессов, а также учитывать их влияние на речной сток и другие компоненты природной среды кроется в том, что отсутствует государственная программа изучения наледей. Наземные наблюдения в России выполняются благодаря грантовым поддержкам [19-21 и др.] отдельными научными группами, результаты которых часто не сопоставимы друг с другом. Необходима комплексная программа изучения наледей на ключевых объектах в разных регионах, которая будет учитывать особенность климата, геологического строения, рельефа, геокриологических характеристик, ландшафтов на формирование наледей. Данные спутникового мониторинга позволят дать ретроспективную оценку изменения их размеров ^{[5, 22]}. В прошлом веке на наледях рр. Самокит, Верина, Анмангында и др. проводились комплексные наблюдения за их водным режимом и ролью наледей в формировании речного стока. Поэтому именно эти объекты и предлагается рассмотреть в первую очередь в качестве ключевых объектов сети наледеметрических пунктов.

В бассейне р. Анмангында (бассейн р. Колымы) был получен длительный ряд данных (1962–1991 гг.) о площади, объеме, толщине гигантской наледи в разные сезоны года, а также гидрологическом режиме реки в створе ниже наледи (1962–1987) и на крупном притоке руч. Ольчан (1962–1978 гг.). На основе этих материалов были получены интегральные кривые уменьшения морфометрических характеристик наледи, которые применялись в расчетах на других объектах ^{[10, 23, 24]}.

В данной работе на примере бассейна р. Анмангында предлагается изучить гидрологический режим, а также, используя длительный ряд данных с помощью методов математического моделирования, разработать метод учета наледного питания.

Объект исследования

Площадь водосбора р. Анмангында (впадает в р. Детрин, приток р. Колыма), для которого проводится расчет, составляет 376 км². Территория исследования относится к зоне сплошного распространения многолетней мерзлоты, ее мощность достигает 300 м, уменьшается до 80–100 м в долинах рек и прерывается в таликовых зонах ^[25]. Геологическое строение представлено осадочными породами верхнего триаса, юры и мела. По данным ^[26] глубина нулевых теплооборотов составляет 11,5 и 13,5 м на высоте 712 и 1027 м соответственно.

Высота водосбора р. Анмангында составляет от 700 до 1850 м. Долина реки на участке образования наледи (рис. 1) имеет трапецеидальную форму с пологими склонами. Протяженность наледной поляны составляет около 6—7 км, с абсолютной высотой 700—770 м. Склоны покрыты лиственничным лесом на мохово-лишайниковой подстилке, местами кедровый стланик и пятна открытой тундры. Почвы горнолесные подзолистые, тундровые и горнотундровые.

Рис. 1. Анмангындинская наледь.

Дно долины сложено аллювием, в верхней части он слабо прикрыт кустарничком и моховым покровом. На участках, где лед лежит продолжительное время, образуется кочкарник, а также густой ивняк. Наледь образуется ежегодно, в современный период ее площадь может достигать 5,4 км2 (2000 г.), а исторический максимум составил 6,8 км² (1978 г.). Средний объем наледи за период 1963-1991 гг. был 8,5 млн м³, что составляет 22,6 мм годового слоя стока. Современные оценки показывают значительное уменьшение общих размеров Анмангындинской наледи ^[27], так в 2000-2022 гг. максимальный объем по расчетным оценкам не превышал 8,2 млн м³, а в среднем составляет 7 млн м³.

Водный режим реки Анмангында характеризуется данными гидрологического поста, расположенного ниже наледи. На участке полыньи выше наледи и в средней части наледной поляны наблюдается круглогодичный сток.

Климат исследуемого района резко континентальной, средняя годовая температура на станции Усть-Омчуг (576 м над ур. моря) за период 1967-2022 гг. составила –9,8ºС (в июле +14,5ºС, в январе -33,4ºС). В самые холодные месяцы (декабрь-январь) температура опускается до -57℃, а в теплые поднимается до +33,6℃ (июль). Среднегодовая сумма осадков составляет 341,2 мм (м/с Усть-Омчуг).

Снежный покров в горах устанавливается в конце сентября, а по всему бассейну в середине октября. По данным полевых исследований 2021-2023 гг. высота снега в разных частях водосбора на конец марта составляет от 48 до 84 см, тогда как на м/с Усть-Омчуг — 25 см.

Материалы и методы исследования

Исторические данные

Материалы Колымского управления гидрометеослужбы были собраны, оцифрованы и проанализированы за период 1962-1991 гг. Они включали в себя информацию о динамике площади, объема, толщины льда в разные сезоны года. Подробно анализ и результаты описаны в работах ^{[27, 28]}. Отмечено, что за последние 60 лет максимальные размеры наледи сократились на 15 и 17% для площади и объема соответственно (в среднем с 5,5 до 4,7 км², с 8,5 до 7,0 млн м³).

Информация о гидрологическом режиме получена для нижнего створа реки Анмангында (1962-1987) ^[29].

Современный ряд наблюдений

Оценка современной площади наледи выполнялась на основе анализа космических снимков спутников Landsat (датчики TM, ETM+ и OLI) и Sentinel-2 (датчик MSI) за период 1999-2022 гг. и 2016-2022 гг., соответственно. Точность определения площади льда по снимкам зависит от его разрешения, отклонение от фактического размера наледи может составлять до 15%.

Данные космических снимков не позволяют оценить объем льда, поэтому используется расчетная формула (1):

W=0,364*S^1,167, (1)

где W — объем льда (тыс. м³), S — площадь льда, тыс. м², коэффициенты a и n были предложены в работе ^[27]. Для исторического периода (1963–1990 гг.) формула показала удовлетворительные результаты. Отклонение расчетных от наблюденных величин составило 0,01—2,09 млн м³ или 0—5 мм стока.

С помощью беспилотного летательного аппарата в 2021 и 2022 гг. была выполнена съемка наледной поляны. По результатам работы были построены ортофотопланы, определена площадь и объем наледи ^[30], а также изменение толщины льда. Точность определения объема зависит от многих факторов — стабилизации камеры, количества контрольных точек, освещения и др. Сравнение наземных наблюдений с данными БПЛА показало, что ошибка толщины льда в теплый период года составила от 0,04 до 1,70 м. В случае грубых ошибок — материалы в анализе не использовались. Так, точность определения объема варьировалась от 0,02 до 0,46 млн м³ (5—13%), что составляет от 0 до 1,5 мм стока для бассейна р. Анмангында.

Модель «Гидрограф»

Математическая модель «Гидрограф» с распределёнными параметрами комбинирует в себе как физически обоснованные, так и концептуальные методы, что позволяет достичь баланса между сложностью расчетов и доступностью параметризации модели ^[31]. Она применялась для расчета стока рек бассейна р. Колымы, где распространены многолетнемерзлые породы ^{[32, 33]}, и результаты были приняты удовлетворительными. В расчетах не проводилась оценка вклада наледного питания в сток рек, хотя на территории Северо-Востока страны он может составлять 0,2-11,0 % в период весеннего половодья (май-июнь, бассейн р. Колыма) ^[8].

В ранее разработанной версии модели «Гидрограф» при расчете процессов формирования стока в зоне распространения многолетнемерзлых пород не учитывалась доля вклада наледного питания в сток рек. Авторами работы был предложен расчетный блок «Наледь», описывающий процессы разрушения наледей и учета наледного питания в формировании речного стока. В рамках адаптации алгоритма и программы модели были модифицированы как расчетные схемы, так и интерфейс.

Для использования модели «Гидрограф» проводится схематизация бассейна реки. Используются две системы — 1) набор репрезентативных (расчетных) точек (РТ), представляющих собой центры ячеек гексагональной сетки, и 2) система стокоформирующих комплексов (СФК), определяемых типом ландшафта, почвы, растительности, высотной поясностью.

Блок «Наледь» представляет собой алгоритм расчета для массива наледей, приуроченного к конкретной РТ. Для каждой наледи в пределах РТ-площади задается информация о ее характеристиках и местоположении – географические координаты (широта, долгота, высота), максимальная площадь (S_max) и объем (W_max) наледи (определяются по данным полевых исследований или с помощью расчетных формул), средняя толщина льда (D_max), как частное от объема и площади льда. Входные метеорологические данные, а именно эффективные температура (Temp_ef) и дефицит влажности (Def_ef) воздуха, рассчитываются с учетом абсолютной высоты наледи.

За основу блока «Наледь» был взят метод Б. Л. Соколова, который основан на динамике сокращения площади наледи в течение теплого периода, начало которого связывается с датой перехода средней суточной температуры воздуха через нуль в сторону положительных значений ^[10]. Объем льда меняется в зависимости от коэффициента таяния. При разработке метода Б. Л. Соколов использовал ряд данных о динамике наледных явлений по Анмангындинской наледи, которые в настоящий момент не удовлетворяют современным климатическим условиям. Поэтому, по данным спутникового мониторинга за период (2000-2022 гг.) была уточнена региональная кривая разрушения относительной площади Анмангындинской наледи.

Для расчета наледного стока учитывается два фактора, влияющих на разрушение наледи – действие солнечной радиации (Abl_i) и термоэрозионные процессы, когда под воздействием стока речных и подземных вод блоки льда откалываются от основного тела наледи и быстро (в предложенном алгоритме – мгновенно) разрушаются в потоке воды. Такое предложено впервые, ранее блоковое разрушение либо не учитывалось, либо предлагалось решение нескольких уравнений ^[10]. Коэффициент таяния (k_t) был рассчитан по данным полевых исследований 2020-2023 гг. на Анмангындинской наледи, установлено, что для этого объекта он в среднем составляет 1 мм/℃·сутки и варьируется от 0,1 до 4,2 мм/℃·сутки в зависимости от месяца. В модели он задается единым на весь расчетный интервал времени. Также вводится коэффициент испарения со льда (м/мбар·с). Минимальные размеры наледи, для которой проводится расчет, составляют 0,6 км².

Расчетный интервал времени (РИВ) в модели «Гидрограф» может быть любым (но не более суток) в зависимости от наличия данных. В этой работе предлагается интервал сутки, так на каждые сутки определяется как изменились площадь и объем наледи в результате действия солнечной радиации и термоэрозионного разрушения с учетом испарения с поверхности льда. Толщина льда уменьшается равномерно по всей площади наледи в зависимости от эффективной температуры воздуха и коэффициента таяния льда (блок солнечная радиация). Объем отколовшегося льда рассчитывается, как изменение площади, умноженное на среднюю толщину наледи (блок термомеханическое разрушение наледи). Суммарная потеря объема наледи на РИВ пересчитывается в слой и гидрограф наледного стока. Также на каждый РИВ определяются площадь, объем и средняя толщина наледи, а также слой, объем и расход наледного стока, слой испарения со льда, которые используются для дальнейшего расчета. Ниже приведены формулы, используемые в блоке «Наледь» (рис. 2).

Рис. 2. Блок-схема «Наледь».

1) Солнечная радиация — рассчитывается изменение средней толщины наледи за РИВ в результате прихода солнечной энергии.

Если (если слой таяния под действием солнечной радиации превышает толщину льда).

2) Термомеханическое разрушение наледи — рассчитывается изменение площади и объёма наледи за счет блокового разрушения.

3) Баланс — рассчитывается суммарный объем потери льда в i-й РИВ за счет испарения, блокового разрушения и таяния под действием солнечной радиации.

Апробация метода расчета площади и объема наледи на каждый РИВ (рис. 3) показала, что отклонение полученных значений в среднем не превышает 10% от максимальных характеристик льда, а наибольшее отклонение для конкретного дня составило 24% (табл. 1), что принято удовлетворительным. Максимальная ошибка расчета наблюдается в период июль—август, когда размеры наледи составляют в среднем менее 10% от значений площади и объема до начала таяния.

Рис. 3. Апробация метода расчета площади и объема наледи по данным наземных наблюдений (1) и материалам беспилотного летательного аппарата (3). Кривая (2) — результат применения блока «Наледь».

Таблица 1. Отклонение расчетных значений характеристик наледи от наблюденных за каждый год.

© 1998 – 2024 Nota Bene. Publishing Technologies. NB-Media Ltd.