Каждое индивидуально растущее дерево может стать индикатором той окружающей среды, которая вокруг него располагается. Поэтому уровень загрязненности придорожной территории можно оценить по поведению, например, березы повислой, растущей около автомобильной дороги. Причем учетные листья березы реагируют не только на загрязнения от выхлопов автомобилей, но при своем росте и развитии за вегетационный период учитывают комплекс факторов воздуха, воды и почвы.
Цель статьи: выявление закономерностей влияния координат березы повислой, расположенной у кромки автомобильной дороги, на продолжительность вегетационного периода и ширину учетных листьев около дороги в городской среде.
В экологических технологиях с использованием листьев березы повислой приходит понимание о моделировании взаимных связей между параметрами структуры листьев растений методом идентификации [1, 2]. Априори ясно, что именно погода влияет на ход развития и роста (онтогенеза) растений. А на многолетние растения погода влияет через ежегодный онтогенез листвы. Кванты поведения листьев, например, березы повислой [3], распространенной на Северном полушарии, четко зависят от квантов (асимметричных вейвлетов) поведения температуры воздуха и относительной влажности. Метеорологические условия являются сильными факторами активности биологических объектов, и для этого в статье [4] оценены влияние температуры, осадков, атмосферного давления и влажности.
Рост растений – сложный процесс, в его основе лежат такие фундаментальные явления, как ритмичность, полярность, дифференциация, раздражимость, корреляция. Эти процессы являются общими для онтогенеза живых организмов [5, 6]. Онтогенез – индивидуальное развитие организма от зиготы (или вегетативного зачатка) до природной смерти. Благодаря активной деятельности меристем и фотосинтетической активности листьев зеленое растение приобретает ряд черт, которые характеризуют его рост. В процессе онтогенеза растения рост наблюдается на протяжении основных этапов его жизненного цикла [7].
Формирование и отмирание листьев в цикле онтогенеза разделяются на такие этапы: распускание почек, рост и развитие листьев, расцвечивание отмирающих листьев, опадение листьев. Мы предлагаем два этапа онтогенеза – роста до максимума и спада до опадения.
В Берлине [8] по 252 дереву липы на кернах по расстояниям от центра к периферии были выявлены изменения приращения толщины деревьев за 50–100 лет. Семь пород деревьев в Европе [9] показали нелинейные зависимости развития растений от температуры и микроклимата. В биометео-
рологии пока наблюдается неопределенность на будущее. В [10] показано влияние городской растительности на тепловой комфорт человека, равный 22,2 °С. Рациональный фитоклимат позволит противодействовать глобальному потеплению.
Вегетационный период становится одним из важных экосистемных процессов, так как развитие листьев очень чувствительно к температуре воздуха. Поэтому будущее климата в наблюдениях за листьями. Повышение влажности воздуха снижает температуру и накопление биомассы в молодых березках, в особенности восприимчивы листья [11].
Цель статьи – повышение точности индикации качества окружающей листья березы локальной среды сбоку на высоте 1,5–2,0 м со стороны преобладающих ветров по времени вегетации в 2018 г. от первого мая до максимума ширины как среднего из 10 учетных листьев по патенту на изобретение 2606189 в зависимости от расстояния от края дороги, высоты расположения центра зоны 10 учетных листьев над почвой и приведенных географических координат расположения центра комля деревьев.
Материалы и методы исследования
Колебания записываются волновой формулой [1, 2] вида
(1)
где y – показатель (зависимый фактор), i – номер составляющей модели (1), m – количество членов в модели (1), х – объясняющая переменная (влияющий фактор), a1...a8 – параметры модели (1), принимающие числовые значения в ходе структурно-параметрической идентификации в програм-
мной среде CurveExpert-1.40 (URL: http://www.curveexpert.net/), Ai – амплитуда (половина) вейвлета (ось y), pi – полупериод колебания (ось x).
В табл. 1 даны расстояния от дороги и почвы, а также географические и биометрические параметры 10 берез в г. Йошкар-Ола по средней ширине у 10 учетных листьев.
На каждой березе сбоку кроны со стороны преобладающих ветров (примерно северо-запад) выделялась зона диаметром примерно 0,5 м, где выделяли 10 учетных листьев, к корешкам их привязывали белые нитки с бирками, на которых были указаны номера учетных листьев. Приведенная система координат вычислялась следующим образом. С помощью сотового телефона вначале записывали координаты зоны учетных листьев (в градусах, минутах и секундах). Затем эти значения приводили к градусам по формуле ПК = градус + минута / 60 + секунда /3600. Из полученных результатов вычитали: для северной широты начало координат –56,60; для восточной долготы начало координат –47,80. Однако для программной среды Curve Expert 1.40 необходимо умножить полученные значения приведенных координат на 100.
Таблица 1
Параметры расположения берез с 10 учетными листами на каждом дереве
Место учетной березы |
Расстояния L и h, м |
Координаты березы |
Приведенные координаты 10-2, 0 |
Время вегетации ?, сутки |
Макс. ширина листа, , мм |
|||
дороги |
почвы |
с.ш. |
в.д. |
с.ш. ? |
в.д. ? |
|||
1. Ул. Лебедева |
7,0 |
1,60 |
56? 37? 6? |
47? 56? 48? |
1,833 |
14,667 |
110 |
45,92 |
2. Воскресенский пр. |
2,9 |
1,46 |
56?38?6? |
47?54?51? |
3,500 |
11,417 |
110 |
48,66 |
3. Ул. Эшкинина |
5,2 |
1,40 |
56? 37? 58? |
47? 54? 47? |
3,278 |
11,306 |
92 |
39,77 |
4. Пр-т Ленина |
4,9 |
1,60 |
56? 38? 8? |
47? 53? 12? |
3,556 |
8,667 |
99 |
40,11 |
5. Б-р Чавайна |
2,1 |
1,60 |
56? 37? 52? |
47? 54? 47? |
3,111 |
11,306 |
106 |
44,71 |
6. Ул. Панфилова |
4,1 |
1,30 |
56? 37? 40? |
47? 52? 55? |
2,778 |
8,194 |
103 |
39,56 |
7. Ул. Карла Маркса |
3,2 |
1,60 |
56? 36? 47? |
47? 52? 45? |
1,306 |
7,917 |
99 |
44,39 |
8. Ул. Строителей |
4,3 |
1,45 |
56? 36? 54? |
47? 52? 5? |
1,500 |
6,806 |
103 |
42,25 |
9. Ботанический сад |
2,4 |
1,55 |
56? 38? 42? |
47? 52? 50? |
4,500 |
8,056 |
103 |
40,71 |
10. Ул. Осипенко |
2,3 |
1,45 |
56? 38? 42? |
47? 52? 50? |
4,500 |
8,056 |
99 |
44,09 |
Время вегетации до максимального роста |
Средняя максимальная ширина листа |
Рис. 1. Графики рангового распределения биометрических параметров листьев березы повислой
Результаты исследования
и их обсуждение
В табл. 2 даны результаты факторного анализа по данным табл. 1. Здесь факторы L, h, ? и ? являются только ка влияющие переменные, так как они не могут одновременно стать зависимыми показателями. А факторы ? и становятся одновременно влияющими и зависимыми параметрами, поэтому их можно ранжировать по вектору предпорядка предпочтительности «лучше > хуже». Для обоих биометрических параметров ? и принимается вектор «чем больше, тем лучше».
Тогда ранговое распределение будет подчиняться закону экспоненциальной гибели при рангах R = 0, 1, 2, 3,... (рис. 1) по формулам:
– для рангового распределения времени вегетации учетных листьев у 10 берез (0,9912)
; (2)
– для средней ширины у 10 учетных листьев для каждой из 10 берез (0.9852)
(3)
В табл. 2 приведены результаты факторного анализа с использованием тренда, как частного случая формулы (1) при бесконечном периоде колебания.
Коэффициент коррелятивной вариации, как мера адекватности всей системы параметров объекта исследования, определяется как 6,4060 / (6?2) = 0,5338. Первое место среди влияющих переменных занимает время вегетации до конца этапа роста, а второе – средняя ширина учетных листьев при максимуме роста учетных листьев. Среди зависимых показателей как оценочных критериев занимает средняя ширина учетных листьев в конце этапа роста.
Выберем корреляционную матрицу (табл. 3) со средними и сильными бинарными отношениями с коэффициентом корреляции не менее 0,5. Наибольший коэффициент корреляции 0,7574 имеет влияние расстояния от примерного центра зоны расположения 10 учетных листьев на боковой поверхности кроны березы повислой со стороны преобладающих ветров (северо-запад) до края автомобильной дороги.
В табл. 4 расположим модели трендов по убыванию адекватности.
На рис. 2 приведены графики влияния расстояния от дороги до учетных листьев.
Из точечного графика остатков видно, что возможно колебание (табл. 5, рис. 3). Из-за малого количества точек коэффициент корреляции трехчленной модели стал равным 1.0000.
Средняя максимальная ширина учетных листьев березы повислой может стать великолепным индикатором загрязнения придорожной территории. Для идентификации волновых закономерностей необходимо в разных местах принять более 20 берез.
Наибольшая турбулентность наблюдается до 3 м от края дороги. Здесь же по тренду наблюдается максимум средней ширины. А затем от 3 м до 5 м наблюдается снижение ширины листа, то есть в этом интервале расстояния от дороги наблюдается сильное подавление вегетации листьев. После 5 м листья березы начинают расти с возрастающей шириной листьев.
Остальные закономерности из табл. 4 показаны на рис. 4.
Таблица 2
Корреляционная матрица факторного анализа и рейтинг факторов
после идентификации по закономерности тренда
Влияющие факторы (параметры x) |
Зависимые факторы (показатели y) |
Сумма ?r |
Место Ix |
|
?, сутки |
, мм |
|||
Расстояние группы листьев от дороги L, м |
0,0021 |
0,7574 |
0,7595 |
4 |
Высота группы листьев от почвы h, м |
0,2360 |
0,4640 |
0,7000 |
5 |
Приведенная северная широта ?, 0 |
0,1582 |
0,1656 |
0,3238 |
6 |
Приведенная восточная долгота ?, 0 |
0,5326 |
0,7540 |
1,2866 |
3 |
Время вегетации учетных листьев ?, сутки |
0,9912 |
0,6804 |
1,6716 |
1 |
Средняя ширина учетных листьев , мм |
0,6793 |
0,9852 |
1,6645 |
2 |
Сумма коэффиц. корреляции ?r |
2,5994 |
3,8066 |
6,4060 |
– |
Место Iy |
2 |
1 |
– |
0,5338 |
Таблица 3
Корреляционная матрица с коэффициентом корреляции не менее 0,5
Влияющие факторы (параметры x) |
Показатели y |
|
?, сутки |
, мм |
|
Расстояние группы листьев от дороги L, м |
0,7574 |
|
Приведенная восточная долгота ?, 0 |
0,5326 |
0,7540 |
Время вегетации учетных листьев ?, сутки |
0,6804 |
|
Средняя ширина учетных листьев , мм |
0,6793 |
Таблица 4
Параметры тренда бинарных отношений данных по табл. 1
Переменная x |
Показатель y |
Тренд y = aexp(–bxc) + dxeexp(–fxg) |
Коэф. кор. r |
||||||
Экспоненциальный закон |
Биотехнический закон |
||||||||
a |
b |
c |
d |
e |
f |
g |
|||
L, м |
, мм |
0,27696 |
–4,90586 |
0,041252 |
–2,50694e7 |
25,06937 |
27,20913 |
0,43559 |
0,7574 |
?, 0 |
47,77710 |
6,90355e-5 |
3,64290 |
8,14379e-140 |
185,8460 |
6,73438 |
1,21378 |
0,7514 |
|
?, сутки |
0 |
0 |
0 |
0,65415 |
0,90434 |
0 |
0 |
0,6804 |
|
, мм |
?, сутки |
0 |
0 |
0 |
14,42463 |
0,52119 |
0 |
0 |
0,6793 |
?, 0 |
204,11047 |
0,29180 |
1 |
30,74222 |
0,46543 |
0 |
0 |
0,5325 |
Двухчленный тренд |
Остатки от двухчленного тренда |
Рис. 2. Графики влияния расстояния от дороги на среднюю максимальную ширину группы листьев
Таблица 5
Модели (1) влияния расстояния от дороги на максимум средней ширины листа
Номер i |
Вейвлет |
Коэф. кор. r |
|||||||
Амплитуда (половина) колебания |
Полупериод колебания |
Сдвиг |
|||||||
a1i |
a2i |
a3i |
a4i |
a5i |
a6i |
a7i |
a8i |
||
1 |
5,77155 |
0 |
–1,96504 |
0,041252 |
0 |
0 |
0 |
0 |
1,0000 |
2 |
–5,59056e9 |
25,07757 |
30,15966 |
0,43558 |
0 |
0 |
0 |
0 |
|
3 |
–3,34826e-20 |
174,83658 |
46,85684 |
1,02318 |
0,81421 |
4,65539e–5 |
4,22651 |
1,29874 |
|
4 |
2308,83874 |
0 |
1,47595 |
1,73980 |
0,022770 |
0,021648 |
1,01222 |
3,29490 |
0,9523 |
Двухчленный тренд |
Первое колебание |
Двухчленный тренд и одно колебание |
Второе колебание |
Рис. 3. Графики влияния расстояния от дороги на среднюю максимальную ширину группы листьев
Прямое и обратное влияния биометрических параметров листьев по табл. 4 происходят по показательному закону. Подробнее рассмотрим влияние восточной долготы.
По графику на рис. 4 снижение времени вегетации происходит в центральном микрорайоне г. Йошкар-Ола. Наиболее сложным образом изменяется график влияния долготы на среднюю максимальную ширину учетных листьев.
На рис. 5 показана карта-схема расположения 10 учетных берез. Видно, что главный трансект получается вдоль р. Малая Кокшага. Максимальная ширина листьев снижается от запада к востоку на правой стороне города до реки, а на левой стороне реки приращение ширины возрастает, снова снижается в местности с высотными домами.
Влияние долготы на максимальную ширину |
Время вегетации на максимальную ширину |
Максимальная ширина на время вегетации |
Влияние долготы на время вегетации |
Рис. 4. Графики бинарных отношений по данным табл. 4
Рис. 5. Карта-схема расположения 10 учетных деревьев березы
Заключение
Средняя максимальная ширина учетных листьев березы повислой может стать великолепным индикатором загрязнения придорожной территории. Для идентификации волновых закономерностей необходимо в разных местах принять более 20 берез.
На каждой березе сбоку кроны со стороны преобладающих ветров выделялась зона, где выделяли 10 учетных листьев, а к корешкам их привязывали белые нитки с бирками и номерами. Сотовым телефоном записывали координаты зоны учетных листьев. Затем вычисляли приведенные координаты. Для программной среды Curve Expert 1.40 необходимо умножить полученные значения приведенных координат на 100.