CC BY
51
Вычислительные технологии
Том 11, часть 2, Специальный выпуск, 2006
ГИДРОЛОГИЧЕСКИЙ РЕЖИМ ВОСТОЧНО-ЕВРОПЕЙСКОЙ РАВНИНЫ И ЗАПАДНО-СИБИРСКОЙ НИЗМЕННОСТИ В КОНТРАСТНЫХ КЛИМАТИЧЕСКИХ УСЛОВИЯХ
П. А. Торопов Московский государственный университет,, Россия c-mail: tormetOinbox.ru
This article addressees an estimation of the reaction of the hydrologic regime of the East European and West Siberian plains on the long-term global climatic changes which took place during last 20 000 years. The calculations of the annual runoff were based on results of 16 models of global atmospheric circulation, taking part in the international project PMIP.
Введение
Одна из наиболее важных прикладных задач географии — оценка реакции речного стока на крупномасштабные климатические изменения. Кроме того, речной сток, являющийся интегральным показателем условий увлажнения определенного региона, удобно использовать для верификации результатов численного моделирования. Для подобных задач естественно использовать контрастные климатические события последних 20 тыс. лет.
В данной работе выполнена оценка реакции гидрологического режима Восточно-Европейской равнины и Западно-Сибирской низменности на глобальные климатические изменения, происходившие во время позднеплейстоценового похолодания (18-21 тыс. календарных лет назад) и в середине голоцена (5-6 тыс. календарных лет назад), а также проанализировано качество современных гидрологических моделей. Для расчетов использовались результаты 16 моделей общей циркуляции атмосферы (МОЦА), участвующих в международном проекте PMIP (Paleoclimate Modeling Intercomparison Project) [1].
1. Оценка качества моделирования современного климата и гидрологического режима
Проект PMIP направлен, во-первых, на выявление генезиса долгопериодных колебаний климата на примере событий прошлого, а во-вторых — на изучение чувствительности МОЦА различного пространственного разрешения на возмущения граничных условий.
© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2006.
В рамках проекта проводится три серии численных экспериментов. В первую очередь, осуществлялось моделирование современного климата, результаты которого использовались как “эталон” для сравнения с контрастными климатическими условиями. Вслед за этим воспроизводились условия холодной эпохи позднеплейстоценового криохрона (21 тыс. календарных лет назад) и теплого оптимума голоцена (6 тыс. календарных лет назад).
Описание моделей общей циркуляции атмосферы, участвовавших в РМ1Р, можно найти в [1]. Подробное тестирование этих моделей на качество воспроизведения современного климата выполнено в работах [2, 3]. Показано, что результатам численного моделирования можно доверять только при их осреднении по территории, площадь которой не менее 1.3 млн км2 (128 ячеек площадью 2.5x2.5°). И даже в пределах такой большой территории (соответствующей площади бассейна Волги) сезонный ход метеорологических характеристик воспроизводится моделями данного пространственного разрешения (в среднем 3x4°) со значительными погрешностями. Так, средняя ошибка годовой амплитуды температуры составляет 5(±4) °С, осадков — 20(±15) % от современной годовой суммы (в скобках указаны величины межмодельного разброса).
При тестировании МОЦА сталкиваются с еще одной проблемой, связанной с качеством климатических архивов, с которыми обычно сопоставляются результаты численного моделирования. В [2, 4] выполнена верификация некоторых климатических архивов осадков и испарения. Результаты этой верификации представлены в табл. 1. Хорошо видно, что расхождения между годовыми суммами осадков, испарения и слоя стока, осредненных по территории Восточно-Европейской равнины для перечисленных выше архивов, различаются на 20-25%. Это означает, что “климатического эталона”, по крайней мере для этих величин, фактически не существует.
Учитывая эти проблемы, в качестве “эталона” для оценки результатов численного моделирования представляется целесообразным выбрать средний многолетний годовой сток с крупных водосборов Восточно-Европейской равнины и Западно-Сибирской низменности (Волги, Оби, и т.д.), который является надежно измеряемой характеристикой [6].
Предлагается величину “гидрометрического стока” сопоставлять с “климатическим стоком”, полученным по результатам моделирования, который рассчитывается как разность
Таблица 1. Средние многолетние составляющие годового водного баланса в условиях современного климата, осредненные по волжскому бассейну, вычисленные по нескольким климатическим архивам осадков и испарения
Архив Осадки, мм Испарение, мм Сток, мм
NCEP/NCAR 590 586 4
ECMW 507 431 76
GPCP 543 - -
JAEGER 677 - -
LEEMENS 474 437 37
LEGATES 538 - -
СРЕДНЕЕ1 555 489 66
НАБЛЮДЕНИЯ2 179
Примечания.
1 Среднее арифметическое значение между рассмотренными архивами.
2
(ст. “Лебяжья”) [5].
“осадки (Р, мм) минус испарение (Е, мм)”. Полагается, что если ошибка моделирования лежит в пределах -20 % < АУ < 20 %, то МОЦА удовлетворительно воспроизводят годовой сток. Дело в том, что ошибка стока АУ, не превышающая по модулю 20 %, не выходит за пределы экстремальной межгодовой изменчивости, характеризуемой величиной 2аУ, (аУ — межгодовое среднеквадратичное отклонение стока). Подробно методика оценки годового стока с Восточно-Европейской равнины и Западно-Сибирской низменности описана в [2].
Результаты моделирования годового стока с обеих территорий, представленные в табл. 2, позволяют сделать два важных вывода. Во-первых, ни одна из моделей общей циркуляции атмосферы (кроме СС5Н) удовлетворительно не воспроизвела средний многолетний сток со всех выбранных водосборов. Соответственно, успешное воспроизведение
Таблица 2. Современный средний многолетний годовой сток, воспроизведенный 16 МОЦА РМ1Р, с водосборов Черного, Каспийского, Балтийского морей и реки Оби и его относительная ошибка (в % от наблюденного годового гидрометрического стока) [5, 7]
Модель общей циркуляции атмосферы Черное море Каспийское море Балтийское море Река Обь
Сток, Q КМ Откло- нение, % Сток, Q КМ Откло- нение, % Сток, Q КМ Откло- нение, % Сток, Q КМ Откло- нение, %
BMRC 88.9 -71 105 -62 377.0 -20 849 115
ССС2 214.4 -31 266.2 -3 518.9 10 616 56
ССМЗ 349.9 12 330.1 20 535.6 14 691 75
nnQTt ООйГъ 372.8 19 238.5 -13 491.1 4 364 -8
CNRM1 268.3 -14 225.6 -18 568.6 21 337 -15
CSIRO 260.9 -16 248.4 -10 372.7 -21 356 -10
ЕСНАМЗ 201.0 -35 210.7 -23 488.7 4 369 -7
GEN2 274.6 -12 139.6 -49 402.9 -14 293 -26
GFDL 260.6 -16 107.6 -61 416.4 -12 146 -63
GISS 110.6 -65 131.1 -52 308.1 -35 226 -43
LMD5 673.8 116 208.3 -24 330.2 -30 162 -59
MRI2 255.4 -18 165.9 -40 425.4 -10 289 -27
UGAM 376.1 20 264.8 -3 450.1 -4 452 14
икмо 350.6 12 160.6 -41 496.3 5 357 -10
UIUC 521.5 68 228.0 -17 434.0 -8 274 -31
YONU 586.2 87 301.0 10 414.8 -12 452 14
АНС1 354(±116) Ц (±38) 207(±67) -25(±32) 439(±71) -7 (±15) 390±189 -2±46
ЛУЧШ2 308(±37) -1(±12) 263(±34) -4(±13) 454(±49) -4(±11) 361±118 —9±11
НАБЛ3 312 - 274 - 471 - 397 -
Примечания.
1 АНС — среднее значение по ансамблю всех МОЦА. В скобках указаны значения межмодельного разброса величин стока (±ст).
2 ЛУЧШ — среднее значение по ансамблю тех МОЦА, которые наиболее точно воспроизводят современные объемы годового стока (выделены жирным шрифтом и курсивом).
3 НАБЛ — средние значения годового стока, полученные по результатам наблюдений. В скобках указаны значения межгодового среднеквадратичного разброса стока [5, 7].
стока с тех или иных территорий в значительной степени можно считать случайным. Это означает, что нет гарантии, что при возмущении граничных условий МОЦА верно воспроизведет гидрологический режим исследуемой территории. Во-вторых, осреднение результатов по всем 16 моделям дает худший результат, чем осреднение по “лучшим” моделям общей циркуляции атмосферы (не выходящим за ошибку 20%). Исключение составляет лишь водосбор реки Оби, но здесь осреднение по всем моделям дает слишком большую величину межмодельного разброса.
2. Климат и речной сток в контрастных климатических условиях
В экспериментах по воспроизведению климата позднеплейстоценового криохрона, выполненных в рамках PMIP, заданы отрицательные аномалии температуры поверхности океана (в соответствии с данными реконструкций CLIMAP [8]), а также распределение ледяных щитов (с соответствующими вариациями альбедо и высоты над уровнем моря) согласно данным реконструкций, представленным в работе [9]. Содержание С02 в атмосфере взято равным 200 ppm [10]. Остальные граничные условия, в том числе солярный климат, мало отличались от современного.
Откликом на изменения факторов климатообразования стало формирование в Северном полушарии значительно более холодного климата по сравнению с современным. Распределение аномалий температуры обусловлено в первую очередь наличием скандинавского ледникового щита. На Восточно-Европейской равнине средняя годовая температура была в среднем на 10-12 °С ниже современной, а в Западной Сибири — на 8-10 °С. Межмодельный разброс среднегодовой аномалии температуры воздуха на обеих равнинах °
статистически достоверной. На территории восточно-европейского и западно-сибирского регионов 21 тыс. лет назад осадков выпадало в среднем на 30-40% меньше, чем в настоящее время. Подробнее с результатами моделирования климата позднеплейстоценового криохрона можно ознакомиться в [2, 11, 12].
В табл. 3 приведены результаты расчета стока с водосборов Каспийского, Черного и Балтийского морей, а также с бассейна реки Оби. Хорошо видно, что в целом годовой сток в перечисленные выше водосборы снизился примерно вдвое. Согласно данным палеореконструкций, суммарный сток, например, в Каспийское море около 20 тыс. календарных лет назад мог составлять 20-25% от современного [13]. В [14] показано, что сток рек южного склона Восточно-Европейской равнины и сток с Западно-Сибирской низменности в эпохи последнего оледенения уменьшался на 20-25%. Таким образом, результаты моделирования лишь качественно соответствуют данным реконструкции, количественно различаясь в 1.5-2 раза. Тем не менее они позволяют выявить генезис резкого уменьшения речного стока с равнинных территорий Северной Евразии. Причина этого процесса — в существенном уменьшении годовых сумм осадков, которое, в свою очередь, является следствием глобального понижения температуры атмосферы и поверхности океана. В [2, 15] также оценен уровень регрессионного состояния Каспийского моря, величина которого по данным моделирования достигала 50 м, что неплохо соответствует данным палеореконструкций [10, 13, 16, 17].
Перейдем к рассмотрению результатов численного моделирования климата и гидрологического режима теплой эпохи оптимума голоцена. На фоне граничных условий, мало
отличавшихся от современных, происходили существенные сезонно-широтные изменения инсоляции на внешней границе атмосферы (эффект Миланковича), оказавшиеся основным фактором, управляющим изменениями климата Земли в течение этих последних
Таблица 3. Воспроизведение моделями годового стока с водосборов Черного и Каспийского морей, а также бассейна реки Оби 21 ООО лет назад
Модель общей циркуляции атмосферы Черное море Каспийское море Река Обь
Сток, Q КМ Отклонение, % Сток, Q КМ Отклонение, % Сток, Q КМ Отклонение, %
ССС2 218.7 20 78.0 -68 335 -46
CCSR 138.1 -60 105.2 -45 12 -97
ЕСНА 246.5 53 157.1 -17 184 -50
GEN2 134.4 -60 54.6 -71 90 -69
LMD5 575.2 -12 209.5 11 297 83
MRI2 108.1 -54 117.2 -38 253 -13
UGAMP 204.4 -39 83.0 -64 116 -74
АНС 232.2 -22 114.9 -42 184 -38
ЛУЧШ 171 -45 80.5 -66 95 -72
Примечание. Дано отклонение от современного объема стока, воспроизведенного моделями.
Таблица 4. Воспроизведение моделями годового стока с водосборов европейской Арктики, Черного и Каспийского морей, а также бассейна реки Оби 6 ООО лет назад
Модель общей циркуляции атмосферы Черное море Каспийское море Балтийское море Река Обь
Сток, Q КМ Откло- нение, % Сток, Q КМ Откло- нение, % Сток, Q КМ Откло- нение, % Сток, Q КМ Откло- нение, %
BMRC 346.7 290 336.1 220 564.5 224 564.5 224
ССС2 221.9 3 278.6 5 425.7 -3 425.7 -3
ССМЗ 382.7 9 391.7 19 419.1 -1 419.1 -1
nnQTt ООйГъ 317.1 -15 207.2 -13 311.7 1 311.7 1
CNRM1 202.3 -25 195.2 -14 404.9 -1 404.9 -1
CSIRO 76.5 -71 121.1 -51 217.5 О 1 217.5 О 1
ЕСНАМЗ 145.0 СО 1 171.7 -18 445.8 -4 445.8 -4
GEN2 296.2 8 136.1 -3 193.9 -13 193.9 -13
GFDL 263.8 1 111.8 4 249.8 0 249.8 0
GISS 93.2 -16 154.1 18 315.3 14 315.3 14
LMD5 1088.4 61 414.8 100 315.3 56 315.3 56
MRI2 255.1 0 152.1 -8 238.5 -23 238.5 -23
UGAM 367 -2 256.2 -3 322.8 -21 322.8 -21
икмо 379.2 8 165.9 3 293.7 -14 293.7 -14
UIUC 526.7 1 156.4 со 1 243.8 -8 243.8 -8
YONU 594.9 1 297.3 1 258 -2 258 -2
АНС 354 (±116) 14 (±38) 220(±47) 14 (±30) 326(±49) 12 (±29) 326(±49) 12 (±29)
ЛУЧШ 320(±30) -2 (±6) 255(±37) -5 (±8) 404(±18) -6 (±4) 404(±18) -6 (±4)
10 тыс. лет. 6 тыс. календарных лет назад они составили +5 % летом и —5% зимой [2] (в среднем по Северному полушарию).
Согласно результатам численного моделирования, ощутимым отличием климата атлантического оптимума от современного является более теплое лето (па 2-3 °С теплее современного) при практически не отличающейся от современной зиме. Подробно о причинах несимметричной сезонной реакции температуры речь идет в [11]. Результаты моделирования показали, что осадков на территории Восточно-Европейской равнины и ЗападноСибирской низменности в эпоху атлантического оптимума выпадало примерно на 5 % меньше, чем в настоящее время.
В табл. 4 приведены результаты расчетов речного стока с главных водосборов ВосточноЕвропейской равнины и Западно-Сибирской низменности. Хорошо видно, что осреднение между всеми 16 моделями дает статистически незначимое завышение речного стока, а осреднение результатов нескольких “лучших МОЦА”, набор которых определяется для каждого бассейна в отдельности, — слабое занижение стока со всех водосборов. В целом можно сделать вывод о том, что по результатам численного моделирования речной сток 6 тыс. календарных лет назад мало отличался от современного. Получившиеся по результатам модельных расчетов аномалии стока не выходят за пределы межмодельного разброса и естественной изменчивости, составляющей 10-20%. Полученные данные соответствуют принципиальной оценке, основанной на результатах палеореконструкций. Результаты большинства эмпирических палеореконструкций сводятся к тому, что, например, волжский сток 6 тыс. календарных лет назад колебался в пределах ±10 % от современного, т. е. мало отличался от наблюдаемого в современных условиях [16-20]. Данные реконструкций стока реки Оби 6 тыс. календарных лет назад для сравнения с результатами моделирования практически отсутствуют. Однако, учитывая общую направленность климатических изменений на территории Восточно-Европейской равнины и Западно-Сибирской низменности в эпоху оптимума голоцена, нет оснований предполагать, что средний многолетний сток Оби существенно отличался от современного. Правда, в некоторых работах, например в [18], предполагается, что годовая сумма осадков на территории Западно-Сибирской равнины и Восточно-Европейской низменности была на 20-30 % выше, чем в современных условиях.
Заключение
Анализ полученных результатов позволяет сделать следующие выводы.
Оценка качества моделирования показала, что современные модели общей циркуляции атмосферы способны более или менее успешно воспроизводить лишь средний многолетний гидрологический режим тех рек, площадь водосбора которых не меньше волжского (1.3 млн км2). При этом использование полного ансамбля МОЦА нецелесообразно — по-видимому, следует отбирать те модели, которые наилучшим образом воспроизводят годовой сток. То же самое относится и к климатическим архивам осадков и испарения.
Показано, что резкое снижение речного стока с Восточно-Европейской равнины и Западно-Сибирской низменности, а также уровня Каспийского и Черного морей 21 тыс. календарных лет назад имеет климатическую природу, т. е. связано с резким уменьшением годовых сумм осадков. Полученные данные неплохо согласуются с результатами палеореконструкций.
Нарастание температуры не всегда сопровождается ростом годовых сумм осадков и речного стока. Такой вывод можно сделать, основываясь на результатах численного моделирования теплого события голоцена 6 тыс. календарных лет назад. Так, увеличение среднегодовой температуры происходило на фоне практически неизменных условий увлажнения.
Список литературы
[1] JOUSSAUME S., Taylor К. Status of the paleoclimate modeling intercomparison project (PMIP) // Proc. of the First Intern. AMIP Sci. Conf., Monterey, California, May 15-19. 1995. / W.L. Gates (Ed.). WCRP 92. WMO/Technical Document 732. Geneva, 1995. P. 425-430.
[2] ТОРОПОВ П.А. Температурный режим и условия увлажнения Восточно-Европейской равнины в контрастных климатических условиях: Автореф. дис. ... канд. геогр. наук. М.: Изд-во МГУ, 2006. 28 с.
[3] Berger A. Long-term variations of daily insolation and Quaternary climatic changes // J. Atmos. Sci. 1978. Vol. 35. P. 390-403.
[4] ТОРОПОВ П.А. Верификация результатов численного моделирования климата ВосточноЕвропейской равнины // Метеорология и климатология. 2005. № 5. С. 19-30.
[5] Косарев, Тужилкин, Даниялова, Архипкин. Гидрология и экология Черного и Каспийского морей // География, общество, окружающая среда. Т. VI: Динамика и взаимодействие атмосферы и гидросферы. М.: Городец, 2004. С. 218-282.
[6] Евстигнеев В.М., Евсеева Л.С. Климатическая модель изменчивости годового стока // Глобальные изменения природной среды (климат и водный режим). М.: Наука, 2000. С. 211-219.
[7] Иванов В.В., Муждаеа О.В., Соловьева З.С. Многолетние изменения годового и сезонного притока речных вод в арктические моря // VI Всерос. гидролог, съезд: Тез. докл. СПб.: Гидрометеоиздат, 2004. С. 39-41.
[8] CLIMAP project members // Science. 1976. Vol. 191. P. 1131-1137.
[9] Peltier W.R. Ice age paleotopography // Science. 1994. Vol. 265(5169). P. 195-201.
[10] Raynaud D., Jousel J., Barnola J.M. et al. The ice record of greenhouse gases // Science. 1993. Vol. 259(5097). P. 926-934.
[11] КИСЛОВ А.В. Климат в прошлом, настоящем и будущем. М.: Наука/Интерпериодика, 2001. 352 с.
[12] Кислов А.В., ТОРОПОВ П.А. Речной сток с Восточно-Европейской равнины при различных состояниях климата // Водные ресурсы. 2006. № 5. С. 28-41.
[13] Варущенко С.И., Варущенко А.И., Клиге Р.К. Изменение режима Каспийского моря и бессточных водоемов в палеовремени. М.: Наука, 1987. 255 с.
[14] Палеогеография Северной Евразии в позднем плейстоцене-голоцене и проблемы географического прогноза. М.: Наука, 1978. 76 с.
[15] Kislov A., Toropov P. Paleo sea-level changes in the Black and Caspian Seas: Links to river runoff and global climate change // PAGES News. 2006. Vol. 14, N 1. P. 36-37.
[16] Рычагов Г.И. Уровенный режим Каспийского моря за последние 10 000 лет // Вест. МГУ. Сер. География. 1993. № 2. С. 38-49.
[17] Свиточ А.А. Морской плейстоцен побережий России. М.: ГЕОС, 2003. 362 с.
[18] Борзенкова И.И. Эмпирическая палеоклиматология: состояние проблемы и методы исследования // Изменения климата и их последствия: Сб. науч. тр. СПб.: Наука, 2002. С. 75-92.
[19] Борзенкова И.И., Лемешко Н.А. Водный баланс бассейна Волги в начале XXI века (на основе палеоклиматических сценариев) // Метеорология и гидрология. № 7. С. 54-60.
[20] Tshepalyga A. Rapid inundations of the Ponto-Caspian shelves, their origin and impact on early man settlement // The Fourth Intern. Meeting on Global Continental Palaeohydrology, Moscow, 2000. P. 67-70.
[21] Величко А.А., Климанов B.A., Беляев А.В. Реконструкция стока Волги и водного баланса Каспия в оптимумы микулинского межледниковья и голоцена // Изв. РАН. Сер. География. 1988. № 1. С. 27-37.
[22] Тарасов П.Е., Гунова B.C., Успенская О.Н. Уровни озер бассейна Волги как индикатор изменения климата в голоцене // Вест. МГУ. Сер. География. 1997. № 3. С. 14-24.
[23] Harrison S., Digerfeldt G. European lakes as paleohydrological and paleoclimatic indicators // Quat. Sci. Rev. 1993. N 12. P. 233-248.
Поступила в редакцию 9 ноября 2006 г.
