Лесоведение, 2022, № 3, стр. 321-334

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

В. О. Корниенко a*, В. Н. Калаев b

a Донецкий национальный университет,
83001 Донецк, ул. Университетская, 24, Украина

b Воронежский государственный университет
394018 Воронеж, Университетская пл., 1, Россия

* E-mail: kornienkovo@mail.ru

Поступила в редакцию 13.07.2020
После доработки 07.12.2020
Принята к публикации 01.12.2021

Полный текст (PDF)

Аннотация

В работе установлена связь между температурным фактором и биомеханическим параметром (модулем упругости) древесных тканей березы повислой (Betula pendula Roth). При оттаивании древесины модуль упругости в среднем снижается в 2–2.5 раза. Подобное снижение происходит неравномерно, ступенчато. Наибольшее влияние на механическую устойчивость дерева имеет скорость изменения модуля упругости от скорости оттаивания образца/ствола древесного растения. При оттаивании показатели параметров механической устойчивости снижаются в среднем на 45%. Такие изменения отражаются на углах отклонения ствола от вертикали, углах отхождения его скелетных ветвей, сопротивлении изгибу и устойчивости всего дерева к ветровым и гравитационным нагрузкам. В условиях низкого антропогенного загрязнения показатель жесткости на изгиб для растений в возрасте 40–45 лет на 22 ± 2% выше, чем в условиях повышенной антропогенной нагрузки города. В условиях антропогенной нагрузки под влиянием ветровых нагрузок, действия температуры и других погодных факторов, у березы повислой происходит трансформация архитектоники кроны, которая в условиях снежной бури и ледяного дождя приводит к необратимой деформации или обломам стволов. В результате на экспериментальных участках в г. Донецке, в период с 2014 по 2020 гг., выпало 63 дерева, а 168 растений подверглись необратимой деформации и имели высокую аварийность. Обломы были зафиксированы у растений на высоте 2–4 м (35%), 5–6 м (17%), 7–8 м (52%) и 9 м (6%). В процессе корреляционного анализа выявлена сильная положительная связь между морфометрическим коэффициентом d : l (отношение диаметра к длине ствола) и механической устойчивостью (R = 0.87), а также аварийностью (R = 0.79) деревьев березы повислой. В связи с этим представляется возможным использование коэффициента d : l как морфометрического маркера механической устойчивости деревьев березы повислой в условиях юга Восточно-Европейской равнины (Донецкий кряж).

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

Экологическая адаптация древесных растений в городской среде к действию природно-климатических факторов является одной из важнейших задач функциональной экологии (Fournier et al., 2013; Dahle et al., 2017). Активно проводятся исследования по оценке механической устойчивости и аварийности деревьев под влиянием температуры окружающей среды (Green et al., 1999; Szmutku et al., 2011; Корниенко, Нецветов, 2013; Корниенко и др., 2018), ветровых нагрузок (James et al., 2006; Dahle, Grabosky, 2010; James et al., 2014; Dahle et al., 2017; Jelonek et al., 2019), снеговых и ледяных бурь (Nock et al., 2016). Многие авторы также отмечают важность изучения влияния снежных и ледяных бурь, а также ветровых нагрузок на структуру и функционирование экосистем (Fahey et al., 2020; Klein et al., 2020), а также непосредственно оценку и частоту нарушений, вызванных такими природно-климатическими факторами (Curtis, Gough, 2018).

Влажность древесины оказывает влияние на механическую устойчивость растений, их модуль упругости и плотность древесины (Green et al., 1999; Szmutku et al., 2011; Корниенко, Нецветов, 2013). В большинстве работ исследования проведены на образцах с влажностью 12% (техническая древесина) (Zelinka et al., 2007; Wood handbook …, 1999, 2010; Virot et al., 2016). Лишь в отдельных работах влажность исследуемых образцов доводили до 145% (Mishiru, Asano, 1984a,b; Green et al., 1999, Shmytku et al., 2012; Spatz, Pfisterer, 2013; Nocetti et al., 2015), что соответствует особым условиям хранения и эксплуатации древесных материалов. Влажность живой древесины варьирует в пределах 25(30)–225%, в зависимости от вида (Niklas, Spatz, 2010). Сезонная динамика факторов окружающей среды и физиологических показателей растений (переход из фазы покоя в фазу вегетации) отражается на содержании влаги в растительных тканях. Амплитуда колебаний влажности древесины связана со скоростью потери/сорбции влаги, которая уменьшается с увеличением плотности тканей и содержания в них экстрактивных веществ (Sell, 1989). С влажностью древесины связаны также показатели, которые определяют устойчивость целого дерева или его частей к природно-климатическим факторам. Вследствие этого полученные для технической древесины значения некорректно применять в расчетах устойчивости живого дерева (Раздорский, 1955; Burgert et al., 2001). Некоторые работы (Virot et al., 2016), в которых применяют данные по сухой (технической) древесине для моделирования поведения древостоев к действию природно-климатических факторов, сталкиваются с серьезной критикой со стороны научного сообщества (Albrecht et al., 2016) и в итоге сводятся к идеализированной модели, которая не может быть применена в прогнозировании живых древесных насаждений.

