Сбалансированная модель глобального биогеохимического круговорота углерода

ИНСТИТУТ СИСТЕМНЫХ ИССЛЕДОВАНИИ ЛЕСА

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

Для разработки долгосрочного прогноза изменения природной пожарной опасности для зоны тайги и смешанных лесов европейской части России планируется использовать разработанный ранее комплекс программ FORRUS-S (Forest of Russia - STAND), предназначенный для имитационного моделирования и анализа динамических процессов, протекающих в лесных массивах. Исходной информацией для модели служат повыдель-ные данные таксационных описаний и планы лесных насаждений. Размеры моделируемых FORRUS-S лесных массивов составляют до нескольких десятков тысяч гектаров; шаг моделирования - 5 лет. В ходе моделирования прогнозируется изменение средних таксационных характеристик древостоя (высоты, диаметра, возраста, запаса и др.) по ярусам, изменение породного и возрастного составов лесного выдела - первичной учетной единицы лесоустройства [4, 5].

На основе предложенной концептуальной модели будет разработана математическая и имитационная модель динамики природной пожарной опасности, которая в совокупности с моделью FORRUS-S, позволит получить прогноз изменения природной

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

Библиографический список

1. Классификация природной пожарной опасности лесов и классификация пожарной опасности в лесах по условиям погоды, а также требований к мерам пожарной безопасности в лесах в зависимости от целевого назначения лесов, показателей природной пожарной опасности лесов и показателей пожарной опасности в лесах по условиям погоды. Приказ Министерства сельского хозяйства Российской Федерации (Минсельхоз России) от 16 декабря 2008 г., № 532, г. Москва.

2. Чумаченко, С.И. Модель долгосрочного прогнозирования динамики показателей природной пожарной опасности для зоны тайги и смешанных лесов европейской части России. Основные параметры модели / С.И. Чумаченко, Д.Н. Маюк // Науч. труды МГУЛ. Электронный журнал, 2012. - Вып. 10.

- http://www.msfu.ru/joumal.

3. Смирнова, О.В. Концептуальная модель динамики напочвенного покрова / О.В. Смирнова, С.И. Чумаченко // Вестник МГУЛ-Лесной вестник, 2012.

- № 9.

4. Чумаченко, С.И. Моделирование динамики дре-востоев с учетом лесохозяйственного воздействия / С.И. Чумаченко, В.В. Сысуев, М.М. Паленова, М.А. Бредихин, В.Н. Коротков / Тр. VII ежегодной конференция МАИБЛ. Устойчивое развитие боре-альных лесов. - М., 1997. - С. 184-190.

5. Чумаченко, С.И. Имитационное моделирование влияния лесохозяйственных воздействий на лесные экосистемы. Мониторинг биологического разнообразия лесов России: методология и методы / С.И. Чумаченко, М.М. Паленова, В.Н. Коротков, С.В. Починков / отв. ред. А.С. Исаев, ЦЭПЛ РАН.

- М.: Наука, 2008. - С. 314-328.

СБАЛАНСИРОВАННАЯ МОДЕЛЬ ГЛОБАЛЬНОГО

биогеохимического круговорота углерода

В.Д. БУРКОВ, проф., д-р техн. наук, МГУЛ,

В.Ф. КРАПИВИН, проф., д-р физ.-мат. наук, Фрязинское отделение Института радиотехники и электроники им. В.А. Котельникова РАН,

В.С. ШАЛАЕВ, проф., директор ИСИЛ, д-р техн. наук

В соответствии с Киотским Протоколом, ориентированным на предотвращение потепления климата Земли вследствие парникового воздействия антропогенных выбро-

burkov@mgul.ac.ru

сов углекислого газа (СО2), возникает настоятельная необходимость уточнения моделей циркуляции углерода в окружающей среде с целью получения более объективных оценок

86

ЛЕСНОЙ ВЕСТНИК 9/2012

ИНСТИТУТ СИСТЕМНЫХ ИССЛЕДОВАНИИ ЛЕСА

степени и характера данного воздействия. Это вызвано тем, что представления о роли атмосферного СО2 в наблюдаемых изменениях климата основаны на весьма неточных и неполных схемах, описывающих потоки углерода между его основными резервуарами на планете. Анализу возникающих здесь проблем посвящены многие исследования [1-27]. Одной из них является разработка сбалансированной схемы глобальных потоков углерода. Ведь примерно 27 % стока углерода из атмосферы в имеющихся схемах остаются необъясненными.

