Вычислительные технологии
Том 11, часть 2, Специальный выпуск, 2006
РЕЧНОЙ СТОК В БАССЕЙНЕ РЕКИ ЛЕНЫ В УСЛОВИЯХ ВЕРОЯТНОГО ГЛОБАЛЬНОГО ПОТЕПЛЕНИЯ КЛИМАТА
А. Г. Георгиади, И. П. Милюкова Институт географии РАН, Москва, Россия e-mail: [email protected]
Results of estimation of river runoff deviations within Lena river basin possible at global climate warming in XXI century caused by growth of the green house gases of anthropogenous origin in atmosphere are submitted. As scenarios of climatic changes results of experiments on global models of the general circulation of atmosphere and ocean, included in the program of the Intergovernmental Panel on Climate Change (IPCC) which primary goal consists in comparison of various scenarios of global climate warming have been used. For estimation of hydrological changes in Lena river basin the macroscale model of monthly water budget developed in Institute of Geography of the Russian Academy of Science [1, 2] is used.
1. Модель месячного водного баланса
Начальная версия модели, примененная к крупнейшим речным бассейнам Русской равнины, и ее последующие модификации подробно рассматриваются в [1, 2]. В модели описываются основные процессы гидрологического цикла: инфильтрация и аккумуляция влаги в почве, испарение (на основе модифицированного метода Торнтвайта |3|), аккумуляция воды в снежном покрове п снеготаяние (на основе метода В.д. Комарова [4]), формирование поверхностного и подповерхностного подземного стока в реки и полного речного стока. Модель позволяет учесть макромасштабные неоднородности гидрометеорологических полей и других характеристик территории. Такой подход обеспечивает необходимую степень надежности моделирования изменений речного стока. В модели месячного водного баланса изменения речного стока и других элементов водного баланса в речных бассейнах оцениваются в узлах регулярной сетки, что облегчает стыковку с данными экспериментов на моделях климата.
При разработке модели применительно к условиям многолетнемерзлых почвогрунтов предполагалось, что процесс формирования стока в общем виде можно представить следующим образом. В холодное время года осадки выпадают в твердом виде (снег) и накапливаются. Верхние слои почвогрунтов при отрицательных температурах воздуха промерзают, образуя практически водонепроницаемый слой. При положительных температурах начинаются снеготаяние и сезонное протаивание почвогрунтов. Осадки и талые
© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2006.
воды, попадая на поверхность водосбора, частично испаряются с поверхности, остальная часть просачивается в верхний талый слой почвогрунтов, мощность которого постепенно увеличивается по мере оттаивания грунта. Часть воды просачивается в нижележащие горизонты подземных вод. Если происходит насыщение почвогрунтов активного слоя до наименьшей влагоемкости, то из излишков влаги формируется быстрый сток (главным образом поверхностный сток и сток из активного слоя надмерзлотной зоны). Влага, просочившаяся в подземные водоносные горизонты, перераспределяясь, идет на формирование стока из них и фильтруется в более глубокие водоносные горизонты вне зоны активного водообмена.
Модель основана на уравнении среднемноголетнего месячного водного баланса речных водосборов, которое в самом общем виде для каждого расчетного узла регулярной сетки может быть записано следующем образом:
<.(Ь) + <дг(Ь) = Р(Ь) - Е(Ь) - /¿(Ь) - ёШ/йЬ,
где <3(Ь) — сумма поверхностного и подповерхностного (из сезонно-талого слоя) стока, мм; <дг(Ь) — сток из подземного горизонта, мм. Сумма <<.3^) и <дг(Ь) составляет полный речной сток; Р(Ь) — атмосферные осадки, мм; Е(Ь) — испарение, мм; /¿(Ь) — фильтрация воды в глубокие горизонты подземных вод вне зоны активного водообмена, мм; йШ/йЬ — изменение запаса воды в зоне активного водообмена речного бассейна на момент времени Ь
Для расчетов водного баланса в бассейне роки Лены ряд блоков модели был модифицирован применительно к условиям многолетнемерзлых грунтов, а также подготовлены блоки для расчетов сезонного протаивания (промерзания) почвогрунтов, В этой версии модели используются две принципиально различные схемы: одна из них позволяет оценить динамику сезонного оттаивания-промерзания активного слоя многолетнемерзлых грунтов, и, соответственно, в другой схеме, используемой для районов, в которых отсутствуют многолетнемерзлые грунты, рассчитываются динамика глубины сезоннопромерзающего слоя почвогрунтов и ее влияние на формирование стока.
Для расчета перемещения границы оттаивания (промерзания) использовано уравнение [5, 6]
= БЬ' + {(БЬ' + Zj-1)2 + 2А*Тр (ДЬ/рЕ Ь )}1/2,
где Zj — глубина границы протаивания в ^'-е сутки; Ь — удельная теплота плавления льда; А
нии) почвогрунтов; Ас — коэффициент теплопроводности снега; Н и Тр — соответственно высота снежного покрова и температура подстилающей поверхности за рассматриваемый интервал времени; Ак — параметр; Wz — объемная влажность почвогрунта на фронте протаивания (в долях от 1); рв — плотность воды.
При расчетах перемещения границы промерзания почвогрунтов БЬ' = (А/АС)Нпри расчетах перемещения границы оттаивания почвы БЬ' = А/Ак.
При отсутствии снежного покрова задавалась зависимость между температурой воздуха и температурой почвы с учетом фактора вечной мерзлоты. Исходя из имеющихся данных наблюдений можно задать, что
Т/ = Т'"кг, к = 1.0 - 1.3,
Т
Коэффициент теплопроводности талых и мерзлых пород рассчитывался из соотношения [5]
Ат = + 0.1Ш - 1.1) - 0.Щ
ГД6 — удельный вес скелета почвогрунта, г/см; кР = 1.5 ... 1.7 для песков, кР = 1.4 ... 1.5
для супесей, кр = 1.3 ... 1.4 для суглинков и глин (первая цифра для талых почвогрунтов, вторая — для мерзлых); Ш — объемная влажность почвы, %; Ас = 4.32 + 1.64рс + 518.4р^, рс — плотность снега, изменяющаяся во времени: рс. = рС0 + рс._1 рскс, рС0 — начальная
кс
Фильтрация в подземный водный горизонт активного водообмена (1Р) за ;'-е сутки определялась из следующих соотношений: при отсутствии вечной мерзлоты
Ц = kf (Р* + М* - Б*) // ;
при наличии вечной мерзлоты
I* = /к/ (Р* + М* - Б*) /Ш/с) ,
где М — суточное снеготаяние, мм; Ш/с — наименьшая влагоемкость иочвогрунов, мм; к/ — коэффициент вертикальной фильтрации в подземные воды, параметр модели; / — параметр, зависящий от типа вечной мерзлоты в данном узле расчетной сетки. Подземный сток рассчитывался следующим образом:
Qgr = кдг ( Ндт /Итах),
где Нтах — максимальная емкость подземного водоносного горизонта; Идг — емкость подземного горизонта в ;'-е сутки; кдг — коэффициент подземного стока (параметр модели);
иИдг = Идг + 1р - Фдг - ^.
Фильтрация влаги в глубокие горизонты подземных вод вне зоны активного водообмена 1(1 принималась равной постоянной величине.
Изменение мощности подземного водоносного горизонта активного водообмена для ;'-х суток рассчитывалось из соотношения
А И * = И * — И
дг дг дг
Фильтрация воды в подземный водный горизонт активного водообмена 1Р и Qgr для разных типов многолетнемерзлых почвогрунтов (сплошная, прерывистая или островная вечная мерзлота и в их отсутствие) рассчитывалась с эмпирическими весовыми коэффициентами.
2. Общая схема расчетов
В качестве исходной информации для модельных расчетов использовались среднемесячные среднемноголетние данные о температурах воздуха и суммах атмосферных осадков, а также данные о влагозапасах и физических характеристиках почвогрунтов (наименьшей влагоемкости почвогрунтов, их объемном и удельном весах), о характере распространения вечной мерзлоты (сплошная, прерывистая, островная, отсутствие многолетнемерзлых
грунтов), максимальной мощности сезонно-талого слоя почвогрунтов, данные о стоке в узлах регулярной сетки.
Расчеты проводились с суточным шагом по времени. Поскольку для расчетов использовались многолетние среднемесячные данные наблюдений, проводилась интерполяция данных наблюдений по суткам. При задании внутримесячной изменчивости температуры воздуха предполагалось, что среднемесячная температура воздуха приходится на 15-е число каждого месяца, а температура на 1-е число всех месяцев с января по июль принималась па Тс^ меньше, тогда как температура па последний день этих месяцев задавалась на Т°А больше среднемесячной, С августа по декабрь делалась обратная процедура. Внутри месячного интервала проводилась линейная интерполяция между заданными температурами воздуха на 1-, 15-е и последнее число каждого месяца. Величина отклонения температуры воздуха на начало и конец месяца от среднемесячной (±Т^) подбиралась эмпирически.
Расчеты проводились с 1 января по 31 декабря. Однако для расчетов кроме исходной информации о метеоданных и почвенных характеристиках необходимо иметь данные о запасах влаги в почве, запасах воды в снеге, мощности подземного водоносного горизонта на дату начала расчетов (1 января). Для приближенного задания начальных условий проводился расчет этих величин за предыдущие четыре месяца (с сентября по декабрь), в дальнейшем данные расчета использовались в качестве начальных условий на 1 января,
3. Методы калибровки параметров
Модель включает 15 параметров. Семь из них (Т°р, Ьт , кг, Ь/, Ндгшах, Т°Д, рс.) подбирались эмпирически. Для идентификации остальных параметров модели использовалась процедура оптимизации X, Розенброка [7], которая представляет собой метод покоординатного спуска с попеременной заменой осей.
Оптимизация модели проводилась на основе гридированных среднемноголетних данных о среднемесячной температуре воздуха, атмосферных осадках, испарении и испаряемости из архива Международного института прикладного системного анализа (Австрия) [8] и среднемесячного речного стока бассейна реки Лены, полученного на основе интерполяции данных о стоке репрезентативных речных бассейнов в пределах бассейна реки Лены, Для оптимизации параметров использованы данные по 50 % узлов регулярной сетки, которые равномерно распределены по территории речного бассейна. Оставшиеся 50 % данных служили основой для верификации модели.
Оптимизация проводилась в несколько этапов отдельно для разных блоков модели. На первом этапе, используя блоки расчета испаряемости и испарения по многолетним данным о среднемесячных фактических величинах испарения, оптимизировались параметры а, а^ и в- В качестве критерия качества принималась сумма среднеквадратичных отклонений расчетных значений месячного испарения от фактических. На втором этапе оптимизировались параметры для расчетов сезонного оттаивания (промерзания) почвы Ак
кС
значений максимального сезонного оттаивания почвы от фактических. Далее в первом приближении оптимизировались коэффициент стаивания к, коэффициент вертикальной фильтрации в грунтовые воды к/ и коэффициент грунтового стока кдг по данным о стоке, И на четвертом этапе вновь уточнялись параметры Ак и кС по данным о сезонном оттаивании почвы. На последнем этапе по данным о стоке уточнялись эти коэффициенты. Кроме
того, на всех этапах проводился эмпирический подбор и уточнялись остальные параметры модели. Подбирались наиболее устойчивые параметры, величину которых можно задать приблизительно без ущерба для точности расчетов.
Такая многоступенчатая процедура идентификации параметров, как показала практика использования гидрологических моделей, несмотря на свою громоздкость, вполне себя оправдывает. Хотя наиболее достоверной гидрологической информацией обычно являются данные о стоке, для первого приближения бывает удобно использовать и другую гидрометеорологическую информацию, что обычно повышает надежность определения параметров,
4. Верификация модели
Модель верифицировалась по критерию близости рассчитанных и наблюденных данных
0 годовом стоке и стоке за период половодья. Мерой для верификации модели служило соотношение So/a, используемое в России в практике гидрологических прогнозов и расчетов. Здесь So = [Е (Robs - W] /П * = [Е (Robs - Rmeajf] /(n - 1), где Robs - наблюденное значение стока; Rcaic — рассчитанное значение стока; Rmeaj — среднее значение стока, число членов в выборке.
Верификация модели проводилась на основе сравнения данных инструментального периода и результатов математического моделирования. Для этого использована независимая выборка гридированных среднемноголетних наблюденных и рассчитанных по модели месячного водного баланса данных о годовом стоке и стоке периода половодья. Выборка сформирована из оставшихся (независимых) данных по 50% узлов регулярной сетки, равномерно распределенных по территории речного бассейна.
Как в случае с годовым стоком (So/a = 0.68), так и со стоком весеннего половодья (So/о = 0.56) значения этого критерия оказались меньше 0,75, что говорит об удовлетворительной надежности результатов, полученных при использовании модели месячного водного баланса,
5. Описание использованных данных
Данные о современном климате. Использованы среднемноголетние данные о среднемесячных значениях температуры воздуха и атмосферных осадках в узлах регулярной сетки
1 х 1° из дополненного глобального архива, подготовленного в Международном институте прикладного системного анализа [8]. Эти данные, так же как данные по стоку и отклонениям климатических элементов, сглажены.
Данные о речном стоке. Для подготовки данных о среднемноголетних месячных значениях речного стока в узлах регулярной сетки 1 х 1° использована информация о сред-немноголетнем месячном стоке средних рек бассейна реки Лены, Данные о стоке относились к средневзвешенному центру бассейнов, а затем интерполировались в узлы регулярной сетки, В расчетах использовались данные о стоке, осредненные за периоды времени, превышающие 30 лет и относящиеся в большинстве своем к временному отрезку между 30-ми — 40-ми и 80-ми годами 20 в.
Данные об отклонениях климатических элементов от их современных значений. В качестве сценарных оценок измененных в будущем климатических условий использованы результаты численных экспериментов на двух моделях общей циркуляции атмосферы и
океана, участвующих в программе 1РСС: модели Метеорологического института им, Макса Планка, Германия (МР1М ЕСНАМ4/ОР¥3), и модели Геофизической лаборатории гидродинамики Припстопского университета, США (СЕБЬ-1130),
Использованы результаты модельных экспериментов, исходящих из семейства А2 сценариев глобальных социально-экономических изменений в 21 в, из последней, улучшенной серии таких сценариев SR.ES, принятой в программе 1РСС, Использованы два периода (2010-2039, 2040-2069 гг.) для расчета отклонений средпемпоголетпих среднемесячных сумм атмосферных осадков и температуры воздуха от современных значений (рассчитанных для периода 1961-1990 гг.).
6. Обсуждение результатов моделирования
Климатические условия. Согласно рассматриваемым сценариям в равнинной, центральной части бассейна реки Л е н ы ожидается существенное потепление климата, более интенсивное в середине 21 в. На рисунках приведены отклонения среднемесячных значений температуры воздуха (ДТ, °С), атмосферных осадков (ДР, мм) и речного стока (Дф, мм) для условий потепления климата от современных их значений, обобщенных по узлам сетки, покрывающей равнинную часть бассейна р. Лепа в 2010-2039 гг. и 2040-2069 гг. На рис. 1 и 2, а — по сценарию Геофизической лаборатории гидродинамики Припстопского университета, США; б — по сценарию Метеорологического института им. Макса Плап-
Рые. 1. Отклонения среднемесячных значений в 2010 2039 гг.
Рис. 2. Отклонения среднемесячных значений в 2040 2069 гг.
ка, Германия; в — среднемесячные значения температуры воздуха (Т, °С), атмосферных осадков (Р, мм), речного стока мм) для условий потепления современного климата.
Согласно сценарию Метеорологического института им. Макса Планка, потепление может быть более ощутимым. Оба сценария прогнозируют более заметный рост температуры воздуха в холодную часть года, что приведет к сокращению ее продолжительности. Характер изменений атмосферных осадков более сложен. Лишь к середине этого столетия согласно обоим сценариям будет заметна тенденция повышения атмосферного увлажпе-
Следует заметить, что характер впутригодового распределения сценарных изменений среднемесячной температуры воздуха и атмосферных осадков весьма схож с соответствующими изменениями, наблюденными в этом регионе в последние десятилетия в период современного потепления климата, по отличается от современных масштабом изменений.
Гидрологические условия. Согласно климатическим сценариям Геофизической лаборатории гидродинамики Припстопского университета и Метеорологического института им. Макса Планка, в равнинной части бассейна реки Лены в течение первой половины теку того столетия не следует ожидать сколько-нибудь ощутимого повышения как годового стока, так и стока за период половодья (три месяца с наибольшим месячным стоком). Лишь при реализации сценария Геофизической лаборатории гидродинамики Припстоп-
ского университета средний по территории сток половодья может возрасти более чем на 10 мм. Следует заметить, что распределение возможных изменений речного стока по территории этой части бассейна характеризуется существенной неоднородностью, которая заметно снижается к середине столетия.
Согласно обоим сценариям могут произойти значительные изменения внутригодово-го распределения стока. Основные различия в его изменениях между двумя сценариями в первой трети текущего века приурочены к части теплого периода после прохождения основной волны половодья, тогда как прогнозируемые по обеим моделям сценарные изменения структуры половодья очень близки между собой (см. рис. 1). К середине века оба сценария дают по существу одинаковую картину изменений внутригодового распределения стока — волна половодья компактно (практически без изменения формы и объема воды за половодье) сдвигается на один месяц раньше по сравнению с современной ситуацией. Заметим, кстати, что характер изменений внутригодового распределения стока в бассейне рок и Лены сходен с аналогичными изменениями при вероятном глобальном потеплении климата в текущем столетии в бассейне реки Волги.
Заключение
1. Согласно рассмотренным сценариям в равнинной (центральной) части бассейна реки Лены ожидается существенное потепление климата, более интенсивное в середине 21 в.
2. Характер внутригодового распределения сценарных изменений среднемесячной температуры воздуха и атмосферных осадков весьма схож с соответствующими изменениями, наблюденными в этом регионе в последние десятилетия в период современного потепления климата, но отличается от современных масштабом изменений.
3. В течение первой половины текущего столетия не следует ожидать сколько-нибудь ощутимого повышения как годового стока, так и стока за период половодья.
4. Могут произойти значительные изменения внутригодового распределения стока. К середине века оба сценария дают по существу одинаковую картину изменений внутригодового распределения стока — волна половодья компактно (практически без изменения формы и объема воды за половодье) сдвигается на один месяц раньше по сравнению с современной ситуацией.
Авторы признательны Р. Лимансу за предоставление глобальной базы данных о современном климате, подготовленной в Международном институте прикладного системного анализа, г. Лаксенбург, Австрия.
Список литературы
[1] Георгиади А.Г., Милюкова И.П. Масштабы гидрологических изменений в бассейне реки Волги при антропогенном потеплении климата // Метеорология и гидрология. 2002. № 2. С. 72-79.
[2] Георгиади А.Г., Милюкова И.П. Возможные изменения речного стока в бассейнах крупнейших рек Русской равнины в XXI веке // Водное хозяйство России. 2006. Jfa 1. С. 62-77.
[3] Willmott C.J., Rowe С.М., Mintz Y. Climatology of the terrestrial seasonal water cycle // J. of Climatology. 1985. Vol. 5. P. 589-606.
[4] Руководство по гидрологическим прогнозам. Выи. 1. Долгосрочные прогнозы элементов водного режима рек, озер и водохранилищ. Л.: Гидрометеоиздат, 1989. 358 с.
[51 Павлов А.В. Теплофизика ландшафтов. Новосибирск: Наука, Сиб. овд-ние, 1979. 278 с.
[6] Бельчиков В.А., Корень В.И. Опыт использования модели формирования талого и дождевого стока рек лесной зоны европейской территории СССР // Тр. ГМЦ. Вып. 246. Л.: Гидрометеоиздат, 1983. С. 3-20.
[7] Розенброк X., Стори С. Вычислительные методы для инженеров-химиков: Пер. с англ. М.: Мир, 1968. 440 с.
[8] Cramer W.P., Leemans R. Global 30-Year Mean Monthly Climatology, 1930-1960. Vol. 2.1 (Cramer and Leemans). Available online (http://www.daac.ornl.gov/) from the ORNL Distributed Active Archive Center. Oak Ridge National Laboratory. Oak Ridge. Tennessee. U.S.A, 1999.
Поступила в редакцию 9 ноября 2006 г.