Scientific journal
Advances in current natural sciences
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,775

REGULARITIES OF THE DYNAMICS OF THE AVERAGE WIDTH OF THE LEAVES OF SILVER BIRCH NEAR THE MAXIMUM GROWTH IN THE VEGETATION PERIOD

Mazurkin P.M. 1 Kudryashova A.I. 1
1 Volga State University of Technology
It is shown that the achievement of the maximum of the average width of the birch leaves is a great bioindicator of environmental quality. Moreover, the quality assessment can be carried out without measuring the concentration of pollution. According to patent 2606189 for the invention on the side at an altitude of 1.5-2.0 m from the prevailing winds on the dynamics from July 27 to September 23, 2018, the average width of 10 leaves near the road was measured. The parameters of the model with two components showed that the half-period of oscillation at the beginning of the growing season for 10 birches varies from 1.32 to 21.35 days. All models of dynamics are the same in design, but have a different character. The growing season near the maximum growth affects with a correlation coefficient of more than 0.999, so the behavior of 10 leaves has a high quantum certainty. To assess the behavior of groups of leaves, a new criterion was introduced – the coefficient of dynamism equal to the ratio of the difference between the fact and the trend to the fact of the average width of the leaves of the birch. The maximum width of the leaves and gave the highest values of the coefficient of dynamism. This criterion allows to rank the birch in descending order of the influence of vibrations from affecting the surrounding of the birch environment.
vegetation
prevailing winds
10 birch
10 account of leaves
width
dynamics
oscillatory adaptation
the dynamic factor
patterns
1. Zharkova V. The solar magnet field and the terrestrial climate. [Electronic resource]. URL: https://watchers.news/2018/11/11/valentina-zharkova-solar-magnet-field-and-terrestrial-climate-presentation/ (date of access: 01.03.2019).
2. Mazurkin P.M., Kudryashova A.I. Method of measurement of dynamics of growth of leaves of the tree in clean ecological conditions. 18th International Multidisciplinary Scientific GeoConference SGEM 2018. 2018. Vol. 18. Issue 5.1. Р. 517–524. DOI: 10.5593/sgem2018/5.1/S20.067.
3. Foster G.L., Royer D.L., Lunt D.J. Future climate forcing potentially without precedent in the last 420 million years. Nature communications. 2017. 8:14845. DOI: 10.1038/ncomms14845.
4. Mdenpp M. Impacts of temperature and ozone on carbon retention processes of birch and aspen Dissertations in Forestry and Natural Sciences. Finland, Joensuu, on June, 08, 2012. 54 р.
5. Crowther T.W., Glick H.B., Covey K.R., Bettigole C., et al. Mapping tree density at a global scale. Nature. 2015. Vol. 525. Р. 201–205. DOI:10.1038/nature14967.
6. Büntgen U., Wacker L., Galván J.D., Arnold S., et al. Tree rings reveal globally coherent signature of cosmogenic radiocarbon events in 774 and 993 CE. Nature communications. 2018. 9:3605. DOI: 10.1038/s41467-018-06036-0.
7. Mazurkin P.M., Kudryashova A.I. Quantum fitometeorologia: quanta behavior of Meteorological parameters on a three hour Measurements in the vegetation period birch. Journal of Modern Technology and Engineering. 2019. Vol. 4. № 2. P. 72–88.
8. Mazurkin P.M., Kudryashova A.I. Factor analysis of annual global carbon dynamics (according to Global Carbon Budget 2017. Materials of the International Conference «Research transfer» – Reports in English (part 2). November 28, 2018. Beijing, PRC, 2018. P. 192–224.
9. Zhang Y., Bielory L., Georgopoulos P. Climate change effect on Betula (birch) and Quercus (oak) pollen seasons in US. International Journal of Biometeorology. 2014. Vol. 58. № 5. Р. 909–919. DOI: 10.1007/s00484-013-0674-7.
10. Polgar C.A., Primack R.B. Leaf-out phenology of temperate woody plants: from trees to ecosystems. New Phytologist. 2011. Vol. 191. Р. 926–941. DOI: 10.1111/j.1469-8137.2011.03803.x.
11. Rousi M., Pusenius J. Variations in phenology and growth of European white birch (Betula Pendula) clones. Tree Physiology. 2005. Vol. 25 (2). Р. 201–210. DOI: 10.1093/treephys/25.2.201.
12. Mazurkin P.M. Method of identification. 14th International Multidisciplinary Scientific GeoConference & EXPO SGEM2014 (SGEM). Albena, Bulgaria, 2014. Vol. 1 (6). Р. 427–434.
13. Wasserstein R.L., Schirm A.L., Lazar N.A. Moving to a World Beyond «p < 0.05». The American Statistician. 2019. Vol. 73. Issue sup1. Р. 1–19. DOI: 10.1080/00031305.2019.1583913.

Эволюция климата Земли на геологические сроки в значительной степени обусловлена вариациями суммарной солнечной радиации [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] вида

mazur01.wmf,

mazur02.wmf,

mazur03.wmf, (*)

где 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

Вейвлет mazur04.wmf

Коэф.

корр.

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

 

mazur1a.wmf mazur1b.wmf

ул. Лебедева Воскресенский пр-т

mazur1c.wmf mazur1d.wmf

ул. Эшкинина пр. Ленина

Рис. 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

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

mazur2a.wmf mazur2b.wmf

б-р Чавайна ул. Панфилова

mazur2c.wmf mazur2d.wmf

ул. Карла Маркса ул. Строителей

mazur2e.wmf mazur2k.wmf

Ботанический сад ул. Осипенко

Рис. 2. Графики динамики средней ширины группы из 10 учетных листьев березы повислой

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

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

Заключение

По патенту 2606189 на изобретение сбоку на высоте 1,5–2,0 м со стороны преобладающих ветров по динамике с 15 июля по 23 сентября 2018 г. измерялись средняя ширина 10 учетных листьев около автомобильной дороги. Параметры модели с двумя составляющими показали, что полупериод колебания в начале вегетационного периода для 10 берез изменяется от 1,32 до 21,35 суток. Все модели динамики одинаковы по конструкции, но имеют разный характер. Время вегетации около максимума роста влияет с коэффициентом корреляции более 0.999, поэтому поведение 10 листьев имеет квантовую определенность.

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

Средняя ширина 10 учетных листьев березы повислой показала очень малую погрешность, поэтому количество учетных листьев с одной локальной зоны можно уменьшить до 5. Измерение ширины учетных листьев до достижения максимума роста почти в два раза сокращает промежуток времени для проведения экспериментов. Это позволяет увеличить количество учетных деревьев, произрастающих в разных экологических условиях, а также уменьшить, например до двух суток, промежутки между измерениями. Время измерений в каждый год можно начать с 01 июля и продолжать до 20 августа.