Эволюция климата Земли на геологические сроки в значительной степени обусловлена вариациями суммарной солнечной радиации [1] и изменениями в содержании парниковых газов в атмосфере. Изменения климата влияют на растительность [2].
Расширение наземных растений вместе с океаном позволило обеспечить долговременную обитаемость на планете. Человечество использует ископаемые виды топлива, к середине XXI в. значения CО2 будет одинаковы с ранним эоценом (50 млн лет назад). Если концентрация CО2 продолжит расти и в XXIII в., земля будет реагировать как за полтора миллиарда лет [3].
Понимание реакции растений на факторы окружающей среды необходимо для прогнозирования изменений темпов поглощения углерода в лесах при различных условиях окружающей среды. В Финляндии имеется представление о реакции молодых лиственных деревьев на повышение температуры воздуха во взаимодействии с тропосферным озоном. Эти знания повышают шансы на разработку моделей для включения параметров, которые описывают лесную систему в изменяющихся климатических условиях [4]. Можно утверждать, что будущее сдерживания количества СО2 находится в увеличении площади лесов [5, 6].
Динамику температуры [7] можно описать множеством асимметричных вейвлетов (множеством квантов поведения) вплоть до ошибки измерений. Аналогично хорошо квантуется содержание СО2 в атмосфере. По вейвлетам универсальной конструкции изменяется динамика углерода в Европе [8].
Рост растений – сложный процесс, в его основе лежат такие фундаментальные явления, как ритмичность, полярность, дифференциация, раздражимость, корреляция. Эти процессы являются общими для онтогенеза живых организмов. Онтогенез – индивидуальное развитие организма от зиготы (или вегетативного зачатка) до природной смерти. Благодаря активной деятельности меристем и фотосинтетической активности листьев зеленое растение приобретает ряд черт, которые характеризуют его рост. В процессе онтогенеза растения рост наблюдается на протяжении основных этапов его жизненного цикла [9–11]. Поэтому в дальнейших исследованиях можно выявлять закономерности влияния метеорологических параметров на динамику вегетативных органов растений.
В экологических технологиях с использованием листьев березы повислой постепенно приходит понимание необходимости моделирования взаимных связей между параметрами структуры листьев растений методом идентификации [12]. Это изобретение относится к инженерной биологии и биоиндикации окружающей среды измерениями качества ростовых органов различных видов растений, преимущественно древесных растений, например проб в виде листьев берез без их срезания с простой и небольшой листовой пластинкой.
Цель статьи: повышение точности индикации качества окружающей листья березы локальной среды по патенту 2606189 на изобретение сбоку на высоте 1,5–2,0 м со стороны преобладающих ветров по динамике с 15 июля по 23 сентября 2018 г. средней ширины 10 учетных листьев около автомобильной дороги с интенсивным движением.
Материалы и методы исследования
В табл. 1 приведены средние ширины для 10 учетных листьев, принятых без срезания и измеренные с 27 июля по 23 сентября 2018 г. Всего получилось 18 строк. Первая точка 01 мая в 2018 г. оказалась началом вегетационного периода, поэтому ставим 0.
Колебания, в частности и тренды, записываются волновой формулой [12] вида
,
,
, (*)
где y – показатель (зависимый фактор), i – номер составляющей модели (*), m – количество членов в модели (*), x – объясняющая переменная (влияющий фактор), a1...a8 – параметры модели (*), принимающие числовые значения в ходе структурно-параметрической идентификации в программной среде CurveExpert-1.40, Ai – амплитуда (половина) вейвлета (ось y), pi – полупериод колебания (ось x).
Мы исключаем меру значимости 0,05 % [13] и оцениваем адекватность формул по коэффициенту корреляции. Все 10 статистических моделей получились с адекватностью выше 0,999. Поэтому динамика роста и спада листьев имеет высокую квантовую определенность.
Результаты исследования и их обсуждение
В табл. 2 приведены параметры модели (*) с двумя составляющими, первый из которых является трендом в виде биотехнического закона [12]. Цикличность колебания на 01.05.2018 г. в виде полупериода показывает параметр a5i. Полупериод колебания в начале вегетационного периода для 10 берез изменяется от 1,32 (ул. Панфилова) до 21,35 суток (ул. Лебедева).
На рис. 1 и 2 приведены графики двухчленных формул по данным табл. 2. Все они одинаковы по конструкции, но имеют разный характер изменения во времени. Тогда можно утверждать, что динамика параметров листьев березы повислой в период вегетации около максимума приращения является очень высокоадекватной с коэффициентом корреляции более 0,999, то есть среднее поведение 10 учетных листьев имеет квантовую определенность. Но метеопараметры за период вегетации березы повислой имеют сильную квантовую запутанность. Тогда в дальнейшем необходимо понять, как же листья березы за 180 млн лет эволюции получили такую высокую колебательную адаптивность? Такое чувство, будто бы они понимают и ведут своими вегетативными органами осознанно к изменениям в окружающей группу учетных листьев среде.
Таблица 1
Средняя ширина 10 учетных листьев березы повислой для г. Йошкар-Ола
Дата |
Время t, сутки |
ул. Лебедева |
Воскресенский пр-т |
ул. Эшкинина |
Пр. Ленина |
б-р Чавайна |
ул. Панфилова |
ул. Карла Маркса |
ул. Строителей |
Ботсад |
ул. Осипенко |
01.05 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
15.07 |
75 |
36,20 |
36,27 |
35,80 |
34,46 |
37,75 |
35,50 |
38,70 |
35,94 |
39,42 |
43,17 |
22.07 |
82 |
37,83 |
38,03 |
37,81 |
35,92 |
38,78 |
36,15 |
40,95 |
37,95 |
39,48 |
43,30 |
25.07 |
85 |
39,81 |
39,46 |
38,93 |
36,66 |
39,89 |
37,63 |
42,18 |
39,23 |
40,09 |
43,32 |
29.07 |
89 |
40,95 |
40,26 |
39,64 |
37,34 |
40,52 |
38,12 |
42,46 |
40,16 |
40,13 |
43,66 |
1.08 |
92 |
42,71 |
41,57 |
39,77 |
38,16 |
41,45 |
38,89 |
42,72 |
40,55 |
40,26 |
43,74 |
5.08 |
96 |
43,11 |
43,18 |
39,38 |
38,99 |
42,38 |
39,21 |
44,00 |
41,44 |
40,39 |
43,77 |
8.08 |
99 |
43,81 |
44,18 |
38,99 |
40,11 |
43,14 |
39,10 |
44,39 |
41,69 |
40,53 |
44,09 |
12.08 |
103 |
44,77 |
45,97 |
38,24 |
40,00 |
44,45 |
39,56 |
44,05 |
42,25 |
40,71 |
43,93 |
15.08 |
106 |
45,38 |
47,87 |
36,71 |
39,61 |
44,71 |
38,61 |
42,49 |
41,17 |
40,70 |
43,83 |
19.08 |
110 |
45,92 |
48,66 |
36,30 |
39,29 |
43,75 |
38,29 |
42,26 |
40,86 |
40,47 |
43,78 |
22.08 |
113 |
45,81 |
47,38 |
36,03 |
39,16 |
42,90 |
38,22 |
42,11 |
40,34 |
40,30 |
43,50 |
26.08 |
117 |
44,88 |
46,07 |
35,96 |
38,74 |
42,43 |
38,17 |
42,02 |
40,32 |
40,19 |
43,32 |
29.08 |
120 |
44,35 |
45,09 |
35,95 |
38,32 |
42,16 |
38,10 |
41,90 |
40,20 |
40,09 |
43,16 |
2.09 |
124 |
44,08 |
44,47 |
35,86 |
38,21 |
41,12 |
38,01 |
41,71 |
40,17 |
40,02 |
42,93 |
9.09 |
131 |
42,94 |
44,09 |
35,83 |
38,09 |
40,73 |
37,91 |
41,71 |
40,17 |
39,98 |
42,86 |
16.09 |
138 |
42,03 |
43,28 |
35,78 |
37,72 |
39,67 |
37,81 |
41,64 |
40,17 |
39,89 |
42,82 |
23.09 |
145 |
41,02 |
41,89 |
35,40 |
37,47 |
38,33 |
37,62 |
41,50 |
40,06 |
39,70 |
42,54 |
Примечание. Полужирным шрифтом показаны максимальные значения средней ширины.
Таблица 2
Модели (*) динамики средней ширины учетных листьев березы повислой
Номер i |
Вейвлет |
Коэф. корр. r |
|||||||
Амплитуда (половина) колебания |
Полупериод колебания |
Сдвиг |
|||||||
a1i |
a2i |
a3i |
a4i |
a5i |
a6i |
a7i |
a8i |
||
ул. Лебедева |
|||||||||
1 |
1,21685e–5 |
9,17015 |
6,05278 |
0,32571 |
0 |
0 |
0 |
0 |
0,9995 |
2 |
0,44598 |
0,48964 |
0,048008 |
0,79134 |
21,25481 |
0 |
0 |
–2,85512 |
|
Воскресенский пр-т |
|||||||||
1 |
0,0040275 |
2,23873 |
0,00027388 |
1,77914 |
0 |
0 |
0 |
0 |
0,9994 |
2 |
–1,06631e–22 |
13,60735 |
0,074901 |
1,09307 |
18,25100 |
–4,10921e–5 |
1,86324 |
–3,10286 |
|
ул. Эшкинина |
|||||||||
1 |
4,37213e7 |
6,24840 |
21,67608 |
0,14738 |
0 |
0 |
0 |
0 |
0,9998 |
2 |
–0,27508e–17 |
11,50074 |
0,27033 |
0,85753 |
9,26474 |
0,0010104 |
1,76082 |
–4,63677 |
|
пр. Ленина |
|||||||||
1 |
0,021506 |
2,08161 |
0,035573 |
0,88346 |
0 |
0 |
0 |
0 |
0,9999 |
2 |
–2,56548e–24 |
14,15017 |
0,17309 |
0,90399 |
1,55507 |
0,0029436 |
1,41184 |
4,23812 |
|
б-р Чавайна |
|||||||||
1 |
0,079922 |
1,54159 |
0,00039496 |
1,65732 |
0 |
0 |
0 |
0 |
0,9997 |
2 |
–4,74722e–30 |
17,78349 |
0,055937 |
1,20232 |
14,57925 |
0,00056749 |
1,68957 |
–1,52615 |
|
ул. Панфилова |
|||||||||
1 |
3,66251e6 |
7,38171 |
21,73815 |
0,16016 |
0 |
0 |
0 |
0 |
0,9998 |
2 |
–1,03801 |
1,08360 |
1,48447 |
0,27784 |
1,31728 |
0,023483 |
1,07558 |
0,91434 |
|
ул. Карла Маркса |
|||||||||
1 |
1,31774e6 |
7,87565 |
21,57557 |
0,16721 |
0 |
0 |
0 |
0 |
0,9994 |
2 |
–5,06517e–6 |
4,54130 |
0,83355 |
0,50876 |
15,76053 |
0,00020116 |
1,98175 |
–4,85485 |
|
ул. Строителей |
|||||||||
1 |
1,86986e6 |
9,10854 |
24,13993 |
0,16944 |
0 |
0 |
0 |
0 |
0,9998 |
2 |
–7,26191e–19 |
11,79845 |
0,26301 |
0,83912 |
13,57641 |
0,0017335 |
1,61783 |
–3,64787 |
|
Ботанический сад |
|||||||||
1 |
6,59123 |
0,52516 |
0,013763 |
0,82177 |
0 |
0 |
0 |
0 |
1,0000 |
2 |
–1,28495e–21 |
12,70989 |
0,13276 |
0,97603 |
20,55353 |
0,00012933 |
1,60894 |
–0,35495 |
|
ул. Осипенко |
|||||||||
1 |
8,05960e8 |
3,85974 |
20,58087 |
0,11220 |
0 |
0 |
0 |
0 |
1,0000 |
2 |
–0,049369 |
0,83696 |
1,13220 |
0,15655 |
9,57983 |
0,10275 |
0,82053 |
–5,58996 |
ул. Лебедева Воскресенский пр-т
ул. Эшкинина пр. Ленина
Рис. 1. Графики динамики средней ширины группы из 10 учетных листьев березы повислой
Для оценки поведения групп листьев введем новый критерий – коэффициент динамичности, равный отношению разницы между фактом и трендом к факту средней ширины (табл. 3). Максимумы ширины листьев дали и наибольшие значения коэффициента динамичности.
Этот критерий позволяет ранжировать березы по убыванию влияния колебаний.
По табл. 3 места по ухудшению условий распределяются так: 1) 0,0072 – ул. Осипенко, на которой почти нет автомашин; 2) 0,0073 – ботанический сад, очень малое влияние забора вокруг объекта; 3) 0,0189 – ул. Панфилова; 4) 0,0198 – ул. Лебедева; 5) 0,0270 – проспект Ленина; 6) 0,0275 – ул. Строителей; 7) 0,0327 – проспект Карла Маркса; 8) 0,0438 – бульвар Чавайна; 9) 0,0431 – ул. Эшкинина; 10) 0,0487 – Воскресенский проспект.
Так что же лучше: спокойная жизнь при низком коэффициенте динамичности или же наибольшее приращение средней ширины учетных листьев березы повислой? По средней максимальной ширине листа из табл. 1 березы по ухудшению продуктивности расположились так: 1) 48,66 мм – Воскресенский проспект; 2) 45,92 мм – ул. Лебедева; 3) 44,71 мм – бульвар Чавайна; 4) 44,39 мм – ул. Карла Маркса; 5) 44,09 мм – ул. Осипенко; 6) 42,25 мм – ул. Строителей; 7) 40,71 мм – ботанический сад; 8) 40,11 мм – проспект Ленина; 9) 39,77 мм – ул. Эшкинина; 10) 39,56 мм – ул. Панфилова. Конечно же, на скорость приращения и достижение средней максимальной ширины учетных листьев влияет и собственное физиологическое состояние дерева. Например, береза с продуктивными листьями по табл. 1 за четверо суток с 15,08 по 19,08 изменила на Воскресенском проспекте среднюю ширину листьев с 47,87 до 48,66 мм или имела скорость приращения (48,66 – 47,87) / 4 = = 0,198 мм. А на ул. Панфилова средняя скорость роста ширины листьев равна (39,56 – 39,10) / 4 = 0,115 мм.
Таблица 3
Коэффициент динамичности роста листьев березы повислой
Дата |
Время t, сутки |
ул. Лебедева |
Воскресенский пр-т |
ул. Эшкинина |
Пр. Ленина |
б-р Чавайна |
ул. Панфилова |
ул. Карла Маркса |
ул. Строителей |
Ботсад |
ул. Осипенко |
01.05 |
0 |
– |
– |
– |
– |
– |
– |
– |
– |
– |
– |
15.07 |
75 |
0,0127 |
0,0331 |
–0,0371 |
0,0050 |
0,0078 |
–0,0098 |
–0,0281 |
–0,0151 |
–0,0007 |
0,0007 |
22.07 |
82 |
–0,0336 |
–0,0168 |
0,0028 |
–0,0069 |
–0,0224 |
–0,0259 |
–0,0064 |
–0,0090 |
–0,0097 |
–0,0046 |
25.07 |
85 |
–0,0123 |
–0,0143 |
0,0274 |
–0,0050 |
–0,0134 |
0,0035 |
0,0118 |
0,0079 |
0,0022 |
–0,0065 |
29.07 |
89 |
–0,0176 |
–0,0344 |
0,0414 |
–0,0076 |
–0,0194 |
0,0043 |
0,0064 |
0,0133 |
–0,0005 |
–0,0009 |
1.08 |
92 |
0,0041 |
–0,0275 |
0,0431 |
0,0009 |
–0,0098 |
0,0169 |
0,0054 |
0,0119 |
0,0005 |
–0,0002 |
5.08 |
96 |
–0,0084 |
–0,0174 |
0,0333 |
0,0078 |
–0,0013 |
0,0177 |
0,0276 |
0,0217 |
0,0015 |
–0,0001 |
8.08 |
99 |
–0,0048 |
–0,0115 |
0,0244 |
0,0270 |
0,0089 |
0,0110 |
0,0327 |
0,0209 |
0,0038 |
0,0072 |
12.08 |
103 |
0,0047 |
0,0106 |
0,0076 |
0,0155 |
0,0317 |
0,0189 |
0,0224 |
0,0275 |
0,0073 |
0,0042 |
15.08 |
106 |
0,0122 |
0,0408 |
–0,0311 |
0,0009 |
0,0348 |
–0,0067 |
–0,0142 |
–0,0014 |
0,0068 |
0,0028 |
19.08 |
110 |
0,0198 |
0,0487 |
–0,0381 |
–0,0115 |
0,0134 |
–0,0155 |
–0,0191 |
–0,0115 |
0,0013 |
0,0033 |
22.08 |
113 |
0,0170 |
0,0203 |
–0,0417 |
–0,0165 |
–0,0041 |
–0,0167 |
–0,0211 |
–0,0250 |
–0,0024 |
–0,0015 |
26.08 |
117 |
–0,0007 |
–0,0068 |
–0,0372 |
–0,0276 |
–0,0095 |
–0,0158 |
–0,0197 |
–0,0244 |
–0,0039 |
–0,0032 |
29.08 |
120 |
–0,0084 |
–0,0249 |
–0,0320 |
–0,0375 |
–0,0096 |
–0,0151 |
–0,0190 |
–0,0254 |
–0,0052 |
–0,0047 |
2.09 |
124 |
–0,0061 |
–0,0300 |
–0,0266 |
–0,0368 |
–0,0238 |
–0,0131 |
–0,0176 |
–0,0221 |
–0,0050 |
–0,0068 |
9.09 |
131 |
–0,0108 |
–0,0121 |
–0,0119 |
–0,0289 |
–0,0068 |
–0,0054 |
–0,0044 |
–0,0114 |
–0,0015 |
–0,0019 |
16.09 |
138 |
–0,0031 |
0,0082 |
0,0039 |
–0,0226 |
0,0015 |
0,0050 |
0,0101 |
0,0031 |
0,0020 |
0,0045 |
23.09 |
145 |
0,0083 |
0,0262 |
0,0120 |
–0,0085 |
0,0096 |
0,0154 |
0,0253 |
0,0179 |
0,0040 |
0,0061 |
Примечание. Полужирным шрифтом показаны максимальные значения коэффициента динамичности.
б-р Чавайна ул. Панфилова
ул. Карла Маркса ул. Строителей
Ботанический сад ул. Осипенко
Рис. 2. Графики динамики средней ширины группы из 10 учетных листьев березы повислой
Как известно, мониторинг окружающей среды выполняют по множеству видов загрязнений, содержание которых в воздухе, воде и почве не должно превышать предельно допустимых концентраций. Однако косвенный мониторинг через поведение учетных листьев березы повислой, в вегетационный период позволяет без измерения концентраций загрязнений провести экологическую оценку территории места произрастания данной березы. В итоге листья березы реагируют на все изменения в окружающей березу среде, причем учитывают и те воздействия, которые еще неизвестны человеку. В этом проявляется сущность наших способов индикации окружающей среды, и патент 2606189 становится пионерным решением.
Для принятия за критерии экологической оценки скорости приращения необходимо учитывать еще и физиологические параметры дерева, например классы санитарного состояния, возраст и другое. Таким образом, для индикации среды вокруг березы пока рекомендуется коэффициент динамичности при максимуме средних размеров учетных листьев.
Заключение
По патенту 2606189 на изобретение сбоку на высоте 1,5–2,0 м со стороны преобладающих ветров по динамике с 15 июля по 23 сентября 2018 г. измерялись средняя ширина 10 учетных листьев около автомобильной дороги. Параметры модели с двумя составляющими показали, что полупериод колебания в начале вегетационного периода для 10 берез изменяется от 1,32 до 21,35 суток. Все модели динамики одинаковы по конструкции, но имеют разный характер. Время вегетации около максимума роста влияет с коэффициентом корреляции более 0.999, поэтому поведение 10 листьев имеет квантовую определенность.
Для оценки поведения групп листьев введен новый критерий – коэффициент динамичности, равный отношению разницы между фактом и трендом к факту средней ширины. Максимумы ширины листьев дали и наибольшие значения коэффициента динамичности. Этот критерий позволяет ранжировать березы по убыванию влияния колебаний от влияния среды.
Средняя ширина 10 учетных листьев березы повислой показала очень малую погрешность, поэтому количество учетных листьев с одной локальной зоны можно уменьшить до 5. Измерение ширины учетных листьев до достижения максимума роста почти в два раза сокращает промежуток времени для проведения экспериментов. Это позволяет увеличить количество учетных деревьев, произрастающих в разных экологических условиях, а также уменьшить, например до двух суток, промежутки между измерениями. Время измерений в каждый год можно начать с 01 июля и продолжать до 20 августа.