Scientific journal

Advances in current natural sciences

ISSN 1681-7494

"Перечень" ВАК

ИФ РИНЦ = 0,775

Климатические изменения, связанные с увеличением температуры воздуха, наиболее заметны на территории Арктической зоны Российской Федерации. Увеличение температуры воздуха за холодный сезон года способствовало перераспределению вкладов отдельных факторов на формирование тех или иных фаз ледового режима, что в свою очередь привело к неопределенности в вопросах оценки и прогнозирования характеристик ледового режима. Неблагоприятные проявления ледового режима нередко приводят к возникновению опасных гидрологических явлений. Это не только наводнения в период ледохода или затора льда, но и продолжительные осенние и весенние переходные периоды, когда невозможна организация ледовых или паромных переправ и иная хозяйственная деятельность, функционирование которой зависит от гидрологической обстановки. Минимизация рисков и урона, наносимого опасными гидрологическими явлениями, зависит от оценки существующих рисков и их изменений во времени, а также развития методов прогнозирования и оповещения населения.

Целью исследования является комплексное описание ледового режима реки Печоры от села Усть-Цильма до села Оксино и развитие методических подходов к прогнозированию высших уровней воды.

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

Исследуемая область относится к Арктической зоне Российской Федерации. В пределах исследуемой области на реке Печоре расположено три гидрологических поста: с. Усть-Цильма, с. Ермицы и с. Оксино, а также две метеорологических станции, расположенные непосредственно в бассейне реки – ГМС Нарьян-Мар и ГМС Усть-Цильма.

Период наблюдений на данных гидрологических постах и метеорологических станциях является достаточным для расчета основных гидрологических и метеорологических характеристик. Так как наиболее полные наблюдения за характеристиками ледового режима начались с 1950 г., то общий период непрерывных наблюдений составляет 70 лет, с 1950–1951 гидрологического года по 2019–2020 гидрологический год (расчётное гидрологическое полугодие принято с октября одного года по май следующего года). Разбиение всего периода наблюдений на базовый и контрольный выполнено с учетом нарушения стационарности рядов наблюдений, таким образом, базовый период выбран с 1950 по 1980 г., контрольный с 1981 по 2020 г.

Однородность рядов оценена при помощи критериев Стьюдента и Фишера. Для оценки линейных трендов использовался критерий значимости выборочного коэффициента корреляции (R). Гипотеза об отсутствии тренда не опровергалась при уровне значимости 2α = 10 %, если выполнялось условие |R| < 2σR. Критическое значение коэффициента корреляции при данной длительности рядов наблюдений составляет Rкрит ≈ 0,23. Приведение рядов метеорологической и гидрологической информации к многолетнему периоду выполнено методом построения регрессионной зависимости за совместный период наблюдений с привлечением данных наблюдений пунктов-аналогов в однородном районе при условии R > 0,70 [1].

Результаты исследования и их обсуждение

Выявлено статистически значимое увеличение средней за холодный сезон температуры воздуха на всех рассмотренных метеостанциях по отношению к базовому периоду. Наибольший вклад в общее потепление климата вносят октябрь, март и май, так как именно за эти месяцы отмечаются статистически значимые тренды и отмечается нарушение однородности рядов по критерию Стьюдента. Среднее увеличение температуры воздуха за холодный период составляет +1,3 °С, что на 0,2 °С больше, чем в среднем по Двинско-Печорскому бассейновому округу [2]. Средняя многолетняя температура воздуха за холодный сезон составляет -8,9 °С (табл. 1).

Условия ледообразования и формирования основных фаз ледового режима на Печоре являются типичными для крупных незарегулированных водотоков равнинного типа. Для них характерно ледообразование по всей глубине потока. Замерзание сопровождается образованием шуги и движением ее вниз по течению. Ледовый покров образуется при остановке, сплочении и смерзании льда и шуги [3].

Комплексный анализ сумм температур воздуха, уровней воды и типа ледовых явлений позволил выявить, что первичные ледовые явления независимо от уровня воды образуются в диапазоне сумм отрицательных температур воздуха от 0 до -31 °С, устойчивый ледостав образуется в диапазоне отрицательных температур воздуха от -43 до -186 °С и имеет относительно тесную обратную зависимость с уровнем воды (R = -0,62). Деградация и разрушение ледового покрова начинается с переходом температуры воздуха через 0 °С, высший уровень воды за период весеннего ледохода формируется в диапазоне сумм положительных температур воздуха от 32 до 180 °С, связи с уровнем воды за этот период выражены слабо.

Анализ линейных трендов показал статистически значимое уменьшение длительности ледостава за счет смещения сроков замерзания на более поздние и сроков вскрытия на более ранние. Общее снижение длительности ледостава составляет от 12 до 4 дней в селах Усть-Цильма и Оксино соответственно. Похожие оценки были получены ранее в работах [4, 5], но за более короткий временной промежуток.

Несмотря на высокую длительность весеннего ледохода, на больших реках волна половодья имеет сложный распластанный во времени вид, что определяет формирование высших уровней за год как при ледоходе, так и в условиях свободного русла. По данным наблюдений высшие уровни за год формировались за период с ледовыми явлениями лишь в 19 % случаев на участке реки в пределах Усть-Цильмы, в 47 % случаев на участке реки в пределах с. Ермицы и 71 % случаев на участке реки в пределах с. Оксино. Таким образом, с увеличением площади водосбора в пределах рассмотренной территории, вероятность формирования высшего уровня воды за год при наличии ледовых явлений повышается, что объясняется более поздними сроками окончания ледохода в устьевых участках.

