438 Северная Атлантика проведено в работе [23] на основе модели ИВМиМГ СО РАН. Кроме межгодовой изменчивости, в гидроло- гических характеристиках сибирских рек наблю- даются положительные тренды, вызывающие в последние годы интерес климатологов, так как это является отражением трендов в общей клима- тической системе. Так, в последние десятилетия в Арктике наблюдается рост речного стока [24—27]. По данным [28] с 1936 по 1999 гг. сток основных шести сибирских рек вырос на 7%. Этот рост хо- рошо коррелирует с увеличением индекса МАО [28]. Это также может быть связано с положитель- ным трендом МАМ (Noг! Аппшаг Моде), являю- щейся частью Арктической осцилляции [29]. При этом весьма интересным представляется вопрос: будет ли этот процесс продолжаться и в дальней- шем в ХХ столетии, продуцируя обратные связи в климатической системе. Все эти факты свиде- тельствуют о том, что влияние изменчивости ба- ланса пресной воды в Арктике на климатическую систему в будущем веке может быть существен- ным и требует изучения на основе климатических моделей речного стока [30—33], моделей циркуля- ции отдельных океанов [34—36], а также совмест- ных моделей климатической системы проекта СМИР-5 (см. напр., [37]. Метод математического моделирования для исследования процессов в ат- мосфере и океане, развитый в свое время Гурием Ивановичем Марчуком [38], является одним из наиболее эффективных инструментов исследова- ния климатической системы Земли. В настоящей статье рассматриваются результаты моделирова- ния межгодовых изменений пресноводного ба- ланса, которые могут происходить в Северном Ледовитом океане в ХХ веке при вариациях реч- ного стока сибирских рек, рассчитанных по моде- ли ИВМИМГ СО РАН на основе результатов рас- четов модели ИВМ РАН [39] сценария КСР 8.5 проекта СМ1Р-5. 1. РАСЧЕТЫ РЕЧНОГО СТОКА ДЛЯ СИБИРСКОГО РЕГИОНА 1.1. Климатическая модель речного стока Разрабатываемая модель является линейной резервуарной моделью. Модель составлена из ли- нейных резервуаров или каскада резервуаров в ячейках сетки. Скорость изменения стока из ячейки в простейшем варианте модели Калини- на—-Милюкова [24, 30] определяется на основе ре- шения интеграла свертки (Дюамеля). В конкретной реализации модели использова- лась структура, предложенная в институте Макса Планка в ГТамбургского университета [30]. В этом подходе поток воды на суше разделяется на три составляющие: поверхностный сток, грунтовый сток, речной сток. Значения коэффициентов за- ИЗВЕСТИЯ РАН. ФИЗИКА АТМОСФЕРЫ И ОКЕАНА КУЗИН и др. держки для поверхностного и речного стоков определяются по формулам, устанавливающим за- висимость от наклона ячейки или от перепада вы- сот между ячейками, отнесенного к расстоянию между их центрами. Поверхностный и грунтовый стоки рассчитываются для единичных резервуа- ров, а речной сток для каскада резервуаров. Коли- чество каскадов в руслах рек определяется на ос- нове полинома регрессии второго порядка, опре- деляемого по величине коэффициента задержки. Значение коэффициента задержки грунтового стока для ячейки принимается постоянным. Ин- фильтрация воды в грунт проводилась в каждой ячейке на основе характеристик поверхности и процента среднего грунтового влагосодержания на основе гидрологических данных для каждого бассейна. В каждой ячейке производится учет процентного содержания болот и озер. 1.2. Данные для моделирования речного стока и постановка эксперимента Рассматриваемая область покрывает террито- рию 40° No-80° No, 50° Е—170° Е по пространству. При проведении численных экспериментов по климатической модели речного стока было вы- брано разрешение, составляющее 1/3 градуса по широте и долготе соответственно. Для выбранного варианта модели стока необхо- димо задание следующих входных параметров: осадки, испарение, переходы из жидкой в твердую фазы и обратно, инфильтрация в почву. При вери- фикации модели эти данные были взяты из данных реанализа МЕККА для периода 1980—2011 гг. Для сравнения данных расчетов и данных наблюдений на гидрологических постах были использованы сле- дующие ежемесячные данные Гидрометеослужбы [15, 16] и данные К-АгсисМЕТ [Бир://умууг-атсис- пе!.5гопб.ейи/У4.0/таех.б0т||. Для проведения рас- четов в ХХ веке использовались данные расчетов по модели ИВМ РАН сценария КСР 8.5 проекта СМПР-51РСС. Расчеты по модели речного стока состояли из двух экспериментов. Первый эксперимент пред- назначен для верификации модели стока. Для этого были проведены численные эксперименты по моделированию стока сибирских рек по дан- ным реанализа МЕККА и построены среднекли- матические годовые гидрографы, которые срав- нивались с данными измерений. Во втором экс- перименте проведено моделирование речного стока в ХХ веке на основе данных расчетов по модели ИВМ РАН [39]. 1.3. Результаты моделирования речного стока В результате первого этапа были рассчитаны речные стоки на основе данных реанализа МЕККА на период 1980—2011 гг. сшагом по времени шесть No 4 том 51 2015