Physical-mathematical modeling of multiyear dynamics of water-balance and snow-reserve components in Ob-Irtysh river basin

Full Text

ВВЕДЕНИЕ

Настоящая статья — продолжение серии статей, посвященных моделированию формирования составляющих водного баланса (прежде всего, стока) на территориях бассейнов северных рек РФ [2, 5, 6, 8, 9, 12, 21, 22]. Актуальность данной тематики обоснована в указанной серии статей и связана как с созданием более точного инструментария решения гидрологических задач, так и с разработкой методики сценарных прогнозов изменения составляющих водного баланса речных бассейнов северного склона РФ [1, 7, 10].

Территория северных регионов РФ, согласно большинству сценариев изменения климата, подвергнется самым большим изменениям. Этот регион привлекает все больше внимания, поскольку как в самом Северном Ледовитом океане, так и на его шельфе имеются большие запасы ископаемого топлива. Кроме того, указанный регион характеризуется огромным разнообразием водных ресурсов, включая самые большие реки мира. Поэтому вопрос о гидрологическом режиме бассейнов северных рек в современных условиях и о его изменении в связи с возможными изменениями климата крайне важен и требует разработки более точных методов оценки многолетней динамики водных ресурсов северных рек РФ в современный период и в предстоящие ближайшие десятилетия.

В указанных выше работах авторов настоящей статьи используется методика моделирования динамики составляющих водных ресурсов северных рек (как в настоящее время, так и в будущем), основанная на разработанной авторами физико-математической модели взаимодействия подстилающей поверхности суши с атмосферой SWAP (Soil−Water−Atmosphere−Plants) [2–4]. Данная методика продемонстрировала хорошее качество воспроизведения стока северных рек РФ и других показателей гидрологического режима их бассейнов для современного периода, что позволяет рассчитывать на адекватное воспроизведение динамики составляющих водного баланса в ближайшие десятилетия при наличии прогнозов динамики метеорологических элементов.

Цель настоящей работы — проверить возможность воспроизведения динамики составляющих водного баланса бассейна р. Оби — одной из пяти величайших рек земного шара, одной из крупнейших рек бассейна Северного Ледовитого океана и всей России, а также динамики снегозапасов в пределах территории бассейна на основе модели взаимодействия поверхности суши с атмосферой SWAP, которая уже использована для гидрологических расчетов для ряда крупных рек северного склона РФ с перспективой сценарного прогнозирования изменения составляющих водного баланса бассейнов северных рек в связи с возможными изменениями климата.

ФИЗИКО-ГЕОГРАФИЧЕСКОЕ ОПИСАНИЕ БАССЕЙНА р. ОБИ

В основу описания физико-географических характеристик бассейна р. Оби положены данные, приведенные в [15, 16].

Обь берет свое начало на Алтае, длина ее от слияния рек Бии и Катуни составляет 3660 км, а от истока Иртыша (ее основного притока) — 5410 км. Площадь Обь-Иртышского бассейна (рис. 1), самого крупного речного бассейна России — 2990 тыс. км². Среднемноголетний объем речного стока р. Оби — 403 км³/год, что соответствует слою стока в 134 мм/год. Река пересекает с юга на север территорию Западной Сибири и на всем своем протяжении, кроме истока, Обь представляет собой типично равнинную реку с малыми уклонами и широкой заболоченной долиной, достигающей местами ширины несколько десятков километров. Впадая в Обскую губу Карского моря, река образует дельту с многочисленными рукавами и островами.

Главным притоком р. Оби является р. Иртыш. Длина р. Иртыша составляет 4248 км, что превышает длину самой Оби. Иртыш — самая длинная река-приток в мире. Протекает по территории Китая, Казахстана и России. Истоки Иртыша находятся на границе Монголии и Китая, на восточных склонах хребта Монгольский Алтай.

 

Рис. 1. Территория Обь-Иртышского бассейна.

 

Обь с Иртышом пересекают несколько природно-климатических зон: от полупустыни на юге до тундры на севере. Большая часть бассейна покрыта тайгой и болотами. В зоне тундры распространена многолетняя мерзлота.

