Влияние взаимного расположения деревьев на ветроустойчивость одноярусных древесных насаждений Текст научной статьи по специальности «Сельское хозяйство, лесное хозяйство, рыбное хозяйство»

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

© Н.В. Филенкова, Ю.В. Захаров, В.Г Суховольский УДК 630*181.23

ГОУ ВПО «Сибирский государственный технологический университет» г. Красноярск, Россия

Работа поддержана РФФИ грант №06-04-48303.

На основе проведенных натурных исследований рассмотрено влияние взаимного расположения деревьев на ветроустойчивость одноярусных древесных насаждений. Исследования проводились на двух пробных площадях: в горных лесах Саян (Восточная Сибирь) и на территории Инской ленты Минусинских ленточных боров (юг Красноярского края). В качестве показателей, характеризующих взаимодействие деревьев в насаждении, использовались параметры: число связей с соседними деревьями n и относительная высота дерева Hrei. Введены парциальные функции экологической полезности для оценки ветроустойчивости одноярусных насаждений в зависимости от взаимного расположения деревьев. Вычислены оптимальные значения числа связей с соседними деревьями и относительной высоты над пологом леса, характеризующие максимальную устойчивость дерева к ветровому воздействию.

The influence of mutual position of trees on even-ages tree stands wind stability is considered. Field studies were conducted in mountain forests of the Sayan (East of Siberia) and on the territory of Inskoy stripe of Minusinsk pine forest (South of Krasnoyarsk region). The nearby surroundings of the tree n and the relative height of tree under the stand canopy Hrd are used as parameters which characterized trees interaction on the stand. The functions of partial ecological utility for the estimation of the even-ages tree stands wind stability depending on mutual position of trees are introduced. The nearby surroundings of the tree and the relative height of tree under the stand canopy which characterized the maximum wind stability of tree under the wind impact are calculated.

Ветер оказывает сильное воздействие на лесные насаждения, а хозяйственная деятельность в лесу и климатические изменения повышают риск вывала дерева под воздействием ветра. Для оценки ветроустойчивости отдельного дерева ранее была введена модель, описывающая реакцию дерева на ветровое воздействие с учетом влияния структуры корневой системы [1]. Однако ветроустойчивость насаждения характеризуется не только суммой показателей ветроустойчивости отдельных деревьев в лесу, но и вследствие кооперативных эффектов взаимодействия

деревьев в насаждении такая устойчивость определяется также пространственной структурой насаждения и ландшафтно-экологическими условиями произрастания.

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

ветроустойчивость одноярусных древесных насаждений.

Исследования проводились на двух пробных площадях. Пробная площадь № 1 размером 3,7 га располагается в горных лесах Саян (Восточная Сибирь) и находится на высоте около 880 м над уровнем моря. В состав насаждения входили пихта сибирская Abies sibirica L. и кедр сибирский (сосна кедровая) Pinus sibirica L., породный состав насаждения 9П1К, преобладают деревья 2 класса роста. С 1979 г. с интервалом в несколько лет производилось картирование. В

июле 2003 г. проводилось картирование местоположения деревьев и определялись их координаты на пробной площади, высота, измерялся диаметр на уровне 1,3 м, размеры проекции кроны по сторонам света.. В этом же году после картирования под воздействием ветра произошел вывал и излом около 10% деревьев. Новое картирование деревьев на пробной площади, включая картирование деревьев, вываленных ветром, проводилось в июле 2004 г.

Пробная площадь № 2 размером 4,1 га находится на территории Инской ленты Минусинских ленточных боров (юг Красноярского края). Насаждение представляло собой чистый сосняк (Ріпш ,?у!ує,?ґгі,? Ь.) в возрасте 100 - 120 лет. Картирование

местоположения деревьев было проведено в 1999 г. и в 2000 г. после прошедшего урагана, когда произошел вывал и излом около 50% деревьев на пробной площади.

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

Пересекающиеся прямоугольники описывают перекрытие крон деревьев.

Одним из основных факторов, влияющих на сибирской характерен поверхностный тип

вероятность вывала дерева, является тип его корневой системы, а для кедра сибирского и

корневой системы. Известно, что для пихты сосны обыкновенной - стержневой [2].

20

10

+

20

□ □ а

10

X

—I---------------------------------------------1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1—і—I—і—г-

+

+

+

і

+

-=М

-К

-1--------1----------1-1-----------1-1-------1-------1---------1-------1-1-----------1-1-----------1-Г"

10 20 0 Ю 20

Рисунок 1 - Фрагменты плана насаждения: а - пробная площадь №1; Ь - пробная площадь №2

В качестве показателей, характеризующих взаимодействие деревьев в насаждении,

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

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

теоретическое уравнение неизвестно, была построена аппроксимирующая кривая в виде

Р = -

а-.

а2 + п

Функция выбиралась дробно-

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

связей с соседними вероятности вывала

деревьями

при

Рисунок 2 - Зависимость вероятности вывала дерева от числа связей п с соседними деревьями (столбик - стандартная ошибка данных): а - пробная площадь №1; Ь - пробная площадь №2

