Результат интеллектуальной деятельности: СПОСОБ ИЗМЕРЕНИЯ ФЛУКТУИРУЮЩЕЙ АСИММЕТРИИ ЛИСТЬЕВ БЕРЕЗЫ

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

Известен способ измерения листьев у древесных растений (см. книгу: Федорова А.И., Никольская А.Н. Практикум по экологии и охране окружающей среды: Учеб. пос. М.: Гуманит. изд. центр ВЛАДОС, 2001. 288 с. С.123-126), когда срезают по 20-25 листьев берез и других пород деревьев, растущих в разных экологических условиях, складывают в пакеты, а затем засушивают между листами газетной бумаги в лабораторных условиях и измеряют параметры у каждого листа, а затем обрабатывают измеренные данные по средним значениям.

Недостатком является низкая чувствительность (точность) индикации и высокая трудоемкость практического применения способа из-за смешивания листьев у разных деревьев в одну пробу.

Известен также способ измерения по В.М. Захарову флуктуирующей асимметрии листьев березы (Захаров, В.М. К оценке асимметрии билатеральных признаков как популяционной характеристики / В.М. Захаров, В.В. Зюганов // Экология. - 1980. - №1, с.10-16; Захаров, В.М. Онтогенез и популяция (стабильность развития и популяционная изменчивость) / В.М. Захаров // Экология. 2001. - №3. - С.164-168); Здоровье среды: методика оценки / Захаров В.М. [и др.]. - М.: Центр экологической политики России, 2000. - 68 с.; Здоровье среды: практика оценки / Захаров В.М. [и др.]. - М.: Центр экологической политики России, 2000. - 320 с.), учитывающий физико-биологический эффект, что стрессирующие воздействия различного типа вызывают в живых организмах изменения гомеостаза (стабильности) развития, которые могут быть оценены по нарушению морфогенетических процессов (Шкиль, Ф.Н. Применение методики раннего выявления нарушений состояния зеленых насаждений / Н.Ф. Шкиль, В.М. Захаров // Экология большого города. Альманах. Вып.8. Проблемы содержания зеленых насаждений и городских лесов в условиях Москвы. - М.: Прима-М, 2003. - С.50-54).

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

Никакой специальной обработки и подготовки материала не требуется. Материал может быть обработан сразу после сбора или позднее. Для непродолжительного хранения собранный материал можно хранить в полиэтиленовом пакете на нижней полке холодильника. Для длительного хранения можно зафиксировать материал в 60% растворе этилового спирта или гербаризировать.

Можно использовать систему признаков, разработанную для березы. Для измерения лист помещают пред собой стороной, обращенной к верхушке побега. С каждого листа снимают показатели по пяти промерам с левой и правой сторон листа (фиг.1).

Для измерений требуются измерительный циркуль, линейка и транспортир. Промеры 1-4 снимаются циркулем-измерителем, угол между жилками (признак 5) измеряется транспортиром.

Величина асимметричности оценивается с помощью интегрального показателя - величины среднего относительного различия на признак (средняя арифметическая отношения разности к сумме промеров листа слева и справа, отнесенная к числу признаков).

В первом действии находят относительное различие между значениями признака слева и справа - для каждого признака. Для этого находят разность значений измерений по одному признаку для одного листа, затем находят сумму этих же значений и разность делят на сумму.

Подобные вычисления производят по каждому признаку (от 1 до 5). В результате получается 5 значений ординаты для одного листа. Такие же вычисления производят для каждого листа в отдельности.

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

Подобные вычисления производят для каждого листа.

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

Полученный показатель характеризует степень асимметричности организма.

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

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

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

Основным недостатком известного способа является то, что нарушается принцип индивидуальности флуктуирующей асимметрии, когда берутся в каждой выборке 100 листьев (по 10 листьев с 10 растений). Это выполнено в угоду существующей теории средней арифметической величины. Но такого параметра нет у листьев, поэтому каждый листочек имеет свои значения параметров, которые нужно определять. Существующая классическая статистика не может ничего сказать об индивидуальных свойствах каждого из 100 листьев, она рассматривает эти 100 листьев только как однородные случайные события. Но каждый листок - это не случайность, а закономерность. Это биологам очевидно, поэтому они пошли на поводу у математиков.

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

Этот технический результат достигается тем, что способ измерения флуктуирующей асимметрии листьев березы, включающий взятие листьев от учетных деревьев, растущих в одинаковых экологических условиях местопроизрастания, причем листья с одного дерева хранятся отдельно, чтобы можно было проанализировать полученные результаты индивидуально для каждой березы, а для этого следует собранные с одного дерева листья связывать за черешки, причем все листья, собранные для одной выборки, следует сложить в полиэтиленовый пакет, туда же вложить этикетку, в которой указаны номер выборки, место сбора, делая максимально подробную привязку к местности, дату сбора, причем для непродолжительного хранения собранный материал хранится в полиэтиленовом пакете на нижней полке холодильника, а для длительного хранения можно зафиксировать материал в 60% растворе этилового спирта или гербаризировать, при этом для измерения каждый лист помещают перед собой стороной, обращенной к верхушке побега, с каждого листа снимают показатели по пяти параметрам листа с левой и правой сторон листа, а для измерений применяют измерительный циркуль, линейку и транспортир, причем промеры длин снимаются циркулем-измерителем, а угол между жилками измеряется транспортиром, отличающийся тем, что для экологической оценки антропогенного воздействия на территории принимают не менее трех берез примерно в одинаковых условиях местопроизрастания, с каждой березы принимают не менее пяти листьев разных размеров со стороны оцениваемой территории, затем измерения по пяти параметрам листа проводят с применением геодезического транспортира с ценой деления мерной шкалы 0,1 мм, причем все не менее 15 листьев принимают за популяцию из отдельных самостоятельных особей, поэтому далее составляют таблицу результатов измерений без усреднения измеренных значений, а полученную выборку статистическим моделированием подвергают факторному анализу выявлением бинарных отношений между 10 показателями, причем все 100 биотехнических закономерностей идентифицируют в программной среде формулой вида:

,

где y - показатель или зависимый количественный фактор (10 параметров по пяти показателям с двух половинок листа);

x - объясняющая переменная или влияющий фактор (те же 10 параметров каждого листа);

a ₁…a ₈ - параметры модели, получаемые идентификацией.

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

Сущность изобретения заключается в том, что в основу методики положена теория «стабильности развития» («морфогенетического гомеостаза»), разработанная российскими учеными В.М. Захаровым и др. в процессе исследований последствий радиоактивного заражения, в том числе после Чернобыльской аварии (Захаров, В.М. Последствия Чернобыльской катастрофы: Здоровье среды / В.М. Захаров, Е.Ю. Крысанов. - М.: Центр экологической политики России, 1996. - 170 с.). Растения же как продуценты экосистемы в течение своей жизни привязаны к локальной территории и подвержены влиянию почвенной и воздушной сред, наиболее полно отражающих весь комплекс стрессирующих воздействий на экосистему.

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

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

Флуктуирующая асимметрия крайне широко распространенное явление. Им охвачены практически все билатеральные структуры у самых разных живых существ. Понятно, что невозможно подвергнуть анализу известные признаки всех билатерально-симметричных структур, но у исследованных флуктуирующая асимметрия регистрировалась. Более того, это явление имеет место даже при иных типах асимметрии, в этом случае она представляет собой отклонения не от строгой симметрии, а от определенной средней симметрии.

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

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

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

Сущность изобретения заключается также и в том, что каждый листок получает свои значения измеренных по методу В.М. Захарову пяти параметров.

Новизна технического решения заключается в том, что впервые статистическая обработка проводится не по обезличенному методу измерений 100 листьев березы, а четко по учету флуктуации параметров у каждого листа из минимального количества из 15 листьев. При этом каждый лист принимается за биологическую особь, которая питается от материнского дерева, но живет полностью обособленно от других листьев. На него действует разные силы извне, поэтому форма и размеры листа разные. И тогда никак нельзя множество листьев принять за выборку (в классической статистике выборка всегда предполагает независимость событий, но одинаковость по физическому смыслу). Но при этом один лист мешает другому даже по освещенности. Поэтому 15 листьев являются взаимно зависимыми объектами и поэтому не могут быть обработаны методами классической математики. В итоге нужен наш биотехнический подход.

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

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

На фиг.1 приведена схема измерений каждого листа: 1 - ширина левой b' и правой b'' половинок листа (измерение проводили посередине листовой пластинки), мм; 2 - длина и второй от основания листа жилки второго порядка, мм; 3 - расстояние и между основаниями первой и второй жилок второго порядка, мм; 4 - расстояние и между концами этих жилок, мм; 5 - угол α' и α'' между главной жилкой и второй от основания листа жилкой второго порядка; на фиг.2 показан график результата влияния расстояния между основаниями первой и второй жилок с правой стороны листьев на расстояние между основаниями первой и второй жилок с левой стороны листьев при проведении измерений по предлагаемому способу на 15 листьях у трех берез; на фиг.3 - то же на фиг.2 график результата влияния ширины на правой стороне листьев на длину второй жилки на правой стороне листьев; на фиг.4 - то же на фиг.2 график результата влияния длины второй жилки на правой стороне листьев на ширину на правой стороне листьев; на фиг.5 - то же на фиг.2 график результата влияния расстояния между основаниями первой и второй жилок с левой стороны листьев на расстояние между основаниями первой и второй жилок с правой стороны листьев; на фиг.6 - то же на фиг.2 график результата влияния длины второй жилки на левой стороне листьев на длину второй жилки на правой стороне листьев; на фиг.7 - то же на фиг.2 график результата влияния ширины на левой стороне листьев на длину второй жилки на правой стороне листьев; на фиг.8 - то же на фиг.2 график результата влияния длины второй жилки на правой стороне листьев на ширину на левой стороне листьев; на фиг.9 - то же на фиг.2 график результата влияния длины второй жилки на правой стороне листьев на длину второй жилки на левой стороне листьев; на фиг.10 - то же на фиг.2 график результата влияния ширины на левой стороне листьев на длину второй жилки на левой стороне листьев; на фиг.11 - то же на фиг.2 график результата влияния длины второй жилки на левой стороне листьев на ширину на левой стороне листьев; на фиг.12 - то же на фиг.2 график результата α'→α'' влияния угла между главной и второй от основания листа жилками с левой стороны листьев на угол между главной и второй от основания листа жилками с правой стороны листьев; на фиг.13 - то же на фиг.2 график результата α''→α' влияния угла между главной и второй от основания листа жилками с правой стороны листьев на угол между главной и второй от основания листа жилками с левой стороны листьев; на фиг.14 - то же на фиг.2 график результата влияния расстояния между концами первой и второй жилок с левой стороны листьев на длину второй жилки на левой стороне листьев; на фиг.15 - то же на фиг.2 график результата влияния расстояния между концами первой и второй жилок с левой стороны листьев на ширину на правой стороне листьев; на фиг.16 - то же на фиг.2 график результата влияния расстояния между концами первой и второй жилок с правой стороны листьев на расстояние между концами первой и второй жилок с левой стороны листьев; на фиг.17 - то же на фиг.2 график результата влияния ширины на правой стороне листьев на расстояние между концами первой и второй жилок с левой стороны листьев; на фиг.18 - то же на фиг.2 график результата влияния расстояния между концами первой и второй жилок с левой стороны листьев на расстояние между концами первой и второй жилок с правой стороны листьев; на фиг.19 - то же на фиг.2 график результата влияния длины второй жилки на правой стороне листьев на расстояние между концами первой и второй жилок с левой стороны листьев; на фиг.20 - то же на фиг.2 график результата влияния ширины на правой стороне листьев на ширину на левой стороне листьев; на фиг.21 - то же на фиг.2 график результата влияния ширины на правой стороне листьев на длину второй жилки на левой стороне листьев; на фиг.22 - то же на фиг.2 график результата влияния расстояния между концами первой и второй жилок с правой стороны листьев на ширину на правой стороне листьев; на фиг.23 - то же на фиг.2 график результата влияния угла между главной и второй от основания листа жилкой с левой стороны листьев на расстояние между основаниями первой и второй жилок с левой стороны листьев.

Способ измерения флуктуирующей асимметрии листьев березы включает такие действия.

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

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

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

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

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

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

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

2. Условия произрастания. Листья должны быть собраны с растений, находящихся в одинаковых экологических условиях (уровень освещенности, увлажнения и т.д.). Рекомендуется выбирать растения, растущие на открытых участках (полянах, опушках), поскольку многие виды светолюбивы и условия затенения являются для них стрессовыми и могут существенно снизить стабильность развития.

3. Возрастное состояние растения. Для исследования рекомендуется выбирать растения, достигшие генеративного возрастного состояния.

4. Сбор листьев с растения. Для исследования предлагается использовать лист, как орган, обладающий билатеральной симметрией.

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

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

7. Размер листьев должен быть сходным, средним для растения.

8. Поврежденность листьев. Поврежденные листья могут быть использованы для анализа, если не затронуты участки, с которых будут сниматься измерения. Рекомендуется собирать с растения несколько больше листьев, чем требуется, на тот случай, если часть листьев из-за повреждения не сможет быть использована для анализа.

По нашему научно-техническому решению способ исключает учет выделенных выше отличительных признаков 2, 6 и 7 по прототипу и дополнительно включает следующие действия.

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

,

где y - показатель или зависимый количественный фактор (10 параметров по пяти показателям с двух половинок листа);

x - объясняющая переменная или влияющий фактор (те же 10 параметров от каждого листа);

a ₁…a ₈ - параметры модели, получаемые идентификацией.

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

Примеры

Пробы березовых листьев в Кугесинской школе (по прототипу).

Параметры листа имеют следующие условные обозначения:

1) ширина левой b' и правой b'' половинок листа, мм;

