К геогр н. Наумов М. М

Рассмотрение поставленного вопроса начнем с установленных фактов. В агрометеорологии хорошо известен метод сумм эффективных температур (Шиголев А.А.), который однозначно говорит о том, что скорость прохождения фенофаз зависит от напряженности температурного фактора [1, 2]. Согласно этому методу, чем выше напряженность температурного фактора, тем быстрее наступает очередная фенофаза. То есть с повышением уровня фактора тепла развитие растений идет быстрее. Таким образом, метод сумм эффективных температур является первым методом, который рассчитывает биологическую шкалу времени только время в нем измеряется не секундами, а градусами Цельсия. Для перехода в этом методе к единицам времени достаточно нормировать всю шкалу по отношению к константе сумм эффективных температур, которая необходима для прохождения всего онтогенеза. Тогда сразу же однозначно мы получаем в ответ, что биологическое время имеет свойство растягиваться или сжиматься в зависимости от напряженности температурного фактора. Более поздние исследования, например Улановой Е.С. [3], показывают, что скорость развития растений зависит так же от напряженности фактора влаги. Зависимость скорости развития растений от напряженности фактора влаги аналогична зависимости скорости развития растений от напряженности фактора тепла. Обобщением таких исследований является работа В.П. Дмитренко [4] в которой скорость развития растений зависит от напряженности факторов света, тепла, влаги, питания. Рассматривая наступление фенологических фаз развития растений как наступление определенных, биологических событий, проходящих в организме растения, мы можем заключить, что шкала биологического времени растягивается или сжимается в зависимости от движений факторов света, тепла, влаги, питания, то есть от основных, наиболее подвижных факторов внешней среды. Согласно такому подходу можно прогнозировать, что длительность межфазных периодов будет зависеть так же и от остальных, не учтенных факторов внешней среды при их движении.

Рисунок 1 – Представление элементарного приращения ΔТ_эф.отн. метода сумм эффективных температур в виде поворота реального вектора ΔТ_эф.max. ΔТ_эф.отн.=ΔТ_эф.пр.

Переход к дифференциальной форме выражения оси биологического времени можно осуществить в виде:

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

где - среднее значение напряженности факторов внешней среды от момента времени Т_н до текущего значения Т.

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

Среднее значение уровня напряженности факторов внешней среды для определенного отрезка интегрирования [Т_н; Т] находится как интеграл

Подставляя (6) в (5) мы получаем двойной интеграл со своей областью интегрирования

где G_F – квадратная область интегрирования со своей границей.

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

.

Интегрирование уравнения (3) естественным образом дало двухмерное пространство биологического временного поля.

Здесь следует отметить, что развитие растительного организма начинается с семени и заканчивается опять таки - семенем, проходя все этапы, стадии и другие формы развития. В этом смысле нецелесообразно рассматривать пространство временных переменных T_F и T_R в декартовой плоскости. Структура биологического времени должна иметь замкнутый характер движений в связи с замкнутым процессом развития растений (от семени до семени) и повторением цикла онтогенеза при каждой вегетации Смотрите так же работы [9, 10]. В связи с этим, в уравнении (7) мы перейдем к полярным координатам. То есть, мы будем рассматривать ту же самую декартову временную плоскость, только временные переменные будут выражены в полярных координатах:

. (8)

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

для отрезка биологического времени [Т_н; 1].

Учитывая (9), интеграл (8) принимает вид:

, (10)

где θ – угол поворота вектора биологического времени;

При переходе от уравнения (8) к уравнению (10) мы использовали следующую конструкцию: мы рассматриваем движения биологического времени в координатной плоскости переменных T_F и T_R; при этом считается, что физическое время имеет свойство движения равномерно слева направо; движения биологического времени осуществляются по кругу внутри организма; точнее от -επ до +επ не совершая полного оборота, работа [11]; в этом случае взаимосвязь физического и биологического времен, как переход от (8) к (10) осуществится как интегрирование уравнения (10).

Теперь мы можем перейти непосредственно к построению автономной системы двух дифференциальных уравнений биологического временного поля. Рассмотрим подынтегральное выражение уравнения (10). Это дифференциальное уравнение двух переменных T и θ. Оно имеет вид определителя

(11)

и, согласно работам Л.С. Понтрягина и А.А Андронова [12], являются вторыми производными дифференциальной системы уравнений:

Нам остается найти явный вид системы дифференциальных уравнений (12). Для этого необходимо проинтегрировать один раз систему дифференциальных уравнений в частных производных (11) (см. [12]):

Таким образом, мы нашли автономную систему двух дифференциальных уравнений в пространстве двух переменных: Т и θ. Здесь, в (13) переменная t является физическим временем, которое течет равномерно и прямолинейно. Система (13) определяет структуру биологического временного поля для первой половины онтогенеза. Для второй половины онтогенеза аналогично получим:

Таким образом, мы получили две системы двух дифференциальных уравнений: (13) и (14) для двух отрезков онтогенеза до цветения и после цветения. Переход от одной системы уравнений (13) к другой системе уравнений (14) осуществляется в растении в момент цветения. Отметим так же и то, что полученные системы (13) и (14) согласно работе [12] являются грубыми динамическими системами. Грубые динамические системы заполняют области пространства динамических систем в целом. Грубые системы можно рассматривать как наиболее простые в соответствующем пространстве динамических систем [13]. Понятие грубости системы соответствует понятию надежности растительных систем [14].

Таблица 1 - Совокупность интегральных траекторий системы биологического временного поля растений (13) в сопоставлении с клеточной структурой растений