Наиболее затороопасным участком из рассмотренных является участок реки близ с. Ермицы, где повторяемость заторов льда достигает 53 %. Среднее значение высшего за год уровня воды заторного происхождения 990 см, для случаев, когда высший уровень за год формировался при ледоходе, среднее значение составило 993 см, на открытой воде – 957 см. При этом стоит отметить что высший уровень за год в 16 % случаев формируется при ледоходе, в 31 % случаев при заторе и в 53 % при свободном ото льда русле (табл. 2). Таким образом, можно заключить, что высший уровень за год примерно в половине случаев формируется на открытой воде, а в половине – при тех или иных ледовых явлениях, при этом наибольших значений высшие уровни воды достигают при ледоходе или заторе льда. Повторяемость заторов льда близ ГП Усть-Цильма составляет 7 %, при этом формируются ординарные уровни воды, ГП Оксино – 21 %. Повторяемость заторов льда на ГП Оксино существенно увеличилась в последние годы наблюдений, так, с 2004 г. отмечается 83 % от всех случаев заторообразования, при этом формируются высшие за год уровни воды.

Снижение фактических рисков для населения зависит от своевременности и качества выпускаемых прогнозов. Можно отметить, что в настоящее время для прогнозирования высших уровней воды используют различные эмпирические методики. Качество выпускаемых прогнозов при этом зависит от актуальности эмпирической информации, используемой при их разработке. Принцип выбора предикторов для прогнозирования высшего уровня ледохода на различных участках рек зависит от степени стохастичности условий его формирования. Выявлено, что степень стохастичности условий формирования высших уровней за период весеннего ледохода обуславливается географическими факторами (озерность, заболоченность, залесенность, площадь водосбора, приточность, затороопасность) и может быть оценена по зависимости высшего уровня воды за период весеннего ледохода от высшего предледоходного уровня воды. При этом стоит отметить, что высокая стохастичность формирования высшего уровня воды за период весеннего ледохода предопределяет высокую стохастичность формирования и других характеристик уровенного режима. В общем случае можно выделить три типа участков рек. Первый тип – зарегулированные озерами участки реки, отдельные бесприточные участки, где высшие уровни формируются с большой степенью инертности. Ко второму типу относятся крупные незарегулированные реки со множеством притоков, которые способны значительным образом повлиять на ледовый режим и характер вскрытия. Третий тип является промежуточным, к нему в первую очередь относятся предустьевые области рек. Гидрографические, морфометрические и физико-географические особенности реки Печоры обуславливают исключительную стохастичность формирования ее уровенного режима, и только участок реки близ с. Оксино можно отнести к третьему типу. Это означает, что для участков реки у с. Усть-Цильма и Ермицы необходим поиск предикторов вне данных постов, так как связь между высшим уровнем ледохода и высшим предледоходным уровнем воды выражена слабо. Оценку же качества выпускаемых прогнозов необходимо выполнять по отношению стандартной ошибки прогноза к стандартному отклонению прогнозируемой величины (S/σHср), а не по отношению к стандартному отклонению среднего изменения прогнозируемой величины за период заблаговременности (S/σHср), так как данный критерий наиболее строгий. Так, по данным работы [7] критерий S/σHср почти в два раза больше, критерия S/σ∆H. Для с. Оксино необходимо использовать инерционные методы с добавлением второстепенных предикторов, а оценку качества выпускаемых прогнозов производить по критерию S/σ∆H. Выбор оптимальных предикторов для прогнозирования высшего уровня за период весеннего ледохода был выполнен в работе [2], при этом улучшение качества выпускаемых прогнозов может быть достигнуто при использовании метода обучения искусственных нейронных сетей.

При разработке модели прогнозирования высшего уровня воды за период весеннего ледохода с учетом высших уровней воды заторного генезиса в с. Оксино на реке Печоре в качестве предикторов использованы высший ледоходный уровень воды в селе Ермицы и высший предледоходный уровень на посту в с. Оксино, средняя заблаговременность прогноза составила 4 дня. Рассмотрим абсолютные ошибки результатов прогнозирования высшего уровня воды за период весеннего ледохода при использовании регрессионных и нейросетевых прогностических моделей (рис. 2).

Можно отметить, что лучшим качеством выпускаемых прогнозов характеризуется нейронная сеть MLP 2-8-1, оправдываемость прогнозов которой составила 92 %, оправдываемость прогнозов, данных по уравнению множественно регрессии, составила 79 %. Для MLP 2-8-1 отношение S/σ∆ на независимом материале составило 0,46, что свидетельствует о высоком качестве выпускаемых прогнозов.

Заключение

В работе рассмотрены вопросы изменения климата и его влияния на ледовый режим реки Печоры. Установлено увеличение средней за холодный сезон температуры воздуха в районе 1,3 °С, что привело к снижению длительности ледостава от 4 до 12 дней. По данным многолетних наблюдений повторяемость заторов льда близ ГП Усть-Цильма составляет 7 %, при этом формируются ординарные уровни воды, ГП Оксино – 21 %.

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