2) длина жилки l_ж второго порядка, второй от основания листа, мм;

3) расстояние l_осн между основаниями первой и второй жилок второго порядка, мм;

4) расстояние l_к между концами этих же жилок, мм;

5) угол α между главной жилкой и второй от основания листа жилкой второго порядка, град.

Результаты эксперимента отображены в таблице 2, в котором даны средние значения параметров для каждого дерева.

Таблица 2 Исходные данные по средним значениям промеров (по 10 листьев с 10 растений) № Ширина b, мм Длина lж, мм Расстояние lосн, мм Расстояние lк, мм Угол α, град левой b' правой b'' левой правой левой правой левой правой левой α' правой α'' 1 18 20 32 33 4 4 12 12 46 50 2 20 19 33 33 3 3 14 13 50 49 3 18 18 31 31 2 3 12 11 50 46 4 18 19 30 32 2 3 10 11 49 49 5 20 20 30 33 6 3 13 14 46 53 6 12 14 22 22 4 4 11 9 39 39 7 14 12 26 25 3 3 11 11 34 40 8 13 14 25 23 3 3 10 8 39 42 9 12 14 24 25 5 5 9 9 40 32 10 14 14 25 25 4 4 9 8 32 32

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

Показатель асимметрии, рассчитанный в Кугесинской школе, равен 0,042, что соответствует второму баллу шкалы (табл.1). Это означает, что растения испытывают слабое влияние неблагоприятных факторов.

