Распределение деревьев сосны обыкновенной по диаметру в очагах интенсивного биогенного воздействия

РАСПРЕДЕЛЕНИЕ ДЕРЕВЬЕВ СОСНЫ ОБЫКНОВЕННОЙ ПО ДИАМЕТРУ В ОЧАГАХ ИНТЕНСИВНОГО БИОГЕННОГО ВОЗДЕЙСТВИЯ

© И.Н. Павлов, О. А. Барабанова, А.Г.Миронов УДК 630*5: 630*443.3 ГОУ ВПО «Сибирский государственный технологический университет» г.Красноярск, Россия

Представлены результаты исследования распределения деревьев сосны обыкновенной по диаметру в очагах интенсивного биогенного воздействия. Наряду с корневыми патогенами (Armillaria mellea sensu lato) в исследуемых насаждениях были отмечены высокие уровни зараженности деревьев раком-серянкой (Cronartium flaccidum (Alb. et Schw.) Wint.; Peridermium pini (Willd.) Kleb.) и заселенности Blastophagus minor Hart. В процессе усыхания образуются очаги диаметром 10-30 м, в дальнейшем сливающиеся в полосы. Несмотря на одновозрастность древостоев сосны, в распределении, как живых, так и усохших деревьев имеет место две вершины. Успешная аппроксимация оказалась возможной только с использованием многокомпонентной смеси нормальных распределений. Первая составляющая фср живых = 15,4 см) образована в результате отсутствия своевременного отпада отстающих в росте деревьев (по генотипу или произрастающие на маломощных почвах). Среди них и образуется первая составляющая сухостоя со средним значением диаметра 11 см. Отпад происходит естественно (действие патогенов вторично) в результате угнетения со стороны соседних деревьев. В дальнейшем именно на их основе происходит усиление вирулентности и агрессивности опенка. В итоге, вторая составляющая отпада фср = 21 см) образована действием опенка, как первичного патогенна, что способствует рассечению общего распределения живых деревьев на две составляющие.

The results of investigation of diameter distribution of ordinary pine trees in the focuses of intensive biogenic influence have been set out. Along with root pathogens (Armillaria mellea sensu lato), high levels of tree infectiousness with Cronartium flaccidum (Alb. et Schw.) Wint., Peridermium pini (Willd.) Kleb. and population of Blastophagus minor Hart have been noticed in the investigated plantations. In the process of drying, the focuses of 10 to 30 m in diameter are formed, later on merged into stripes. Despite the pine stands are even-aged, two peaks in the distribution of both live and dried trees take place. Successful approximation turned out to be possible only with the use of multi-component mixture of normal distributions. The first component (Daverage of ilve = 15.4 cm) has been resulted from the lack of opportune perishing of the trees decelerating in growth (according to genotype, and growing in thin-layer soils). The first component of dead-and-down with average diameter value 11 cm is formed from them. The perishing happens naturally (the effect of pathogens is secondary) resulting from the suppression from neighbouring trees. Later on, namely they become a base of Armillaria’s virulence and aggressivity enhancement. As a result, the second component of perishing (Daverage = 21 cm) is formed by the action of Armillaria as a primary pathogen, which contributes to dissection of the total distribution of live trees into two components.

Распределение числа деревьев по диаметру является важным показателем, характеризующим состояние древостоя. На него оказывает влияние большое количество факторов: возраст, условия произрастания, породный состав, полнота, характер размещения, антропогенное и биогенное воздействие, пожары, ветровалы и пр.

Форма кривой распределения диаметров определяется видом и интенсивностью естественного изреживания и является его характеристикой (Свалов, Свалов, 1973).

Определенную закономерность в возрастной динамике рядов распределения диаметров стволов отмечают А.А. Макаренко (1975);

Н.Н. Свалов, С.Н. Свалов (1973);

В.В. Кузьмичев, Ю.Н. Савич (1979);

В.Я. Каплунов (1989); А.И. Бузыкин и др.

(2002).

Согласно современным представлениям

о гомеокинезе биологические системы стремятся не к абсолютной стабилизации своего состояния, а к удержанию его в определенных пределах при непрерывном колебательном изменении параметров вокруг определенных точек равновесия (Свирежев, Логофет, 1978). Одним из механизмов

