Изменения режима рек вследствие климатических, а также нередко и антропогенных причин происходят во многих регионах мира, в частности на территории РФ. У динамики различных гидрологических характеристик водотоков территории России нет единых закономерностей. Так, например, для европейской части страны прослежена тенденция небольшого увеличения годового стока; в то же время установлено снижение максимального стока для рек бассейнов нижней Волги, Дона, Днепра, южной части бассейна Оки, для большинства из них значительное. Сток рек Сибири и Дальнего Востока, впадающих в северные моря, растет, сток рек восточнее р. Колыма, включая реки полуострова Камчатка, рек южных частей водосборов Оби, Енисея и Лены, крупных притоков р. Селенга снижается. Рост годовых расходов воды отмечен для рек южного Приморья [1]. Минимальный сток (зимний и периода открытого русла) значимо увеличился на большинстве створов в бассейнах р. Ангары и Лены [2]. Снизился годовой сток одного из главных притоков трансграничной р. Амур – р. Сунгари. Отмечено снижение после 1985 г. и стока р. Туманная, второй по площади водосбора после р. Амур трансграничной реки на территории Дальнего Востока [3]. Часто наряду с изменениями режима рек отмечается повышение частоты катастрофических гидрологических событий. Предполагается, что в ближайшие десятилетия рост повторяемости максимальных паводков согласно модельным прогнозам можно ожидать на реках Сибири и Дальнего Востока [1].
Исток р. Туманная находится на восточном склоне г. Чанбайшань, общая длина реки составляет 549 км, площадь водосбора 33,8 тыс. км2 [4], основные притоки на китайской стороне – р. Гайя и Хуньчунь. На российской территории площадью 26 км2 находится всего 17 км устьевого участка реки, по фарватеру которой проходит государственная граница между РФ и КНДР. Основная часть бассейна реки (70 %) расположена в Яньбэнь-Корейском автономном округе провинции Цзилинь КНР, 30 % – в провинции Хамген КНДР. Среднегодовой сток реки по данным китайской станции Цюаньхэ (площадь водосбора 31 800 км2), расположенной на расстоянии 40 км от устья, равен 215 м3/с [5].
Российская часть бассейна реки составляет около 0,1 %, но она испытывает серьезные последствия от наводнений, формирующихся выше по течению на территории КНР и КНДР. На пограничном участке реки во время крупных паводков происходят значительные русловые деформации, разрушение российского берега, что приводит к прямым территориальным потерям РФ. На отдельных участках русла с 1950 г. уже произошло смещение берега на 250–300 м. В 2008 г. для предотвращения сдвига фарватера реки были построены берегоукрепительные сооружения длиной более 13 км.
Рис. 1. Российская часть бассейна р. Туманная и прилегающие территории: 1 – граница бассейна реки; 2 – российская часть бассейна; 3 – гидрологические посты наблюдений
В связи с этим актуально проанализировать гидрологический режим реки, его особенности, приводящие к негативному воздействию вод. Цель исследования заключается в оценке многолетних колебаний различных характеристик уровенного режима р. Туманная и их изменений. До настоящего времени подобные оценки для российской части бассейна реки не проводились.
Материалы и методы исследования
Гидрологический режим р. Туманная в устьевой части (рис. 1) проанализирован на основе данных об уровнях воды (среднегодовых, минимальных и максимальных среднесуточных в году, среднемесячных), так как российский пост наблюдений только уровенный. Он находится на западной окраине железнодорожной станции Хасан. На участке поста долина реки резко расширяется, пойма имеет ширину 6 км, затапливается при уровне 700 см, отметка опасного явления – 750 см. Поскольку до 15 октября 1959 г. створ находился в 200 м выше, использованы данные об уровнях воды за 1960–2020 гг.
Ряды данных наблюдений не являются непрерывными, слабо коррелированы между собой, за исключением среднемесячных уровней ноября – декабря и января – марта, автокорреляция среднегодовых значений уровней в многолетнем ряду слабая. Для зимних месяцев пропуски в рядах среднемесячных значений уровней заполнены с использованием метода линейной регрессии. Для остальных месяцев, а также для средних за год значений для этой цели использовались осредненные величины уровней воды, а также интерполяция. Проверка рядов наблюдений на однородность при уровне значимости 0,05 проводилась с использованием широко применяемых статистических критериев Фишера, Стьюдента и Манна – Уитни [6–8], на значимость трендов – критерия Стьюдента по отношению к коэффициенту регрессии. Общий ряд наблюдений, с учетом данных, изложенных в [3], был разбит на подпериоды: 1960–1985 и 1986–2020 гг.
Результаты исследования и их обсуждение
Диапазон колебаний среднегодовых и минимальных в году зимних уровней воды р. Туманная у п. Хасан в целом за 1960–2020 гг. не очень значителен: 212–301 см и 162–265 см соответственно. Наибольшим он был у многолетних колебаний максимальных в году среднесуточных уровней (299–863 см), в среднем они превышали меженные зимние почти в 3 раза (2,8), в отдельные годы это превышение достигало 3,5–4 раз. Наибольшее превышение максимального уровня (863 см, 2016 г.) над уровнем зимней межени (162 см, 1968 г.) за весь анализируемый период составило более 7 м.
Анализ характера чередования многоводных и маловодных лет в многолетней динамике среднегодовых уровней воды р. Туманная, проведенный с помощью коэффициентов аномалий их значений (отношение значения уровня i-го года к среднему за период), показал, что выделяются два периода с четкими различиями в величинах уровней (рис. 2, а). После 1986 г. фиксируется маловодный период, характеризующийся преобладанием значений годовых уровней ниже среднемноголетних величин (71 %). До 1986 г. включительно в ряду годовых уровней преобладали значения выше среднемноголетних, их количество составляет 65 %. Рассчитанные коэффициенты аномалий максимальных в году уровней свидетельствуют, что с 2010 г. наблюдается почти непрерывная группировка лет с их повышенными значениями, превышающими среднемноголетние (553 см) на 12–56 % в зависимости от года. Исключение составляют 2011 и 2014 гг., когда значения максимальных уровней были ниже средних за рассматриваемый период (рис. 2, б).
Визуальный анализ графического отображения многолетней динамики данных характеристик дает основания предположить, что ряды их значений не являются однородными (рис. 3, а и б).
Проверка многолетних рядов годовых значений, среднесуточных минимальных и максимальных в году уровней воды на однородность по критериям Фишера и Стьюдента показала, что неоднородными по дисперсии являются ряды значений среднегодовых и минимальных в году уровней, а неоднородными по среднему являются ряды всех трех характеристик. В 1986–2020 гг. (второй интервал разбиения исходного ряда) размах колебаний среднегодовых уровней незначительно увеличился по отношению к значениям за 1960–1985 гг., минимальных и максимальных в году уровней снизился. Средние значения годовых и минимальных в году уровней снизились с 264 до 243 и с 209 до 188 (на 8 и 10 % соответственно); значения максимальных в году уровней увеличились с 505 до 589 см (на 17 %).
а) среднегодовые уровни
б) максимальные в году уровни
Рис. 2. Коэффициенты аномалий значений уровней воды р. Туманная (ст. Хасан)
Таблица 1
Многолетние характеристики уровней воды по месяцам, см
Месяц |
I |
II |
III |
IV |
V |
VI |
VII |
VIII |
IX |
X |
XI |
XII |
среднее (I) |
217 |
217 |
225 |
245 |
260 |
293 |
301 |
336 |
315 |
258 |
244 |
241 |
среднее (II) |
201 |
200 |
204 |
222 |
251 |
265 |
288 |
309 |
292 |
229 |
214 |
207 |
снижение, см (%) |
16 (7) |
17 (8) |
21 (9) |
23 (9) |
9 (3) |
28 (10) |
13 (4) |
27 (8) |
23 (7) |
29 (11) |
30 (12) |
34 (14) |
SD (I) |
26 |
28 |
25 |
25 |
33 |
60 |
52 |
61 |
59 |
29 |
29 |
34 |
SD (II) |
35 |
30 |
23 |
25 |
44 |
47 |
55 |
74 |
74 |
39 |
45 |
40 |
изменение, см/% |
9 (35) |
2 (7) |
-2 (-8) |
0 (0) |
11 (33) |
-13 (-22) |
3 (6) |
13 (21) |
15 (25) |
10 (34) |
16 (55) |
6 (18) |
Примечание. I, II – подпериоды (1960–1985, 1986–2020 гг.); SD – стандартное отклонение.
Использование критерия Манна – Уитни не подтвердило неоднородности рядов значений только максимальных в году уровней воды. Достоверным по критерию Стьюдента является линейный тренд на снижение для годовых уровней.
Режим уровней внутри года по месяцам. Многолетние колебания среднемесячных уровней значительно различаются. Коэффициенты вариации рядов значений изменяются от 0,11 (март, апрель) до 0,21–0,23 (август, сентябрь). Рассчитанные уравнения линейного тренда показывают, что для всех месяцев, за исключением мая, тренд отрицательный. Наиболее всего среднемесячные уровни воды снизились в октябре – декабре (более чем на 10 %) (табл. 1).
Проверка рядов среднемесячных значений уровней воды по критериям Стьюдента (с учетом разности дисперсий) и Манна – Уитни показала, что неоднородными по среднему значению являются ряды уровней в марте, апреле и октябре – декабре, для этих месяцев результаты применения критериев совпадают. Значимость тренда по критерию Стьюдента подтверждена только для апреля, июня и октября (табл. 2). Наибольшими темпами значения уровней воды в р. Туманная за период с 1960 по 2020 г. снижались в июне и августе, в среднем на 0,82 и 0,96 см в год (рис. 4).
Опасность наводнений. За период с 1960 по 2020 г. вода р. Туманная выходила на пойму 11 раз, при этом отметки опасного явления (750 см) за рассматриваемый период были превышены в 1965, 1986, 2000, 2016, 2019 и 2020 гг. на величину от 28 до 113 см, т.е. половина случаев наводнений на российской части бассейна наблюдалась в период с 2010 по 2020 г. Максимальное за 61-летний период наблюдений превышение имело место в 2016 г. во время тайфуна Лайонрок, когда уровень воды достигал 863 см (обеспеченность 1,64 %) и был размыт берег реки на ширину от 3 до 10 м, построенное берегоукрепление частично разрушено.
а) среднегодовые (1) и минимальные в году (2) уровни воды
б) максимальные в году уровни воды
Рис. 3. Характеристики гидрологического режима р. Туманная, 1960–2020 гг.
Таблица 2
Оценка однородности и значимости трендов рядов наблюдений за уровнями воды по месяцам воды на р. Туманная – п. Хасан
Месяц |
I |
II |
III |
IV |
V |
VI |
VII |
VIII |
IX |
X |
XI |
XII |
Оценка однородности |
||||||||||||
Критерий Фишера (p = 0,05) |
О |
О |
О |
Н |
Н |
О |
Н |
Н |
Н |
Н |
О |
Н |
Критерий Стью-дента (p = 0,05) |
О |
Н |
Н |
Н |
О |
Н |
О |
О |
О |
Н |
Н |
Н |
Критерий Манна – Уитни |
О |
О |
Н |
Н |
О |
О |
О |
О |
О |
Н |
Н |
Н |
Оценка значи-мости трендов |
– |
– |
– |
+ |
– |
+ |
– |
– |
– |
+ |
– |
– |
Примечание. О – однородный; Н – неоднородный; «+» – тренд значимый; «-» – незначимый.
а) июнь, 1960–2020 гг.
б) август, 1960–2020 гг.
Рис. 4. Многолетняя динамика уровней воды р. Туманная – п. Хасан по отдельным месяцам
Заключение
Полученные оценки хорошо согласуются с данными о снижении водности р. Туманная, опубликованными ранее китайскими исследователями, но касающимися только годового стока (по данным станции наблюдений Цюаньхэ за 1959–2011 гг.). В результате проведенной работы выявлены изменения различных характеристик гидрологического режима реки, общей закономерностью которых является отрицательный линейный тренд рядов почти всех проанализированных характеристик уровней с темпами их снижения в целом за 1960–2020 гг. от 0,1 до 1 см в год. Исключение составляют максимальные в году уровни, увеличивавшиеся в течение 61 года согласно рассчитанному уравнению регрессии на 2 см в год, и уровни мая с условным ростом 0,05 см в год. Достоверность тренда подтверждена по критерию Стьюдента для средних уровней за год, а также за апрель и октябрь. Размах межгодовых колебаний характеристик, напротив, увеличился, за исключением максимальных и минимальных в году суточных уровней и средних уровней марта и июня. Половина случаев превышения уровней отметки опасного явления (3 из 6) наблюдалась в 2010–2020 гг., что свидетельствует об увеличении их повторяемости. Исторический максимум уровня наблюдался в 2016 г. во время тайфуна Лайонрок, когда во время наводнения при отметке уровня в 863 см (обеспеченность 1,64 %) берегоукрепление на российском берегу реки было частично разрушено. Таким образом, максимальные в году суточные уровни – это единственная характеристика, значения которой увеличивались с середины 1980-х гг., при этом выросла частота случаев превышения уровней отметки опасного явления. Все это свидетельствует о возможности в дальнейшем усиления катастрофичности опасных гидрологических событий.
Для подобного рода оценок важно определение предположительной «точки перелома», года, после которого изменяется тенденция динамики гидрологического режима. Практически все исследованные характеристики имели локальные пиковые значения в 1986 г., что расходится с данными по ст. Цюаньхэ, возможно потому, что из-за отсутствия в 1985 г. наблюдений на ст. Хасан пропущенные значения были заменены средними.
Возможно, на формирование неоднородности исследуемых характеристик уровенного режима р. Туманная у п. Хасан помимо климатического оказывает влияние и антропогенный фактор, т.е. характер использования территории водосбора реки, а также большое количество построенных на китайской части водохранилищ различной емкости. Зарегулирование притоков реки привело к значительному снижению размаха колебаний минимальных в году (зимних) уровней воды, в меньшей степени оказало влияние на межгодовые вариации максимальных в году их значений из-за небольших размеров водохранилищ и, возможно, отсутствия в них противопаводочной емкости.
Полученные результаты дают определенное понимание закономерностей изменений гидрологического режима р. Туманная и могут быть использованы для прогноза различных характеристик водности реки, а также для обоснования проектов реконструкции защитных берегоукрепительных сооружений.