Далее данные таблицы обрабатывали в программной среде Curve Expert-1.3, результаты обработки представлены в таблице 3.

Асимметрия березовых листьев из Кугесинской школы по предлагаемому способу. Полный факторный анализ таблицы 2 включает в себя 10 факторов и 10²=100 факторных отношений, из которых десять будут показывать влияние фактора самого на себя в виде ранговых распределений (монарные отношения), а остальные 100-10=90 зависимостей будут показывать характер бинарных отношений.

Если общее количество факторов равно m, то получим формулу общей численности биотехнических закономерностей m², из которых количество монарных отношений будет равно m, а количество бинарных факторных отношений (отношений факторов друг с другом) составит m²-m.

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

Общее уравнение тренда (тенденции), то есть детерминированной закономерности без учета волновых возмущений для всех 10²=100 факторных отношений, имеет вид

где y - показатель или зависимый количественный фактор (в нашем примере параметр листа);

x - объясняющая переменная или влияющий фактор;

a ₁…a ₈ - параметры модели, получаемые идентификацией на ПЭВМ.

Корреляционная матрица полного факторного анализа приведена в таблице 3. Она включает бинарные и ранговые распределения показателей.

Сильные факторные связи при коэффициенте R≥0,7 приведены с четырьмя значащими цифрами.

