CC BY
56
Вычислительные технологии
Том 13, Специальный выпуск 3, 2008
Моделирование углеродного баланса болотных экосистем южной тайги при различных сценариях изменения климата
Е, А, Дюкарев, Е, А, Головацкая Институт мониторинга климатических и экологических систем, СО РАН,
Томск, Россия E-mail: egor@imces.ru, golovatskaya@imces.ru
CO2 and CH4 emission fluxes and net primary biological production was estimated for different peatland ecosystems. A local model of the carbon cycle in peatland ecosystems accounting for accumulation of carbon in vegetation and peat deposit was developed. A forecast of development of peatland ecosystems in the south taiga is provided for various scenarios of climate change.
Введение
Наблюдающееся в XX в. глобальное потепление проявляется во всех регионах России. От конца XIX к концу XX в. потепление для территории России в целом составило около +1°С [1]. Существует мнение о возможной связи через парниковый эффект заметного повышения концентрации углекислого газа в атмосфере с глобальными изменениями климата. Основная планетарная функция наземных экосистем заключается в формировании и регулировании потоков углерода в биосфере, обеспечении непрерывной работы углеродного цикла. Функционирование наземных экосистем в качестве источников или стоков углекислоты определяется балансом между двумя мощными биосферными процессами — фотосинтетической продукцией органического углерода и выделением С02 при дыхании и разложении органического вещества. Болота — единственные в наземной биоте экосистемы, которые обеспечивают постоянный долговременный сток в них углерода.
В настоящей работе представлены результаты многолетних (1999-2006 гг.) натурных измерений потоков и запасов углерода олиготрофных болотных экосистем Западной Сибири, на основании которых построена локальная модель углеродного цикла болотной экосистемы и выполнен прогноз развития торфоболотных экосистем южной тайги при различных сценариях изменения климата.
1. Объекты и методы исследования
В качестве объекта исследования выбрано торфяное месторождение Бакчарское, расположенное между реками Икса и Бакчар (Томская область), которое характеризуется
© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2008.
закономерной сменой болотных фитоценозов олиготрофпого тина: сосново-куетарнич-ково-сфагновый фитоценоз — высокий рям; такой же фитоценоз с угнетенным древостоем — низкий рям; открытая осоково-сфагновая топь. Мощность торфяной залежи изменяется от 1 до 3 м при возрасте от 3 до 5 тысяч лет [2].
Измерения потоков эмиссии С02 с поверхности торфяных залежей проводились камерным методом с использованием инфракрасного газоанализатора ОПТОГАЗ 500.4 два раза в месяц па высоком ряме и открытой топи, а па низком ряме — в круглосуточном режиме каждые 3 ч в течение семи дней ежемесячно (с мая по сентябрь).
Для последующей оценки углеродного баланса одновременно проводились наблюдения
2
(па высоте 0.3...'2 м). Сопровождающие метеорологические наблюдения выполнялись па низком ряме в соответствии с требованиями стандартных метеопаблюдепий,
2. Экспериментальные исследования круговорота углерода в болотных экосистемах
Чистая первичная продукция (NPP) является показателем накопления углерода в виде
растительного вещества. Результаты исследования показали, что в среднем олиготроф-
пые болотные экосистемы имеют одинаковые величины па высоком и низком ряме и
22
учета древесного яруса). Основной вклад в продукцию па исследуемых олиготрофпых болотных биогеоценозах вносят корпи растений (53 59%) и сфагновые мхи (22 30%). Динамика продуктивности за исследуемые годы определяется погодными условиями (рис. 1,а).
NPP
фитоценозов, показал, что продукция па высоком ряме и открытой топи зависит от температуры воздуха в течение вегетационного периода, па низком ряме температура играет незначительную роль, в то время как количество осадков имеет сильное влияние па продуктивность.
2
2
J Корни
□ Мок
□ Куттарютлот
□ 'Грата
Высокьш р ям Низкьш р ям Олфьггал того
] "Эмиссия С02 П Эмиссия СН4 □ Сток с юдамя
Рис. 1. Чистая первичная продукция (а) и эмиссия углерода (б) в различных фитоценозах за 1999 2006 гг.
сфагновой топи) (рис, 1,6). Установлена зависимость эмиссии С02 от температуры (г = 0.59) и влажности воздуха (г = — 0.40) воздуха, а также от содержания С02 в приземном слое воздуха. Общий поток углерода из торфоболотных экосистем состоит из эмиссии С02 и СН4 с поверхности торфяной залежи и выноса углерода с болотными
водами. Эмиссия СН4 составляет от 5 % (в высоком ряме) до 13 % (на топи) от эмиссии 2
эмиссии (рис. 1,6).
Сравнение NPP и эмиссии углерода с поверхности торфяной залежи показало, что
2
с поверхности торфяной залежи всех исследуемых биогеоценозов. Таким образом, в настоящее время на территории исследуемого болотного массива преобладает процесс депонирования углерода в виде органического вещества растений и торфа.
3. Прогноз развития торфоболотных экосистем при изменении климата
Круговорот углерода в болотной экосистеме можно представить схематически в виде классической системы блоков и потоков. Блоками являются содержание углерода в атмосфере (А), запасы углерода в живой фитомассе (Live) и мортмассе (Mort). В качестве потоков рассматриваются: нетто-поступление углерода из атмосферы в живую фитомассу — разница между усвоением углерода в процессе фотосинтеза (Pli ) и дыханием растений (Rl); переход углерода живой фитомассы в мортмассу (D) — отмирание,
2 4
ности (Rd). Имеются сведения о значительной роли процессов реассимиляции (геА) в углеродном балансе болот [3]. Однако в связи с трудностью количественных оценок в предложенной модели она не учитывается. Некоторый вклад в формирование баланса вносит вынос углерода с болотными водами (X). В рамках локальной модели считаем, что атмосферная концентрация углерода меняется независимо от потоков углерода между пулами органического вещества. Для изучения долговременных изменений в системе будем рассматривать годовые величины потоков без учета сезонной динамики.
Рассмотрим отношение содержания углерода в фито- и мортмассе к запасам фитомассы в начальный (современный) момент (L0): M = Mort/L0, L = Live/L0, Изменение содержания углерода в живой массе растений и мортмассе будет описываться уравнениями, которые в безразмерном виде выглядят следующим образом:
Здесь ИРР — годичная первичная продукция, которая связана с общими запасами фитомассы линейной зависимостью ИРР = ИРР0Ц 7 — коэффициент, характеризующий отмирание фитомассы; Яр — минерализация мортмассы, X — сток с болотными водами.
Годичная продукция растительности с температурой и количеством С02 в воздухе связана соотношением
где С — концентрация углекислого газа в воздухе; Т — температура воздуха; р0 — коэффициент, связанный с современной продуктивностью растительности; р и р2 —
dL/dt = (1 — Y)pL, p = NPP0/L0, dM/dt = YP L — f — x, f = Rd/Lo,x = X/Lq.
(1)
p = po(1 + pi AT )(1 + Р2ДС ),
(2)
коэффициенты, определяющие увеличение продуктивности при росте температуры и концентрации С02 [4].
Скорость эмиссии углерода из залежи зададим соотношением
/ (Т) = /с (1 + /1ДТ + / ДС), (3)
где /1 и ¡2 — коэффициенты, связанные с изменением эмиссии при увеличении темпе-
2
Для анализа влияния изменений климата па процесс аккумуляции углерода в торфе выполнены расчеты по модели (1) с учетом соотношений (2) и (3), определяющих
зависимости основных потоков углерода от климатических характеристик, Копцептра-2
[6], Прогноз изменения климата в XX в, на территории Западной Сибири взят из данных расчетов но модели НасЮМЗ, выполненных в Центре Хедли (Великобритания) [7], Решение дифференциального уравнения (1) проводилось численно с использованием явного метода интегрирования. Шаг но времени в модели составлял один год. Исходные данные для моделирования приведены в табл. 1.
Т а б л и ц а 1. Параметры мидели для различных фитицеышив
Параметр Высокий рям Рям Тонь
Фитомаееа (Ь0), гС/м2 591.8 539.5 436.9
Мортмаееа (Мо), гС/м2 2036 1940 1290
ИРРо, гС/м2/год 320.4 285.8 259.8
Э миссия, г С / м2 / год 185.4 135.2 111.2
Сток, г С /м2 /год* 15.0 15.0 15.0
Безразмерный параметр
Мо/Ьо 3.44 3.60 2.95
Ро 0.54 0.53 0.59
/о 0.31 0.25 0.25
X 0.03 0.03 0.03
* — оценка по литературным данным
Рис;. 2. Динамика фито- и мортмаееы для экосистем южной тайги при неизменном климате
Рис. 3. Изменение запасов фито- и мортмаееы при различных сценариях изменения климата для ряма
Т а б л и ц а 2. Влияние изменений климата на темпы накопления углерода торфяной залежью. Приведено изменение темпов аккумуляции относительно современных, %
Высокий рям Рям Тонь
Фитомасса
А1 29.0 28.3 32.3
А2 23.4 22.9 26.0
В1 16.2 15.8 17.9
В2 15.7 15.3 17.4
Мортмасса
А1 114.0 100.7 100.9
А2 91.6 79.1 81.0
В1 60.8 52.9 54.0
В2 60.5 52.4 53.4
Варьирование начальных значений эмиссии и продуктивности показывает, что в зависимости от соотношения между параметрами /0 и р0 может происходить как рост, так и деградация торфяной залежи. На основании анализа литературных данных и собственных исследований нами определены следующие параметры, характеризующие чувствительность моделируемых процессов к изменениям климата: /1 = 0.035, /2 = -0.002, р1 = 0.01 Р2 = 0.002, 70 = 0.99. Модельный расчет, в котором климат считался неизменным в течение всего XX в., показал, что происходит рост запасов мортмаееы для всех исследованных экосистем (рис. 2).
При анализе климата в соответствии с различными сценариями установлено нелинейное изменение запасов фитомассы и мортмаееы (рис. 3). Увеличение концентрации С02 в атмосфере и сопутствующее ему потепление климата вызывают рост скорости накопления углерода в торфяной залежи. Наиболее сильно этот эффект проявляется к концу века.
В табл. 2 приведены данные но увеличению темпов накопления углерода в различных экосистемах относительно современных темпов аккумуляции, которое произойдет к концу столетия. К 2100 в южной тайге скорость роста пула мортмаееы возрастет на 52 114% но сравнению с современными темпами. Наибольшие изменения в темпах аккумуляции углерода наблюдаются в экосистеме высокого ряма.
Заключение
Согласно модельным расчетам, в зависимости от типа экосистемы и выбранного сценария изменения климата скорость депонирования углерода растительностью к 2100 г, возрастет на 22±6 %, а скорость накопления в торфяной залежи — на 75±22 % по сравнению с современными темпами накопления.
Список литературы
[1] Груза Г.В., Ранькова Э.Я. Колебания и изменения климата на территории России // Физика атмосферы и океана. 2003. Т. 39, № 2. С. 166-185.
[2] Большое Васюганское болото. Современное состояние и процессы развития / Под ред. М.В. Кабанова. Томск: Изд-во ИОА СО РАН, 2002. 229 с.
[3] Naumov А.V. Carbon budget and emission of greenhouse gases in bog ecosystems of Western Siberia // Eurasian Soil Sci. 2004. Vol. 37. P. 58-64.
[4] Tapko A.M. Антропогенные измерения глобальных биосферных процессов. М.: Физмат-лит, 2005. 232 с.
[5] Головацкая Е.А., Дюкарев Е.А., Ипполитов И.И., Кабанов М.В. Влияние ландшафтных и гидрометеорологических условий па эмиссию СО2 в торфоболотных экосистемах // Доклады Академии наук. 2008. Т. 418, № 4. С. 539-542.
[6] IPCC Climate Change 2001: the Scientific Basis. Cambridge: Cambridge Univ. Press, 2001. 881 p.
[7] Pope V.D., Gallani M.L., Rowntree P.R., Stratton R.A. The impact of new physical parametrizations in the Hadley Centre climate model — HadAM3 // Climate Dynamics. 2000. Vol. 16. P. 123-146.
(http://www.met-oifice.gov.uk/research/hadleycentre/models/HadCM3.html).
Поступила в редакцию 21 февраля 2008 г.