Из рисунка 2 видно, что вероятность вывала дерева уменьшается с увеличением числа соседей. Это наблюдается как в случаях вывала малого числа деревьев на пробной площади №1, так и в случае массового вывала деревьев на пробной площади №2.

Еще одним параметром, определяющим влияние взаимодействия деревьев в насаждении на ветроустойчивость, является относительная высота дерева над пологом насаждения Нге1. Для каждой пробной площади по данным наблюдений была построена зависимость вероятности вывала дерева от относительной высоты дерева над пологом насаждения (рисунок 3). Можно полагать, что в предельном случае при Нге1 <

0 вероятность Р вывала должна равняться нулю, а при Нге1 ^ да Р ^ 1. Кроме того, предполагается, что зависимость между Р и Нге1 монотонна.

Из разнообразных функций,

удовлетворяющих этим условиям, для аппроксимации данных натурных

наблюдений была выбрана дробно-

Н т

рациональная функция Р =---------——, где

А + Н—

коэффициенты А и т определялись из условия минимума квадратов невязки. Величина А характеризует относительную высоту дерева над пологом насаждения Нге1ш при вероятности вывала 'А

Рисунок 3 - Зависимости вероятности вывала дерева от относительной высоты дерева над пологом леса (столбик - стандартная ошибка данных): а - пробная площадь №1; Ь - пробная площадь №2

На рисунке 3 видно, что для обеих площадей вероятность вывала дерева с высотой, меньшей высоты полога леса, близка к нулю. На пробной площади №2 для деревьев с высотой, близкой к высоте полога леса, вероятность вывала достаточно высока (фактически произошел массовый вывал таких деревьев), а для таких же деревьев на пробной площади №1 вероятность вывала мала. Как видно, форма графиков на рис. 3 качественно совпадает, это говорит о независимости этого кооперативного эффекта от места и условий произрастания дерева.

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

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

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

экологической полезности [3]. Парциальная функция экологической полезности

кооперативного эффекта устойчивости дерева в зависимости от относительной высоты Нге1 имеет следующий вид:

г = 1 --

H

гєі

л+нгв—

где А и т - константы, определяемые из натурных данных.

Функция экологической полезности конкуренции деревьев за свет минимальна, когда превышение высоты дерева над соседями стремится к нулю. В качестве функции парциальной экологической полезности конкуренции деревьев за свет, в зависимости от относительной высоты дерева над пологом леса Нге1, можно использовать функцию

г = 1 --

В

В + Не

Парциальная функция экологической полезности кооперативного эффекта устойчивости дерева в зависимости от числа связей с ближайшими соседями п имеет следующий вид

Г3 =1 -

І

І 2 + П

где *\, и к- константы, определяемые из

натурных данных.

Функция экологической полезности конкуренции деревьев за питание максимальна,

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

п

Б + п2 '

Общую функцию экологической полезности Я ( к, п ) можно записать как произведение функций парциальной

полезности: Я{Нге1, п) = г • г2 • г3 ■ г4.

Согласно принципу максимальной

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

Константы для пробной площади № 1: *1 = 2,01, ь = 9,25, Л = 6,3, Б = 600, В = 2.

Оптимальные значения: п = 4, Нге1 =0,96.

Константы для пробной площади № 2: *1 = 33,33, *2 = 49,32, Л = 1,22, Б = 30, В = 1,5.

Оптимальные значения: п = 4, Нге1 =0,73.

Рисунок 4 - Общая функция экологической полезности: а - пробная площадь №1; Ь - пробная площадь №2

т

Сравнивая полученные оптимальные

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

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

характеризующих устойчивость дерева к вывалу, которые лежат в окрестности

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

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

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

можно ввести представление об интервалах таких относительных изменений,

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

пропорционально логарифмической

производной функции полезности по параметру

1 дЯ

п / (п ) =

Я дп

Рисунок 5 - Зависимость / (п) от числа связей с соседними деревьями п:

а - пробная площадь №1 (/’ = -0.019); Ь -пробная площадь №2 (/’ = -0.104)

На рисунке5 видно, что устойчивая конфигурация с числом связей с соседними деревьями, лежащим в интервале 3 < п < 5, для пробной площади №1 имеет относительные изменения величины функции полезности в 5 раз меньшие, чем для пробной площади №2. Чем меньше тангенс угла наклона кривых на рисунке 5, тем больше интервал возможных субоптимальных решений. Если тангенс угла наклона большой, то интервал возможных субоптимальных решений с погрешностью не большей, чем в предыдущем случае, сужается до одного оптимального значения числа связей с соседними деревьями.

Выводы:

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

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

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

максимальную устойчивость дерева к ветровому воздействию.

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

1. Захаров, Ю.В., Суховольский, В.Г. Модели устойчивости деревьев и насаждений к воздействию ветра // Лесоведение.- 2004. - № 2. - С. 61-67.

2. Калинин, М.И. Формирование корневой системы деревьев. М.: Лесная пром-сть, 1983. 152 с.

3. Суховольский, В.Г. Экономика живого: Оптимизационный подход к описанию процессов в экологических сообществах и системах. - Новосибирск: Наука, 2004. - 140 с.