Климат Западной Сибири умеренно континентальный. Зимой и летом в южную часть бассейна поступает арктический воздух. Весной и осенью наблюдается активное проникновение циклонов, приносящих влажный воздух с Атлантики. Континентальность климата возрастает в основном в направлении с севера на юг. Средняя температура января понижается до –15°С в юго-западной части бассейна и до –30°С в его северо-восточной части. Средняя температура июля достигает 5°С на севере бассейна и 20°С в его южной части.

Согласно оценкам, приведенным в указанных выше работах, годовая сумма осадков изменяется от 400 до 800 мм. Наибольшее количество осадков выпадает в верховьях Оби, приуроченных к горным системам Алтая. Второй их максимум — на Полярном Урале. На равнинной части больше всего осадков выпадает в Среднем Приобье (600–650 мм). К северу сумма осадков уменьшается.

В соответствии с распределением по территории бассейна увлажнения и тепловых ресурсов суммарное испарение в пределах водосборной площади Оби в основном подчиняется широтной зональности, за исключением горных областей, где преобладает влияние высотной поясности. В истоках Оби и Иртыша (высокогорная часть бассейна) суммарное испарение составляет 250–300 мм/год. Наибольшие величины испарения (450 мм) приходятся на среднюю часть водосбора, расположенную в зоне мелколиственных лесов и южной тайги. Севернее значения испарения снижаются.

Основную массу воды р. Обь получает в результате таяния снегов. Река характеризуется сравнительно невысоким и растянутым (до 4 мес) весенне-летним половодьем, повышенным стоком летне-осеннего периода и низкой зимней меженью. Средняя продолжительность половодья изменяется от верховьев до устья от 100 до 140 дней. На период весенне-летнего половодья приходится существенная доля годового стока. Питание р. Иртыша смешанное: в верховьях снеговое, ледниковое и, в меньшей степени, дождевое; в нижнем течении — снеговое и дождевое.

Речной сток в бассейне Оби зарегулирован в основном малыми и небольшими водохранилищами, в регионе имеется только одно крупное водохранилище в верхнем течении реки — Новосибирское. Среднемноголетний объем стока воды этого водохранилища в створе гидроузла равен 50–55 км³/год, что составляет 12–14% среднемноголетнего объема речного стока в устье Оби.

Что касается Иртыша, то эта река на казахской территории в своем верхнем течении зарегулирована водохранилищами Иртышского каскада ГЭС: Бухтарминским, Усть-Каменогорским и Шульбинским. Среднемноголетние сбросы через гидростворы Бухтарминской, Усть-Каменогорской и Шульбинской ГЭС составляют соответственно 4.6, 4.8 и 1.7% среднемноголетнего объема стока в устье реки Оби.

Основываясь на приведенных данных, в целом речной сток Обь-Иртышского бассейна можно считать слабо зарегулированным, что и дает возможность попытаться использовать для расчета динамики составляющих водного баланса бассейна модель SWAP без учета режима сбросов указанных выше водохранилищ.

В пределах равнинной части Обь-Иртышского бассейна расположены тундровая, лесотундровая, лесная, лесостепная и степная зоны с характерными для них почвами и растительностью. В северной части территории равнины распространены торфянисто-болотные и торфяно-болотные почвы, в южной части наряду с ними — обычны солонцы, солоди, встречаются также солончаки. Зональные типы почв — тундрово-глеевые, подзолистые, дерново-подзолистые, черноземы и темно-каштановые — приурочены к относительно дренированным территориям.

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

СХЕМАТИЗАЦИЯ ОБЬ-ИРТЫШСКОГО БАССЕЙНА И ЕГО ИНФОРМАЦИОННОЕ ОБЕСПЕЧЕНИЕ

Использование SWAP для проведения модельных расчетов динамики составляющих водного баланса и снежного покрова Обь-Иртышского бассейна предполагает разбиение его территории на ряд расчетных ячеек, как, например, в [17, 26]. Как видно из рис. 2, Обь-Иртышский бассейн до стоковой ст. Салехард представлен совокупностью 441 ячейки регулярной пространственной сетки с разрешением 1×1 по широте и долготе.

В пределах каждой ячейки моделирование ее теплового и водного режима проводится независимо от других ячеек. Величина стока с расчетных ячеек бассейна (с учетом площадей, рассчитанных на основе реальных водосборных территорий, привязанных к ячейке [26]) служит входной информацией для следующего блока модели — блока трансформации стока в речной сети. В настоящей работе за основу была принята схематизация речных бассейнов, разработанная в [31] и подвергнутая авторами небольшой коррекции. В полученной схематизации ячейки бассейнов соединены каналами русловой сети, трансформация стока в которых может рассчитываться с использованием различных подходов. В данной работе для этой цели применена простая линейная модель формирования водного баланса в русле [24, 28].