Ко всему вышесказанному можно добавить, что большинство упомянутых научных работ нацелены на исследования лесных массивов (естественных насаждений), исследований по изучению механической устойчивости насаждений в городской среде (городских посадок) крайне мало. Так, М.В. Нецветов с соавторами (Нецветов и др., 2009; Корнієнко та ін., 2009; Нецветов, Суслова, 2009) в условиях юго-востока Украины изучали механическую устойчивость древесных растений к действию вибрационно-акустических нагрузок техногенной природы. В рамках вибрационной экологии эти исследования являются фундаментальными и основополагающими, однако физико-механические свойства тканей деревьев из коллекции Донецкого ботанического сада изучались только на территории дендрария (зона низкой антропогенной нагрузки). В условиях действия антропогенных факторов, значения биомеханических параметров и, как следствие, механическая устойчивость древесины могут отличаться (Dahle et al., 2017). Таким образом, представляется крайне важным изучение физико-механических свойств тканей древесных растений в условиях урбанизированной среды, а также определение их механической устойчивости и аварийности при сочетанном действии природно-климатических факторов.

Целью работы является оценка влияния природно-климатических факторов и антропогенного загрязнения на деревья березы повислой в условиях юга Восточно-Европейской равнины (Донецкий кряж) на примере г. Донецка, для чего были обобщены результаты исследования городских посадок березы повислой в 2014–2020 гг.

В задачи исследования входили: 1) исследование влияния температуры на модуль упругости древесных тканей березы повислой in vitro (лабораторные исследования); 2) оценивание механической устойчивости березы повислой, произрастающей в условиях техногенной нагрузки, а также рисков при температурных изменениях; 3) изучение влияния природно-климатических факторов на аллометрию березы повислой в городских посадках г. Донецка; 4) оценка влияния сезонных природно-климатических факторов на аварийность березы повислой.

ОБЪЕКТЫ И МЕТОДИКА

В качестве объекта исследования нами была выбрана береза повислая (Betula pendula Roth). Согласно литературным данным, вид обладает высоким адаптивным потенциалом в условиях промышленного загрязнения (Neverova et al., 2013), т.е. обладает комплексом реакций, обеспечивающих приспосабливаемость растений к экстремальным условиям среды. По значению проницаемости мембран как интегрального показателя функционального состояния растительных тканей в условиях антропогенной нагрузки городской среды береза повислая имеет низкие значения, что говорит о высокой устойчивости к атмосферным загрязнениям и лучшей работе систем регуляции и поддержания гомеостаза по сравнению с другими исследованными видами в условиях городской среды (Сарбаева и др., 2013). По уровню устойчивости к высоким температурам, действию “суховея” и водоудерживающей способности вид среднеустойчив (Михеева и др., 2011).

Представленность березы повислой в городских посадках центральной части современного Донецка составляет ~5% (Суслова и др., 2012; Глухов и др., 2016), а в целом по городу ~2% (Поляков, 2009) от всех видов древесных растений. Критический возраст в условиях города составляет 50 лет (Корниенко, Калаев, 2018). Растение быстрорастущее, морозо- и засухоустойчиво, нетребовательно к почвам, в условиях города Донецка используется как в одиночных, так и групповых, аллейных посадках. Вдоль автомагистралей березу повислую используют как в первом, так и во втором ряду.

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

В качестве места проведения исследований было выбрано четыре территории (участки №№ 1–4). Длина исследуемых насаждений березы повислой на участке № 1 составляет ~385 м, ширина – около 50 м (пересечение Ленинского пр. и ул. Одесской). Растения посажены рядами, всего 7 полных рядов. В каждом ряду с расстоянием в 2–3 м высажены 112 ± 3 берез. Возраст деревьев ~45 лет. На участке № 2 (парк “Кованых фигур”, пересечение основных автомагистралей – ул. Артема с ул. Университетской и пр-та Ватутина с пр-том Мира), растения произрастают вдоль пешеходных дорожек, высажены как солитеры. Возраст 7–20 лет. На участке № 3 (пр-т Ильича) деревья произрастают в линейных насаждениях вдоль автомагистрали. Возраст 40–50 лет. На участке № 4 (Дендрарий ГУ “Донецкий ботанический сад”, южный массив, вдоль Макеевского шоссе), возраст деревьев 40–50 лет. Растения произрастают одиночно, в куртинах либо небольшими группами (3–5 ед.).

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

Таблица 1.  

Значения параметров антропогенной нагрузки от автотранспорта и вибрационно-акустического шума

Номер участка Интенсивность движения, ед. ч–1 $\bar {x} \pm {{s}_{x}}$ Вибрационно-акустический шум, дБА
возле дороги (бордюр) 1 ряд деревьев окраина древостоя
$\bar {x} \pm {{s}_{x}}$ MAX $\bar {x} \pm {{s}_{x}}$ MAX $\bar {x} \pm {{s}_{x}}$ MAX
1 1060 ± 54 78 ± 5 92 74 ± 2 80 69 ± 1 73
2 1681 ± 211 79 ± 4 110 76 ± 2 85 62 ± 5 68
3 1334 ± 30 79 ± 1 86 73 ± 1 79 68 ± 1 73
4 45 ± 5*

