Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

MATHEMATICAL MODELING OF THE RELATIONSHIP OF CLIMATIC FACTORS AND POPULATION DYNAMICS CEPHUS PYGMAEUS L. AND COLLYRIA COXATOR VILL. IN DIFFERENT PHASES OF ONTOGENESIS OF WINTER WHEAT IN THE STAVROPOL TERRITORY

Glazunova N.N. 1
1 FSBEI HPE «Stavropol State Agrarian University»
Eight year old body of empirical data on the population dynamics of an ordinary summer wheat stem sawfly and its entomophage calls in different phases of development of winter wheat (booting, earing, flowering, milk ripeness, waxy ripeness) analyzed using statistical methods (regression and correlation analyzes) and using basic computer program " Usage -6", which allowed to mathematically describe the populations of insects , weather and climatic conditions in the Stavropol region. The quantitative relationship between the size of populations and its phytophage yaytseparazitom in the form of plots and empirical equations in different phases of development of winter wheat during the period from 1998-2005. Shows the combined effect of climatic factors on the development of insect populations in the form of surface equations linking populations with temperature and amount of precipitation in different periods of ontogenesis of winter wheat. According to the results of mathematical modeling of the proposed changes the dynamics of population ordinary stem sawfly and its entomophage calls in different periods of ontogenesis of winter wheat and climatic factors in the form of empirical equations of second order for the Stavropol Territory.
and mathematical modeling
phase of development of winter wheat
precipitation
temperature
Collier
Common sawflies

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

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

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

В расчетах были приняты следующие показатели: длительность фазы, , дни; средняя точка продолжительности фазы, дата месяца; средняя температура фазы, °С (х) и средняя величина осадков фазы, мм (у).

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

Полученные данные были обработаны с использованием базовой компьютерной программы «Статистика-6». Результаты хода изменений средних фазовых температур и объема осадков в различные периоды онтогенеза озимой пшеницы обыкновенного хлебного пилильщика и коллирии в 1998-2005 годы.

С помощью программы «Статистика 6» мы рассчитали и аналитические зависимости функции температуры «t = f (τ)» для исследуемого временного интервала для каждого года начиная с 1998 по 2005 гг. включительно.

Рассчитанные программой уравнения хода изменений средних фазовых температур имеют следующий вид:

t, °С(1998) = 17,58 - 0,3904τ + 0,0197τ2 – 0,0002τ3 + 6Е-07τ4, (1)

t, °С(1999) = 13,718 – 0,364τ + 0,0202τ2 – 0,0002τ3 + 6Е-07τ4, (2)

t, °С(2000) = 9,3214 + 0,3583τ - 0,0125τ2 + 0,0003τ3 – 2Е-06τ4, (3)

t, °С(2001) = 13,42 + 0,1911τ - 0,0064τ2 + 0,0002τ3 – -1Е-06τ4, (4)

t, °С(2002) = 9,2325 + 0,6506τ - 0,0262τ2 + 0,0005τ3 – 3Е-06τ4, (5)

t, °С(2003) = 10,25 + 1,1537τ - 0,0504τ2 + 0,0008τ3 – 4Е-06τ4, (6)

t, °С(2004) = 13,558 + 0,029τ - 0,0007τ2 + 3Е-05τ3 – 4Е-07τ4, (7)

t, °С(2005) = 11,167 + 0,7529τ - 0,0318τ2 + 0,0005τ3 – 3Е-06τ4, (8)

Также были получены и аналитические зависимости функция осадков «Vос = f (τ)» для 5 фаз онтогенеза озимой пшеницы за исследуемый временной интервал для каждого года.

Рассчитанные программой уравнения изменения количества осадков в эти же периоды имеют следующий вид:

Vос, мм(1998) = 17,581 – 0,3904τ + 0,0197τ2 – 0,0002τ3 + 6Е-07τ4 (9)

Vос, мм(1999) = 173,03 – 10,061τ + 0,4215τ2 – 0,0072τ3 + 4Е-05τ4 (10)

Vос, мм(2000) = 235,39 – 17,51τ + 0,7605τ2 – 0,012τ3 + 6Е-05τ4 (11)

Vос, мм(2001) = 185,36 – 3,8552τ + 0,13τ2 – 0,0027τ3 + 2Е-05τ4 (12)

Vос, мм(2002) = 9,2325 + 0,6506τ – 0,0262τ2 + 0,0005τ3 – 3Е-06τ4 (13)

Vос, мм(2003) = -48,035 + 6,5934τ –0,2714τ2 + 0,0044τ3 – 2Е-05τ4 (14)

Vос, мм(2004) = 109,25 + 16,32τ –0,317τ2 + 0,0018τ3 (15)

Vос, мм(2005) = 48,91 + 4,9541τ – 0,2368τ2 + 0,0036τ3 – 2Е-05τ4 (16)