Необходимое при проведении гидрологических расчетов информационное обеспечение модели SWAP включало в себя заданную для каждой расчетной ячейки приземную метеорологическую информацию (температуру и влажность воздуха, атмосферное давление, скорость ветра, атмосферные осадки, а также величин приходящих потоков длинно- и коротковолновой радиации) с трехчасовым временным разрешением и параметры подстилающей поверхности, а именно параметры почвы, растительного покрова, топографические данные бассейна. В данной работе метеорологические характеристики были получены из данных наблюдений 115 метеорологических станций (рис. 2), для которых качество информации соответствовало поставленной задаче (при этом величины приходящих потоков коротковолновой и длинноволновой радиации рассчитывались по данным стандартных метеорологических измерений по методике, приведенной в [2, 29]). Интерполяция значений метеорологических элементов в центры расчетных ячеек бассейна проводилась с использованием модифицированной процедуры кригинга [3, 11]. Параметры подстилающей поверхности бассейна взяты из глобальных баз данных, а затем наиболее важные их них уточнены, как и в вышеупомянутых работах авторов [2–4], на основе метода оптимизации с использованием данных по стоку, измеренному на трех стоковых станциях рек Оби и Иртыша (рис. 2). Проверка результатов расчетов проводилась по данным о стоке на тех же стоковых станциях, но за более длительные периоды наблюдений.

 

Рис. 2. Расчетная схема Обь-Иртышского бассейна и его речной сети. Светло-серым цветом отмечен бассейн р. Иртыша до стоковой ст. Ханты-Мансийская (2), темно-серым — часть бассейна р. Оби до стоковой ст. Колпашево (3), не закрашена нижняя часть Обь-Иртышского бассейна до стоковой ст. Салехард (1).

 

При построении базы данных растительного покрова, адаптированной под модель SWAP, за основу была принята глобальная база данных из проекта GSWP-2 (Second Global Soil Wetness Project) [18, 32]. При этом использованы классификация типов подстилающей поверхности IGBP (International Geosphere-Biosphere Programme), включающая в себя 17 типов поверхности суши. База данных скорректирована для устранения грубых несоответствий и ошибок, а также дополнена рядом специфических параметров, необходимых при работе с моделью SWAP. Подробности указанной корректировки изложены в [3, 13]. Там же содержится и перечень параметров растительного покрова, используемых в модели SWAP.

При подготовке адаптированной под модель SWAP базы данных по почвенным параметрам за основу были приняты данные из информационной базы ISLSCP-II) / GSWP-2 (International Satellite Land-Surface Climatology Project, Initiative II), которые, в свою очередь, взяты из IGBP-DIS (Data Information System). Указанные данные также частично скорректированы в целях устранения неправдоподобных значений ряда параметров и согласования их между собой и с параметрами растительного покрова [13].

Для последней группы данных (топографических характеристик, включающих среднюю высоту ячейки расчетной сетки и углы наклона поверхности ячейки в меридиональном и широтном направлениях) высоты взяты из информационной базы Центра данных по системам наблюдений за ресурсами Земли — EDC (Earth Resources Observation Systems Data Center), углы наклона рассчитаны по данным о средних высотах соседних ячеек.

Набор характерных для водосборов северных рек параметров, подлежащих оптимизации, определен и обоснован в [3, 11]. Для Обь-Иртышского бассейна число оптимизируемых параметров равнялось 8. К ним относились, в частности, корректирующие множители (общие для всего бассейна) к значениям коэффициента фильтрации почвы K₀, толщины корнеобитаемого слоя h_r, альбедо растительного покрова в теплый период (в отсутствии снега) alb_Sum, альбедо растительного покрова в зимний период (при наличии снега на элементах растительности) alb_Win, обозначенные как k_K0, k_hr, k_albSum и k_albWin соответственно. Кроме того, принималось, что для каждой ячейки расстояние от поверхности почвы до глубины практически водонепроницаемого слоя h₀ = k_h0 h_r, где k_h0 — множитель, общий для всего бассейна, который также калибровался. Наконец, калибровались такие параметры, как альбедо свежевыпавшего снега alb_sn, коэффициент шероховатости n в формуле Маннинга и эффективная скорость движения воды в русле u_e.