поддержания стабильности в экосистеме является организация отрицательной обратной связи (Одум, 1986). На разных возрастных этапах максимальный запас древостоя формируется различным количеством деревьев. Естественный отпад не только удаляет из древостоя генетически худшие экземпляры, но и создает структуру, обеспечивающую устойчивость ценоза. В свою очередь, загущенность среди деревьев нижних ступеней толщины и, как следствие, высокая положительная асимметрия выполняют роль обратной связи,

определяющей, какое количество деревьев необходимо элиминировать из древостоя.

При определенном сочетании факторов (высокое плодородие почвы в сочетании с малой мощностью; равномерное размещение возобновления и его одновозрастность) происходит нарушение своевременного отпада. При произрастании леса на маломощных дерновых почвах,

подстилаемых на глубине 20-30 см твердыми горными породами, корневая система сосны развивается в верхнем, легкопересыхающем слое (исследуемое нами насаждение). В итоге стержневой корень отсутствует или сильно редуцирован. До определенного времени насаждение развивается как высокопродуктивное. В дальнейшем после превышения эдафического потенциала поддерживать необходимую потребность насаждения в элементах питания и влаге наступает снижение устойчивости фитоценоза в целом. В первую очередь страдают деревья, располагающиеся в условиях жесткой конкуренции, на маломощных почвах и с плохо развитой корневой системой. Особую опасность для них представляет наступление даже кратковременного засушливого периода. При этом отрицательные последствия могут не проявляться длительное время. Как правило, недостаток влаги в почве ведет к потере части ассимиляционного аппарата, тонких корней и как следствие - снижение устойчивости к болезням и вредителям. При этом с момента заражения до гибели может пройти достаточно длительный промежуток времени. Анализ годичных колец 80-110летних погибших хвойных показал, что деревья были инфицированы

дереворазрушающими грибами комплекса Armillaria mellea sensu lato 50 лет назад (Shaw, Kile, 1991).

Внедрение столь сильного патогена безусловно не может не сказаться на строении древостоя. Чаще всего подобное воздействие ведет к усложнению структуры, требующей для своего описания оригинальных функций.

Для описания случайных величин в прикладных исследованиях используется около 30 законов и семейств распределений (Лемешко, 1998). Естественно, что это не покрывает существующего многообразия распределения случайных величин. Простые модели часто не могут описать рост числа случайных величин, возникающих на «хвостах» распределений, которые

определяются другими параметрами (например, почвенная мозаичность; развитие болезней и вредителей, экстремальное климатическое воздействие).

Сложные биологические, технические, социальные системы включают совокупности элементов, принадлежащих, как правило, нескольким классам. Примером такой совокурности является система, состоящая из нормально функционирующих (1-й класс) и ослабленных (2-й класс) объектов (элементов). Каждый класс объектов можно охарактеризовать некоторой функцией плотности распределения, описывающей закон распределения анализируемого признака. Распределение в такой системе может быть унимодальным, бимодальным и полимодальным.

Программа и методика исследований

В нижнем течении р. Караульной (ООПТ зеленая зона г. Красноярска) на прилегающих сопках в последнее десятилетие наблюдается интенсивная куртинная гибель сосны обыкновенной. Усыханию подвержены деревья всех классов Крафта (рисунок 1). Погибают экземпляры, не испытывающие каких-либо признаков угнетения, с хорошо развитой кроной. Гибель наступает очень быстро. Образуются очаги диаметром 10-30 м, в дальнейшем сливающиеся в полосы. Часто наблюдается гибель деревьев по периметру существующих окон в древостое. Куртинное усыхание, наличие плодовых тел АгтШайа те11еа 8еп8и 1а1о, а также ряд других диагностических признаков (Павлов, Миронов, 2004) свидетельствуют об

интенсивном развитии армилляриоза. Исключительно быстрая гибель отдельных деревьев сосны, ослабленных корневым патогеном, вызвана их интенсивным заселением малым сосновым лубоедом. Для усохших деревьев также характерна более высокая частота встречаемости рака-серянки.

Для изучения кольцевых процессов распространения агрессивных видов грибных патогенов на полигоне 10,8 га было заложено 13 круговых ППП общей площадью 1,17 га. Использование круговых постоянных пробных площадей с картированием из центра через измеряемые угол и расстояние (Павлов, 2005) имеет значительные

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

