Scientific journal
Scientific Review. Biological science
ISSN 2500-3399
ПИ №ФС77-57454

ECOLOGICAL METHODS OF ASSESSING THE QUALITY OF THE ENVIRONMENT FOR FLUCTUATING ASYMMETRY OF BIRCH LEAVES

Mazurkin P.M. 1
1 Volga State University of Technology
The methods of environmental assessment of the quality of the surrounding growing birch medium by fluctuating asymmetry of the leaves after they stop growing are briefly presented. In comparison with the average values ??of the measurement parameters for each of the five leaves from at least three birch trees, the accuracy of modeling increases and the complexity of measurements decreases. The fluctuating asymmetry can be captured by statistical modeling from a much smaller volume of measurements. The smallest measurement we recommend is 15 (three trees with five differently sized leaves from different places). This will reduce the volume of measurements 100/15 7 times. The advantage of the proposed methods according to our patents for inventions is the technical simplicity of execution, since the equipment requires only a measuring compass and a geodetic protractor with a scale of division of 0.1 mm. Therefore, the invention can be widely implemented in school ecological circles, school forest districts, and even in kindergartens, as well as in geographical and other expeditions in the study of the quality of the territory on the properties of birch foliage. The methods obtained patents for inventions of the Russian Federation 2549779, 2556980, 2556985, 2556987, 2569748.
three birch trees
five leaves each
10 parameters
factor analysis
regularities

В экологических технологиях с использованием листьев березы повислой [1–3] постепенно приходит понимание необходимости моделирования взаимных связей между параметрами структуры листьев растений методом идентификации [4].

Наши российские изобретения 2549779, 2556980, 2556985, 2556987, 2569748 (соавтор – студентка [5–7]) относятся к инженерной биологии и биоиндикации окружающей среды измерениями качества ростовых органов различных видов растений, преимущественно древесных, например проб в виде листьев берез с простой и небольшой листовой пластинкой.

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

Методика измерений и анализа

На рис. 1 приведена схема измерений каждого листа: 1 – ширина левой b′ и правой b″ половинок листа (измерение проводили посередине листовой пластинки), мм; 2 – длина l′ж и l″ж второй от основания листа жилки второго порядка, мм; 3 – расстояние l′осн и l″осн между основаниями первой и второй жилок второго порядка, мм; 4 – расстояние l′к и l″к между концами этих жилок, мм; 5 – угол α′ и α″ между главной жилкой и второй от основания листа жилкой второго порядка.

Сбор материала следует проводить после остановки роста листьев (в средней полосе начиная с июля).

Для экологической оценки антропогенного воздействия на территории принимают не менее трех берез примерно в одинаковых условиях местопроизрастания. С каждой березы принимают не менее пяти листьев разных размеров со стороны оцениваемой территории. Затем измерения по пяти параметрам листа проводят с применением геодезического транспортира с ценой деления мерной шкалы 0,1 мм, причем все не менее 15 листьев принимают за популяцию из отдельных самостоятельных особей. Поэтому далее составляют таблицу результатов измерений без усреднения измеренных значений.

maz1urk1.tif

Рис. 1. Схема измерений листа

Полученную выборку статистическим моделированием подвергают факторному анализу выявлением бинарных отношений между 10 показателями, причем все 100 закономерностей идентифицируют в программной среде CurveExpert-1.40 двухчленной формулой вида

maz03.wmf, (1)

где y – показатель или зависимый количественный фактор (10 параметров по пяти показателям с двух половинок листа); x – объясняющая переменная или влияющий фактор (те же 10 параметров от каждого листа); maz04.wmf – параметры модели, получаемые идентификацией.

По результатам факторного анализа выявлением бинарных отношений между 10 показателями выполняют экологическую оценку территории по коэффициенту коррелятивной вариации, а затем из 100 биотехнических закономерностей отбирают имеющие коэффициент корреляции не менее 0,7 и рассматривают попарно пять параметров листа, а также экологическую оценку проводят по отличиям между структурой и параметрами уравнений.

Результаты измерений в чистой зоне города Звенигово Республики Марий Эл

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

Таблица 1

Результаты промеров листьев березы, произрастающих в чистой зоне города (по 5 листьев с 3 растений)

№ дерева

№ листа

Ширина b, мм

Длина lж, мм

Расстояние lосн, мм

Расстояние lк, мм

Угол α, град

левой

b′

правой b″

левой l′ж

правой l″ж

левой l′осн

правой l″осн

левой

l′к

правой l″к

левой α′

правой α″

1

1

18

17

32

30

8

8

10

10

30

28

2

19

21

33

36

9

8

12

13

34

39

3

22

19

36

32

6

8

14

14

37

33

4

17

20

31