Исследования Кондратьева и др. [9, 10, 22-24] показали, что учет пространственной неоднородности стоков и источников атмосферного углерода и уточнение параметров его потоков на границах атмосфера-суша и атмосфера-океан повышают точность оценок, но не устраняют ни указанного дисбаланса, ни расхождений между потоками углерода и динамикой изменения парциального давления углекислого газа в земной атмосфере.

Реализуемая в рамках глобального углеродного проекта GCP (Global Carbon Project) программа контроля потоков углерода в условиях развития современной цивилизации, к сожалению, пренебрегает этим противоречием. Это в значительной мере снижает доверие к существующим сценариям ожидающих Землю природных катаклизмов, вызванных таянием арктических и антарктических льдов, и искажает реальную картину антропогенного влияния на климатические процессы.

Наиболее принципиальный недостаток программы GCP состоит в том, что она не учитывает геологическую компоненту современного круговорота углерода [2]. По имеющимся оценкам, весьма приближенным, в породах земной коры и мантии (литосфере) сосредоточено ~ 1023 г углеродсодержащих газов [11, 14]. Находящаяся в них масса углерода в ~ 104 раз превышает его количество, присутствующее сегодня в биосфере (над земной поверхностью). Между биосферой и литосферой происходит постоянный очень интенсивный обмен углеродом, носящий са-моорганизующий характер. По данным ис-

следований [1, 2], в результате действия принципа Ле Шателье содержание подвижного углерода в системе стремится к устойчивому соотношению

п/т. = C = const,

где п и т. - соответственно общее количество и среднее время пребывания подвижного углерода (во всех формах) в i-м резервуаре системы;

С - константа, характеризующая скорость круговорота углерода в системе.

Для современного состояния круговорота углерода в биосфере С = (2,7 ± 0,05) 1017 г СО2/год.

Если условие (1) выполняется, то уход углерода из одного какого-либо резервуара системы вследствие круговорота пополняется его поступлением из других резервуаров. Если условие (1) не выполняется, то в системе возникают некомпенсированные перетоки вещества, которые возвращают ее к равновесию.

В процессе глобального круговорота подвижный углерод многократно пересекает земную поверхность. Над поверхностью Земли, играющей роль «геохимического барьера», он циркулирует преимущественно в окисленном виде (СО2), а под этой поверхностью - главным образом в восстановленной форме (СН4). Это обстоятельство требует совместного изучения глобальных круговоротов СО2 и СН4 [10].

Среди всех факторов, управляющих перераспределением углерода над (и под) поверхностью планеты, наиболее важную роль в современную эпоху играют два процесса, носящих в основном региональный характер:

1) перенос углерода метеогенными водами и

2) хозяйственная деятельность человека. Благодаря первому фактору избыточный углерод на континентах поступает с метеогенными водами под земную поверхность, где, восстанавливаясь до углеводородов (УВ), участвует в процессах современного нефтегазообразования [3]. Вследствие этого промышленные месторождения нефти и газа размещаются, во-первых, в пределах крупных осадочных бассейнов, дренирующих огромные по площади территории и, во-вторых, тяготеют к

ЛЕСНОЙ ВЕСТНИК 9/2012

87

ИНСТИТУТ СИСТЕМНЫХ ИССЛЕДОВАНИИ ЛЕСА

зонам разломов земной коры. Наличие разломов, с одной стороны, облегчает проникновение метеогенных вод под поверхность Земли, а с другой - способствует разгрузке этих вод от транспортируемого ими углерода [2, 4].

Фактические данные свидетельствуют [16], что заметное пополнение месторождений нефтью и газом происходит за 10-50 лет. Столь малое время формирования залежей обусловлено высокой скоростью переноса углерода через земную поверхность водами [2,

3]. По измерениям, вероятно заниженным, эта скорость в масштабе всей планеты составляет ~ 5 1015 гС/год [6, 7]. Данный процесс крайне неравномерен [15]. Основная масса углерода опускается «вниз» в пределах материков и на их окраинах. А поступает «наверх» в зонах срединных океанических хребтов, по континентальным разломам, а также при многочисленных вулканических извержениях на океаническом дне. Количественно эти восходящие и нисходящие потоки углерода на континентах и особенно в океанах изучены, к сожалению, пока еще крайне слабо [7]. Но это говорит о том, что глобальный круговорот углерода нельзя изучать изолированно также и от планетарного круговорота воды [9].