координат. Исследуемый древостой сосны обыкновенной (9С1Ос + Б ед. Лц,) 5 класса возраста имеет полноту 1,0 (2 бонитет).

Изучение изменения в строении

древостоев в условиях интенсивного

биогенного воздействия представляет не только теоретический, но и практический интерес.

Впервые алгоритм разделения совокупности на две нормальные

составляющие методом моментов предложил К. Пирсон. В дальнейшем метод был упрощен, появились работы по разделению смеси биноминальных распределений, Пуассона, Вейбулла (Никитин,

Швиденко, 1978).

Рисунок 1 - Куртинное усыхание сосны обыкновенной

Уравнение функции плотности смеси вероятностных распределений ^х),

описывающее закон распределения анализируемого признака во всей

(объединенной) генеральной совокупности имеет вид:

^х) = ЕР] х Г|

]=1

к --------------------

ЕР; = 1, Р] > 0, ] = 1,к,

н

где р] - удельный вес (доля) _] - го компонента смеси (априорные вероятности появления наблюдений из ] - го класса объектов), £](х) = Ф^а^х) - плотность

распределения анализируемого признака для _]- го класса, а] - векторные параметры, от

значений которых зависит закон распределения _]- го компонента смеси,

] = 1, к - номер класса, к - число классов.

Результаты и обсуждение

Апробировано шесть чаще всего применяемых распределений (нормальное, бета, Вейбулла, логнормальное, Пуассона, экспоненциальное). Насколько хорошо выборка описывается теоретическим

законом, проверяли с использованием различных критериев согласия (Колмогорова, Лемана-Розенблата, Пирсона). Основная идея критериев согласия состоит в измерении расстояния между функциями эмпирического и теоретического распределения

(Орлов, 2002).

В качестве контроля было выбрано сходное по таксационно- лесоводственным показателям насаждение без признаков заболевания, примыкающее к древостою с признаками армилляриоза (рисунок 2).

Асимметрия распределения по диаметру живых деревьев составляет 0,17±0,09. Гипотеза о нормальности представленного распределения отклонена на основании критерия Шапиро-Уилкса (р=0,00008). По критерию Колмогорова-Смирнова гипотеза

была отклонена при несколько меньшем уровне вероятности (р<0,1).

Асимметрия распределения по диаметру сухостоя составляет 0,45. Однако гипотеза о нормальности распределения не может быть отклонена на основании критериев Шапиро-Уилка (р=0,42) и Колмогорова-Смирнова (р>

0.20). Причина противоречия заключается в малом числе наблюдений.

со

ф

_а

т

ф

ф

о

5

90

80

70

60

50

40

30

20

10

0

-10

8 12 16 20 24 28 32 36 40 44 48 52 56 60

2

Диаметр, см

Рисунок 2 - Распределение деревьев сосны обыкновенной в контроле

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

8 3 8 Я Я Я

«I 5§

^ ^ нн

к

к

са

н

о

н

н

о

о

о

X

_ г

нн О

-а Л

со о £