31

5

7

13

11

44

34

5

16

17

29

31

12

11

12

16

30

31

2

б

16

17

28

29

5

6

10

10

34

37

7

17

15

31

29

3

2

10

9

43

37

8

18

18

33

32

5

5

11

9

37

36

9

16

19

28

32

5

5

10

11

38

42

10

23

21

37

40

5

6

12

13

38

37

3

11

18

19

31

33

5

4

13

12

49

49

12

17

16

30

30

5

4

10

11

46

43

13

23

21

40

40

4

3

15

15

48

42

14

16

20

37

35

6

8

13

12

44

44

15

23

22

37

40

6

7

14

15

46

44

Таблица 2

Корреляционная матрица полного факторного анализа пластических признаков листьев березы для оценки состояния среды

Влияющие факторы x

Зависимые факторы (показатели y)

Сумма

Σrx

Место

Ix

b′

b″

l′ж

l″ж

l′осн

l″осн

l′к

l″к

α′

α″

Ширина b′, мм

1

0,6530

0,8105

0,8210

0,1700

0,0905

0,6621

0,6639

0,2740

0,1222

5,2672

6

Ширина b″, мм

0,7179

1

0,7139

0,9185

0,2413

0,4052

0,7404

0,5721

0,4132

0,3913

6,1138

2

Длина l′ж, мм

0,8058

0,6806

1

0,8340

0,2299

0,1746

0,7572

0,5652

0,4036

0,2785

5,7294

5

Длина l″ж, мм

0,8175

0,9144

0,8133

1

0,3181

0,3546

0,7190

0,6315

0,3420

0,4419

6,3523

1

Расстояние l′осн, мм

0,2712

0,4867

0,3173

0,2346

1

0,8960

0,1701

0,4796

0,6639

0,4962

5,0156

8

Расстояние l″осн, мм

0,2776

0,4977

0,0818

0,2226

0,9262

1

0,2068

0,5337

0,6369

0,5953

4,9786

9

Расстояние l′к, мм

0,5290

0,7585

0,7693

0,6791

0,4647

0,3308

1

0,7367

0,5290

0,2888

6,0859

3

Расстояние l″к, мм

0,5851

0,7079

0,4921

0,6560

0,5585

0,4618

0,7553

1

0,2643

0,4204

5,9014

4

Угол α′, град

0,2573

0,3117

0,3686

0,3450

0,7079

0,6413

0,4834

0,1871

1

0,7796

5,0819

7

Угол α», град

0,1405

0,3644

0,2091

0,4233

0,4177

0,5406

0,2787

0,1283

0,7796

1

4,2822

10

Сумма Σry

5,4019

6,3749

5,5759

6,1341

5,0343

4,8954

5,773

5,4981

5,3065

4,8142

54,8083

Место Iy

6

1

4

2

8

9

3

5

7

10

0,5480

Полный факторный анализ включает 10 факторов и 102 = 100 факторных отношений. Для всех них уравнение имеет вид формулы (1).

Корреляционная матрица факторного анализа приведена в табл. 2.

Коэффициент коррелятивной вариации экологического множества из 15 листьев (по 5 листьев с 3 деревьев) равен 54,8083 / 102 = 0,5480. Этот критерий применяется при сравнении различных мест взятия проб березовых листьев.

Анализ бинарных отношений между факторами

Для этого из данных табл. 2 исключим монарные отношения, оставим только бинарные отношения с сильными факторными связями (табл. 3) с коэффициентами корреляции не менее 0,7.

Таблица 3

Корреляционная матрица сильных бинарных отношений пластических признаков листьев березы при условии r ≥ 0,7

Влияющие факторы x

Зависимые факторы (показатели y)

b′

b″

l′ж

l″ж

l′осн

l″осн

l′к

l″к

α′

α″

Ширина b′, мм

   

0,8105

0,8210

           

Ширина b″, мм

0,7179

 

0,7139

0,9185

   

0,7404

     

Длина l′ж, мм

0,8058

   

0,8340

           

Длина l″ж, мм

0,8175

0,9144

0,8133

     

0,7190

     

Расстояние l′осн, мм

         

0,8960

       

Расстояние l″осн, мм

       

0,9262

         

Расстояние l′к, мм

 

0,7585

0,7693

       

0,7367

   

Расстояние l″к, мм

 

0,7079

       

0,7553

     

Угол α′, град

       

0,7079

       

0,7796

Угол α″, град

               

0,7796

 

Остались 22 сильные бинарные зависимости. Наивысшей силой обладает формула биотехнической закономерности l″осн → l′осн. Расстояние между основаниями первой и второй жилок с правой стороны листьев наибольшим образом влияет на расстояние между основаниями первой и второй жилок с левой стороны листьев.