Антропогенный фактор перераспределения углерода между наземными и подземными резервуарами проявляется в добыче и использовании нефти (4 млрд т/год) и природного газа (350 трлн м3/год). Тем самым суммарная масса потребляемого ископаемого углерода составляет ~ 7 1015 гС/год, что сопоставимо с его количеством, поступающим на поверхность за счет естественной циркуляции. Поэтому хозяйственная деятельность людей, неразрывно связанная с этими объемами утилизации УВ, наряду с воздействием на климат способна во многом влиять и на региональные процессы генерации нефти и газа в недрах [4].

В проекте GCP данное обстоятельство также необходимо учитывать. Современная практика транспортировки нефти и газа на многие тысячи километров от мест добычи может привести к существенному перераспределению мировых ресурсов углеводородов, причем в ближайшие десятилетия. Интенсив-

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

С учетом высказанных выше соображений за основу более корректного изучения парникового эффекта может быть взята общая схема круговорота углерода [10], дополненная «геологическим» блоком (рис. 1). Эта схема наряду с описанием основных потоков СО2 в биосфере учитывает также возможность накопления подвижного углерода под земной поверхностью, чем замыкается и становится уравновешенной.

Базовая модель

Схема рис. 1 уравновешивает все потоки углерода, отражая обмены углеродом на границах атмосферы с земными покровами и акваториями морей и океанов и учитывая роль дегазационных процессов и движение метеогенных вод через земную поверхность. Согласно [2] глобальный круговорот углерода включает три основных круга его циркуляции, отличающихся временем обращения углерода между его резервуарами. Для точного расчета количественных оценок всех потоков углерода необходимо определить пространственные и предметные распределения акваторий, земных покровов и геосферы с указанием их собственных характеристик. Современный уровень знаний позволяет решить эту задачу лишь частично и, тем самым, уменьшить интервал неопределенности в оценках будущего изменения климата. Однако используемая при этом модель глобального круговорота углерода должна отражать не только пространственную мозаику его резервуаров, стоков и источников, но и обеспечивать динамический расчет их мощностей. Сделанные до настоящего времени расчеты с помощью подобных моделей не в полной мере учитывали сведения о состоянии и классификации земных покровов, а тем более об изменчивости акваторий Мирового океана. Поэтому формирование схемы пространственного распределения почвенно-растительных формаций и квазиоднородных акваторий Мирового океа-

88

ЛЕСНОЙ ВЕСТНИК 9/2012

ИНСТИТУТ СИСТЕМНЫХ ИССЛЕДОВАНИИ ЛЕСА

Рис. 1. Блок-схема глобального биогеохимического круговорота углерода в системе атмосфе-ра-суша-океан-геосфера. Описание резервуаров и потоков углерода дается в табл. 1

на позволит преодолеть указанные недостатки предшествующих моделей.

Система балансовых уравнений для схемы рис. 1 с учетом обозначений таблице имеет вид

