Портал | Содержание | О нас | Авторам | Новости | Первая десятка | Дискуссионный клуб | Чат Научный форум --> Первая десятка "Русского переплета" Темы дня:

Президенту Путину о создании Института Истории Русского Народа. |Нас посетило 40 млн. человек | Чем занимались русские 4000 лет назад?

| Кому давать гранты или сколько в России молодых ученых?

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

ПРИ ВОЗМОЖНОМ ГЛОБАЛЬНОМ ПОТЕПЛЕНИИ КЛИМАТА

Э. Д. ЕРШОВ

Московский государственный университет

им. М.В. Ломоносова

Динамика мерзлоты определяется как естественноисторическим ходом развития природной обстановки, так и изменениями климата за счет воздействия антропогенного фактора, которое может проявляться в виде глобального потепления за счет, например, увеличения содержания парниковых газов в атмосфере. Анализ публикаций, посвященных реакции мерзлоты на возможное изменение климата, показывает, что часто для этой цели используют качественные методы прогноза (метод аналогий, экспертных оценок). В работах, где применяют количественные методы, не всегда четко сформулированы математические модели, а основной трудностью обычно является задание корректных начальных климатических и геокриологических условий.

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

В настоящее время существуют несколько сценариев возможного изменения климата под влиянием антропогенных факторов [1, 3]. Нами использован сценарий, предложенный сотрудниками Института глобального климата и экологии РАН (ИГКЭ), согласно которому на территории криолитозоны России к 2100 году среднегодовая температура приземного воздуха повысится на 4-8?C в зависимости от географического положения района.

Закономерности деградации многолетнемерзлых пород (ММП) могут быть представлены в виде пространственно-временной изменчивости трех характеристик: время начала оттаивания ММП, повышение их среднегодовой температуры, понижение кровли ММП. Необходимые для таких оценок долгосрочные изменения температуры ММП требуют решения нестационарной задачи теплопроводности (задачи Стефана) с меняющимися во времени граничными условиями. Такая задача точного аналитического решения не имеет, ее реализация возможна только численными методами с использованием ЭВМ. Методика выполненных прогнозных расчетов изложена в [2].

Выбранные для исследования пункты (метеостанции с относительно длинными рядами наблюдений) примерно равномерно расположены на территории криолитозоны. В каждом пункте назначались несколько типов разреза пород. Так, для равнинных районов криолитозоны выбраны пески, суглинки и торф, подстилаемый песками или суглинками. В соответствии со строением разреза, составом и криогенным строением пород назначали необходимые для расчетов теплофизические характеристики. Для возможности сопоставления результатов расчетов свойства каждого типа пород приняты одинаковыми для всей криолитозоны. Для рыхлых пород коэффициент теплопроводности талого песка - 1,6, мерзлого - 1,8, талого суглинка - 1,2, мерзлого - 1,4, талого торфа - 0,6, мерзлого - 0,9 ккал/м.ч. " ?C, объемная теплоемкость талого песка - 660, мерзлого - 490, талого суглинка - 720, мерзлого - 520, талого торфа - 660, мерзлого - 490 ккал/м3 " ?C, затраты тепла на фазовые переходы песка - 16 000, суглинка - 25 000, торфа - 50 000 ккал/м3. Для горных районов приняты теплофизические свойства, характерные для скальных пород выбранных реперных пунктов.

При прогнозном математическом моделировании были приняты следующие допущения. В течение расчетного периода (110 лет) повышение среднегодовой температуры воздуха происходит по линейному закону; теплофизические характеристики растительного и снежного покровов, а также теплофизические свойства верхнего слоя пород остаются неизменными; не учитывались возможная осадка оттаивающих ММП и развитие криогенных процессов. Расчетная область задавалась на 30-50 м глубже современного положения подошвы ММП. В процессе расчетов температурное поле фиксировалось по глубине расчетной области через 5-10 лет. Результаты моделирования, выполненного при таких допущениях, позволили выявить некоторые пространственно-временные закономерности процесса деградации криолитозоны.

Период времени (t) от начала потепления климата до момента, когда среднегодовая температура ММП повысится до 0?C (в дальнейшем начнется прогрессирующее оттаивание ММП), оказался прямо пропорционален среднегодовой температуре ММП на момент начала потепления климата (t0) и обратно пропорционален тренду потепления (Dt) и может быть рассчитан по уравнению t = A| t0 |/ Dt, где A - коэффициент, который зависит от тренда повышения температуры воздуха и может быть получен из графика зависимости времени начала оттаивания ММП от их начальной среднегодовой температуры для каждого темпа потепления климата. На рис. 1 в качестве примера приведены графики для пунктов, различающихся строением разреза (торф, суглинок, песок), ландшафтными и климатическими условиями (использованы прогнозные расчеты по пунктам Нум-То, Надым, Тазовский, Воркута, Амдерма), но имеющих один тренд потепления. Расчеты показали, что при Dt = 0,06?C/год A = 1,13, а при Dt = 0,045?C/год A = 2,48.