Таблица 3 Корреляционная матрица полного факторного анализа пластических признаков листьев березы для опенки состояния среды Влияющие факторы x Зависимые факторы (показатели y) Сумма Σrx Место Ix b' b'' α' α'' Ширина b', мм 0,9977 0,9070 0,9595 0,9644 0,615 0,621 0,8436 0,9064 0,9753 0,9236 8,7135 3 Ширина b'', мм 0,9079 0,9953 0,9511 0,9803 0,8247 0,655 0,7364 0,9326 0,9050 0,8515 8,7398 2 Длина , мм 0,9498 0,8723 0,9845 0,9580 0,255 0,585 0,8211 0,8436 0,9426 0,9161 8,1280 7 Длина , мм 0,9642 0,9772 0,9573 0,9970 0,689 0,379 0,7954 0,8818 0,9402 0,9470 8,5281 4 Расстояние , мм 0,099 0,018 0,319 0,065 1 1,0000 0,022 0,153 0,266 0,127 3,0690 10 Расстояние , мм 0,563 0,315 0,494 0,382 0,498 1 0,526 0,464 0,387 0,658 5,2870 9 Расстояние , мм 0,9023 0,8609 0,8487 0,8505 0,021 0,532 1 0,8539 0,8141 0,9655 7,6489 8 Расстояние , мм 0,9359 0,7951 0,9362 0,8994 0,8443 0,472 0,8566 0,9936 0,7084 0,8374 8,2789 5 Угол α', град 0,9606 0,9911 0,9521 0,9897 0,497 0,506 0,648 0,678 0,9909 0,9477 8,1611 6 Угол α'', град 0,9494 0,9885 0,9516 0,9753 0,8990 0,6999 0,7806 0,676 0,9240 0,9857 8,8300 1 Сумма Σry 8,2298 7,7204 8,354 8,0616 6,143 6,4499 7,0297 7,3829 7,8535 8,1595 75,3843 - Место Iy 2 6 1 4 10 9 8 7 5 3 - 0,7538