Поскольку территория Обь-Иртышского бассейна очень большая, для улучшения качества расчетов она разделена на три части (рис. 2) (как и бассейн р. Лены в [8]), для каждой из которых находился свой оптимальный набор параметров. Для нахождения оптимальных параметров подстилающей поверхности первой части Обь-Иртышского бассейна — бассейна р. Иртыша – использовались наблюдения за суточным стоком на стоковой ст. Ханты-Мансийск (калибровочный период – 1978–1987 гг.), для второй части бассейна — наблюдения за суточным стоком на стоковой ст. Колпашево (калибровочный период – 1968–1977 гг.), а для оставшейся третьей части — наблюдения за суточным стоком на стоковой ст. Салехард (калибровочный период – 1968–1977 гг.) (рис. 2). В последнем случае входной сток для двух ячеек, соседствующих с ячейками, в которых находятся стоковые станции Ханты-Мансийск и Колпашево, принимался равным не нулю (как для всех остальных пограничных ячеек третьей части бассейна), а величине рассчитанного стока соответственно с первой и второй частей Обь-Иртышского бассейна.

Калибровка модельных параметров для выбранных бассейнов проводилась автоматизировано с использованием метода глобального поиска экстремума целевой функции SCE-UA [19]. При этом ставилась задача не только достижения максимума целевой функции, в качестве которой использовалась эффективность расчета суточного стока по Нэшу и Сатклиффу Eff [27], но также и соблюдения условия для абсолютного значения систематической ошибки расчетов: |Bias| ≤5%.

С использованием статистик Eff и Bias для каждого из трех подбассейнов найден свой набор оптимальных значений параметров, что позволяет проводить расчеты составляющих водного баланса для любой части Обь-Иртышского бассейна.

РЕЗУЛЬТАТЫ МНОГОЛЕТНИХ РАСЧЕТОВ СОСТАВЛЯЮЩИХ ВОДНОГО БАЛАНСА ОБЬ-ИРТЫШСКОГО БАССЕЙНА

Оценка полученных при гидрологическом моделировании результатов проводилась по двум статистическим критериям соответствия рассчитанных и измеренных суточных значений стока: систематической ошибке Bias, равной разности между средними рассчитанных и измеренных значений стока, и упомянутой эффективности расчета Eff по Нэшу и Сатклиффу.

На рис. 3 показаны результаты сравнения измеренных и рассчитанных значений суточных слоев стока для упомянутых выше трех стоковых станций Обь-Иртышского бассейна. Полученные результаты свидетельствуют о хорошем соответствии измеренных и рассчитанных данных — для ст. Салехард (расчетный период — 1967–1999 гг.) эффективность расчета суточного стока Eff = 0.85, систематическая ошибка Bias = 5.5%, для ст. Колпашево (расчетный период – 1967–1994 гг.) Eff = 0.82, Bias = 3.4%, для ст. Ханты-Мансийск (расчетный период – 1978–1997 гг.) Eff = 0.84, Bias = –6.2%.

 

Рис. 3. Динамика измеренных (1) и рассчитанных (2) суточных значений слоя стока р. Оби в районе ст. Салехард (а) и ст. Колпашево (б), а также р. Иртыша в районе ст. Ханты-Мансийск

 

Можно заметить на рис. 3б, что рассчитанные пиковые значения стока в районе ст. Колпашево часто завышены, и в то же время часто занижены его меженные значения. Такое несколько худшее соответствие форм измеренных и рассчитанных гидрографов стока (по сравнению с гидрографами стока, измеренного на других станциях) отражается на немного сниженном значении Eff (0.82). Указанный эффект в районе ст. Колпашево может быть связан с недоучетом регулирования стока Новосибирским гидроузлом.

Представляет интерес сравнение полученных результатов для Обь-Иртышского бассейна с результатами аналогичных расчетов, выполненных в работе [30] на основе всемирно известной и широко используемой в гидрологических расчетах модели VIC [25]. В [30], к сожалению, приведена статистика только для месячных значений стока, поэтому и сопоставление будет проводиться только для месячных значений.