ggs((p Xz,t) dt

jeQs meQs

где S - резервуар углерода (A - атмосфера, U - верхний фотический слой океана над термоклином, L- суша, I - промежуточный фотический слой океана под термоклином, D - глубокие слои океана, N - придонный слой океана, G - геосфера, );

ф - широта, X - долгота, z - глубина, t - время;

H.s - поток углерода из j-го резервуара в резервуар S;

HSm - сток углерода из резервуара S в m-й резервуар;

QS - множество резервуаров углерода, граничащих с резервуаром S;

aS - запас углерода в резервуаре S.

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

ЛЕСНОИ ВЕСТНИК 9/2012

89

ИНСТИТУТ СИСТЕМНЫХ ИССЛЕДОВАНИИ ЛЕСА

Таблица

Потоки углерода в форме СО2 в биосфере, учитываемые в имитационной модели глобального биогеохимического круговорота углерода, блок-схема которой представлена на рис. 1

Поток углерода Идентификатор Оценка потока

Выделение при сжигании растительности hv 6,9

Выделение при сжигании ископаемого топлива H GA 3,6

Десорбция HUA HUA 97,08

Сорбция Hau 100

Выветривание горных пород H-al 0,04

Вулканические эманации H2 GA 2,7

Поглощение наземной растительностью ^AL 224,4

Дыхание растений H2la 50-59,3

Дыхание людей H la 0,7

Дыхание животных H4la 4,1

Выделение при разложении гумуса почв H la 139,5

Выделение корнями растений H6la 56,1

Жизнедеятельность населения H OL 0,3

Жизнедеятельность животных H2ol 3,1

Отмирание растительности H OL 31,5-50

Выпадение в донные осадки Hng 0,1-0,2

Растворение морских отложений HgN 0,1

Разложение детрита в фотическом слое океана HOL, H1OI 35

Разложение детрита в глубоких слоях океана Hod 5

Подъем с глубинными водами -^NA HDD HIU 45

Опускание с поверхностными водами и за счет гравитационного оседания TT H H nUP nID nDN 40

Фотосинтез HuO, HIO 69

Подземный сток Hlu 0,5

Поверхностный сток H2lu 0,5-0,6

Дыхание живых организмов в океане H2OU H2OI 25

Геосферный сток углерода из биосферы Hlg 70

Поступление углерода из геосферы в биосферу Hga 17

буют элементы биогеоценотического блока. Для этого всю поверхность суши Z покроем однородной сеткой географических пикселей Z.. = {(ф,А): фг1<ф<ф;; j <A<A.} с границами по широте (ф._1 ,ф.) и долготе (А ,А) и площадью о. Количество пикселей определяется доступной базой данных, т.е. выбором размеров сетки (Аф, АА): i = 1,..., n; n = [180/Аф]; j = 1,.,k; k = [180/АА]. В каждом таком пикселе возможно размещение N типов поверхностей, включая типы почвенно-растительных формаций, акваторий водных систем и антропогенных объектов. Динамика растительного покрова s-го типа подчинена закону

dB /dt = R -M - T,

s s s s’ v '

где Rs - фотосинтез,

M и T - потери биомассы B за счет ее

s s s

отмирания и эвапотранспирации соответственно.

Компоненты правой части уравнения (3) являются функциями характеристик окружающей среды - освещенности, температуры, влажности атмосферы и почвы, концентрации СО2 в воздухе. Существуют различные способы и формы параметризации этих функций. Одной из них является модель Коллатца [21], послужившая базовой для глобальной биосферной модели SiB2 [25]. Температура, влажность и скорость испарения влаги в слое растительности и почве связаны зависимостями с биометрическими параметрами и энергетическими потоками в системе атмосфера_растение_почва. При этом по

90

ЛЕСНОЙ ВЕСТНИК 9/2012

ИНСТИТУТ СИСТЕМНЫХ ИССЛЕДОВАНИИ ЛЕСА

аналогии с электростатикой вводится понятие «сопротивления» и потоки рассчитываются по простой формуле: поток = разность потенциалов/сопротивление. Модель SiB2 учитывает потоки собственного и скрытого тепла через испарение водяного пара для растительности и почвы, а потоки СО2 разделяются по классам растительности С3 и С4, что существенно повышает точность параметризации функций правой части уравнений (3). Согласно [21], три фактора регулируют функцию R : эффективность фотосинтетической ферментной системы, количество фотосинтетически активной радиации (ФАР), поглоще-ниное хролофиллом клетчатки, и способность растительного вида к усвоению и экспорту во внешнюю среду продуктов фотосинтеза. Применение принципа Либиха и учет данных о распределении типов растительных покровов по пикселям {£..}, о парциальных давлениях СО2 и О2, о температуре и плотности атмосферы, об уровне освещенности позволяет рассчитать значения потоков Н в уравнении (2) для всех пикселей суши.

Модель круговорота углерода в системе атмосфера-океан подробно описана в работе [23]. Она базируется на той же сетке географических пикселей, но объединенных по зональному принципу согласно классификации Tarko [26]. Вся толща океанской среды рассматривается как единый биогеоценоз, в котором основным связующим фактором является поток органического вещества, продуцируемого в поверхностных слоях и затем проникающего до максимальных глубин океана. Регулятором потоков углерода в этой среде является карбонатная система, параметрическое описание которой дано в работе [24].

Многие параметры глобальной модели круговорота углерода могут измеряться в режиме спутникового мониторинга, что дает возможность применить адаптивную схему расчета характеристик парникового эффекта [24]. Эта схема позволяет обучать модель в непрерывном режиме коррекции ее структуры и параметров. Действительно, спутниковые измерения в видимом и ближнем инфракрасном диапазонах обеспечивают оперативные оценки ФАР и таких характеристик

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

Результаты моделирования

Таким образом, замыкание потоков углерода в окружающей среде с помощью схемы рис. 1 открывает возможность для дальнейших целенаправленных исследований глобального круговорота углерода и повышения достоверности в прогнозах парникового эффекта.

По имеющимся предварительным расчетам роли наземных экосистем и Мирового океана в поглощении атмосферного СО2 выражаются потоками 0,7±1,0 и 2,4±0,3 ГтС/год соответственно, а сток в геосферный резервуар равен 1,8±0,7 ГтС/год. С ростом атмосферной концентрации СО2 ожидается повышение роли наземной биоты и снижение потока углерода на границе атмосфера-океан. Гео-сферный сток углерода слабо зависит от климатических параметров и проявляет высокую инерционность. Также в последующем необходимо учитывать различие резервуаров углерода по характеристикам его подвижности и принимать во внимание механизмы самоорганизации глобального круговорота углерода. К сожалению, GCP не предусматривает изучение подобных механизмов захоронения углерода, а это означает, что глобальный проект в очередной раз не даст желаемых результа-

ЛЕСНОИ ВЕСТНИК 9/2012

91

ИНСТИТУТ СИСТЕМНЫХ ИССЛЕДОВАНИЙ ЛЕСА

Рис. 2. Усредненные по долготе скорости поглощения атмосферного СО2 наземными и океанскими экосистемами при трех сценариях антропогенной эмиссии углерода: 6,8 ГтС/год (пунктирная кривая, 2004 г), 10 ГтС/год (жирная сплошная кривая, прогноз на 2020 г.) и 20 ГтС/год (тонкая сплошная кривая, прогноз на 2020 г.). Обозначение:

H. = H„ + H2,, + H1,, - Hu - Я\, - Hrj - нгл

A AU AL AL UA LA LA LA

Рис. 3. Годовое распределение потока углерода через границу атмосфера-океан в различных широтных зонах

тов. В то же время предлагаемая здесь схема позволяет максимально учесть все аспекты изменчивости глобального круговорота углерода и отразить его взаимосвязь с круговоротами других химических веществ [2, 10].

Численные эксперименты требуют введения определенных предположений относительно функциональных представлений потоков углерода как функций характеристик окружающей среды и задания сценариев антропогенной активности и климатических трендов. Дополнительно к уже ранее введенным в табл. 1 параметризациям потоков определим HLG = CJTBG , Hga = CGIXGA, где CB и CG - запасы углерода в биосфере (4 1019 гС) и

геосфере (1 1013 гС) соответственно, tbg и tga - время жизни углерода в биосфере (600 лет) и геосфере (1,3 млн лет) соответственно.

Результаты моделирования приведены на рис. 2. Видно, что с возрастанием антропогенной нагрузки на глобальный круговорот углерода стабилизирующая роль биосферы сохраняется даже при увеличении к 2020 г. промышленных выбросов углерода в атмосферу в три раза по сравнению с современным периодом. Детальный анализ соотношения потоков углерода показывает, что эта стабилизация достигается в основном увеличением потоков HLG, HAL2 и HAU. Отсюда следует, что для сохранения равновесия климатической

92

ЛЕСНОЙ ВЕСТНИК 9/2012

ИНСТИТУТ СИСТЕМНЫХ ИССЛЕДОВАНИИ ЛЕСА

системы необходимо предотвратить нарастание процессов обезлесивания и загрязнения Мирового океана, особенно в высоких широтах (рис. 3).

Как уже неоднократно отмечалось, внимание к проблеме углекислого газа и другим малым газовым составляющим атмосферы вызвано в основном обеспокоенностью ученых, экономистов и политиков возможными последствиями из-за ожидаемого изменения климата. Поскольку сжигание ископаемого топлива и обезлесивание поверхности суши являются наиболее значимыми источниками антропогенного возрастания концентрации СО2 и СН4 в атмосфере, то именно эти газы изучаются наиболее внимательно. Оба эти газа поглощают инфракрасное излучение земной поверхности и поэтому воздействуют на планетарный климат. Возникшая в конце ХХ столетия перед политиками задача уменьшения этого воздействия нашла решение в виде Протокола Киото. Реальное и объективное решение этой задачи лежит в глубоком изучении биогеохимических круговоротов всех газов, так или иначе связанных с круговоротом углерода. Именно на таком подходе сосредоточено внимание Проекта GCP, который нацелен на комплексный анализ имеющихся данных о круговороте углерода в окружающей среде и на проведение необходимых измерений. Проект объединяет усилия Международной геосферно-биосферной программы (IGBP), Международной программы по изучению человеческого фактора в глобальном изменении окружающей среды (IHDP) и Всемирной программы изучения климата (WCRP). Круговорот углерода рассматривается как центральный круговорот земной системы, связанный с климатом, круговоротом воды, круговоротами биогенных солей и производством биомассы за счет фотосинтеза на суше и в океанах. Участие человеческого фактора в круговороте углерода определяется факторами производства пищи и формирования условий среды обитания. Однако Проект GCP не учитывает роль малоподвижного углерода и геосферных процессов в регулировании сбалансированного сочетания потоков углерода на Земле, без чего невозможно оце-

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

Библиографический список

1. Баренбаум, А.А. Механизмы самоорганизации при глобальном геохимическом круговороте вещества на Земле / А.А. Баренбаум // Синергетика. - М.: МГУ. - 2000. - № .3. - С. 275-295.

2. Баренбаум, А.А. Галактика, солнечная система, Земля / А.А. Баренбаум. - М.: ГЕОС, 2002. - 392 с.

3. Баренбаум, А.А. Механизм формирования скоплений нефти и газа / А.А. Баренбаум // Докл. АН.

- 2004. - 399. - № 6. - С. 802-805.

4. Баренбаум, А.А. Интенсификация притока глубинных углеводородов / А.А. Баренбаум, С.Н. Закиров, Э.С. Закиров, И.М. Индрупский и др. // Доклады АН. - 2006. - 406. - № 2. - С. 221-224.

5. Барцев, С.И. Глобальные обобщенные модели динамики углекислого газа / С.И. Барцев, А.Г. Дегер-менджи, Д.В. Ерохин // Пробл. окр. среды и природных ресурсов. - 2003. - № 12. - С. 11-28.

6. Войтов, Г.И. Химизм и масштабы современного потока природных газов в различных геоструктурных зонах Земли / Г.И. Войтов // Журнал Всес. хим. о-ва им. Д.И. Менделеева. - 1986. - 31. - № 5.

- С. 533-539.

7. Войтов, Г.И. О холодной дегазации метана в тропосферу Земли // Теоретические и региональные проблемы геодинамики / Г.И. Войтов // Тр. ГИН РАН.- М.: Наука. - 1999. - Вып. 515. -С. 242-251.

8. Дегерменджи, А.Г. Глобальные малоразмерные модели динамики и устойчивости биосферы / А.Г. Дегерменджи, С.И. Барцев // Пробл. окр. среды и природных ресурсов.- 2003. - № 7. - С. 32-49.

9. Кондратьев К.Я. Глобальный круговорот углерода и климат / К.Я. Кондратьев, В.Ф. Крапивин // Исследование Земли из космоса. - 2003. - № 1. -С. 3-15.

10. Кондратьев, К.Я. Моделирование глобального круговорота углерода / К.Я. Кондратьев, В.Ф. Крапивин. - М.: Физматлит. - 2004. - 336 с.

11. Корценштейн, В.Н. Растворенные газы подземной гидросферы Земли / В.Н. Корценштейн. - М.: Недра, 1984. - 232 с.

12. Крапивин, В.Ф. Парниковый эффект и глобальный биогеохимический цикл углекислого газа / В.Ф. Крапивин // Пробл. окр. среды и природных ресурсов. - 2000. - № 2. - С. 2-18.

13. Крапивин, В.Ф. Глобальные изменения окружающей среды / В.Ф. Крапивин, К.Я. Кондратьев. -СПб.: СПб ун-т. - 2002. - 724 с.

ЛЕСНОЙ ВЕСТНИК 9/2012

93