Коэффициент коррелятивной вариации экологического (раздел биологического) множества из 100 листьев (по 10 листьев от 10 деревьев) равен 75,3843/10²=0,7538. Этот критерий применяется при сравнении различных мест взятия проб березовых листьев. Иначе говоря, этот статистический показатель становится мерилом для сопоставления разных мест.

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

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

Таблица 4 Корреляционная матрица сильных бинарных отношений пластических признаков листьев березы при условии r≤0,7 Влияющие факторы x Зависимые факторы (показатели γ) b' b'' α' α'' Ширина b', мм   0,9070 0,9595 0,9644     0,8436 0,9064 0,9753 0,9236 Ширина b'', мм 0,9079   0,9511 0,9803 0,8247   0,7364 0,9326 0,9050 0,8515 Длина , мм 0,9498 0,8723   0,9580     0,8211 0,8436 0,9426 0,9161 Длина , мм 0,9642 0,9772 0,9573       0,7954 0,8818 0,9402 0,9470 Расстояние , мм           1,0000         Расстояние , мм                     Расстояние , мм 0,9023 0,8609 0,8487 0,8505       0,8539 0,8141 0,9655 Расстояние , мм 0,9359 0,7951 0,9362 0,8994 0,8443   0,8566   0,7084 0,8374 Угол α'', град 0,9606 0,9911 0,9521 0,9897           0,9477 Угол α', град 0,9494 0,9885 0,9516 0,9753 0,8990   0,7806   0,9240