Для ст. Салехард при использовании модели SWAP для расчета месячных слоев стока Eff = 0.87, Bias = 5.5%, при использовании VIC Eff = 0.79, Bias = –0.9%. При расчете месячного стока р. Иртыша в районе ст. Ханты-Мансийск для модели SWAPEff = 0.88, Bias = –6.2%. В [30] сток р. Иртыш был рассчитан для стоковой ст. Тобольск, расположенной на ~400 км выше по течению ст. Ханты-Мансийск, при этом Eff  = 0.78, Bias = 2.1%.

Приведенное сравнение показывает, что эффективность расчетов стока с использованием модели SWAP выше эффективности модели VIC. Систематическая же ошибка при расчете стока в районе Обь-Иртышского бассейна с использованием SWAPнесколько выше, чем с использованием VIC. Однако в обоих случаях их оценки близки к точности измерения стока. Поэтому можно сделать вывод, что качество расчетов стока Обь-Иртышского бассейна с применением модели SWAP, по крайней мере, не ниже качества расчетов с использованием VIC.

На рис. 4 показаны климатические (усредненные по расчетным периодам) годовые гидрографы месячных слоев стока для трех использованных в работе стоковых станций. Видно, что соответствие рассчитанных и измеренных гидрографов достаточно хорошее. Статистика соответствия подтверждает этот факт: для трех станций расчета климатического месячного стока 0.92< Eff <0.98, 3.2< Bias <5.5%.

 

Рис. 4. Климатические (усредненные по расчетным периодам) измеренные (1) и рассчитанные (2) годовые гидрографы месячных слоев стока р. Оби в районе станций Салехард (а) и Колпашево (б), а также р. Иртыша в районе ст. Ханты-Мансийск (в).

 

Модель SWAP позволяет рассчитывать все составляющие водного баланса как всего бассейна, так и отдельных его частей. На рис. 5 приведена динамика годовых величин осадков, суммарного испарения и стокообразования для трех выбранных частей Обь-Иртышского бассейна, а также для всего бассейна в целом.

Интерполированные по территории бассейна с помощью операции кригинга значения измеренных годовых осадков, в среднем равные 470 мм, варьируют по трем выделенным подбассейнам и по годам в диапазоне от 350 до 650 мм/ год, что примерно соответствует данным, представленным в [15, 16]. При этом наименьшие величины осадков относятся к бассейну р. Иртыша, наибольшие — к северной части Обь-Иртышского бассейна.

Рассчитанное испарение также варьирует в среднем от 360 (в районе бассейна р. Иртыш) до 260 мм в нижней (северной) части бассейна р. Оби, составляя в среднем по всему Обь-Иртышскому бассейну ~320 мм, что также приблизительно соответствует данным об испарении в [15, 16].

Что касается стока, то наибольшая интенсивность его формирования приходится в основном на бассейн непосредственно р. Оби, меньшая — на бассейн р. Иртыш.

ОЦЕНКА ЕСТЕСТВЕННОЙ НЕОПРЕДЕЛЕННОСТИ ВЕЛИЧИН РЕЧНОГО СТОКА С ОБЬ-ИРТЫШСКОГО БАССЕЙНА

Как видно из приведенных в предыдущем разделе результатов оценок составляющих водного баланса и данных измерений речного стока с Обь-Иртышского бассейна, все указанные переменные (в масштабе как суток, так и лет) обладают естественной изменчивостью. Данная изменчивость, как показано в [14, 20, 23], обусловлена хаотическим характером атмосферных процессов (так называемым “климатическим шумом”), неустойчивостью по отношению к малым ошибкам в задании начальных условий при расчете эволюции метеорологических характеристик. Как показано в [20], формирование речного стока и его естественной неопределенности в рамках конкретного бассейна происходит под влиянием вынуждающей силы (“forcing”), в качестве которой выступают атмосферные процессы, при этом неустойчивость их характеристик приводит к неустойчивости значений речного стока. Причина возникновения атмосферного хаоса — объективная неустойчивость по отношению к начальным условиям: малое изменение начального условия со временем приводит к достаточно большим изменениям динамики системы. Для нелинейной системы с диссипацией невозможно конкретно предсказать ход ее развития, поскольку начальные условия состояния атмосферы никогда не могут быть заданы с абсолютной точностью.

 