maz1urk2.wmf

maz1urk3.wmf

Рис. 2. График результата влияния l″осн → l′осн

Рис. 3. График результата влияния b″ → l″ж

maz1urk4.wmf

maz1urk5.wmf

Рис. 4. График результата влияния l″ж → b″

Рис. 5. График результата влияния l′осн → l″осн

maz1urk6.wmf

maz1urk7.wmf

Рис. 6. График результата влияния l′ж → l″ж

Рис. 7. График результата влияния b′ → l″ж

Для примера на рис. 2 показан график результата влияния l″осн → l′осн расстояния между основаниями первой и второй жилок с правой стороны листьев на расстояние между основаниями первой и второй жилок с левой стороны листьев при проведении измерений по нашим способу на 15 листьях у трех берез. На рис. 3 показан график результата b″ → l″ж влияния ширины на правой стороне листьев на длину второй жилки на правой стороне листьев. На рис. 4 – график результата l″ж → b″ влияния длины второй жилки на правой стороне листьев на ширину на правой стороне листьев. На рис. 5 – график результата l′осн → l″осн влияния расстояния между основаниями первой и второй жилок с левой стороны листьев на расстояние между основаниями первой и второй жилок с правой стороны листьев. На рис. 6 – график результата l′ж → l″ж влияния длины второй жилки на левой стороне листьев на длину второй жилки на правой стороне листьев. На рис. 7 – график результата b′ → l″ж влияния ширины на левой стороне листьев на длину второй жилки на правой стороне листьев.

Сравнение предлагаемых способов с прототипом

Таким образом, сравнение показывает, что предлагаемая информационная технология по обработке значений одних и тех же параметров у отдельных листьев в 57 / 22 = 2,59 раза строже по сравнению с прототипом. При этом коэффициент корреляции также в 0,7538 / 0,5480 = 1,38 раза меньше, что указывает на лучшую изменчивость. Поэтому можно сделать вывод о том, что принятие средних арифметических значений параметров листьев березы до 100 (10 деревьев по 10 листочков) является искусственным, сглаживающим изменчивость флуктуирующей асимметрии. В реальности она намного вариабельнее.

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

Распишем формулы у 22 сильных бинарных факторных связей, ранжируя в табл. 4 по убыванию коэффициента корреляции.

Таблица 4

Параметры закономерностей сильных бинарных связей факторов

п/п

maz07.wmf

maz08.wmf

Коэф.

корр.

Первая составляющая

Вторая составляющая

a1

a2

a3

a4

a5

a6

a7

a8

1

l″осн → l′осн

3,071555

0

0,071176

0

0,037407

0,39827

0,36397

0

0,9262

2

b″ → l″ж

216,90035

0

0,18129

0

0,48931

0,82839

0,079967

0

0,9185

3

l″ж → b″

0,13005

2,26719

0,057208

0

35,24084

0

0

0

0,9144

4

l′осн → l″осн

0,37763

2,21272

0,18243

0

0

0

0

0

0,8960

5

l′ж → l″ж

285,55703

0

0,13188

0

0,61588

0,94539

0,016692

0

0,8340

6

b′ → l″ж

83,922

0

0,10854

0

0,62999

0,88169

0,050474

0

0,8210

7

l″ж → b′

429,19236

0

0,15290

0

0,25493

0,89348

0,028748

0

0,8175

8

l″ж → l′ж

9,96212

0

0,019433

0

0,74876

1,025059

0,00041489

0

0,8133

9

b′ → l′ж

45,87364

0

0,081235

0

0,85588

0,92382

0,029908

0

0,8105

10

l′ж → b′

0,0097932

2,10278

0,0099842

0

7,62307

0

0

0

0,8058

11

α′ → α″

5,52381

0

0,0057460

0

2,68220

0,62964

0,0054722

0

0,7796

12

α″ → α′

1,99219

0

0,039772

0

1,065425

1,035282

0,010853

0

0,7796

13

l′к → l′ж

33,97637

0

0,041107

0

0,39111

0,61879

0,15405

0

0,7693

14

l′к → b″

0,098318

3,010187

0,19554

0

2,58419

0

0

0

0,7585

15

l″к → l′к

0,017059

3036064

0,19256

0

4,97922

0

0

0

0,7553

16

b″ → l′к

110,10051

0

0,40834

0

0,60848

1,01504

0,00041684

0

0,7404

17

l′к → l″к

4,27876

0,78319

0,023133

0

10,46382

0

0

0

0,7367

18

l″ж → l′к

0,087963

2,27808

0,058024

0

25,015367

0

0

0

0,7190

19

b″ → b′

125,22754

0

0,16146

0

0,13403

0,73296

0,12324

