Санкт-Петербургский государственный университет
Сезонные закономерности распределения осадков на территории Западного Саяна и их учет при моделировании стока
Г.В. Пряхина, Е.С. Зелепукина
На основе анализа метеоданных с суточным разрешением за период с 1970 по 1985 г. установлено, что на территории Западного Саяна, слабо обеспеченной данными режимных гидрометеорологических наблюдений, более тесные связи наблюдаются для годовых сумм осадков, а менее тесные -- для сезонных сумм. Получены зависимости изменения осадков теплого и холодного периодов для семигумидных и семиаридных районов исследуемой территории. Установлено, что на южном семиаридном макросклоне наибольшую величину достоверности имеет линейная корреляция: вне зависимости от высотного интервала величина плювиометрического градиента остается постоянной. Для северного семигумидного макросклона изменения сумм осадков с высотой определяются логарифмической зависимостью: для низкогорья (до 500 м) характерно быстрое нарастание количества осадков (градиент теплого периода достигает 100 мм / 100 м), а выше 1300 м градиенты и теплого и холодного периодов практически не превышают 2 мм / 100 м. Достоверность выявленных зависимостей была подтверждена модельными расчетами. Объекты и периоды моделирования выбирались с учетом наличия как метеорологических, так и гидрологических данных (водосборные бассейны рек Амыл и Ус для семигумидного и семиаридного районов соответственно). Было установлено, что применение сезонных плювиометрических градиентов заметно улучшает сходимость рассчитанных и наблюденных суточных гидрографов при математическом моделировании стока. Показано, что для семиаридных районов учет пространственного распределения осадков оказался более значимым. Введение высотных градиентов наиболее существенно улучшает сходимость результатов моделирования в годы с наименьшими значениями критерия Нэша -- Сатклифа и снижает количество лет с неудовлетворительным качеством расчетов, что особенно важно при моделировании стока в условиях недостаточности гидрологической информации.
Ключевые слова: пространственная интерполяция осадков, вертикальные плювио- метрические градиенты, моделирование стока.
Seasonal distribution patterns of precipitation on the territory of West Sayan Ridge and their consideration in runoff modeling
Specific height dependences of the precipitation based on the daily meteodata analysis (for the period 1970-1985) were obtained for Western Sayan, which is under conditions of insufficient data provided by the scarce network of stations. Analysis revealed that the annual precipitation sum is characterized by closer dependence than seasonal precipitation. It has been generally observed that the linear correlation has the highest reliability for the southern macroslope of the ridge. At the same time the dependence of the northern slope is primarily determined by height: distribution of precipitation is defined by logarithmic function. Thus, the low mountain pluviometric gradient (below 500 m) is characterized by a quick increase in precipitation (in the warm period the gradient reaches a value of 100 mm/100 m). As for the mid-mountain parts (above 1300 m), gradients of warm and cold periods do not exceed 2 mm/100 m. Reliability of identified dependences has been confirmed by model calculations. All the modelled objects and periods were selected considering the existence and availability of both hydrological and meteorological data -- the river Us and Amyl basins (right tributaries of Yenisei, Krasnoyarsk district, Russia) which represent typical runoff conditions for the semi-arid and semi-humid areas of the Western Sayan. The investigation also found that the application of seasonal pluviometric gradients during runoff modeling has improved the convergence of calculated and obtained hydrographs. It was shown that taking into account vertical pluviometric gradients turned out to be more significant for semi-arid regions. And furthermore, it has improved the quality of runoff modeling notably for periods characterized by the least value of Nash-Sutcliffe criteria, and reduced the amount of years with improper quality of calculation. This is especially important under conditions of insufficient hydrometeorological data.
Keywords: spatial interpolation of the values of meteorological elements, pluviometric gradients, runoff modeling.
Введение
Горные районы характеризуются большей по сравнению с равнинными территориями интенсивностью протекания гидрологических, геоморфологических и других процессов. Широкий диапазон высот, значительные уклоны, различная экспозиция склонов, формы и ориентация речных долин, а также сложная ландшафтная структура приводят к значительной пространственной неоднородности в распределении метеорологических элементов, что, в свою очередь, определяет особенности режима стока с горных водосборов.
Несмотря на многочисленные исследования по выявлению закономерностей распределения осадков с высотой, проводимые для различных горных стран, при решении разнообразных практических задач проблема определения региональных значений плювиометрических градиентов остается весьма актуальной. Это связано в первую очередь с низкой плотностью существующей сети метеостанций и постов, невысокой репрезентативностью данных из-за расположения станций преимущественно в долинах и наличия пропусков в рядах наблюдений, а также с определяющим влиянием гипсометрического и орографического факторов на пространственную изменчивость поля осадков и пр. Значения плювиометрических градиентов зависят не только от высоты, но и от сезона, вида осадков, экспозиции склонов, длины ряда наблюдений и временных интервалов, по которым определяется градиент (Marchand, 1986; Smadja, 1991; Martinez del Castillo et al., 2012; Sanchez Martin, 1995).
В качестве объекта исследования был выбран среднегорный хребет субширотного простирания Западный Саян, расположенный между более высокими поднятиями Горного Алтая и Восточного Саяна. Несмотря на сравнительно малую высоту и отсутствие развитого гляциально-нивального пояса, хребет представляет собой важный ороклиматический барьер: наветренные северные и северо-западные склоны характеризуются относительно высокой для континентального сектора влагообеспеченностью, а подветренные южные и юго-восточные -- недостаточной.
Выявление высотных зависимостей распределения осадков в Саянах затруднено по причине отсутствия инструментальных наблюдений в верхних частях горных хребтов в отличие от территории Алтая, значительно лучше обеспеченной метеорологическими данными, в том числе и в высокогорье. В этой ситуации можно ориентироваться лишь на имеющиеся данные снегомерных съемок (Геткер и Жданов, 1992) и ряд работ по расчету значений весенних осадков с суточным разрешением (Бураков и Гордеев, 2013). В связи с тем, что при моделировании гидрографов стока необходимы зависимости, учитывающие пространственное изменение осадков во все периоды года, целью настоящего исследования стало выявление региональных зависимостей распределения осадков по высоте на территории среднегорий Западного Саяна для различных сезонов. Проверка корректности полученных закономерностей была выполнена путем расчета суточных гидрографов стока.
Материалы и методы
Методика выявления региональных особенностей изменения осадков с высотой, представленная в работе, основана на корреляционном анализе сетевых данных о количестве осадков с суточным и месячным разрешением за 1970-1985 гг. по 18 метеостанциям и постам, расположенным на высотах примерно от 200 до 1400 м (рис. 1). Выбор указанного временного интервала обусловлен отсутствием пропусков в рядах данных на максимальном количестве имеющихся в регионе станций и постов; кроме того, период 1961-2013 гг. не отличается значимым трендом в межгодовом ходе осадков в Алтае-Саянском регионе (Чередько, 2015).
сезонный осадки моделирование сток
Рис. 1. Расположение гидрологических и метеорологических станций в районе исследования. Зависимости изменения осадков от высоты: I -- для I группы метеостанций, II -- для II группы; а -- холодный период, Ь -- теплый период, с -- год
Для выделения групп станций с более тесными корреляционными связями был применен кластерный анализ. В качестве правила объединения использовалась одиночная связь, в качестве меры расстояния -- величина й = 1 - г, где г -- коэффициент линейной корреляции между суммами осадков на станциях. Поскольку предварительный анализ распределения осадков с высотой на изучаемой территории выявил различия для холодного и теплого периодов (Пряхина и др., 2017), все расчеты проводились отдельно для сумм осадков за теплый (апрель-октябрь), холодный (ноябрь-март) периоды и за год. Из-за пропусков в рядах данных на станциях Нижне-Усинское, Буйба и Арадан расчеты производились по данным шести лет наблюдений.
Для каждого года определялись уравнения аппроксимирующих кривых изменения количества осадков с высотой, по которым вначале было рассчитано количество осадков для стометровых высотных интервалов, а затем -- отношения сумм осадков на определенной высоте к сумме осадков на станциях (Mi). Уравнения аппроксимирующих кривых изменения Mi с высотой характеризуют зависимости распределения осадков в различные периоды года.
Проверка адекватности предложенного подхода к интерполяции осадков осуществлялась путем сравнения наблюденных и смоделированных гидрографов стока. Для этого использовалась моделирующая система с распределенными параметрами «Гидрограф», разработанная Ю.Б.Виноградовым (Виноградов и Виноградова, 2010) и показавшая свое успешное применение в различных физико-географических условиях Hydrograph Model Website. URL: http://www.hydrograph-model.com (дата обращения: 10.09.2018). (Vinogradov et al., 2011).
Моделирование проводилось для двух случаев:
1. Без учета вертикальных градиентов. При этом для интерполяции суточных осадков по площади водосбора использовался метод линейной триангуляции, при котором для каждой расчетной точки (РТ) гексагональной сетки, равномерно покрывающей площадь бассейна, использовались данные трех ближайших метеостанций, причем вес каждой метеостанции учитывался пропорционально расстоянию до РТ.
2. С учетом вертикальных градиентов. В этом случае для расчета сумм суточных осадков в РТ количество наблюденных осадков на метеостанции умножалось на соответствующий поправочный коэффициент:
где a, b -- коэффициенты уравнений, аппроксимирующих кривых изменения Mi с высотой; n -- номер РТ; Hn -- высота расчетной точки, м; Wi -- весовой коэффициент метеостанции, определяемый по уравнению:
где i -- порядковый номер метеостанции; Di,n -- расстояние между соответствующей метеостанцией и РТ, м; m -- число метеостанций в группе.
В качестве критерия сравнения был взят критерий эффективности Нэша -- Сатклифа (NS), используемый Американской ассоциацией гражданских инженеров, а также Всемирной метеорологической организацией (Shamseldin and O'Connor, 2001).
Результаты и обсуждение
В ходе корреляционного анализа выявлено, что на исследуемой территории более тесные связи наблюдаются для годовых сумм осадков, менее тесные -- для сумм за теплый период: статистически значимые коэффициенты корреляции изменяются от 0,95 до 0,60, соответственно. Более высокие корреляционные связи между годовыми осадками по сравнению с сезонными подтверждают известный факт, что с увеличением периода осреднения корреляция возрастает (Каган, 1979). Уменьшение корреляционных связей с увеличением расстояния между станциями особенно заметно в теплый период, когда возрастает вероятность образования локальных облачных систем и осадки носят очаговый характер.
С помощью кластерного анализа в северо-восточном предгорье Западного Са- яна была выделена группа станций с наиболее близкими связями: Минусинск, Ку- рагино, Ермаковское, Каратуз, Казыр. С этой группой были агрегированы станции Верхний Амыл и Оленья речка, так как для них максимальные значения d в теплый период не превышают 0,20. Вторая группа объединяет станции, расположенные южнее главного водораздела Западного Саяна в диапазоне высот от 660 м до 1404 м: Нижне-Усинское, Арадан, Буйба, Оленья речка. Данные некоторых станций (Абакан, Верхний Кужебар, Черепановка, Усть-Шадат и др.) не использовались в расчетах в связи со статистически незначимыми связями рядов осадков с данными остальных станций. Таким образом, группу I составили станции, отражающие закономерности распределения осадков семигумидной области Западного Саяна, группу II -- семиаридной (см. рис. 1). Отметим, что граница, разделяющая области с разными зависимостями распределения осадков, не совпадает с осевой частью хребта, а смещена в подветренную сторону вследствие перераспределения твердых атмосферных осадков за счет метелевого переноса (Ревякин, 1981).
Результаты расчетов показали, что во все годы и для всех периодов для первой группы станций большую величину достоверности аппроксимации имеют уравнения, представляющие логарифмическую зависимость (коэффициент детерминации изменяется от 0,94 до 0,97), а для второй -- линейную (коэффициент детерминации изменяется от 0,93 до 0,99) (7е1ерикта е! а1., 2017), что соотносится с зависимостями, используемыми в прогнозах стока (Балабанова и Заборцева, 1984).
Для семигумидной области быстрое увеличение градиента осадков прослеживается до 500 м (см. рис. 1, I): градиент теплого периода в предгорьях составляет более 100 мм / 100 м, холодного -- около 40 мм / 100 м; более плавное нарастание наблюдается примерно до 1300 м, а выше -- градиенты и теплого и холодного периодов составляют в среднем всего 2 мм / 100 м. Подобная закономерность распределения летних осадков на Алтае (в долине р. Актру) наблюдается до высоты 2500 м, после чего скорости увеличения количества осадков резко снижаются вплоть до достижения нулевой границы градиента (Галахов и др., 1987). Учитывая, что в Саянах выше 1400 м данные режимных наблюдений отсутствуют, полагаем, что полученные зависимости можно обоснованно использовать только до уровня верхней границы леса, поскольку по данным снегомерных съемок (Авдеева и Бураков, 2003) выше 1650 м в подгольцовом поясе имеет место резкое уменьшение снегозапасов вследствие усиленного перевевания и выдувания снега.
В семиаридных частях хребта (см. рис. 1, II) полученные градиенты теплого периода составляют около 100 мм / 100 м, а холодного -- 30 мм / 100 м. Выявленную линейную зависимость можно использовать примерно до верхней границы леса, которая на южном склоне прослеживается до уровня 1800-1850 м (ОаулШпа апё 7е1ерикта, 2017); выше -- в подгольцовом поясе -- скорее всего, происходит снижение значения градиента аналогично отмеченному на северном семигумидном макросклоне.
Следует добавить, что максимальная величина отклонения рассчитанных значений сумм осадков от наблюденных составила для некоторых станций для теплого периода 14 %, а для холодного -- 25 %. Это может быть связано с неодновременностью выпадения осадков в предгорной и горной частях водосбора, что часто приводит к завышению рассчитанных с помощью плювиометрических градиентов сумм осадков в горах при интенсивных ливнях на равнине в период прохождения атмосферных фронтов (Гордеев, 2012). Значительные различия градиентов жидких и твердых атмосферных осадков в сочетании с пространственно-временной дискретностью выпадения осадков в разных частях водосбора определяют целесообразность использования при моделировании стока именно сезонных зависимостей вместо осредненных годовых.