Остались 57 сильных бинарных зависимостей. Тогда доля сильных факторных бинарных связей будет равна 100×57/(100-10)=63,33%. Это очень высокий показатель

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

Таблица 5 Результаты промеров листьев березы, произрастающих в чистой зоне (по 5 листьев с 3 растений) № дерева № листа Ширина b, мм Длина lж, мм Расстоян. lосн, мм Расстояние lк, мм Угол α, град левой b' правой b'' левой правой левой правой левой правой левой α' правой α'' 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 6 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

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

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

Таблица 6 Корреляционная матрица полного факторного анализа пластических признаков листьев березы для оценки состояния среды Влияющие факторы x Зависимые факторы (показатели γ) Сумма Σrx Место Ix b' b'' α' α'' Ширина 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 Длина , мм 0,8058 0,6806 1 0,8340 0,2299 0,1746 0,7572 0,5652 0,4036 0,2785 5,7294 5 Длина , мм 0,8175 0,9144 0,8133 1 0,3181 0,3546 0,7190 0,6315 0,3420 0,4419 6,3523 1 Расстояние , мм 0,2712 0,4867 0,3173 0,2346 1 0,8960 0,1701 0,4796 0,6639 0,4962 5,0156 8 Расстояние , мм 0,2776 0,4977 0,0818 0,2226 0,9262 1 0,2068 0,5337 0,6369 0,5953 4,9786 9 Расстояние , мм 0,5290 0,7585 0,7693 0,6791 0,4647 0,3308 1 0,7367 0,5290 0,2888 6,0859 3 Расстояние , мм 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

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

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

Таблица 7 Корреляционная матрица сильных бинарных отношений пластических признаков листьев березы при условии r≥0,7 Влияющие факторы x Зависимые факторы (показатели y) b' b'' α' α'' Ширина b', мм     0,8105 0,8210             Ширина b'', мм 0,7179   0,7139 0,9185     0,7404       Длина , мм 0,8058     0,8340             Длина , мм 0,8175 0,9144 0,8133       0,7190       Расстояние , мм           0,8960         Расстояние , мм         0,9262           Расстояние , мм   0,7585 0,7693         0,7367     Расстояние , мм   0,7079         0,7553       Угол α', град         0,7079         0,7796 Угол α'', град                 0,7796

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

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

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

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

Таблица 8 Параметры закономерностей сильных бинарных связей факторов № п/п x→y Коэфф. корр. Первая составляющая Вторая составляющая a1 a2 a3 a4 a5 a6 a7 a8 1 3,071555 0 0,071176 0 0,037407 0,39827 0,36397 0 0,9262 2 216,90035 0 0,18129 0 0,48931 0,82839 0,079967 0 0,9185 3 0,13005 2,26719 0,057208 0 35,24084 0 0 0 0,9144 4 0,37763 2,21272 0,18243 0 0 0 0 0 0,8960 5 285,55703 0 0,13188 0 0,61588 0,94539 0,016692 0 0,8340 6 83,922 0 0,10854 0 0,62999 0,88169 0,050474 0 0,8210 7 429,19236 0 0,15290 0 0,25493 0,89348 0,028748 0 0,8175 8 9,96212 0 0,019433 0 0,74876 1,025059 0.00041489 0 0,8133 9 45,87364 0 0,081235 0 0,85588 0,92382 0,029908 0 0,8105 10 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 33,97637 0 0,041107 0 0,39111 0,61879 0,15405 0 0,7693 14 0,098318 3,010187 0,19554 0 2,58419 0 0 0 0,7585