* Естественный шум в насаждениях. Примечание. На рис. 1–4 $\bar {x} \pm {{s}_{x}}$ – среднее значение ± стандартное отклонение.

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

Таблица 2.  

Показатели токсических веществ в почвах и атмосферном воздухе на территории проведения исследований

Номер участка Загрязнение почвы, $\bar {x} \pm {{s}_{x}}$ Загрязнение атмосферного воздуха, $\bar {x} \pm {{s}_{x}}$
Zn, мг кг–1 Cu, мг кг–1 Cr, мг кг–1 CO, мг м–3 H2S, мг м–3 NH3, мг м–3
1 33.0 ± 4.4 42.5 ± 7.4 1.0 ± 0.1 4.5 ± 0.5
2 <0.01 202.3 ± 217.6 1.2 ± 0.5 4.3 ± 0.5 0.025 ± 0.005 0.11 ± 0.1
3 33.0 ± 3 42.5 ± 4.4 1.0 ± 0.1 4.3 ± 0.5 0.025 ± 0.005 0.11 ± 0.1
4 0.005 ± 0.001 0.2 ± 0.1 3.5 ± 0.5 0.006 ± 0.001
ПДК 23* 3* 6* 5** 0.008** 0.2**

 * Подвижная форма. ** Максимальная разовая концентрация. Примечание. Обозначения см. табл 1; “–” – не определилось.

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

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

Температурную зависимость модуля упругости проводили на побегах древесных растений, отобранных на исследуемых участках. Возраст отобранных побегов 3–5 лет, их срезали при отрицательных температурах в период покоя с нижней части кроны, не затеняемой в течение светового дня. Для контроля температуры образцов в них со стороны среза помещали термопару (модель TZ-A/BL алюминий) на глубину 1.5 см. Срезы покрывали герметичным материалом.

Модуль упругости древесных волокон определяли по величине изгиба цилиндра, обрезанной ветви, горизонтально защемленной в тисках, в ответ на приложение силы на ее свободном конце (методика детально описана в работах Корниенко, Нецветова (2014) и Корниенко и др. (2018)). Модуль упругости при температуре 15°С (температура воздуха в лаборатории) отличался, поэтому для сравнения его температурных зависимостей использовали значения, приведенные к MOЕ (modulus of elasticity=модуль упругости) при T = 15°С. Отбор образцов и исследования влияния температуры на модуль упругости древесины (MOE) проводились в феврале – начале марта каждого года (2014–2020 гг.).

