close

Вход

Забыли?

вход по аккаунту

?

Бархатов Ю.В., Дегерменджи А.Г

код для вставкиСкачать
Новый подход к
моделированию динамики
продуктивности и газообмена
бореального леса
Бархатов Ю.В., Дегерменджи А.Г.
Институт биофизики СО РАН, Красноярск,
Академгородок, 660036
Регистрируемые климатические изменения
планетарного масштаба, как естественные, так и
антропогенные, обусловливают новые режимы
функционирования природных экосистем.
Изменяются не только физические и
биофизические характеристики взаимодействия
составляющих литосферы, атмосферы,
гидросферы и биосферы, но и скорости и
направленности круговоротов основных
биогенных элементов, возникают новые,
усиливаются или ослабляются прямые и обратные
связи, ускоряются или замедляются первичные
процессы сукцессии экосистем.
доля статей по
ключевым словам
"global change" в
общем числе статей
3,E-04
3,E-04
2,E-04
доля статей по
ключевым словам
"Model(ling)" и
"Biospere" в общем
числе статей
2,E-04
1,E-04
5,E-05
0,E+00
1970
1980
1990
2000
2010
Динамика количества статей по глобальной
биосферной тематике за 1975-2003 годы в
“импактовых” журналах (по данным ISI).
Природные экосистемы Сибири играют
существенную роль в стабилизации атмосферы и
гидросферы в условиях современного климата,
особенно в связи с аккумулированием тепличных
газов из атмосферы. Лесные экосистемы Средней
Сибири играют важную роль в поддержании
структурного баланса углерода и баланса
минеральных веществ, обладают высокой
способностью к саморегуляции и установлению
равновесия процессов обмена вещества с
окружающей природной средой. В настоящее
время неизвестно, как будут функционировать и
изменяться существующие природные экосистемы
на территории Сибири в условиях меняющегося
глобального и особенно регионального климата.
1,E-04
1,E-04
"сток углерода"
"моделирование
бореальных лесов"
8,E-05
6,E-05
4,E-05
2,E-05
0,E+00
1990
1995
2000
Динамика роста доли публикаций по
тематикам «сток углерода» и
«моделирование бореальных лесов»
(по базе данных ISI)
2005
Задачи моделирования бореальных лесов Сибири:
- определение углеродного баланса в экосистеме
Северной Евразии. Возможность расширения модели
на экосистемы соответствующей группы (Евразия,
Северная Америка).
- изучение направленности и интенсивности
круговоротов лесных экосистем в меняющемся
климате, количественная оценка круговоротов с
определением, какие из них и в каких условиях
становятся
источниками
или
остаются
аккумуляторами тепличных газов
- оцененка буферной роли природных экосистем по
основным биогенным элементам для обширных
территорий Сибири и способность поглощения
углерода сибирскими лесами и болотами.
АТМОСФЕРНЫЙ БЛОК
Экосистема
бореального
леса
Данные
спутникового
зондирования
Данные
наземных
наблюдений
ГИДРОБЛОК
Блоки модели
Блоки контроля
O2 CO2
Е
O2 CO2
животные
бактерии
дерево
O2 CO2
подстилка
останки
животных
травяной покров
O2 CO2
биогенные
элементы
Е
почва
бактерии
Блок-схема математической модели ценоза бореального
леса Каждый из блоков модели в свою очередь может
представлять собой сложную структуру – субмодель.
атмосфера
лес
почва
животные
ручьи, реки
седименты
Схема потоков углерода в экосистеме
бореального леса
Используются дифференциальные уравнения вида
dX
n
dt
где -
F ( X 1 ,..., X n ,..., X m )
X n одна из переменных экосистемы
Зелен ая ф и том асса
dX
1
dt
X
1
Х 1 ( Х1 )
Х1
A A
1
Y X1
1
1
1
2
R
(
at
bt
c
)
1 K T ( T T opt ) 1 K W ( W W opt ) 1 K E ( E E opt )
0
t – врем я, п рош ед ш ее с н ачала вегетати вн ого сезон а, Т – тем п ература, W – влаж н ость, E освещ ен н ость
Н еф отоси н тези рую щ ая ф и том асса
dX
2
dt
Х 1 Х 2
X 2
ˆ X N
2
N KN
Д ревеси н а
0
dX 3 dt
X 2 t t end
Т рава
dG
dt
G
G ( G G ) A A
2
YG
1
1
1
2
R
(
a
t
b
t
c
)
G
1
1
1
1 K T 1 ( T T opt 1 ) 1 K W 1 ( W W opt 1 ) 1 K E 1 ( E E opt 1 )
0
П од сти л к а
B2 B 2
dP
dt
YP
Г етеротроф ы
dA
A( A A )
dt
A A1 A 2
A 1 1 ˆ A 1
1 1
A
X1
X1 K
X1
X 1 K X1
,
A 2 2 ˆ A 2
X1
, 2 1
G
A G KG
G
G KG
,
A X1
X 1 K X1
О стан к и гетеротроф ов
B
dD
A A D 1
dt
YD
Б ак тери и 1
dB 1
B 1 ( B1 B1 )
dt
D
B1 ˆ B1
D KD
Б ак тери и 2
dB 2
B 2 ( B2 dt
P
B 2 ˆ B 2
P KP
А зот
dN
dt
X 1 X 2
YN
B2
)
* 1
* 1
B 1 1 D B 1 B 2 1 P B 2
YD
YD
G
G KG
Традиционным подходом для глобальных вегетационных
моделей является так называемый "клеточный" (cell grid)
метод. В представляемой модели используется более
простой в техническом плане, но не менее действенный
метод.
Этот
метод
основан
на
разделении
рассматриваемого региона на компартменты, однородные
по своим структурно-функциональным характеристикам
(биоценозы), и каждый из этих компартментов будет
рассматриваться как единое целое. Моделирование
компартмента-биоценоза
происходит
при
помощи
обыкновенных дифференциальных уравнений. Число
компартментов и границы между ними определяются
путем анализа баз данных по исследуемому региону спутниковых и других. Подобный подход позволяет
совместить как простоту описательной части, так и
возможность
глубокой
проработки
исследуемых
процессов.
Данный метод предпочтительнее еще и потому, что
имеет дело с реальным биологическим объектом –
биоценозом, тогда как традиционные методы создают
искусственную градацию.
Биоценоз как условно независимая часть
биосферы
Модельное представление биоценоза
Параметр Х
расстояние
Центр биоценоза
границы биоценоза
Схема компартментов включает границы биоценозов, а также их
условные центры. Модельный центр не обязательно физический
центр компартмента-биоценоза, это точка, в которой значения
определяющего компонента биоценоза равно среднему по всему
биоценозу. Эта точка считается центром координат при
определении координат всех других точек биоценоза (показано
пунктирными линиями). Параметры, влияющие на динамику
биомассы определяющего компонента биоценоза изменятся по
отношению к центральной точке, т.е. имеется их градиент. Форма
записи зависимости параметров может быть разной, в данном
случае выбрана простая линейная (векторная) зависимость
(x, y) .
Имея два изменяющихся параметра (направление их векторов
показано на рисунке), мы рассчитываем отклонение фитомассы в
точках А и Б от средней. Таким образом можно рассчитать
динамику любого компонента системы в любой интересующей
нас точке.
Расчет значений параметров экосистемы в любой точке
компартмента
B
Пример расчета изменения
фитомассы (условные еденицы)
при условии линейного
изменения градиентов внешних
условий для точек
A
A (x1,y1), B(x2,y2)
Градиент
Температуры
Градиент
влажности
Универсальная система для расчета экологических моделей (на
основе MatLab)
Вывод
результатов в
виде
зависимости
одной
переменной от
другой
Главное меню и
используемые
файлы
Вывод
результатов в
табличном виде
формата Excel
Ввод информации о модели –
формулы, коэффициенты, параметры,
начальные условия
Вывод результатов в
виде зависимости от
времени
60
ф и т о м а с с а , т / га
50
о б щ а я ф и то м а с с а
40
з е л е н а я ф и то м а с с а
30
20
10
0
0
500
1000
в р е м я , с у тки
Численный расчет модели –
условный лиственный лес
1500
Численный расчет модели на
больших временах влияние лесных
пожаров
53
fito m ass kg 103 h a -1
52,5
52
no c hanging
tem perature
51,5
51
s ingle inc reas e of
tem perature
50,5
50
perm anetly inc reas ing
tem perature
49,5
49
48,5
48
0
500
1000
tim e , d a ys
Численный расчет модели –
влияние изменения температуры
1500
C O 2 flu x velo sit y, g 10 3 h a -1 d ay -1
300
250
200
150
100
50
0
-50 0
200
400
600
800
-100
-150
-200
tim e , d a ys
Численный расчет модели –
скорость потока CO2
1000
1200
1400
1600
Документ
Категория
Презентации
Просмотров
9
Размер файла
433 Кб
Теги
1/--страниц
Пожаловаться на содержимое документа