о Л (3 «

я й

(н Я

й я

и и

а ч

и и « ч: о <и £ & Л

Л

О

X

ж

о

Критерий согласия Показатели Число компонентов нормального распределения

живые деревья сухостой

1 2 1 2

Колмогорова (О) Орас 0,0378 0,0206 0,1468 0,0373

а 0,200 0,200 0,200 0,200

Окр(а) 0,0440 0,0440 0,0879 0,0879

Вывод: гипотеза не отклоняется не отклоняется отклоняется не отклоняется

£ 15 ^ 8 (Кю2)р3с 0,236 0,047 0,635 0,033

а 0,500 0,500 0,500 0,500

0,118 0,118 0,118 0,118

Вывод: гипотеза отклоняется не отклоняется отклоняется не отклоняется

О. СЗ V с8 ^ Храс 42,555 20,458 68,634 15,458

число параметров, с 2 5 2 5

Число степеней свободы 15 12 15 12

а 0,050 0,058 0,050 0,217

X к Р 2 24,996 20,512 24,996 15,465

Вывод: гипотеза отклоняется не не отклоняется отклоняется отклоняется

Для нарушенных древостоев, при проверке согласия (по критериям Колмогорова, Лемана-Розенблата, Пирсона) по значениям достигнутого уровня значимости, ни одна из функций не может быть принята.

Поэтому для аппроксимации было решено воспользоваться

многокомпонентной смесью нормальных распределений:

1 f (х - )2

ґ|(х) =

ехр

л

28-

I = 1,к

V Ч у 2 (6) где XсрJ , SJ - оценки математического ожидания и дисперсии j - й компоненты к компонентной функции плотности смеси нормальных распределений древостой сосны обыкновенной по диаметру х.

В таблице 2 представлены результаты проверки согласованности эмпирической и теоретических функций смесей

распределений по критериям согласия Колмогорова (Б), Лемана - Розенблатта (ю2), Пирсона (х2). Учитывая рекомендации (Орлов, 2002; Справочник по ..., 1989, 1990) проверку согласованности функций

распределения проводили при уровнях значимости а1 = 0,2 - для критерия

Колмогорова, а2 = 0,5 - для критерия Лемана - Розенблатта, а3 = 0,05 - для критерия Пирсона.

На рисунке 3 представлено распределение по диаметру деревьев, произрастающих в очагах заболевания, с выделением живых и сухостойных.

Погибшие деревья образуют ряд распределения с двумя вершинами. Данная закономерность имеет биологическую интерпретацию. Экземпляры, образующие первое нормальное распределение со средним значением диаметра 11 см, усыхают естественно (действие патогенов вторично) в результате угнетения со стороны соседних деревьев. Доля первой составляющей (р) несколько меньше второй (0,42). В дальнейшем усиление вирулентности и агрессивности опенка связано с его

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

образована действием опенка, как

первичного патогена.

Асимметрия распределения по диаметру составляет 1,09±0,19. Гипотеза о

нормальности представленного

распределения отклонена на основании критериев Шапиро-Уилкса (р=0,0000) и Колмогорова-Смирнова (р<0,01).

Диаметр, см

------1 ------2

Рисунок 3 - Распределение числа деревьев по диаметру

Живые деревья, произрастающие на круговых площадках в очагах усыхания, также образуют ряд распределения с двумя вершинами. Первая составляющая (Оср = 15,4 см) образована угнетенными деревьями (доля числа наблюдений - 14,5 %), вторая (Оср = 30 см) - состоит из особей входящих в основной полог. Деление на два класса вызвано патогенным действием опенка. Вторая вершина сухостоя (О ср=21 см) свидетельствует о направлении и объектах воздействия биотрофов, что способствует рассечению общего распределения живых деревьев. Асимметрия распределения живых деревьев по диаметру составляет 0,01. Несмотря на столь близкое значение к нормальному распределению, критерии Шапиро-Уилкса (р=0,0003) и Колмогорова-Смирнова (р<0,01) не подтвердили гипотезу о нормальности представленного

распределения.

Заключение

Наличие двухвершинного (и более) распределения является закономерным в разновозрастных многоярусных

насаждениях. В описываемых случаях речь идет об одновозрастном насаждении сосны обыкновенной с примесью березы и осины. Для одновозрастного древостоя

светолюбивой сосны это не является типичной формой. И действие патогенов является следствием избыточной густоты среди деревьев из низших ступеней толщины и направлено на формирование более симметричного распределения.

Е.Л. Маслаков (1984) указывает, что ранговая структура популяций с возрастом насаждений меняется мало. Относительное сходство в строении древостоев было выявлено еще работами Фекете, Шиффеля, Вейзе и других лесоводов. Н.В. Третьяковым (1927) сформулирован закон единства в строении насаждений. Хотя в ряде работ (Дыренков, 1969; Никитин, 1972;

Лебков, 1973)

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

деревьев-лидеров (30%), средних -

спутников

(Маслаков, 1984). Данное соотношение

(близкое к нормальному) являясь результатом сложных взаимодействий

генотипа, ценотических условий,

конкуренции, эдафической мозаичности,

будет наилучшим для оптимального

использования ресурсов среды (солнечного света, влаги, элементов питания и пр.). Судя по распределениям, показанным на рисунке

3, у насаждений, подверженных воздействию дереворазрушающих грибов биотрофного комплекса, отсутствуют оптимальные

пропорции, которые отмечает Е.Л. Маслаков. Отсутствие своевременной активной

дифференциации среди отстающих деревьев, с последующим отпадом наименее

приспособленных сопровождается общим ослаблением насаждения и, при наличии инфекции грибов биотрофного комплекса, распадом. В случае интенсивного

воздействия АгтШайа те11еа 8еи8и Мо не происходит естественного и положительного для древостоя отпада только наиболее мелких деревьев. Усыхание происходит

куртинами до 10 -15 деревьев. При этом погибают деревья всех классов Крафта.

Библиографический список

1. Бузыкин, А.И., Пшеничникова, Л.С.,

Суховольский, В. Г. Густота и

продуктивность древесных ценозов. -

Новосибирск: Наука, 2002. -152 с.

2. Дыренков, С.А. Закономерности

дифференциации деревьев по толщине и их математическая интерпретация //ЭВМ и математические методы в лесном хозяйстве. -Л.: ЛенНИИЛХ, 1969. -С. 63-75.

3. Каплунов, В.Я. Прогнозирование строения древостоев по диаметру // Лесоведение. -1989. -№5. -С. 68-73.

4. Кузьмичев, В.В., Савич, ЮН. Влияние густоты посадки на рост сосновых культур // Лесоведение. -1979. -№6. -С. 56-63.

5. Лебков, В.Ф. Теория строения

древостоев и ее роль в оценке продуктивности лесов //Вопросы

лесоведения. -Красноярск: Наука, 1973. -С. 90-104.

6. Лемешко, Б.Ю. Асимптотически оптимальное группирование наблюдений в критериях согласия //Заводская лаборатория. -1998. -Т. 64. - №1. - С.56-64.

7. Макаренко, А.А. О свойствах рядов распределения деревьев в древостоях // Лесоведение. -1975. -№6. -С. 42-49.

8. Никитин, К. Е. Применение ЭВМ в лесной таксации. -М.: Урожай, 1972. - 133 с.

9. Никитин, К.Е., Швиденко, А.З.

Методы и техника обработки лесоводственной информации. -М.: Лесн.

пром-ть, 1972. -272 с.

10. Одум, Ю. Экология. Т.1. -М.: Мир, 1986. -328 с.

11. Орлов, А.И. Эконометрика. - М.: Экзамен, 2002. - 576 с.

12. Павлов, И.Н. К методике

картографирования куртинного усыхания деревьев с применением

геоинформационных технологий // Лесное хозяйство - 2005. - № 4 - С.38-39.

13. Павлов, И.Н., Миронов, А.Г. Куртинное усыхание хвойных насаждений в измененных условиях среды (вопросы диагностики) //Непрерывное экологическое

образование и экологические проблемы: Сборник статей по материалам Всероссийской научно - практической конференции. Т.1. -Красноярск: СибГТУ,

2004. -С. 173-185.

14. Свалов, Н.Н., Свалов, С.Н. О динамике рядов распределения диаметров стволов в одновозрастных сосняках // Лесоведение. -1973. - №5. -С. 58-62.

15. Свирежев, Ю.М., Логофет, Д.О. Устойчивость биологических сообществ. -М.: Наука, 1978 . -352 с.

16. Справочник по прикладной

статистики. Т.1.: Пер. с англ./Под ред. Э. Ллойда, У. Ледермана, Ю.Н. Тюрина. - М.: Финансы и статистика, 1989, - 510 с.

17. Справочник по прикладной

статистики. Т.2.: Пер. с англ./Под ред. Э. Ллойда, У. Ледермана, С.А. Айвазяна, Ю.Н. Тюрина. - М.: Финансы и статистика, 1990. -526 с.

18. Ушанов С.В. Применение

многомерных статистических методов при принятии решений. - Красноярск: СибГТУ, 2003. - 239 с.

19. Феллер, В. Введение в теорию

вероятностей и ее приложения. Т. 2. -М.: Мир, 1984. - 738 с.

20. Феллер, В. Введение в теорию

вероятностей и ее приложения. Т.1. - М.: Мир, 1984. - 528 с.

21. Shaw C. G., Kile G. A. Armillaria root

disease. Agricultural handbook №691. Forest service US Department of Agriculture. -Washington D.C., 1991. -233 р.