РЕАКЦИЯ ДРЕВЕСНОЙ РАСТИТЕЛЬНОСТИ НА ВОЗДЕЙСТВИЕ ЭКОЛОГИЧЕСКИХ ФАКТОРОВ

 

Чжан C.А., Пузанова О.А., Чжан Л.А.

(ФГБОУ ВПО «БрГУ», г. Братск, РФ)

 

In this paper, the predictability of declining productivity of plantations in the region under the influence of  prolonged exposure to industrial emissions

 

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

Распределение насаждений  по группам и ко­личественная оценка условий местопроизрастания - это организация структурной основы математической модели. Ее содержание - взаи­мосвязи между таксационными показателями. Законченный вид мо­дель приобретает после выявления всех необходимых закономерно­стей.

В ходе исследований состояния насаждений по данным пробных площадей были выявлены 3 зоны влияния промышленных выбросов  в районе города Братска и установлены характеристики этих зон: Зона I или зона экстремального загрязнения древостоев.  Имеет вытянутую конфигурацию  по направлению преобладающих ветров, с Ю-З на С-В. Содержания фтора в хвое свыше 8 мг/кг, средний балл категории состояния 2,3. В Ю-В направлении зона заканчивается и граничит с санитарной зоной алюминиевого завода. Средний радиус зоны составляет 5-6 км. Зона II - зона сильного загрязнения древостоев. В эту зону входят древостои, интенсивно накопившие твердые и растворимые загрязнители,  содержания фтора в хвое от 2-8 мг/кг, средний балл категории состояния 1,8. Эта зона имеет сильно вытянутую конфигурацию и значительную площадь, т.е. можно отметить, что почти вся обследованная территория может быть отнесена ко второй зоне. Максимальная протяженность зоны с ЮВ на СЗ составляет  около 60-70 км. Ширина зоны составляет около 30 км. Зона III - зона слабого загрязнения древостоев. Не имеет четко выраженной конфигурации, т.к. не найдены пробные площади, которые можно было бы отнести к условно чистым или фоновым зонам.

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

Приближенную оценку закона распределения показателя условий местопроизрастания по каждой зоне можно по­лучить, используя коэффициенты асимметрии и эксцесса. Расчеты этих показателей  деревьев  обработка и статистический анализ проведены согласно методике, изложенной в работе Г. Ф. Лакина [1].

Величина классового ин­тервала рассчитывалась по формуле

i = (x_max - x_min)/(1+3.32*lg n)

Результаты вычислений коэффициентов асимметрии и эксцесса сведены в таблицу 1.

 Расчет производился по данным показателя условий местопроизрастания, полученным на основе уравнения рег­рессии вида d_1.3 = f(h).

 

Таблица 1 - Коэффициенты асимметрии и эксцесса рядов распределения показателей условий местопроизрастания  деревьев по зонам загрязнения

	Зоны загрязнения
	Экстремального загрязнения	Сильного загрязнения	Слабого загрязнения
Коэффициент ассиметрии	0.49	0.10	0.11
Коэффициент эксцесса	0.27	0.58	1.17

 

Оценка полученных корреляционных отношений: значение корреляционного отношения между высотой и возрастом изменяется от 0.731 до 0.946. Это количественно больше, чем для высоты и диаметра, где этот показатель колеблется от 0.621 до 0.894. Численные величины выборочных ошибок корреляционных отно­шений в целом меньше для показателей h и d, чем h и А. Их абсолют­ные колебания составляют соответственно 0.005-0.057 и 0.240-0.212.

На рисунках 1-2 представлены зависимости между таксационными показателями по зонам загрязнения.

Рисунок 1 – Зависимость  между диаметром и возрастом древостоев по зонам загрязнения

Рисунок 2– Зависимость  между диаметром и процентом распределения деловой древесины древостоев по зонам загрязнения

 

Для видов взаимосвязей h = ƒ(d_1.3), h = f(A) были получены уравне­ния регрессии (с аппроксимацией по типу показательной функции) и остаточные дисперсии (табл. 2).

 

Таблица 2 - Зависимости между h и A, h и d; по зонам загрязнения

Зоны загрязнения	Вид связи	Сосна
		Уравнение регрессии	σ2ост, см2
I	h = ƒ(A) h = ƒ(d1.3)	h = -0,473А2+17,181А-93,215 h = -0,007d2+0,758d+2,773	2.14 0.62
II	h = ƒ(A) h = ƒ(d1.3)	h = 0,6142A2-16,484А+163,34 h = 0,0018d2+0,42d +6,28	2.16 0.88
III	h = ƒ(A) h = ƒ(d1.3)	h = -0,0712A2+6,338А-24,941 h = 0,0029d2+0,264d+9,67	3,61 1.2

 

Математические зависимости между диаметром и возрастом древостоев можно представить следующими уравнениями:

1зона:            y = 0,2768x² - 10,768x + 154,75;  R² = 0,9304

2 зона:           y = 0,2485x² - 8,8244x + 131,38; R² = 0,7089

3 зона:           y = 0,0774x² - 1,8333x + 66,75;   R² = 0,9759

Математические зависимости между диаметром и процентом распределения деловой древесины древостоев можно представить следующими уравнениями:

1 зона:      y = -0,0015x² + 0,2381x + 69,5;  R² = 0,9524

2 зона:    y = 0,0863x² - 2,5536x + 91,5;   R² = 0,944

3 зона:       y = 0,0045x² + 0,2589x + 76,875; R² = 0,8854

На основании полученных математических зависимостей можно сделать следующие выводы: на основании многолетних наблюдений за древостоями по материалам пробных площадей накоплен большой материал, позволяющий говорить о возможности прогнозирования снижения продуктивности региона; под влиянием длительного воздействия промышленных выбросов нарушается корреляция между высотой и диаметром сосновых насаждений. Высота насаждений подверженных воздействию промвыбросов меньше, чем высота неповрежденных насаждений (на 1,5-2 м). Исключение составляют деревья диаметром 40 см и более. Здесь явно преимущество в высоте модельных деревьев перед табличными данными.

Литература

1. Лакин Г.Ф. Биометрия: Учебное пособие для биол. спец. вузов. - Изд. -4-е, перераб. и доп. -М: Высшая школа, 1990. -352 с.