Рис. 5. Динамика годовых значений осадков (1), суммарного испарения (2) и стокообразования (3) для трех выбранных частей Обь-Иртышского бассейна: бассейна р. Иртыша до ст. Ханты-Мансийск (а), верхней части бассейна р. Оби до ст. Колпашево (б), нижней части бассейна р. Оби до ст. Салехард (в), а также для Обь-Иртышского бассейна в целом (г).

 

В этом разделе рассмотрим оценку масштаба указанной неопределенности стока для р. Оби. В качестве меры неопределенности Un той или иной положительно определенной геофизической характеристики Х (в данном случае речного стока), имеющей случайную составляющую, будем использовать следующий показатель [14]:

Un = (x_0.975 – x_0.025) / M(X), (1)

где x_0.975 и x_0.025 — квантили 97.5- и 2.5%-й обеспеченности величины X, M(X ) — ее математическое ожидание. Указанная характеристика представляет собой отношение диапазона, в который может с 95%-й обеспеченностью попасть случайная величина, к ее среднему значению. Саму разность x_0.975 – x_0.025 назовем абсолютной неопределенностью величины Х.

На рис. 6а показаны средние за 33 года значения рассчитанного годового стока и диапазон его естественной изменчивости (абсолютной неопределенности), а также измеренные и рассчитанные конкретные реализации стока за 1967–1999 гг. Само значение неопределенности Un годового стока для р. Оби (в предположении, что функция распределения вероятности значений стока может быть аппроксимирована логнормальным распределением) оценено на основе годовых значений стока в районе ст. Салехард и оказалось равным 0.57.

На рис. 6б показан климатический (осредненный за 1967–1999 гг.) гидрограф рассчитанного месячного стока р. Оби, диапазон его абсолютной неопределенности (оцененной с принятием указанного выше предположения) и конкретные реализации измеренных месячных значений стока за указанные годы. На рис. 6в представлен годовой ход значений неопределенности месячного стока.

 