Для оценки прочности и механической устойчивости березы повислой в условиях города использовали следующие параметры: Pcr и mcr — предельно допустимая нагрузка и масса, при действии которых ствол начинает деформироваться или обламывается; Hcr – критическая высота ствола, при достижении которой действие собственного веса привело бы к необратимой деформации или облому; RRB – относительное сопротивление изгибу; EI – сопротивление изгибу (методика детально описана в работе Корниенко, Калаева (2018).

Результаты визуального осмотра березы повислой фиксировались с помощью фотоаппарата Nikon Coolpix S2600, обработку и анализ изображений проводили в программе AxioVision Rel. 4.8. В результате выполнения работы для изучения архитектоники кроны, углов отхождения ствола и скелетных ветвей от нормали было обработано свыше 1500 электронных фотографий. В программе AxioVision Rel. 4.8. углы измерялись с помощью функций Measure → Angle с точностью до 1°. Среднее значение параметра не демонстрировало полной картины поведения древесного растения при изгибе. В связи с этим мы выделили три участка (А, Б, В), разделяющих ствол растения по основным линиям изгиба, с указанием места произрастания дерева (рис. 1).

Рис. 1.

Углы отклонения стволов березы повислой на экспериментальном участке № 1 (общая схема). А – 10.43°. Б – 55.85°. В – 113.17°.

В условиях антропогенной нагрузки аллометрические исследования стволов березы повислой проводили на экспериментальном участке № 1 в период 2016–2017 и 2019–2020 гг. Такие временные рамки связаны с погодными условиями, при которых получены основные результаты по исследованию архитектоники кроны и аварийности березы повислой. Так, в эти годы на территории проведения эксперимента отмечались следующие природно-климатические условия: после длительного потепления в течение зимнего сезона температура воздуха резко опускалась до значений –20°C, выпадали сильные осадки (мокрый снег, снег), метель (порывы ветра 15–20 м с–1), налипание мокрого снега, наледи на побегах/стволах.

Диаметр ствола измеряли мерной вилкой, высоту деревьев – с помощью электронного высотомера HEC Haglof (Швеция). Жизнеспособность исследованных деревьев определяли по 8-бальной шкале Савельевой (Савельева, 1975), где 8 баллов это здоровое растение, сухие ветви в кроне отсутствуют, ствол не имеет повреждений, а 0 баллов это полное усыхание всего дерева. Возрастное состояние дерева определяли по Чистяковой (Чистякова и др., 1989). В работе имеются обозначения возрастных классов: I – 5–15 лет; II – 16–25 лет; III – 26–55 лет; IV – 56–75 лет. Аварийность каждого дерева определяли по методике, предложенной в работе (Корниенко, Приходько, 2018). При оценке древесных насаждений нами были использованы следующие методы: визуальный, ретроспективный, инструментальный (а – взятие керна прирастным буром для определения доли гнили в древесине, %; б – биомеханические испытания “Pullingtest”). Плотность (объемный вес) древесных волокон определяли методом погружения и взвешивания.

Для статистической обработки данных использовали программы “Statistica 8" (StatSoft Inc.) и “Excel 2010” (Microsoft Corporation). При построении графиков выводили уравнения регрессии и оценивали величину достоверности аппроксимации (R2). Достоверность отличий средних значений полученных данных между коэффициентом d : l и состоянием дерева (norm, deform, crushed) определяли с использованием t-критерия Стьюдента. Для выявления достоверности между коэффициентом d : l и механической устойчивостью использовали коэффициент корреляции Пирсона (r), который считали отличным от нуля при p < 0.05.

РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ

Биомеханические исследования березы повислой. Плотность древесных тканей в условиях г. Донецка составляет 980 ± 7 кг м–3 вне зависимости от места произрастания.

Модуль упругости в условиях антропогенной нагрузки (экспериментальные участки №№ 1–3) города составляет 5.03 ± 0.77 ГН м–2, в условиях Донецкого ботанического сада (участок № 4 – дендрарий) – 4.30 ± 0.46 ГН м–2. Различия показателей биомеханических параметров достоверны при p < 0.05.

Влияние температуры на MOE древесных тканей березы повислой in vitro (лабораторные исследования). Зависимость MOE от температуры в общем виде для всех образцов имеет тенденцию к нелинейному падению при повышении температуры от 255 до 317 K (рис. 2).

Рис. 2.

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

При оттаивании изменение модуля упругости древесины у всех образцов происходит неравномерно, ступенчато. Число изломов кривой MOE (T) варьирует от 1 до 3, а температуры, при которых изменяется угол наклона MOE (T) – в пределах 2°–4°. Наибольшая крутизна наклона отмечается на участке от –16 до 0°C. Наши экспериментальные данные согласуются с данными, полученными в работах японских ученых A. Mishiro и I. Asano (1984а, б). Стоит отметить, что значения точек изгибов могут варьировать от видовых особенностей растений (например, анатомических и биохимических (Корниенко, Нецветов, 2014)). В наших экспериментах температура образцов изменялась со временем, как и изменялась скорость ее роста. Крутизна кривой модуля упругости от температуры снижается по мере падения скорости роста температуры образца. Видовые отличия в данном процессе не столь изучены и могут иметь место, например, при наличии ассимилянтов в живице, которое влияет на скорость льдообразования в сосудах.

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

Рис. 3.

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

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

Для березы повислой отношение EI к d : l имеет степенную зависимость (R2 = 0.83). При оттаивании EI уменьшается на ~45%, что при сочетанном действии с природно-климатическими факторами окружающей среды может вызвать потерю механической устойчивости растения и привести к необратимой деформации стволов (рис. 4а).

Рис. 4.

Зависимость параметров механической устойчивости от отношения d/l (а–в) и возрастных классов (г). 1 – значение показателя в зимний период при отрицательной температуре, 2 – значение показателя при оттаивании, 3 – значение показателя в летний период при положительной температуре, 4 – значения, полученные на территории ГУ “Донецкий ботанический сад” (контроль), 5 – значения, полученные в условиях антропогенной нагрузки г. Донецка.

Отношение RRB к d : l имеет линейную зависимость (R2 = 0.87), а Pcr к d : l степенную зависимость (R2 = 0.97). Параметры тесно связаны с показателем модуля упругости тканей растения и зависят от отношения диаметра ствола к его высоте. Наибольшему риску подвержены молодые растения (до 5–7 лет) и угнетенные вследствие роста в плотном древостое и борьбы за свет взрослые растения (имеют тонкие и высокие стволы d : l < 0.01, рис. 4б).

Предельно допустимая нагрузка дает точные значения критической массы при действии которой произойдут необратимые изменения/облом ствола (рис. 4в). Отношение d : l можно использовать также следуя зависимостям на рис. 3в, как морфометрический маркер устойчивости древесного растения. При значении d : l в диапазоне от 0.004 до 0.015 критическая масса варьирует в пределах от 11 до 815 кг, при этом ствол находится в замороженном состоянии, а с учетом температурных изменений (рис. 2 и 3) механическая устойчивость дерева снижается на ~45%, значения Pcr лежат в диапазоне от 6.2 до 458.9 кг.

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

Полученные температурные зависимости MOE и параметров механической устойчивости объясняют эффекты деформации и обломов стволов березы повислой в зимний период.

Исследование архитектоники кроны деревьев березы повислой, произрастающих в условиях антропогенного загрязнения г. Донецка (2017–2020 гг.). В условиях низкой антропогенной нагрузки на территории дендрария ГУ “Донецкий ботанический сад” деревья березы повислой обладают высоким баллом жизнеспособности (в хорошем состоянии – 5–7 баллов по Савельевой), они не подвергались крушению и необратимым изгибам при ветровалах (2017–2020 г.), снеговых и ледяных штормах 2017 и 2020 г. Это объясняется зависимостью log EI от возраста (рис. 4г). В условиях относительного контроля показатель жесткости на изгиб на 22 ± ± 2% выше, чем в условиях антропогенной нагрузки города (p < 0.05). При этом углы отхождения от вертикали составляют не более 5°, а отношение d : l > 0.02.

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

Таблица 3.  

Влияние места произрастания в древостое на угол отхождения ствола дерева березы повислой (2017 г.)

Место произрастания Число,
шт.
Угол наклона ствола от вертикали, °
$\bar {x} \pm {{s}_{x}}$ доля от общего числа, % MAX доля от общего числа, %
Окраина древостоя со стороны жилого дома 91 16 ± 5 60 43 40
Окраина древостоя со стороны трассы 220 15 ± 6 40 26 60
Середина древостоя 352 4 ± 2 90 12 ± 3 10

Примечания. Обозначения см. табл. 1.

В результате обследований в 2020 г. выявлено, что угол наклона стволов растений изменился и трансформировался (табл. 4). Среднее значение параметра уже не отражало и не описывало поведения древесного растения при изгибе. В связи с этим мы выделили три участка, разделяющих ствол растения по основным линиям изгиба: участок “А” от основания до первого изгиба ствола (~3 м), участок “Б” – середина ствола растения, третий участок “В” с линией изгиба до вершины кроны.

Таблица 4.  

Влияние места произрастания в древостое на угол отхождения ствола дерева березы повислой (2020 г.)

Место произрастания Угол наклона ствола от вертикали, град $\bar {x} \pm {{s}_{x}}$
участок А участок Б участок В
Окраина древостоя со стороны жилого дома 18 ± 3 57 ± 7 118 ± 12
Окраина древостоя со стороны трассы 14 ± 3 49 ± 8 96 ± 21
Середина древостоя 7 ± 5 36 ± 26 80 ± 44

Примечания. Обозначения см. табл. 1.

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

Влияние сезонных природно-климатических факторов на аварийность деревьев березы повислой, произрастающих в условиях антропогенной нагрузки. В результате ветровала 2017 г. в первом ряду от жилого дома, там, где углы отклонения от вертикали имели значения 40°–45°, выпало несколько растений. Исследовав их аллометрические параметры, экспериментально установили, что коэффициент d : l у таких деревьев ≤0.01.

В 2020 г. в условиях снежной бури и ледяного дождя выпало 63 дерева. Обломы были зафиксированы у растений на высоте 2–4 м (35%), 5–6 м (17%), 7–8 м (52%) и 9 м (6%). В основном это растения с относительно тонкими стволами, d : l у которых составляет 0.01. Для таких растений были получены зависимости сопротивления изгибу от d : l с учетом их состояния: нормально растущие (norm, имеющие вертикальное направление 90° ± 5°), деформированные (deform) или подвергшиеся облому ствола (crushed) в результате действия климатических факторов (снежных бурь и ледяного дождя) (рис. 5). Отношение EI к d : l для всех групп деревьев имеет степенную зависимость и высокий коэффициент детерминации R2 ~ 0.89–0.99. Для вертикально растущих деревьев отношение d : l имело диапазон 0.03–0.05.

Рис. 5.

Зависимость сопротивления изгибу ствола дерева от отношения d : l. 1 – стволы имеют строго вертикальное направление роста, 2 – деформированные, изогнутые стволы, 3 – подверглись облому ствола.

Деревья, произрастающие в условиях антропогенной нагрузки, обладают меньшей устойчивостью, что связано, во-первых, со значениями дендрометрических параметров стволов растений (рис. 6), во-вторых, с увеличением высоты расположения биомассы дерева, и как следствие – приложением нагрузок в виде силы ветра, а также дополнительной массы при выпадении осадков в верхней трети от высоты дерева. Поэтому нами и были зафиксированы обломы стволов в основном двух вариантов: первый – на высоте 2–4 м (35%), т.е. у основания, когда при прогреве ствола образуется неоднородная структура и на высоте – 2–4 м создается концентрация напряжения механически нагруженного ствола; и второй на высоте 7–9 м (58% от общего числа обрушившихся растений), т.е. в верхней трети, что является классической экологической стратегией выживания дерева в условиях действия критических статических природно-климатических нагрузок (налипание снега, оледенение).

Рис. 6.

Зависимость диаметра ствола дерева от d : l как морфометрического маркера механической устойчивости. 1 – стволы имеют строго вертикальное направление роста, 2 – деформированные, изогнутые стволы, 3– подверглись облому ствола.

Использование коэффициента d : l как морфометрического маркера механической устойчивости деревьев березы повислой на юге Восточно-Европейской равнины (Донецкий кряж). Актуальность исследований, посвященных изучению связи между морфометрическими параметрами ствола d и h(l) и механической устойчивостью деревьев к ветровым и гравитационным нагрузкам отмечал еще Г.Б. Кофман в своей работе “Рост и форма деревьев” (1986). Он с целью нахождения механически устойчивого состояния дерева применял соотношение h(l) : d ~ 102 (т.е. отношение высота : длина ствола дерева к его диаметру у основания). Теоретически было установлено, что при h > 110r (h – высота, r – радиус у основания) деревьям угрожает потеря упругой механической устойчивости. Однако автор при разработке математической модели, описывающей устойчивость живых деревьев, использовал физико-механические параметры (E – модуль упругости и ρ плотность древесных тканей) технической (сухой) древесины хвойных пород, чего, по мнению других исследователей (Burgert et al., 2001), делать нельзя. Так, сухая древесина, обладающая высокой жесткостью и хрупкостью, существенно отличается от зеленой, живой древесины, которая является более упругой (имеется в виду значение модуля упругости, которое коррелирует со значением влажности (w)), и менее хрупкой (Burgert et al., 2001). Также значимым фактором механической устойчивости к действию природно-климатических факторов является индивидуальная адаптация древесных видов к условиям среды произрастания (Thomas, 2011). Поэтому расчеты, выполненные на растениях одного вида, без учета условий произрастания и особенностей физико-механических свойств древесины живых растений, некорректно использовать в качестве универсальной модели механического поведения дерева при нагрузках, а также вывода морфометрического коэффициента устойчивости/аварийности (Niklas, Spatz, 2012; Albrecht et al., 2016). В связи с этим актуальной задачей является нахождение для каждого вида древесных растений связи состояния (аварийности) и механической устойчивости с коэффициентом d : l, а также диапазона вариации этого соотношения с учетом условий произрастания.

В нашем исследовании при статистической обработке данных методом корреляционного анализа была выявлена сильная положительная связь между морфометрическим коэффициентом d : l и механической устойчивостью (r = 0.87), а также аварийностью (r = 0.79) деревьев березы повислой (табл. 5 ).

Таблица 5.  

Связь состояния дерева (аварийность) и его механической устойчивости с коэффициентом d : l

Состояние дерева Отношение диаметра к длине ствола, d : l Жесткость на изгиб, Н/м–2
Norm 0.28–0.43 1.59 × 106–8.93 × 106
Deform 0.08–0.015 1.01 × 104–10.49 × 104
Crushed 0.05–0.017 1.01 × 104–8.19 × 105

Примечание. Norm – стволы имеют строго вертикальное направление роста, Deform – деформированные, изогнутые, Crushed – подверглись облому ствола.

Отношение d : l для деревьев березы повислой, произрастающих в условиях антропогенной нагрузки г. Донецка, в среднем составляло 0.01 ± 0.003. При таких значениях наблюдались необратимые деформации или обломы стволов/скелетных ветвей. Это может быть объяснено тем, что в условиях антропогенной нагрузки происходит нарушение физиолого-биохимических, клеточных, молекулярных и генетических процессов в растении (Поляков, 2009; Alonso-Serra et al., 2020). Внешним результатом действия загрязнения являются изменения морфометрических параметров стволов и кроны древесных растений (Поляков, 2009). Считается, что морфометрические параметры кроны и ствола древесного растения оказывают больший вклад в вариацию механической устойчивости и аварийности деревьев (Sellier, Fourcaud, 2009). В связи с этим представляется возможным использование коэффициента d : l как морфометрического маркера механической устойчивости деревьев березы повислой в условиях юга Восточно-Европейской равнины (Донецкий кряж). Возможно, для других видов растений связь диапазона варьирования d : l в пределах 0.01 и состояния дерева будет также прослеживаться, этот показатель может быть универсальным, однако это требует дальнейших исследований.

ЗАКЛЮЧЕНИЕ

Установлено влияние температуры на модуль упругости древесных тканей березы повислой. Показатель достоверно снижается при оттаивании в среднем в 2–2.5 раза, снижение MOE происходит неравномерно, ступенчато. Наибольшее влияние на механическую устойчивость дерева имеет скорость изменения модуля упругости (dMOE:dt) от скорости оттаивания образца/ствола древесного растения (dT:dt).

При оттаивании показатели параметров механической устойчивости снижаются в среднем на 45%. Такие изменения отражаются на углах отклонения ствола от вертикали, углах отхождения его скелетных ветвей, сопротивлении изгибу и устойчивости всего дерева к ветровым и гравитационным нагрузкам. В условиях низкого антропогенного загрязнения показатель механической устойчивости logEI для растений в возрасте 40–45 лет на 22 ± 2% выше, чем в условиях повышенной антропогенной нагрузки города (p < 0.01). В условиях антропогенной нагрузки под влиянием ветровых нагрузок, температуры, снежных и ледяных бурь происходит трансформация архитектоники кроны березы повислой. Такие изменения приводят к необратимой деформации или обломам стволов, в условиях снежной бури и ледяного дождя. В результате 168 растений подверглись необратимой деформации и имели высокую аварийность, выпало 63 дерева. Обломы были зафиксированы у растений на высоте 2–4 м (35%), 5–6 м (17%), 7–8 м (52%) и 9 м (6%).

В процессе корреляционного анализа была выявлена сильная положительная связь между морфометрическим коэффициентом d : l и механической устойчивостью (r = 0.87), а также аварийностью (r = 0.79) деревьев березы повислой. В связи с этим представляется возможным использование коэффициента d : l как морфометрического маркера механической устойчивости деревьев березы повислой в условиях юга Восточно-Европейской равнины (Донецкий кряж).

Список литературы

  1. Глухов А.З., Хархота Л.В., Пастернак Г.А., Лихацкая Е.Н. Современное состояние дендрофлоры города Донецка // Самарский научный вестник. 2016. Т. 2(15). С. 20–24.

  2. Корниенко В.О., Калаев В.Н. Механическая устойчивость древесных пород и рекомендации по предотвращению их аварийности в городских насаждениях. Воронеж: Роза Ветров, 2018. 92 с.

  3. Корниенко В.О., Калаев В.Н., Елизаров А.О. Влияние температуры на биомеханические свойства древесных растений в условиях закрытого и открытого грунта // Сибирский лесной журн. 2018. № 6. С. 91–102.

  4. Корниенко В.О., Нецветов М.В. Влияние отрицательных температур на механическую устойчивость дуба красного (Quercus rubra L.). // Промышленная ботаника. 2013. Вып. 13. С. 180–186.

  5. Корниенко В.О., Нецветов М.В. Криоскопия влаги и температурная зависимость модуля упругости древесины // Вісті Біосферного заповідника “Асканія-Нова”. 2014. Т. 16. С. 88–94.

  6. Корнієнко В.О., Нецветов М.В., Нікуліна В.М., Суслова О.П. Дослідження стійкості дерев до вібрацій // Вісник Львівського університету. Серія фізична. 2009. Вип. 44. С. 185–193.

  7. Кофман Г.Б. Рост и форма деревьев. Новосибирск: Наука, 1986. 211 с.

  8. Михеева М.А., Федорова А.И. Влияние высоких температур на устойчивость древесных растений в городской среде // Вестник Воронежского государственного университета. Серия: География. Геоэкология. 2011. № 2. С. 166–175.

  9. Нецветов М.В., Суслова Е.П., Никулина В.М., Корниенко В.О. Механическая устойчивость деревьев к антропогенным вибрационным нагрузкам // Фізичні методи в екології, біології та медицині. Програма і збірник тез ІІ міжнародної конференції. Львів–Ворохта, Україна, 2-6 вересня 2009 р. Львів: Видавничий центр ЛНУ ім. І. Франка, 2009. С. 22–23.

  10. Нецветов М.В., Суслова Е.П. Механическая устойчивость деревьев и кустарников к вибрационным нагрузкам // Промышленная ботаника. 2009. № 9. С. 60–67.

  11. Поляков А.К. Интродукция древесных растений в условиях техногенной среды. Донецк: Ноулидж, 2009. 268 с.

  12. Раздорский В.Ф. Архитектоника растений. М.: Советская наука, 1955. 432 с.

  13. Савельева Л.С. Устойчивость деревьев и кустарников в защитных лесных насаждениях. М.: Лесн. пром-сть, 1975. 168 с.

  14. Сарбаева Е.В., Воскресенская О.Л., Воскресенский В.С. Оценка устойчивости древесно-кустарниковых растений в урбанизированной среде // Современные проблемы науки и образования. 2013. № 2. URL: http://science-education.ru/ru/article/view?id=9011 (дата обращения: 27.04.2020).

  15. Суслова О.П., Поляков О.К., Нецвєтов М.В., Дацько О.М., Лихацька О.М. Життєздатність деревних рослин у міських вуличних насадженнях на південному сході України // Промышленная ботаника. 2012. Вып. 12. С. 12–18.

  16. Чистякова А.А., Заугольнова Л.Б., Полтинкина И.В., Кутьина И.С., Лшцинский Н.Н. Диагнозы и ключи возрастных состояний лесных растений. Деревья и кустарники: методические разработки для студентов биологических специальностей. М.: Прометей, 1989. 106 с.

  17. Albrecht A., Badel E., Bonnesoeur V., Brunet Y., Constant T., Defossez P., Langre E., Dupont S., Fournier M., Gardiner B., Mitchell S. J., Moore J. R., Moulia B., Nicoll B. C., Niklas K. J., Schelhaas M., Spatz H.-Ch., Telewski F. W. Comment on “Critical wind speed at which trees break” // Physical review. 2016. E 94. P. 067001-1–067001-2.

  18. Alonso-Serra J., Shi X., Peaucelle A., Rastas P., Bourdon M., Immanen J., Takahashi J., Koivula H., Eswaran G., Muranen S., Help H., Smolander O., Su Ch., Safronov O., Gerber L., Salojarvi J., Hagqvist R., Mähönen A., Nieminen K., Helariutta Y. Elimäki locus is required for vertical proprioceptive response in birch trees // Current Biology. 2020. V. 30. Is. 4. P. 589–599.

  19. Burgert I., Bernasconi A., Niklas K.J., Eckstein D. The influence of rays on the transverse elastic anisotropy in green wood of deciduous trees // Holzforschung. 2001. V. 55. Is. 5. P. 449-454.

  20. Curtis P.S., Gough C.M. Forest aging, disturbance and the carbon cycle // New Phytol. 2018. V. 219. P. 1188–1193.

  21. Dahle G.A., James K.R., Kane B., Grabosky J.C., Detter A. A review of factors that affect the static loadbearing capacity of urban trees // Arboriculture & Urban Forestry. 2017. V. 43. Is. 3. P. 89–106.

  22. Dahle G.A., Grabosky J.C. Variation in modulus of elasticity (E) along Acer platanoides L. (Aceraceae) branches // Urban Forestry & Urban Greening. 2010. V. 9. P. 227–233.

  23. Fahey R.T., Atkins J.W., Campbell J.L., Rustad L.E., Duffy M., Driscoll Ch.T., Fahey T.J., Schaberg P. G. Effects of an experimental ice storm on forest canopy structure // Canadian J. Forest Research. 2020. V. 50. Is. 2. P. 136–145.

  24. Fournier M., Dlouhá J., Jaouen G., Almeras T. Integrative biomechanics for tree ecology: beyond wood density and strength // J. Experimental Botany. 2013. V. 64. № 15. P. 4793–4815.

  25. Green D.W., Evans J.W., Logan J.D., Nelson W.J. Adjusting modulus of elasticity of lumber forchanges in temperature // Forest Products J. 1999. V. 49(10). P. 82–94.

  26. James K.R., Dahle G.A., Grabosky J., Kane B., Detter A. Tree biomechanics literature review: dynamics // Arboriculture & Urban Forestry. 2014. V. 40. Is. 1. P. 1–15.

  27. James K.R., Haritos N., Ades P.K. Mechanical stability of trees under dynamic loads // American Journal of Botany. 2006. V. 93. Is. 10. P. 1522–1530.

  28. Jelonek T., Tomczak A., Karaszewski Z., Jakubowski M., Arasimowicz-Jelonek M., Grzywiński W., Kopaczyk J., Klimek K. The biomechanical formation of trees // Drewno 2019. V. 62. № 204. P. 5–22.

  29. Klein R.W., Koeser A.K., Kane B., Landry S.M., Shields H., Lloyd S., Hansen G. Evaluating the likelihood of tree failure in Naples, Florida (United States) following hurricane Irma // Forests. 2020. V. 11. Is. 485. P. 1–10.

  30. Mishiro A., Asano I. Mechanical properties of wood at low temperatures: effect of moisture content and temperature on bending properties of wood. Part I. Moisture content below the fiber saturation point // J. Japan Wood Research Society. 1984a. V. 30(3). P. 207–213.

  31. Mishiro A., Asano I. Mechanical properties of wood at low temperatures: effect of moisture content and temperature on bending properties of wood. Part II. Moisture content beyond the fiber saturation point // Ibid. 1984b. V. 30(4). P. 277–286.

  32. Neverova O.A., Legoshchina O.M., Bykov A.A. Anatomy of leaves of Betula pendula (Roth.) affected by air emissions in industrial area of Kemerovo city // Middle East J. Scientific Research. 2013. V. 17. Is. 3. P. 354–358.

  33. Niklas K.J., Spatz H.-C. Worldwide correlations of mechanical properties and green wood density // American J. Botany. 2010. V. 97. Is. 10. P. 1587–1594.

  34. Niklas K.J., Spatz H.-Ch. Plant Physics. Chicago: University of Chicago Press, 2012. 426 p.

  35. Nocetti M., Brunetti M., Bacher M. Effect of moisture content on the flexural properties and dynamic modulus of elasticity of dimension chestnut timber // European J. Wood and Wood Products. 2015. V. 73. P. 51–60.

  36. Nock C.A., Lecigne B., Taugourdeau O., Greene D.F., Dauzat J., Delagrange S., Messier Ch. Linking ice accretion and crown structure: towards a model of the effect of freezing rain on tree canopies // Annals of Botany. 2016. V. 117. Is. 7. P. 1163–1173.

  37. Sell J. Eigenschaften und Kenngrößen von Holzarten // Lignum, Baufachverlag AG Zürich. 1989. 87 p.

  38. Sellier D., Fourcaud T. Crown structure and wood properties: influence on tree sway and response to high winds // American J. Botany. 2009. V. 96(5). P. 885–896. https://doi.org/10.3732/ajb.0800226

  39. Spatz H.-C., Pfisterer J. Mechanical properties of green wood and tree risk assessment // Arboriculture & Urban Forestry. 2013. V. 39(5). P. 218–225.

  40. Szmutku M.B., Campean M., Laurenzi W. Influence of cyclic freezing and thawing upon spruce wood properties // Pro Ligno. 2012. V. 8(1). P. 35–43.

  41. Szmutku M.B., Campean M., Sandu A.V. Microstructure modifications induced in spruce wood by freezing // Ibid. 2011. V. 7. Is. 4. P. 26–31.

  42. Thomas S.C. Size- and age-related changes in tree structure and function // Tree Physiology. 2011. Ch. 2. P. 33–64.

  43. Virot E., Ponomarenko A., Dehandschoewercker E., Quere D., Clanet C. Critical wind speed at which trees break // Physical Review. 2016. E 93. P. 023001-1–023001-7.

  44. Wood handbookWood as an engineering material. Madison: WI, Forest Products Laboratory, 2010. 508 p.

  45. Wood handbookWood as an engineering material. Madison: WI, U.S. Department of Agriculture, Forest Service, Forest Products Laboratory, 1999. 463 p.

  46. Zelinka S.L., Stone D.S., Rammer D.R. Equivalent circuit modeling of wood at 12% moisture content // Wood Fiber Science: Journal of the Society of Wood Science and Technology. 2007. V. 39(4). P. 556–565.

Дополнительные материалы отсутствуют.