CC BY
148
SEARCH OF THE INDICATOR-PRECURSORS FOR THE DISTURBANCE OF NATURAL PROCESSES IN THE
FOREST, AGROFOREST AND URBOECOSYSTEMS Krapivin V.F., Prof. Institute of Radio Engineering and Electronics. VA Kotelnikov RAS, Dr. Sci. (Tech.); Shalaev V.S., Prof. MSFU, Dr. Sci. (Tech.); Burkov V.D., Prof. MSFU, Dr. Sci. (Tech.); Soldatov V. Yu., IRE RAS, Ph.D (Tech.)
[email protected], [email protected] Moscow State Forest University (MSFU), 1st Institutskaya st., 1, 141005, Mytischi, Moscow reg., Russia Institute of Radioengineering and Electronics (IRE) of RAS, Mokhovaya 11-7, Moscow, 125009, Russia
The problem of the numerical index choice for the state offorest, agroforest or urban ecosystem is analyzed. An algorithm for the biocomplexity indicator calculation as the state index of such ecosystems is proposed. The algorithm is based on the conceptual model of the interaction between the vegetation community and the environment. The model formalizes the energetic exchange process in this interaction by means of introduction of the (V, W)-exchange equation with V as the expendable energy, W as the consumable resources. Such formalization of the two environments' interaction allows the (V, W)-exchange equations, the solution ofwhich provides the stable trend in the ecosystem evolution. The calculation of the biocomplexity indicator is considered as an example for the global system of soil-plantformations with spatial distribution of the 30 types within the 4°*5° geographic scale. In this case, the anthropogenic impact on the forest ecosystems can significantly change the biosphere biological complexity. For example, if the forest area is decreased or increased by 10 percents to 2050, the biosphere biocomplexity will reduce by 12-15 percents or grow by 18-21 percents respectively. The algorithm proposed for the biocomplexity indicator assessment allows reaching the balanced combination between municipal parks and buildings in the urban ecosystem.
Keywords: indicator, monitoring, forest fire, algorithm, stability, information, biocomplexity, biosphere, scenario.
МОДЕЛИРОВАНИЕ ГЛОБАЛЬНЫХ ЦИКЛОВ УГЛЕРОДА И МЕТАНА
В.Ф. КРАПИВИН, проф., ИРЭ им. В.А.Котельникова РАН, д-р физ.-мат. наук, В.С. ШАЛАЕВ, проф., директор ИСИЛМГУЛ, д-р техн. наук, В.Д. БУРКОВ, проф. МГУЛ, д-р техн. наук
[email protected], [email protected] ФГБОУ ВПО «Московский государственный университет леса» 141005, Московская обл., г. Мытищи-5, ул. 1-я Институтская, д. 1, МГУЛ ФГБУН Институт радиотехники и электроники им. В.А. Котельникова РАН
125009, Москва, ул. Моховая 11, кор. 7
Проблема парникового эффекта по причине антропогенного воздействия на глобальный круговорот парниковых газов остается в стадии обсуждения в научной литературе. В данной работе предложена математическая модель совместного биогеохимического круговорота углерода и метана при рассмотрении пространственной структуры распределения их источников и стоков. Все резервуары и потоки углерода и метана, учитываемые в модели, разделены на категории с различными временными масштабами. Пространственный масштаб для источников и стоков углерода и метана на суше выбран 4° по широте и 5° по долготе. Мировой океан параметризуется точечной моделью с выделением четырех уровней по глубине. Взаимодействие между атмосферой и резервуарами углерода на суше и в океане включает процессы фотосинтеза, разложения, дыхания и горения, а также процессы дегазации земной коры. Глобальный круговорот метана в системе атмосфера-гидросфера-суша описывается схемой, отражающей основные его потоки, радиационный потенциал которых является функцией времени. Структура модели совместного круговорота углерода и метана включает 12 блоков, реализующих базовые процедуры расчета их потоков и согласования разнородной информации. Приведены результаты имитационных экспериментов, рассматривающих сценарии изменения площади лесов и демонстрирующих их значительную роль в изменении климата.
Ключевые слова: климат, двуокись углерода, метан, модель, параметризация, поток, лес, парниковый эффект.
Глобальные циклы углерода и метана явля- точниках и стоках. Согласно существующим
ются объектом многочисленных междуна- данным, основными источниками СО2 явля-
родных и национальных научных программ, ются цементная промышленность, сжигание
нацеленных на параметризацию и понимание ископаемых топлив и использование земель-
обратных связей, существующих в системе ных ресурсов. В 2013г. скорость нарастания
биосфера-климат-общество (СБКО). Между- концентрации парниковых газов в атмосфере
народный проект GCP (Global Carbon Project) составила 2,1 % при стоках в океан 28 % и на
является одним из них. Этот проект нацелен суше 23 %. При этом замкнутость, например,
на реализацию различных процедур для акку- круговорота углерода по стокам и источникам
муляции знаний о парниковых газах и их ис- составила 80-85 %. Это означает, что неопре-
деленность сбалансированности глобального цикла углерода составляет 15-20 % . Аналогичная ситуация касается и других парниковых газов.
Многие авторы пытаются создавать модели для описания глобальных циклов углерода и метана [1-10]. Эти попытки стимулируются растущим интересом к проблеме глобального изменения климата и спорами о причинах этих изменений. Дискуссии на эту тему детально проанализированы в работе [6] в связи с противоречиями по проблеме глобальной экодинамики. В связи с этим СгаскпеП et а1. [3] детально классифицировали угрозы устойчивости окружающей среды, включая жизнь человека на Земле. Отмечено, что только одна причина изменения климата из двенадцати имеет антропогенный характер. Другие причины имеют природный характер и только косвенно связаны с деятельностью человека.
Процесс моделирования циклов углерода и метана характеризуется серией моделей, которые имеют иерархию пространственных шкал и элементов окружающей среды. Многообразие таких моделей объясняется сложностью циклов углерода и метана, а также недостатком данных и знаний об их источниках и стоках. Вообще надежность оценки роли парниковых газов в изменении климата зависит от детальности рассмотрения глобальных биогеохимических циклов многих химических веществ, которые так или иначе взаимодействуют с парниковыми газами. Отсюда и следует сложность глобальной модели СКБО, которая бы могла оценить последствия от реализации всего многообразия антропогенных сценариев изменения окружающей среды.
Кгаршп, Varotsos [6] показали, что существует некоторая оптимальная структура глобальной модели СБКО, когда имеющиеся базы данных позволяют реализовать имитационные эксперименты с оценкой глобальной экодинамики. Degermendzhi [4] развил альтернативный подход к моделированию системы биосфера-климат в условиях некомплектности глобальных баз данных. Предложенный метод основан на так называемых маломасштабных
моделях, когда вся биосфера делится на отдельные компартменты. Такой подход экономно использует имеющиеся базы данных, предъявляя к ним минимальные требования.
Однако проблема парникового эффекта за счет СО2 и СН4 остается в состоянии неопределенности. Эта статья развивает комбинированную модель совместного биогеохимического цикла углерода и метана, принимая к рассмотрению пространственное распределение почвенно-растительных формаций (ПРФ) и океанических экосистем. Концептуальная диаграмма глобального цикла углерода включает множество природных и антропогенных процессов, которые определяют его динамику и имеют различные временные масштабы от десятков до сотен и тысяч лет. Это обстоятельство ставит задачу изучения характерных времен полного перемешивания атмосферы. Хорошо известно, что концентрация СО2 в атмосфере значительно изменяется в течение года. Различие между максимальными и минимальными оценками концентрации СО2 варьирует от 10 млн-1 на Южном полюсе до 15 млн-1 на Северном полюсе. Потоки углерода формируются в результате многих природных и антропогенных процессов, указанных в табл. 1. Блок-схема этих процессов указана на рис. 1.
Как следует из табл. 1, все резервуары и потоки углерода делятся на категории, различающиеся по временным масштабам. Модель биогеохимического цикла углерода должна это отражать и описывать его потоки между резервуарами. Очевидно, что существующие базы данных и знаний не могут поставить детальную входную информацию в эту модель о роли каждого дерева, животного, микроба, листа, озера, реки, морской акватории, ландшафта и т.д. Поэтому построение глобальной модели цикла углерода ведется с постоянным усложнением и детализацией [5]. Рис. 1 определяет некоторый достигнутый на этом пути уровень сложности биогеохимического цикла углерода и дает возможность описать его на ограниченном уровне пространственной и предметной детализации, фокусируя внимание на процессах, которые важны на региональном и глобальном уровне. К сожалению, ряд потоков углерода имеют очень приближенные
Рис. 1. Блок-схема глобального цикла углерода в системе атмосфера-океан-суша. Резервуары и
потоки CO2 перечислены в табл. 1 Fig. 1. Block diagram of the global carbon cycle in the atmosphere-ocean-land. Tanks and CO2 fluxes are listed in Table. 1
оценки, включая, например, потоки F3 и F26 В соответствии с этим для модели выбран пространственный масштаб 4°х5°, для которого имеется необходимая информация.
Пространственная экология предлагает различные фрагментации (ПРФ) и океанических акваторий, но в рамках схемы рис. 1 пока невозможно найти необходимых данных. В региональном масштабе, например, географические базы данных Европы и Канады имеют масштаб 1:1000000. Аналогичные базы данных многих стран создаются с ориентацией на экономические аспекты сельского хозяйства и лесной промышленности. База данных ФАО практически покрывает всю сушу в рамках масштаба 1:5000000. К сожалению, эти базы данных из-за их узкой предметной ориентации не позволяют в предлагаемой модели перейти к более детальной пространственной классификации элементов суши и гидросферы. Поэтому здесь, согласно [2, 5, 6], для наземных экосистем принимается классификация на 30 типов ПРФ, а мировой океан описывается моделью Тарко [8].
Как видно из рис.1, роль Мирового океана в глобальном цикле углерода параметризуется через процессы обмена на границе атмосфера-океан и перемешивания между различными слоями воды с участием соответствующих экосистем [6].
Потоки углерода, представленные на рис. 1 и в табл. 1, используются для балансовых уравнений
дt "" 0Ф ь дХ ' й } ,
где С^ - резервуар углерода в я-м пикселе пространственной структуры на суше и в океане, / и Js - источники и стоки углерода соответственно, ф - широта; X - долгота; ? - время;
V^У^) - поле ветра в я-м пикселе.
Взаимодействие между атмосферой и резервуарами углерода на суше и в океане выражается в потоках углерода, формируемых экологическими, геофизическими и биогеохимическими процессами, включая фотосинтез, дыхание, разложение, горение,
Таблица 1
Резервуары и потоки углерода в форме СО2 в биосфере, учитываемые в схеме рис. 1 Reservoirs and fluxes of carbon in the form of CO2 in the biosphere recognized in the scheme of Fig. 1
Средние оценки
Резервуары и потоки CO2 Идентификатор резервуаров (ГтC) и потоков (rrC/год)
Атмосфера CA 650-750
Фотический слой выше термоклина Cu 380-520
Океан Промежуточный фотический слой (под термоклинном) Океанические растения (водоросли) Живые элементы Ci Сф CL 280-610 3 1,5
Детрит Глубокие слои океана C dt cd 29 28500-33890
Придонный слой океана Cb 5000
Мертвое органическое вещество ( в почве и торфе) cs 1500-3000
Суша Растительность Cv 550-610
Вечная мерзлота cp 1672
Живые организмы clo 0,5
Ископаемые топлива Уголь Нефть cc co 4500 500
Газ Cg 5000
Эмиссия за Растительность f8 6,9
счет сжигания Ископаемые топлива и производство цемента F28 4000
Выветривание горных пород F4 0,04
Извержение вулканов F5 2,7
Поглощение наземной растительностью F6 224,4
Растения F7 50-59,3
Дыхание Население f10 0,7
Животные F11 4,1
Разложение гумуса почв f9 139,5
Эмиссия Корни растений f15 56,1
Таяние вечной мерзлоты f2 190
Жизненные Население F12 0,3
функции Животные F13 3,1
Разложение растительных остатков F14 31,5-50
Осаждение в донные отложения f.6 0,1-0,2
Растворение морских осадков f,7 0,1
Разложение Фотический слой f22 35
детрита Глубокие слои океана F.s 5
Подъем глубоких вод f.9 45
Сток с поверхностными водами и гравитационное осаждение f20 40
Фотосинтез F„ 69
Подземный сток f23 0,5
Поверхностный сток f24 0,5-0,6
Дыхание живых организмов в океане f25 25
Процессы дегазации F. 21,16
Сток в недра Земли f26 1,3
Десорбция f27 97,08
Сорбция F3 100,0
извержение вулканов, выветривание горных пород и т.д. Детальное описание этих процессов и их параметризация выполнены в работах [5, 6, 9]. Некоторые из них требуют уточнения.
Пусть Рк(ф, к, 0 есть продукция фотосинтеза для растительности тика к в пикселе Е.. с географическими координатами (ф., к) в
У 1 .
момент времени Функция Рк(ф, к, 0 зависит от осадков и температуры атмосферы. Тогда
Таблица 2
Список блоков модель глобальных циклов углерода и метана (МГЦУМ), схемы которой
представлены на рис. 1 и 2 List of model units of the global cycles of carbon and methane, circuits which are shown in Fig. 1 and 2
Блок Описание блока
СПР Согласование пространственного распределения пикселей Е. и акваторий ОЛ с доступными базами данных и знаний
УМА Управление моделями и алгоритмами для параметризации потоков Е (г = 1,., 28) и С. ( = 1,., 13)
РП Расчет потоков Е (г = 1,., 28) и С. ( = 1,., 13) с учетом пространственной структуры системы атмосфера-суша-океан
КП Корректировка потоков Е (г = 1,., 28) и С. ( = 1,., 13) за счет рассмотрения циркуляции атмосферы, гидрологических процессов и течений океана
АОБ Алгоритм обновления базы данных и знаний
АСРИ Алгоритм согласования разнородной информации
ФС Формирование сценариев изменения пространственной структуры почвенно-растительных формаций, потоков Е3, Е27 и Е28 вследствие загрязнения поверхности океана и изменения технологий использования минеральных ресурсов
ОИИК Оценка индикаторов изменения климата в результате реализации сценариев.
УИП Учет изменений потоков Е (г = 1,., 28) и С. ( = 1,., 13) в результате вариаций климатических параметров
СДОС Слежение за динамикой обратных связей в окружающей среде при изменении их параметров
ПИРВ Пересчет индексов радиационного воздействия в зависимости от динамики Е/С. в соответствующих пикселях и акваториях
ИИК Информационный интерфейс и контроль потоков данных между блоками МГЦУМ
поток углерода между атмосферой и живой биомассой на наземном пикселе Е.. может
У
быть описан простой моделью
^(фЛО = С П(Ф, \ ъ
где С6 - коэффициент, отражающий эффективность механизма фотосинтетического отклика.
Цикл углерода в океане включает потоки, которые формируются термохалинной циркуляцией, обменом между атмосферой и поверхностным слоем океана, фотосинтезом, окислительными процессами, осаждением и поглощением СО2 из атмосферы. Описание моделей этих процессов дано в работах [5, 6, 9]. Описывающие физико-химические процессы газового обмена на границе воздух-вода потоки F27 и F3 были детально изучены, описаны в [1] и, как результат, поток F3 может быть представлен формулой
^=11/(71)^/(1+0,5М
где у(Ть) - показатель влияния температуры на растворимость СО2 в воде; р5 - соленость воды; Т - температура воды; ра - атмосферное давление. Потоки F27 и F3 являются функционально симметричными в зависимости от
кислотности воды: при рН<8,11 океан поглощает С02, а при рН>8,11 океан выделяет СО2.
Совместная модель глобальных циклов углерода и метана (МГЦУМ) реализует схемы рис. 1 и 2, и ее блоки представлены в табл. 2.
Глобальные циклы углерода и метана, главным образом, определяют динамику климата. Источниками метана служат заболоченные территории (22 %), добыча угля, нефти и газа (19 %), процессы брожения (16 %), выращивание риса (12 %), сжигание биомассы (8 %), мусорные свалки (6 %), обработка сточных вод (5 %), отходы животноводства (5 %), термиты (4 %) и океан (3 %). Средняя скорость формирования СН4 в каждом пикселе е.. на суше и в акватории Ок океана изменяется от 0,4х106 т/г в лесистой местности до 280х106 т/г на рисовых плантациях и от 0,1 х106 т/г на океаническом шельфе до 6,7х106 т/г в открытом океане [5, 6].
Знание этих потоков метана при их пересчете на радиационный потенциал СО2 обеспечивает реализацию схемы рис. 1 с уточненными характеристиками потоков Е (г = 1,..., 28) в пикселях Е.. и акваториях О
sK>
Рис. 2. Глобальный цикл метана в системе атмосфера-гидросфера-суша Fig. 2. The global cycle of methane in the atmosphere-hydrosphere-land
где источники и стоки метана существуют. Обычно принимается, что отношение радиационных потенциалов С02:МН4 составляет 1:72, 1:25 и 1:7,6 в течение 20, 100 и 500 лет соответственно.
Площади многих источников и стоков СН4 могут занимать лишь часть пикселя 4°х5°. Каждый пиксель е.. имеет природные и антропогенные участки, которые влияют по-разному на потоки метана. Поэтому МГЦУМ позволяет использовать такую информацию, и тем самым пространственная структура модели уточняется.
Потоки метана С. (г = 1-13) являются функциями температуры, географических координат и других характеристик окружающей среды. Например, Panikov и Dedysh [7] разработали модель эмиссии метана в атмосферу для снегом покрытых болот в Западной Сибири, что позволяет рассчитывать потоки СН4 в пикселях с вечной мерзлотой. В этом случае поток С можно описать соотношением
dC/dt = цХ(1-У)/У, где ц - характерная скорость роста микробной популяции с биомассой X;
У - прирост биомассы на единицу потребленного субстрата.
При этом микробная биомасса X формируется по закону
dX/dt=yX,
где ц = ^ SR/(S + KS) S - концентрация ката-болического субстрата;
R - функция физиологического состояния;
KS - константа насыщения, численно равная той концентрации субстрата, при которой ц = 0.5ц .
г г г max
Природные ветланды и рисовые плантации поставляют в атмосферу более 30 % метана из всех источников. Поток C3 опишем уравнением
C3 = Я/l(ГWз(PHf4(rP, где Hr - гетеротрофное дыхание,
T - температура почвы,
h - положение уровня водной поверхности,
r - окислительно-восстановительный по-
p
тенциал,_/\ (/=1-4) - функции, параметризующие скорости эмиссии метана. Потоки C1 и C2 характеризуют основные стоки атмосферного метана и описываются его реакцией с OH [10].
Описанная здесь глобальная модель цикла двух парниковых газов позволяет изучать роль почвенно-растительных формаций
Таблица 3
Модельные оценки изменений запасов углерода при условиях, когда все хвойные
леса в Северном полушарии будут сожжены. Обозначения: CX=C+CD+CB Model estimates of carbon stock changes under conditions where all the coniferous forests in the northern hemisphere will be burned. Designations: CX=C+CD+CB
Время после воздействия, годы Отклонение в запасах углерода, Гт Время после воздействия, годы Отклонение в запасах углерода, Гт
АС, AQ АСх АС, АС, АС, АСх
0 150,2 -5,7 15,7 0,1 60 25,1 -30,4 8,5 13,1
10 114,6 -34,1 30,8 3,1 70 16,7 -25,7 6,7 13,5
20 84,1 -44,2 22,5 7,4 80 11,2 -21,5 5,1 14,0
30 64,4 -43,9 19,4 8,5 90 5,9 -17,9 4,3 14,3
40 48,2 -40,7 14,6 10,4 100 2,1 -14,8 2,8 13,9
50 35,2 -33,8 12,1 11,7 200 -9,1 -3,4 -0,9 13,5
Таблица 4
Модельные оценки изменений запасов углерода при условиях, когда все леса в Северном полушарии будут сожжены. Обозначения: CX=C+CD+CB Model estimates of carbon stock changes under conditions where all the forests in the northern hemisphere will
be burned. Legend: CX=C+CD+CB
Время после Отклонение в запасах углерода, Гт Время после Отклонение в запасах углерода, Гт
воздействия, АС АС, АС АС воздействия, АС АС, АС АС
годы годы
0 238,1 -7,9 24,9 0,1 60 44,2 -49,1 14,3 20,7
10 174,2 -31,6 47,9 4,9 70 33,1 -41,6 10,3 21,9
20 138,9 -67,6 39,2 10,0 80 20,7 -35,0 7,5 22,8
30 107,9 -90,3 32,2 13,8 90 12,9 -30,1 5,5 23,3
40 82,0 -64,3 24,1 16,8 100 7,3 -24,7 3,7 23,6
50 60,9 -56,9 18,4 19,1 200 -12,7 -5,9 -1,7 21,7
Таблица 5
Модельные оценки изменений запасов углерода при условиях, когда все тропические леса будут сожжены. Обозначения: CX=C+CD+CB Model estimates of carbon stock changes under conditions where all tropical forests are burned.
Legend: CX=C+CB+CB
Время после воздействия, годы Отклонение в запасах углерода, Гт Время после воздействия, годы Отклонение в запасах углерода, Гт
АС, АС, АС, АСх АС, АС, АС, АСх
0 406,2 -20,0 42,2 0,2 60 2,9 -12,4 3,0 22,8
10 264,4 -93,7 74,1 8,0 70 -5,8 -7,5 0,5 22,8
20 162,2 -84,8 48,0 14,9 80 -11,6 -4,2 -0,9 22,6
30 90,6 -38,6 27,9 19,1 90 -13,2 -2,6 -1,7 22,5
40 45,4 -36,5 15,0 21,3 100 -14,5 -1,9 -2,1 21,8
50 18,3 -21,6 7,5 22,4 200 -13,2 -2,3 -1,9 17,7
и океанских экосистем в изменениях климата путем рассмотрения сценариев антропогенного воздействия на них. Примеры гипотетических катастрофических изменений лесных территорий приведены в табл. 3-5 и на рис. 3. Видно, что наземная биота играет роль стока атмосферного СО2 в динамике основных ре-
зервуаров углерода. Аналогичные имитационные эксперименты дают возможность оценить роль всех типов ПРФ в поглощении СО2. Например, на рис. 3 показана роль сокращения лесных покровов в динамике концентрации СО2 в атмосфере. Как следует из результатов, приведенных на рис. 3, при сокращении гло-
Рис. 3. Динамика концентрации CO2 для различных сценариев изменения площадей лесов: 1 - скорость изменения площади лесов остается на уровне 2010 г.; 2 - все леса к 2050 г. будут сведены; 3 - площадь лесов к 2050 г. будет сокращена на 50 %; 4 - то же на 10 %; 5 - к 2050 г площадь лесов увеличится на 50 %; 6 - то же на 10 % Fig. 3. Dynamics of CO2 concentration for different scenarios of change in forest area: 1 - the rate of change in forest area remains at 2010 .; 2 - all forests by 2050 will be reduced; 3 - forest area in 2050 will be reduced by 50%; 4 - the same at 10%; 5 - 2050 forest area will increase by 50%; 6 - the same 10%
Таблица 6
Динамика отношения интегральных скоростей (C6) поглощения СО2 растительным покровом в процессе его замены Dynamics ratio of the integrated absorption (C6) rate of CO2 vegetation in the process of replacing it
Сценарий замены почвенно-растительных формаций в пикселях 4°х5° F F6 измененная ПРФ/ б первоначальная ПРФ
Первоначальная ПРФ Измененная ПРФ 2020 г. 2030г. 2050 г. 2100 г.
Арктические пустыни и тундры Лесотундра 2,54 2,04 1,97 2,15
Тундры Лесотундра 0,96 0,94 0,95 1,03
Горная тундра Лесотундра 1,42 1,15 1,01 1,04
Северотаежные леса Среднетаежные леса 1,63 1,44 1,11 1,12
Южнотаежные леса Среднетаежные леса 1,98 1,66 1,45 1,33
Широколиственные хвойные леса Среднетаежные леса 3,91 3,69 2,86 1,95
Субтропические широколиственные и Влажные вечнозеленые
хвойные леса тропические леса 3,07 2,55 2,44 1,77
Ксерофитные редколесья Влажные вечнозеленые
и кустарники тропические леса 21,49 19,46 17,95 18,34
Степи умеренно засушливые Влажные вечнозеленые
и засушливые тропические леса 22,15 18,33 16,32 14,77
Пампы и травяные саванны Влажные вечнозеленые
тропические леса 99,23 77,12 68,54 70,09
Субтропические пустыни Влажные вечнозеленые
тропические леса 187,65 153,23 138,39 140,67
Альпийские и субальпийские луга Влажные вечнозеленые
тропические леса 790,04 766,42 751,26 767,33
Переменно-влажные листопадные Влажные вечнозеленые
тропические леса тропические леса 1,42 1,33 1,23 1,25
Тропические ксерофитные открытые Влажные вечнозеленые
редколесья тропические леса 67,56 60,19 56,59 57,34
Тропические саванны Влажные вечнозеленые
тропические леса 5,98 4,98 4,67 5,08
Тропические пустыни Влажные вечнозеленые
тропические леса 25,99 24,76 23,57 22,65
Субтропические и тропические травя- Влажные вечнозеленые
но-древесные заросли тугайкого типа тропические леса 17,11 15,78 14,92 13,99
Субтропические полупустыни Влажные вечнозеленые
тропические леса 0,95 1,21 0,97 1,09
бальной площади лесов на 19 % к 2050 г. по сравнению с 2010 г. (4,033 млн га) концентрация СО2 к концу 21-го столетия возрастет на 53 %. Наоборот, увеличение площади лесов на 10 % в тот же период уменьшит концентрацию СО2 в атмосфере на 12 %. Другие сценарии, результаты реализации которых представлены на рис. 3, показывают катастрофические последствия в изменении климата. Следовательно, человечество находится в ситуации, когда сокращение площади лесов стимулирует негативные изменения окружающей среды. Более того, процесс изъятия наземного углерода, как это следует из табл. 3-5 и рис. 3, существенно зависит от стратегии управления лесным хозяйством. Поэтому даже рассмотренный здесь вариант пространственного распределения почвенно-растительных формаций и выбор пиксельной структуры наземных экосистем МГЦУМ помогает получить оценки последствий реализации стратегий управления растительным покровом как в глобальном и региональном масштабах, так и в каждом пикселе 4°х5°. Например, в табл. 6 показано, как поток С6 изменяется в зависимости от смены растительного покрова.
В заключение отметим, что структура МГЦУМ не изменится при рассмотрении другой схемы пространственного распределения почвенно-растительных формаций. К сожалению, географическая сетка 4°х5° достаточно груба для выявления более тонких эффектов управления растительными покровами. Имеющиеся базы данных для пространственных шкал 1°х1°; 0,5°х0,5° и 1кмх1км позволяют уточнить некоторые блоки модели, но для этого необходимо дополнить их список с описанием биопродукционных процессов как функций вегетационных индексов типа NDVI. Однако многие пиксели в этом случае оказываются в неопределенности. Так, например, в зонах вечной мерзлоты важными являются геофизические, геохимические и экологические процессы, которые связаны с изменением климата и еще недостаточно изучены. Известно, что процессы таяния вечной мерзлоты значительно изменяют земную поверхность, и поэтому совершенствование МГЦУМ неизбежно требует развития ГИС для этой зоны.
Конечно, число нерешенных задач глобального моделирования парникового эффекта значительно. Достаточно схематичны сведения о связи земного климата с космической пылью. Известно, что 100-300 т космической пыли поступает в земную атмосферу ежедневно. С учетом этого предстоит расширить блоковую структуру МГЦУМ.
Данная работа была частично поддержана Российским фондом фундаментальных исследований в рамках гранта № 13-0100023.
Библиографический список / References
1. Алексеев, В.В. Физическое и математическое моделирование экосистем / В.В. Алексеев, И.И. Крышев, Т.Г. Сазы-кина. - Санкт-Петербург: Гидрометеоиздат, 1992.- 367 с. Alexeev V.V., Kryshev I.I., Sazykina T.G. Fizicheskoe i matematicheskoe modelirovanie ekosistem [Physical and mathematical modeling of the ecosystems]. St-Petersburg: Hydrometeoizdat, 1992. 367 p.
2. Базилевич, Н.И. Картосхемы продуктивности и биологического круговорота главнейших типов растительности суши / Н.И. Базилевич, Л.Е. Родин // Известия Всесоюзного географического об-ва, 1967. -T. 99. - № 3. - С. 190-194.
Bazilevich N.I., Rodin L.E. Kartoskhemy produktivnosti i biologicheskogo krugovorota glavneishikh tipov rastitel'nosti sushi [Schematic map of productivity and biological cycle of the main types of the land vegetation]. Proceedings of Russian Geographical Society, 1967. V. 99. № 3. pp. 190-194.
3. Cracknell A.P., Krapivin V.F., Varotsos C.A. Global Climatology and Ecodynamics: Anthropogenic Changes to Planet Earth. Chichester, U.K.: Springer/Praxis, 2009. 518 p.
4. Degermendzhi A.G. New directions in biophysical ecology. In: A.P. Cracknell, V.F. Krapivin, and C.A. Varotsos (eds.). Global Climatology and Ecodynamics. Chichester, U.K.: Springer/Praxis, 2009. pp. 379-396.
5. Kondratyev K.Ya., Krapivin V.F., Varotsos C.A. Global Carbon Cycle and Climate Change. Chichester, U.K.: Springer/PRAXIS, 2003. 372 p.
6. Krapivin, V.F., Varotsos C.A. Biogeochemical cycles in globalization and sustainable development. Chichester, U.K.: Springer/Praxis, 2008. 562 p.
7. Panikov N.S., Dedysh S.N. Cold season CH4 and CO2 emission from boreal peat bogs (West Siberia): winter fluxes and thow activation dynamics. Global Biogeochemical Cycles, 2000. V. 14. № 4. pp. 1071-1080.
8. Tarko A.M. Analysis of Global and Regional Changes in Biogeochemical Carbon Cycle: A Spatially Distributed Model. Laxenburg, Austria: Interim Report, IR-03-041, IIASA, 2003. 28 p.
9. Williams, R.G. Ocean Dynamics and the Carbon Cycle. / R.G. Williams and M.J. Follows. Massachusetts: Massachusetts Institute of Technology, 2011. 434 p.
10. Xu, S. A process-based model for methane emission from flooded rice paddy systems / S. Xu, P. Jaffe, D.L. Mauzerall. Ecological Modelling, 2007. V.205. pp. 475-491.