Для установления количественной взаимосвязи численности популяций обыкновенного хлебного пилильщика и коллирии во время лета имаго в разные годы (за 1998-2005 годы) исследований озимой пшеницы в виде эмпирических уравнений. Был проведен анализ их численности с использованием базовой компьютерной программы «Статистика-6» для четырех фаз развития озимой пшеницы: трубкование, колошение, цветение, молочная спелость.

В расчетах были приняты следующие показатели: длительность фазы, , дни; численность популяции, экз/м2 (Рi). Взаимосвязь между численностью популяций обыкновенного хлебного пилильщика и его энтомофагом была установлена в виде аналитических зависимостей функции

Рi = f (τ), (17)

где Рi – численность популяции, экз/м2;

τ – длительность фаз, дни.

Расчеты значений Рi были выполнены по ежегодным эмпирическим данным, полученным за восемь лет исследований (1998-2005 гг.) по динамике численности обыкновенного хлебного пилильщика и коллирии.

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

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

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

Рисунок 1 – Динамика лета обыкновенного хлебного пилильщика в различные фазы онтогенеза за 1998-2005 гг.

Рисунок 2 – Динамика лета коллирии в различные фазы онтогенеза за 1998-2005 г.

Проведенный анализ показал, что аналитические зависимости численности популяции обыкновенного хлебного пилильщика и коллирии в различные фазы онтогенеза за период 1998-2005 года имеют вид уравнений линейной регрессии третьей степени Рк = а+вτ+сτ2+dτ3. Данные уравнения описываю лёт фитофага и его энтомофага, который в Центральном предкавказье начинается в конце фазы трубкования озимой пшеницы в первые особи фитофага нами были отмечены с 10.05. и продолжается до фазы начала молочной спелости последние особи попадались до середины июня 15.06.

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

Полученные аналитические зависимости изменения численности лета обыкновенного хлебного пилильщика и коллирии в различные периоды онтогенеза озимой пшеницы за 1998-2005 гг. были разбиты по исследуемым фазам (трубкование, колошение, цветение, молочная спелость) и проанализированы по годам с использованием базовой компьютерной программы «Статистика-6».

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

у = а + вх + сх2 + dх3 + fx4, (18)

где у – численность популяции, Р, экз/м2;

х – среднее значение фазы развития, τ, дни.

Принимая во внимание тот факт, что численность популяций насекомых и паразитов зависит от двух факторов: температуры, t, °С и количества осадков, Vос, мм, то эти зависимости можно представить в виде поверхности, определяемой уравнениями типа

Р, экз/м2 = а + в*t + с*Vос (20)

Р, экз/м2 = а - в*t + с*Vос+d*t2-e*t*Vос - f*Vос2 , (21)

На рисунках 3-6 представлены графические и аналитические зависимости «Р = f(t; Vос)» в виде поверхности второго порядка для обыкновенного хлебного пилильщика и коллирии в период колошения и цветения за 1998-2005 годы, полученные согласно базовой компьютерной программе «Статистика-6».

Изменение динамики лета обыкновенного хлебного пилильщика и коллирии в фазу трубкования в зависимости от температуры и количества осадков происходят по следующим уравнениям:

Рп= 4,1648+0,0567*t -0,1437*Voc -0,0009*t2-0,0008*t*Voc +1,8017Е-5*Voc2 (22)

Рк= 1,2616+0,0751*t +0,0034*Voc -0,0022*t2-0,0004*t*Voc +1,043Е-5*Voc2 (23)

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

Рп= -128,9742-11,9895*t -0,2229*Voc +0,2726*t2+0,0136*t*Voc -0,0002*Voc2 (24)

Рк=77,3991-7,2005*t -1,1094*Voc +0,1643*t2+0,0071*t*Voc -0,0002*Voc2 (25)

Рисунок 3 – Изменение динамики лета обыкновенного хлебного пилильщика в период колошения в зависимости от температуры и количества осадков (экз/м2).

Рисунок 4 – Изменение динамики лета коллирии в период колошения в зависимости от температуры и количества осадков (экз/м2).

Рисунок 5 – Изменение динамики лета обыкновенного хлебного пилильщика в период цветения в зависимости от температуры и количества осадков (экз/м2).

Рисунок 6 – Изменение динамики лета коллирии в период цветения в зависимости от температуры и количества осадков (экз/м2).

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

Рецензенты:

Ченикалова Е.В., д.б.н., профессор кафедры химии и защиты растений ФГБОУ ВПО «Ставропольский государственный аграрный университет», г. Ставрополь.

Стародубцева Г.П., д.с.-х.н., профессор, заведующая кафедрой физики ФГБОУ ВПО «Ставропольский государственный аграрный университет», г. Ставрополь.