0

0,7179

20

b″ → l′ж

138,011555

0

0,16128

0

0,75209

0,91987

0,043949

0

0,7139

21

l″к → b″

0,14397

3,050669

0,22213

0

0

0

0

0

0,7079

22

α′ → l′осн

33,68506

0

0,044599

0

0

0

0

0

0,7079

Матричное представление модели (1) компактное, но для наглядности распишем три первых бинарные связи по отдельности в виде формул:

– влияние расстояния между основаниями первой и второй жилок с правой стороны листьев на такое расстояние (рис. 2) с левой стороны листьев

maz17.wmf; (2)

– влияние ширины на правой стороне листьев на длину второй жилки (рис. 3) на правой стороне листьев

maz18.wmf; (3)

– влияние длины второй жилки на правой стороне листьев на ширину на правой (рис. 4) стороне листьев

maz19.wmf; (4)

– влияние расстояния между основаниями первой и второй жилок с левой стороны листьев на расстояние между основаниями первой и второй жилок с правой стороны листьев

maz20.wmf; (5)

– влияние длины второй жилки на левой стороне листьев на длину второй жилки на правой стороне листьев

maz21.wmf; (6)

– влияние ширины на левой стороне листьев на длину второй жилки на правой стороне этих же листьев

maz22.wmf; (7)

– влияние длины второй жилки на правой стороне на ширину на левой стороне листьев

maz23.wmf; (8)

– влияние длины второй жилки на правой стороне листьев на длину второй жилки на левой стороне листьев

maz24.wmf; (9)

– влияние ширины на левой стороне на длину второй жилки на левой стороне листьев

maz25.wmf; (10)

– влияние длины второй жилки на левой стороне на ширину на левой стороне листьев

maz26.wmf; (11)

– влияние угла между главной и второй от основания листа жилками с левой стороны листьев на угол между главной и второй от основания листа жилками с правой стороны листьев

maz27.wmf; (12)

– влияние угла между главной и второй от основания листа жилками с правой стороны листьев на угол между главной и второй от основания листа жилками с левой стороны листьев

maz28.wmf; (13)

– влияние расстояния между концами первой и второй жилок с левой стороны листьев на длину второй жилки на левой стороне листьев

maz29.wmf; (14)

– влияние расстояния между концами первой и второй жилок с левой стороны листьев на ширину на правой стороне листьев

maz30.wmf; (15)

– влияние расстояния между концами первой и второй жилок с правой стороны листьев на расстояние между концами первой и второй жилок с левой стороны листьев

maz31.wmf; (16)

– влияние ширины на правой стороне листьев на расстояние между концами первой и второй жилок с левой стороны листьев

maz32.wmf; (17)

– влияние расстояния между концами первой и второй жилок с левой стороны листьев на расстояние между концами первой и второй жилок с правой стороны листьев

maz33.wmf; (18)

– влияние длины второй жилки на правой стороне листьев на расстояние между концами первой и второй жилок с левой стороны листьев

maz34.wmf; (19)

– влияние ширины на правой стороне листьев на ширину на левой стороне листьев

maz35.wmf; (20)

– влияние ширины на правой стороне на длину второй жилки на левой стороне листьев:

maz36.wmf; (21)

– влияние расстояния между концами первой и второй жилок с правой стороны листьев на ширину на правой стороне листьев

maz37.wmf; (22)

– влияние угла между главной и второй от основания листа жилкой с левой стороны листьев на расстояние между основаниями первой и второй жилок с левой стороны листьев

maz38.wmf. (23)

Заключение

Флюктуирующую асимметрию можно уловить статистическим моделированием из гораздо меньшего объема измерений. Наименьший объем измерений мы рекомендуем 15 (три дерева по пять разных по крупности листьев из разных мест). Это позволит снизить объем измерений в 100 / 15 ≈ 7 раз. Однако при этом повышается точность анализа флуктуирующей асимметрии.

Все парные закономерности из табл. 2 имеют общую конструкцию по формуле (1). В частных случаях образуются формулы (2–23) при условиях получения некоторыми параметрами модели значений 0 или 1. Например, формула (22) получается при условиях a = 1 и a5 = 0. А уравнение (23) получилось при условиях a2 = 0, a4 = 1 и a5 = 0. Тогда можно сделать вывод о том, что все 10 параметров листьев на момент проведения измерений взаимно изменяются по общей трендовой модели (1). При этом конкретные количественные значения параметров этой модели зависят от факторов окружающей данный лист воздушной среды. В динамике онтогенеза параметры листьев изменяются по волновым закономерностям [4]. Тогда появляется возможность перехода на индивидуальные измерения каждого учетного листа, но для этого необходимы более совершенные измерительные приборы.