Рис. 6. Значения неопределенности слоя стока р. Оби в районе стоковой ст. Салехард: (а) — абсолютная неопределенность годового стока (черная жирная прямая линия — среднее за расчетный период значение рассчитанного годового стока, пунктирные прямые линии показывают диапазон абсолютной неопределенности годового стока, черная и серые тонкие линии показывают конкретные реализации годового стока — измеренного и рассчитанного соответственно); (б) — абсолютная неопределенность месячного стока (черная жирная и пунктирные линии показывают средние за расчетный период значения рассчитанных месячных слоев стока и диапазон их неопределенности соответственно, серые точки — конкретные реализации в разные годы рассчитанных месячных значений слоя стока; (в) — внутригодовая динамика неопределенности месячных значений слоя стока.

 

Как видно из приведенных рис. 6а, 6б, относительные значения неопределенностей для месячного стока, естественно, выше, чем для годового. При этом пик значений неопределенности месячного стока приходится на апрель–май, т.е. на период начала весеннего половодья. Это объясняется тем, что в этот период действуют две причины, обусловливающие неопределенность стока: неопределенность метеорологических характеристик, от которых зависит интенсивность стокообразования, и неопределенность момента начала формирования стока, который меняется от года к году.

РАСЧЕТЫ ФОРМИРОВАНИЯ СНЕЖНОГО ПОКРОВА В ОБЬ-ИРТЫШСКОМ БАССЕЙНЕ

Поскольку, как сказано выше, основная часть стока р. Оби формируется за счет снеготаяния, представляет интерес оценить качество моделирования с использованием SWAP-динамики снегозапасов в Обь-Иртышском бассейне.

На рис. 7а показана многолетняя динамика величин снегозапасов, осредненных по всему бассейну. Средние по бассейну снегозапасы определялись на основе данных измерений 128 маршрутных снегосъемок на 115 снегомерных станциях, расположенных непосредственно в бассейне или рядом с ним (рис. 2). Поскольку не всегда имелись данные снегосъемок за конкретные сутки для каждой станции, в целях повышения качества оценки средних по бассейну измеренных снегозапасов осреднения данных снегомерных наблюдений проводились только для тех суток, для которых количество маршрутов превышало 85. Рассчитанные средние по бассейну снегозапасы находились путем осреднения их величин, полученных для каждой ячейки на основе модели SWAP с учетом площадей ячеек. Среднеквадратическое отклонение измеренных и рассчитанных средних по бассейну снегозапасов оказалось ~10 мм, эффективность расчетов по Нэшу и Сатклиффу — 0.80. Учитывая точечный характер процедуры измерения снегозапасов, а также редкость снегомерной сети для территории Обь-Иртышского бассейна (одна снегомерная станция на ~25 тыс. км²), соответствие измеренных и рассчитанных средних по бассейну снегозапасов можно признать достаточно хорошим.

Хорошее качество расчета снежного покрова подтверждает интегральная функция его распределения P (обеспеченность) по территории бассейна. На рис. 7б приведен пример сопоставления интегральных функций распределения климатических величин снегозапасов (усредненных по расчетному периоду) по бассейну на 31 марта, полученных на основе маршрутных снегомерных наблюдений и данных о снегозапасах в расчетных ячейках Обь-Иртышского бассейна, рассчитанных на основе модели SWAP.

 

Рис. 7. Сопоставление характеристик снегозапасов в районе Обь-Иртышского бассейна, полученных по данным снегомерных измерений и рассчитанных на основе модели SWAP: (а) — динамика усредненных по бассейну р. Оби измеренных (кружки) и рассчитанных (линия) снегозапасов; (б) — интегральная функция распределения P по территории бассейна климатических (усредненных по расчетному периоду) значений снегозапасов, полученных на основе снегомерных наблюдений (1) и рассчитанных с помощью модели SWAP (2) на 31 марта.

 

ВЫВОДЫ

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

Проведенные с использованием модели тепловлагообмена подстилающей поверхности суши с атмосферой SWAP расчеты составляющих водного баланса Обь-Иртышского бассейна показали, что модель SWAP в сочетании с информационным обеспечением, основанным на глобальных базах данных по параметрам подстилающей поверхности и на данных наблюдений метеорологических станций, при проведении указанной выше оптимизации нескольких наиболее важных модельных параметров для трех частей Обь-Иртышского бассейна позволяет достаточно хорошо воспроизводить гидрографы суточного и месячного стока в различных створах р. Оби и ее притоков, а также динамику и других составляющих водного баланса в бассейне.

Проведена оценка естественной неопределенности годовых и месячных величин слоя стока в замыкающем створе р. Оби в районе ст. Салехард. Данная неопределенность обусловлена климатическим шумом метеорологических элементов, которые формируют динамику гидрологического цикла территории Обь-Иртышского бассейна. Показано, что пик значений неопределенности месячного стока р. Оби приходится на апрель–май, т.е. на период начала весеннего половодья.

Воспроизведение моделью SWAP динамики средних по Обь-Иртышскому бассейну снегозапасов также можно признать вполне удовлетворительным, что свидетельствует об адекватной формализации процессов формирования снежного покрова моделью SWAP.

Благодарности

Авторы выражают признательность сотрудникам Мирового центра данных по речному стоку (The Global Runoff Data Centre, D − 56068 Koblenz, Germany) за предоставленные данные по стоку рек Оби и Иртыша.

 

Работа выполнена при финансовой поддержке РФФИ (проект № 17-05-00555) ― раздел “Оценка естественной неопределенности значений речного стока с Обь-Иртышского бассейна” и РНФ (проект № 16-17-10039) ― разделы “Схематизация Обь-Иртышского бассейна и его информационное обеспечение”, “Результаты многолетних расчетов составляющих водного баланса Обь-Иртышского бассейна”, “Расчеты формирования снежного покрова в Обь-Иртышском бассейне”.

About the authors

E. M. Gusev

Water Problems Institute of the Russian Academy of Sciences

Author for correspondence.
Email: sowaso@yandex.ru
Russian Federation, Moscow

O. N. Nasonova

Water Problems Institute of the Russian Academy of Sciences

Email: sowaso@yandex.ru
Russian Federation, Moscow

E. A. Shkurko

Water Problems Institute of the Russian Academy of Sciences

Email: sowaso@yandex.ru
Russian Federation, Moscow

L. Ya. Dzhogan

Water Problems Institute of the Russian Academy of Sciences

Email: sowaso@yandex.ru
Russian Federation, Moscow

G. V. Ayzel

Water Problems Institute of the Russian Academy of Sciences

Email: sowaso@yandex.ru
Russian Federation, Moscow

References

Supplementary files