ВЛИЯНИЕ МЕЖГОДОВОЙ ИЗМЕНЧИВОСТИ СТОКА СИБИРСКИХ РЕК пространяются далее в СЛО [20, 21] в системе по- верхностных течений. Указанные выше значения годового стока для Енисея, Лены и Оби являются только некоторыми осредненными значениями за период в несколько десятилетий. По данным Гидрометслужбы [22, 23] суммарный годовой расход Енисея, Лены и Оби за рассматриваемый период с 1936 по 1990 гг. претер- певает существенные межгодовые вариации. Изме- нения составляют от 25 до 30%. Важную роль в этом играет межгодовая климатическая изменчивость атмосферной циркуляции и характеристик поверх- ности [24, 25]. Основным механизмом, влияющим на процессы в Арктическом регионе, являются Арктические колебания, или более точно — Кольце- вая мода Северного полушария (КМС) [26, 27]. В то же время корреляция между вариациями климати- ческих характеристик и речным стоком не столь вы- сока, как можно было ожидать. Для годового стока только 31—55% вариаций объясняются поверхност- ными и атмосферными условиями [28]. Кроме межгодовой изменчивости, в гидрологи- ческих характеристиках сибирских рек наблюдаются устойчивые тренды, вызывающие в последние годы повышенный интерес исследователей, так как это является отражением трендов в общей климатиче- ской системе и может продуцировать обратные свя- зи. Так, в последние десятилетия в Арктике наблюда- ется рост речного стока в Арктический океан [18, 20, 29, 30]. Во всей Евразии с 1936 по 1999 гг произошло увеличение речного стока на 128 км3 [18] или 10 км` в год [31]. Существующие гипотезы, объясняющие наблю- даемый тренд, различны, но в настоящее время нет полного подтверждения ни одной из них [18, 28, 30]. В работе [28] наличие трендов в последние десяти- летия связывается с увеличением индексов Северо- Атлантического колебания и Кольцевой моды Се- верного полушария. Анализ связи стока сибирских рек с характеристиками атмосферной циркуляции с использованием эмпирических и модельных дан- ных представлен в работе [32]. Влияние изменчивости речного стока на клима- тическую систему может быть существенным и тре- бует изучения на основе климатических моделей [33], поскольку метод математического моделирова- ния в настоящее время является одним из наиболее эффективных инструментов исследования климати- ческой системы Земли [34, 35]. На основе совмест- ных моделей крупномасштабной циркуляции океана и морского льда были рассчитаны пути распростра- нения речных вод в Арктическом бассейне [36—39]. Как правило, траектория движения вод определя- лась на основе решения задачи о распространении трассера от точечных источников, расположенных в устьях крупных рек, впадающих в Северный Ледо- витый океан. В этих работах, а также в многочис- & ИЗВЕСТИЯ РАН. ФИЗИКА АТМОСФЕРЫ И ОКЕАНА 833 ленных исследованиях, проводимых в рамках меж- дународного проекта АОМПР (Атсйс Осеап Моде! ГШмегсотрапзопо — Рго)ес!, — БИр://ЛуУУУМПО!.еаи/ра- ве.й0?р14=29900), в численных моделях динамики океана условия на границах для притока пресной воды рек задаются согласно имеющимся данным о среднеклиматическом сезонном изменении стока основных 13 рек Арктического региона. Целью дан- ной работы является исследование влияния межго- довой изменчивости стока сибирских рек на рас- пределение потоков пресной воды, поступающей из Арктики в Северную Атлантику. На основе чис- ленных расчетов исследуется чувствительность мо- дельного распределения к малым изменениям внешнего воздействия, связанным с вариациями стока основных сибирских рек. ПОСТАНОВКА ЧИСЛЕННОГО ЭКСПЕРИМЕНТА Структура численной модели Для проведения исследования использовалась совместная численная модель океан—лед, адаптиро- ванная к области Арктики и Северной Атлантики. Для описания океанических процессов использо- валась региональная численная модель динамики океана, разработанная в ИВМиИМГ СО РАН. Основ- ные положения последней версии модели представ- лены в работах [40, 41]. Уравнения модели выписаны в криволинейных ортогональных координатах с ис- пользованием традиционных приближений гидро- статики и Буссинеска. На поверхности океана ис- пользуется приближение “твердой крышки”. Урав- нения динамики решаются методом выделения баротропной и бароклинной мод, при этом баро- тропная часть представляется в виде решения урав- нения относительно функции тока. При решении за- дачи по времени используется комбинация явного и полунеявного подходов с использованием метода расщепления по физическим процессам и геометри- ческим направлениям, основы которого были зало- жены в работах Г.И. Марчука [42, 43]. Уравнение, описывающее перенос течениями, аппроксимирует- ся на основе применения схемы ООТСКЕЗТ [44]. Используемая модель морского льда, известная под названием эластическая вязкопластическая мо- дель (е1а50с у5сои5-р1азис), является модификацией стандартной вязкопластической модели динамики льда [45]. Данная модель подробно документирова- на в работе [46]. Толщина снежного покрова и льда вычисляется на основе термодинамической модели [47] для каждой категории льда. Горизонтальный перенос льда осуществляется с помощью адвектив- ной схемы полулагранжева типа [48]. Для моделирования процессов взаимодействия Арктического бассейна и Северной Атлантики рас- сматривается область Северного Ледовитого океана No 6 том 46 2010