CC BY
60
БИОЛОГИЧЕСКИЕ НАУКИ И ФИЗИКА
УДК 574 ББК 28.080.3
Володченкова Людмила Александровна
^искатель
кафедра социально-экономической географии и туризма Омский государственный педагогический университет
г. Омск
Гуц Александр Константинович
доктор физико-математических наук, профессор кафедра кибернетики Омский государственный университет г. Омск
Volodchenkova Lyudmila Alexandrovna
Applicant for a Degree Chair of Social and Economic Geography and Tourism
Omsk State Pedagogical University Omsk
Guts Alexander Konstantinovich
Doctor of Mathematics and Physics,
Professor Chair of Cybernetics Omsk State University Omsk
Прогнозирование состояний соснового леса при нарастании засухи или пожарной опасности Predictive Modeling of Pine Forest State in Case of Drought or Fire Danger
Дается описание равновесных состояний соснового леса в процессе развития засухи или пожарной опасности на основе математической теории катастроф. Показано как можно прогнозировать будущие состояния леса до пожара и после пожара.
The description of equilibrium states of pine forest in case of drought or fire danger on the basis of the mathematical theory of catastrophes is given in the article. It is shown how to predict the future states of forest before and after a fire.
Ключевые слова: экология, лесной биоценоз, засуха, пожар, прогнозирование.
Key words: ecology, forest biogenesis, drought, fire, predictive modeling.
Целью статьи является изучение и прогнозирование равновесных состояний леса до и после пожара засухи (пожара). Для теоретического описания наблюдающегося экологического кризиса мы используем математическую теорию катастроф Рене Тома [5].
В качестве возможных внешних факторов, влияющих на лесную экосистему, взяты влажность почвы * , мозаичность леса т , наличие конкуренции к и степень пожарной опасности (равно - степень засухи) а .
В работах [2,3] было показано, что четырехярусному лесу можно поставить в соответствие так называемую потенциальную функцию вида
а
V(х,к,т,а,м>) = — (х - х0)6 + к(х - х0)4 + т(х - х0)3 + а(х - х0)2 + *(х - х0),
6
где х - продуктивность фитомассы, Первый член а(х - х0) /6 (а > 0) определяется наличием только четырех ярусов леса. Учет каждого нового яруса увели-
х на единицу Коэффициент а = а\а2аза4 где а1
чивает показатель степени х на единицу. Коэффициент _ «\«2«3«4, где доля фитомассы ]-го яруса в фитомассе всего леса.
Далее
52
к = -ск(С1 - С10),т = ст(-1),а =-с^(КП - КП0),* = ^(Ж -
М
где
С1 - индекс конкуренции [\];
5 2 / М
- коэффициент дисперсии, являющий показателем равномерности
5 2 / М
распределения деревьев в пространстве; если близко к нулю, то распре-
деление регулярное, к единице - случайное, а чем более единицы, - тем контагиознее (пятнистее), т.е. мозаичнее;
КП - комплексный показатель пожарной опасности В.Г.Нестерова текущего дня, рассчитываемый по формуле:
КП = 1 Т (Т1 - т;), і=1
т тг
где 1 - температура воздуха 1 -го дня; 1 - температура точки росы 1 -го дня;
1 =1 - день, когда был дождь, 1 =1’’"’п - дни без осадков, п - число дней после последнего дождя [4];
Ж
влажность почвы;
постоянные коэффициенты.
Величины С10, КП 0,Ж() - это критические значения факторов, обозначающие границы экологической устойчивость фитоценоза. Через х° обозначим минимально допустимую для изучаемого типа леса продукцию фитомассы.
Увязывание антропогенного фактора с комплексным показателем В.Г.Нестерова объясняется тем, что возникновению пожара по вине людей способствует засушливая погода. Комплексный показатель В.Г.Нестерова как раз и вычисляется с целью оценивания степени пожарной опасности в лесу по условиям погоды. Очевидно, что одновременно наша модель хорошо описывает засуху.
Рассмотрим сосновые насаждения в южных экорегионах Западной Сибири, который характеризуется следующими значениями параметров до катастрофы (пожара, засухи):
а = 90 • 5 •3 •2; С1 =12, С10 - 6 (сосна, [\]); 5 /М =1 (случайное распределе-
х — 8 5
ние деревьев); 0 ’ т/га за год (согласно [6] для сосновых насаждений Ом-
ской области); Ж = 28% Ж° = 27% .
Таким образом, мы зафиксировали все факторы, кроме показателя пожарной опасности КП.
«Общероссийская шкала имеет пять классов пожарной опасности в лесу по условиям погоды (табл.\). Для отдельных регионов разработаны региональные шкалы пожарной опасности в лесу по условиям погоды, учитывающие местные особенности и те, в которых значения комплексного показателя по классам отличаются от значений общероссийской шкалы» [4].
Таблица 1.
Шкала пожарной опасности в лесу по условиям погоды
Класс пожарной опасности в лесу по условиям погоды Значение КП Степень пожарной опасности
I До 300 —
II От 30\ до \000 Малая
III От \00\ до 4000 Средняя
IV От 400\ до \0000 Высокая
У От \0000 Чрезвычайная
Исходя их этой таблицы, принимаем, что КП0 =10003,6 .
Для установления равновесий, возникающих при нарастании пожарной опасности для соснового насаждения, будем увеличивать комплексный показатель пожарной опасности КП, начиная с IV класса. Иначе говоря, моделируем процесс пожарной опасности в лесу, меняя показатель КП от значения 4003,6 до \2003,6.
Компьютерные эксперименты показали, что можно принять следующие
^ ск = 472,5, ст — \, са — \000, с„ — 2 • \03,
значения для коэффициентов: к т а * (каждый
в соответствующих единицах измерения).
Потенциал соснового леса имеет вид
V(х,к,т,а,*) — 450(х - 8,5)6 - 2835(х - 8,5)4 -- \000(КП - КП0)(х - 8,5)2 + 2 • \03(Ж - Ж0)(х - 8,5) (\)
С точки зрения математики равновестия лесной экосистемы - это либо точка минимума, либо точка максимума, либо так называемая точка перегиба
~ д, V — V,кта*)(х) — V(х,к,т,а,*) ТГ
потенциальной функции (к4 ' 4 7. Как правило, устойчивые
равновесия системы - это минимумы функции V — V(x, к, ^ a, *). На рисунках
графика функции V — к,т,а,м,) (х) они изображаются ямками (соответственно максимумы - неустойчивые равновесия - изображаются вершинами горок).
Для конкретного набора внешних факторов у функции V к,т,а,№) (х) мо-
х х
жет быть несколько точек минимумов. Скажем, это точки ’• •’ . В каком из
этих равновесий находится система? Как узнать это? Для этого используют различные правила, т.е. способы, по которым выбираются равновесия. Рассмотрим два основных правила.
Правило Максвелла предписывает взять в качестве равновесия системы та-
(к ,т,а, м)
(х)
достигает наименьшего
кую точку минимума, в которой функция значения.
Правило максимального промедления предписывает состоянию оставаться в минимуме при изменении факторов до тех пор, пока он существует.
В момент, когда происходит исчезновение старого минимума и становится необходимым переход в новый, происходит экологическая катастрофа.
Окончательный результат компьютерных экспериментов представлен в виде графиков потенциала леса на рис.1-6.
КП =4003.6.
Ж = Жп +1%.
КП =9998.6,
Ж = Ж +1%.
Продукция х=8,5.
Продукция х=8,5.
Рис.3. Класс пожарной опасности IV.
КП =10000,
Ж = Ж0 +1%.
КП =10006,6,
Ж = Жп +1%.
Продукция х=8,5.
Продукция х=6,3.
Рис.5. Класс пожарной опасности V. Рис.6. Класс пожарной опасности V.
КП =10203,6.
Ж = Ж0 +1%.
КП =12003,6,
Ж = Ж0 +1%.
Продукция х=4,6.
Продукция х=1,8.
Мы видим, что засуха (пожар) имела место при V классе пожарной опасности, поскольку резко, скачком, уменьшилась продуктивность леса (рис.4). По мере дальнейшего увеличения комплексного показателя пожарной опасности КП , т.е. по мере развития засухи (пожара), продуктивность падает с 8.5 до 1.8. Это следует расценивать как итог засухи (или лесного пожара). При этом экосистема никак не могла оказаться в правой ямке (рис.5), хотя это и устойчивое равновесие. На рис.6 ямки имеют уже одинаковую глубину. Если считать, что, например, пожар миновал, имеет место гарь, и воспользоваться правилом Максвелла, то правая ямка - это возможное альтернативное состояние леса. Поскольку оно соответствует более высокой продуктивности, то можно говорить, что переход в правое состояние - это восстановление леса после пожара.
Потенциалом леса (1) можно пользоваться для прогноза. Состояние леса -это точка на кривой графика потенциала с соответствующими значениями
С1,КП,Ж,х . По ней можно судить, является ли состояние леса устойчивым, далеко ли до устойчивого, и к какому устойчивому состоянию оно устремится в ходе сукцессии после засухи или послепожарной сукцессии.
Библиографический список
1. Вайс, А. А. Классификация деревьев и горизонтальная структура ценозов [Электронный ресурс] / А. А. Вайс. // Научный журнал КубГАУ. - 2007. - №.31(7). - С.1-13. - Режим доступа: http://ej.kubagro.ru/2007/07/pdf/14.pdf.
2. Володченкова, Л.А. Описание и моделирование катастроф лесных биоценозов [Текст] / Л.А. Володченкова, Н.А. Калиненко // Вестник Челябинского государственного педагогического университета. - 2009. - № 12. - C.298-309.
3. Володченкова, Л.А. Катастрофы типа «бабочка» в лесных экосистемах [Текст] / Л.А. Володченкова, А.К. Гуц // Математические структуры и моделирование. - 2009. -Вып.19. - С.48-68.
4. Воробьев, Ю.Л. Лесные пожары на территории России: Состояние и проблемы [Текст]/ Ю.Л. Воробьев, В.А. Акимов, Ю.И. Соколов; Под общ. ред. Ю.Л.Воробьева; МЧС России. - М.: ДЭКС-ПРЕСС, 2004. - 312 с.
5. Постон, Т. Теория катастроф и ее приложения [Текст]/ Т. Постон, И. Стюарт. - М.: Мир, 1980. - 570 с.
6. Швиденко, А. Модели и таблицы биологической продуктивности [Электронный ресурс]/ А. Швиденко, Д. Щепащенко, С. Нильссон, Ю. Булуй // Леса и лесное хозяйство России. - Режим доступа: http://www.iiasa.ac.at/Research/FOR/forest_cdrom/home_ru.html.
Bibliography
1. Poston, T. Theory of Catastrophes and Its Applications [Text] / T. Poston, I. Stewart I. - М: Mir, 1980. - 570 p.
2. Shvidenko, А. Models and Tables of Biological Productivity [Electronic Resource] / A. Shvidenko, D. Shchepashchenko, S. Nilsson, Yu. Buluy // Forests and Forestry of Russia. - Access Mode: http://www.iiasa.ac.at/Research/FOR/forest_cdrom.
3. Volodchenkova, L.A. The Description and Modelling of Forest Biocenosis Disasters [Text] / L.A. Volodchenkova, N.A. Kalinenko // Bulletin of Chelyabinsk State Pedagogical University. -2009. - № 12. - P.298-309.
4. Volodchenkova, L.A. “Butterfly” Catastrophes in Forest Ecosystems [Text] / L.A. Volodchenkova, A.K. Guts // Mathematical Structures and Modelling. - 2009. - №.19. - P.48-68.
5. Vorobyev, Yu.L. Forest Fire on the Territory of Russia: State and Problems [Text] / Yu.L.Vorobyev, V.A. Akimov, Yu.I. Sokolov. - M.: DEKS-PRESS, 2004.
6. Weiss, A.A. Classification of Trees and Horizontal Cenosis Structure [Electronic Resource] / A.A. Weiss // Scientific Journal of KubSAU. - 2007. - №.31 (7). - P.1-13. - Access mode: http://ej.kubagro.ru/2007/07/pdf/14.pdf.