Кроме того, мы можем отметить, что динамическая система (13), а так же (14) разбивает область G_R на конечное число связных ячеек – растительных клеток, заполненных обыкновенными временными траекториями. Такие клетки разбиваются на два класса: на класс клеток, примыкающих к циклу без контакта, и на класс внутренних клеток. Каждая такая клетка имеет в составе границы одну особую траекторию, являющуюся «элементом притяжения», или стоком, и одну особую траекторию, являющуюся «элементом отталкивания», или, источником. Каждая клетка имеет в составе границы одну особую траекторию, являющуюся стоком [12]. Будем считать, что для систем (13) и (14) существует взаимно однозначное и взаимно непрерывное временное преобразование области G_R в себя. Тогда качественная структура разбиения этих временных областей на траектории определяется особыми траекториями, направлением движения по замкнутым траекториям и характером устойчивости предельных особых траекторий [12]. Это утверждение означает, что предельный цикл, отвечающий клеточному ядру, определяет управляющую структуру растительного организма. В то же время существуют и другие траектории, отвечающие другим растительным структурам (другим структурам клетки). К теоретическим результатам можно отнести так же и следующее. Временная структура по типу седла сопоставляется с хроматином клетки. Тогда, на основании [12], сепаратрисы (усы седла) в нормальном состоянии клетки не идут из седла в седло, рис. 2

а б в

Это означает, что клетка не делится и хроматин находится в свободном состоянии. Как только (с движением биологического времени) усы седла начинают идти из седла в седло происходит разрыв грубости системы - движение хроматина (образуется веретено), клетка делится, рис. 3.

dt=Δt=1/5000 – шаг разрешения (15 мин), весь период времени от 0 до 5000 соответствует периоду времени 0.01≤Т≤1 или Т_н≤Т≤Т_цв (53-е суток); О=θ – угол поворота вектора времени.

Рисунок 4 – Фазовый портрет динамической системы биологического временного поля (13) для двух случаев: А – уровень напряженности факторов внешней среды U_sr=0.9; В – уровень напряженности факторов внешней среды U_sr=0.5.

траекторий биологического временного поля. Интегральные траектории заполняют всю область пространства G_R. Интегральные траектории биологического временного поля имеют всевозможные направления и пронизывают всю область определения, которая соответствуют целостности растительного организма. Константы интегрирования с₁ и с₂ соответствуют культуре подсолнечника и составляют значения: с₁=0.65 отн. ед. биологического времени и с₂=0.28 отн. ед. биологического времени. Эти константы могут быть записаны в виде: и , где N₁ и N₂ числа, которые дают в сумме N₁+N₂ число генов растения, то есть констант биологического времени.

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

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

Литература:

1. 
   Шиголев А.А. Температура как количественный агрометеорологический показатель скорости развития растений и некоторых элементов их продуктивности. // Труды ЦИП, 1957, -Вып.53.-С.75-81.
   
2. 
   Шиголев А.А. Методика составления фенологических прогнозов. В кн.: Сборник методических указаний по анализу и оценке сложившихся и ожидаемых агрометеорологических условий. - Ленинград: Гидрометеоиздат,1957.-С.5-18.
   
3. 
   Уланова Е.С. агрометеорологические условия и урожайность озимой пшеницы. – Ленинград: Гидрометеоиздат, -1975. – 302 С.
   
4. 
   Дмитренко В.П. О совместном учете влияния факторов внешней среды на продолжительность межфазных периодов растений. // Труды УкрНИГМИ, - 1987, Вып. 223, С. 3 – 23.
   
5. 
   Наумов М.М. Рост растений и биологическое время. // Вісник ОДЕКУ, 2005, Вип. 1, С. 72-78.
   
6. 
   Наумов М.М. Рост растений и биологическое время. Результаты расчетов. // Вісник ОДЕКУ, 2006, Вип. 2, С. 101-107.
   
7. 
   Наумов М.М. Онтогенетические колебания биологического времени растений: яровые и озимые формы. // Український гідрометеорологічний журнал. – 2011. № 8, С. 106 – 115.
   
8. 
   Наумов М.М. о существовании грубой динамической системы биологического времени у растений на примере культуры подсолнечника. // Метеорологія, кліматологія та гідрологія. – 2001, Вип. 43, С. 136 – 146.
   
9. 
   Наумов М.М. Векторный характер биологического времени растений. // Метеорологія, кліматологія та гідрологія. – 2004, Вип. 48, С. 226 – 234.
   
10. 
    Наумов М.М. Векторный характер биологического времени растений. Циркуляция биологического времени. // Метеорологія, кліматологія та гідрологія. – 2005. Вип. 49, С. 328 – 339.
    
11. 
    Наумов М.М. Ротор биологического временного поля растений и электромагнитная теория Максвелла. // Український гідрометеорологічний журнал. 2009. № 5. – С. 147 – 157.
    
12. 
    Понтрягин Л.С. Избранные труды. М., «Наука», 1988, т. 2, С. 30-34.
    
13. 
    Баутин Н.Н., Леонтович Е.А. Методы и приемы качественного исследования динамических систем на плоскости. М., «Наука», 1990, -486 С.
    
14. 
    Гродзинский Д.М. Надежность растительных систем. К., «Наукова думка», 1983, -367 С.
    
15. 
    Атабекова А.И., Устинова Е.И. Цитология растений. М., ВО «Агропромиздат», 1987, -246 С.
    
16. 
    Ризниченко Г.Ю. Лекции по математическим моделям в биологии. М-Ижевск, Изд. РХД, 2002, -236 С.