Изменение среднегодовой температуры ММП в будущем существенно зависит от их начальной температуры (рис. 2). При относительно высоких начальных температурах ММП большая часть поступающего в грунт тепла идет на фазовые переходы при опускании кровли ММП, поэтому повышение температуры ММП относительно медленное. При низких исходных температурах ММП интенсивность повышения их температур увеличивается. Моделирование также показало, что при тренде потепления 0,06?C/год скорость многолетнего оттаивания мерзлого торфа составит около 6 см/год, суглинка - 13, песка - 20 см/год (рис. 3).

Выявленные закономерности позволяют установить количественные оценки интенсивности деградации ММП при их разной исходной среднегодовой температуре на любой момент времени от начала повышения температуры воздуха для разных пород. Результаты такой оценки на 50-й и 110-й год от начала потепления климата для некоторых типов грунтов представлены в табл. 1.

Через 50 лет от начала потепления в южной части криолитозоны, где современные температуры ММП не опускаются ниже - 2?C, начнется повсеместное разобщение на глубине сезонных и многолетнемерзлых толщ (формирование неслитых ММП). Слитые толщи мерзлых пород (ММП сливающегося типа) сохранятся на территориях с их современной среднегодовой температурой ниже - 3?C. Граница залегающих с поверхности ММП сместится на север приблизительно на 350-400 км, то есть примерно до широты Полярного круга. Южнее этой границы возникнет зона со сложным чередованием слитых и неслитых ММП: при их исходной температуре -1 ╙ - 2?C кровля ММП начнет опускаться, при начальной температуре ниже - 2?C сохранятся слитые ММП.

Состояние криолитозоны Западной Сибири через 110 лет от начала потепления климата отражено на рис. 4. При составлении карты использованы границы и градации среднегодовых температур ММП, показанные на геокриологической карте России масштаба 1 : 2 500 000. На карте выделены четыре области с различными прогнозируемыми геокриологическими условиями.

1. Область, где ныне существующие ММП за расчетный период полностью протают. Оттаивание ММП здесь начнется через 5-10 лет от начала потепления. Среднегодовые температуры пород в этой области будут по всему разрезу положительные. Северная граница области будет приблизительно соответствовать современной границе распространения высокотемпературных (не ниже - 0,5?C) ММП.

2. Область, где практически повсеместно ММП оттают с поверхности. В зависимости от интенсивности повышения температуры воздуха оттаивание ММП начнется через 10-50 лет от начала потепления воздуха. В разрезе выделятся два горизонта пород: талые со среднегодовой температурой от 0 до 2?C и подстилающие их ММП, температура которых меняется по глубине от 0 до - 0,5?C. Эта область соответствует современному распространению ММП со среднегодовой температурой от - 0,5 до - 5?C.

3. Область, где будут преобладать ММП со среднегодовой температурой до - 2,5?C и встречаться участки, где, как и во второй области, произойдет оттаивание ММП с поверхности. Оттаивание ММП на таких участках начнется не ранее чем через 80 лет от начала потепления. Южная граница этой области является прогнозной границей распространения ММП на 110-й год от начала потепления. Граница с поверхности ММП сместится за этот период на 700-800 км к северу.

4. Область сплошного распространения ММП. Она соответствует современному распространению ММП со среднегодовой температурой ниже - 7?C. Среднегодовая температура ММП этой области, за исключением таликовых зон под крупными реками и озерами, останется отрицательной, но на 4-6?C выше, чем в современных условиях.

Во всех областях талые, протаивающие и мерзлые породы будут характеризоваться нестационарным температурным режимом. Их температура по разрезу, глубины сезонного промерзания, сезонного и многолетнего оттаивания будет меняться во времени на протяжении всего прогнозного периода. Итогом проведенных исследований явилась "Прогнозная геокриологическая карта России на 110-й год_" масштаба 1 : 25 000 000 (рис. 5). Масштаб карты позволяет отразить закономерности изменения геокриологической обстановки в генерализованном виде.

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

ЛИТЕРАТУРА

1. Антропогенные изменения климата / Под ред. М.И. Будыко, Ю.А. Израэля. Л.: Гидрометеоиздат, 1987. 406 с.

2. Ершов Э.Д., Максимова Л.Н. и др. Реакция мерзлоты на глобальные zизменения климата // Геоэкология. 1994. ╧ 5.

3. Парниковый эффект, изменения климата и экосистемы. Л.: Гидрометеоиздат, 1989.

* * *

Эдуард Дмитриевич Ершов, доктор геолого-минералогических наук, профессор, зав. кафедрой геокриологии геологического факультета МГУ, заслуженный деятель науки РСФСР, лауреат Государственной премии Российской Федерации в области науки и техники и премии им. М.В. Ломоносова первой степени. Автор более 300 научных статей, 15 монографий и восьми учебников и учебных пособий в области теплофизики, физикохимии и петрографии мерзлых пород, региональной, динамической, литогенетической и исторической геокриологии.