№ п/п x→y Коэфф. корр. Первая составляющая Вторая составляющая a1 a2 a3 a4 a5 a6 a7 a8 15 0,017059 3036064 0,19256 0 4,97922 0 0 0 0,7553 16 110,10051 0 0,40834 0 0,60848 1,01504 0,00041684 0 0,7404 17 4,27876 0,78319 0,023133 0 10,46382 0 0 0 0,7367 18 0,087963 2,27808 0,058024 0 25,015367 0 0 0 0,7190 19 125,22754 0 0,16146 0 0,13403 0,73296 0,12324 0 0,7179 20 138,011555 0 0,16128 0 0,75209 0,91987 0,043949 0 0,7139 21 0,14397 3,050669 0,22213 0 0 0 0 0 0,7079 22 33,68506 0 0,044599 0 0 0 0 0 0,7079

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

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

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

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

- влияние расстояния между основаниями первой и второй жилок с левой стороны листьев на расстояние между основаниями первой и второй жилок (фиг.5) с правой стороны листьев

- влияние длины второй жилки на левой стороне листьев на длину второй жилки (фиг.6) на правой стороне листьев

- влияние ширины на левой стороне листьев на длину второй жилки на правой (фиг.7) стороне листьев

- влияние длины второй жилки на правой стороне листьев на ширину на левой (фиг.8) стороне листьев

- влияние длины второй жилки на правой стороне листьев на длину второй жилки (фиг.9) на левой стороне листьев

- влияние ширины на левой стороне листьев (фиг.10) на длину второй жилки на левой стороне листьев

- влияние длины второй жилки на левой стороне листьев (фиг.11) на ширину на левой стороне листьев

- влияние угла между главной и второй от основания листа (фиг.12) жилками с левой стороны листьев на угол между главной и второй от основания листа жилками с правой стороны листьев

- влияние угла между главной и второй от основания листа жилками с правой стороны листьев (фиг.13) на угол между главной и второй от основания листа жилками с левой стороны листьев

- влияние расстояния между концами первой и второй жилок с левой стороны листьев (фиг.14) на длину второй жилки на левой стороне листьев

- влияние расстояния между концами первой и второй жилок с левой стороны листьев (фиг.15) на ширину на правой стороне листьев

- влияние расстояния между концами первой и второй жилок с правой стороны листьев (фиг.16) на расстояние между концами первой и второй жилок с левой стороны листьев

- влияние ширины на правой стороне листьев на расстояние между концами первой и второй жилок (фиг.17) с левой стороны листьев

- влияние расстояния между концами первой и второй жилок (фиг.18) с левой стороны листьев на расстояние между концами первой и второй жилок с правой стороны листьев

- влияние длины второй жилки на правой стороне листьев (фиг.19) на расстояние между концами первой и второй жилок с левой стороны листьев

- влияние ширины па правой стороне листьев (фиг.20) на ширину на левой стороне листьев:

- влияние ширины на правой стороне листьев (фиг.21) на длину второй жилки на левой стороне листьев

- влияние расстояния (фиг.22) между концами первой и второй жилок с правой стороны листьев на ширину на правой стороне листьев

- влияние угла между главной и второй от основания листа жилкой (фиг.23) с левой стороны листьев на расстояние между основаниями первой и второй жилок с левой стороны листьев

Из 22 формул с сильными факторными связями попарно группируем те, которые в общей математической модели поведения популяции из 15 листьев образуют циклы внутри динамики флуктуирующей асимметрии и дают рейтинг параметров листьев.

Первое место - 3 (номер из фиг.1) - расстояние и между основаниями первой и второй жилок второго порядка, мм:

Здесь видно, что асимметричность наблюдается и по структуре статистической модели: у второго уравнения нет стрессового возбуждения.

Второе место - 2 (фиг.1) - длина и второй от основания листа жилки второго порядка, мм:

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

Третье место - 5 (фиг.1) - угол α' и α'' между главной жилкой и второй от основания листа жилкой второго порядка, град:

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

Четвертое место - 4 (фиг.1) - расстояние и между концами этих жилок, мм:

Здесь структура упростилась, и при этом постоянный член получил разный знак, то есть структура модели (1) все же изменилась.

Пятое место - 1 (фиг.1) - ширина левой b' и правой b'' половинок листа (измерение проводили посередине листовой пластинки), мм:

Пара этого уравнения (20) находится (см. табл.6) в зоне средней тесноты связи и получает коэффициент корреляции 0,6530.

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

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

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