Введение. Инструментарий дистанционного зондирования Земли (ДЗЗ) позволяет увеличить объемы поступающей информации и повысить оперативность ее обработки [1–3]. Для оценки состояния, динамики развития сельхозкультур и прогнозирования их продуктивности все чаще используют методы удаленного мониторинга и контроля. Один из таких показателей представляет собой нормализованный разностный вегетационный индекс (normalized difference vegetation index, NDVI) [4–6].
При анализе сезонной динамики NDVI, проведенном рядом авторов [7, 8], выявлены более высокие значения данного показателя в летние месяцы, что объясняется увеличением биомассы и количественного состава хлорофилла в вегетирующих растениях. При этом понижается яркость растений в красной зоне и возрастает в инфракрасной области спектра [9], увеличивая расчетное значение NDVI.
Высокая степень взаимосвязи вегетационного индекса и количества накопленного хлорофилла в вегетирующих растениях подтверждается результатами исследований [10–12]. Особо подчеркивается наличие корреляционной зависимости между интегрированными максимальными значениями NDVI и их продуктивностью [13–17]. Изучением динамики изменения вегетационного индекса на посевах озимой пшеницы, возделываемой по интенсивной технологии, удалось показать, что наибольшие значения NDVI соответствуют моменту начала фазы колошения, достигая пика при 0,80–0,88 [18].
В то же время характер изменения сезонных значений вегетационного индекса и его отражение агробиологических показателей растений существенно зависят от внешних воздействий, особенностей возделываемых культур, региональных условий и т. д. [7, 19, 20]. Это не дает возможности предложить единый подход к оценке условий развития культуры и применению единой прогностической модели.
Кроме того, наращивание разработки новых средств защиты и стимулирования растений предполагает соответствующие решения в части оценки их эффективности на основе применения инновационных технологий, в т. ч. средствами ДЗЗ.
Таким образом, предлагаемая работа, направленная на моделирование во времени урожайности сельскохозяйственных культур под воздействием новых систем защиты (препарат «Нигор++») с применением данных NDVI является актуальной и современной.
Теоретическая значимость исследования состоит в определении роли динамики изменения показателей вегетативного индекса NDVI при оценке эффективности новых систем защиты сельскохозяйственных культур, а также получении прогностических моделей урожайности культур на основе полиномиальных (второй степени) функций.
Практическая ценность работы заключается в получении корреляционных зависимостей между урожайностью и пиковыми значениями вегетативного индекса по данным многолетних наблюдений, необходимых для формализации прогноза урожайности; использовании полученных математических моделей при достоверном планировании производственно-экономических задач конкретного хозяйства.
В опытно-производственном хозяйстве Орловского ГАУ в течение ряда лет проводится последовательная работа по внедрению информационно-цифровых технологий в производство [21–23]. В 2021 г. университетом выполнялась комплексная НИР [24], одна из содержательных частей которой заключалась в получении экспериментальных данных результатов исследования эффективности нового биологически активного средства «Нигор++».
Цель исследования – оценка эффективности применения нового биопестицида и разработка прогностической модели урожайности зерновых культур на основе данных ДЗЗ.
Задачи: получение, обработка и анализ ретроспективных данных и текущих сезонных значений NDVI по яровому ячменю и озимой пшенице; сравнительная оценка действия нового биологического средства защиты растений (СЗР) защитно-стимулирующего действия «Нигор++»; численная оценка многолетнего корреляционного отношения максимальных значений NDVI и урожайности культур и вывод о возможности использования вегетационного индекса как показателя урожайности культур с разработкой адекватного уравнения регрессии в условиях опытно-производственного хозяйства Орловского ГАУ.
Объекты и методы. Полевые исследования сезона 2021 г. выполнялись на опытно-производственных участках Орловского ГАУ. Тип почвы: серая лесная почва. Культуры – озимая пшеница Московская 39 (площадь посевов 48,1 га), яровой ячмень Раушан (площадь посевов – 17,4 га). Посев озимой пшеницы провели 10.09.2020, ярового ячменя – 07.05.2021 обработанными фунгицидным протравителем «Скарлет» семенами.
Ранее проведенными исследованиями показано [25–27], что биологический препарат «Нигор++» (Патент РФ № 2463759) увеличивает ассимиляционную поверхность, общую биомассу, корневую емкость, накопление зеленых пигментов, устойчивость к болезням, урожайность зерна за счет увеличения сохранности растений на делянке, размера колоса, массы семян, а также значительно улучшает качество зерна (количество протеина, клейковины и седиментации).
На опытных посевах озимой пшеницы препарат «Нигор++» вносили 25 мая 2021 г. в фазу «начала выхода в трубку» и повторно 16 июня 2021 г. в фазу «колошения» по флаговому листу. Посевы ярового ячменя обрабатывались 10 июня и 28 июня 2021 г.
Ретроспектива данных NDVI обеспечивалась архивными спутниковыми снимками за 2016–2020 гг. Для получения многолетних статистических данных по динамике изменения вегетационного индекса NDVI использовали архивы геопортала «КосмосАгро» [28]. Приложение ScanEx Geomixer [29] обеспечивало оперативную аналитическую обработку поступающей информации.
В течение всего вегетативного сезона 2021 г. выполнялась аэрофотосъемка и расчет индексов NDVI опытных участков; строилась сезонная динамика вегетационного индекса. Для получения текущих данных использовался беспилотный летательный аппарат (БПЛА) «Агрофлай Квадро».
Индекс вегетации NDVI определяли по формуле от первоисточника [4]; оценку точности моделей – показателем средней абсолютной ошибки MAPE (Mean Absolute Percentage Error) [30].
Результаты и их обсуждение. Исследованиями в полевых условиях сезона 2021 г. установлены следующие условия, влияющие на сроки развития растений и активность бактерий биостимулирующего препарата (яровой ячмень/ озимая пшеница):
- Фаза «3 лист/кущение». Сумма осадков за период – 14/35 мм; средняя температура воздуха – 20,2/10,5; средний уровень ультрафиолета – 16,7/10,1 Вт · м-2.
- Фаза «Выход в трубку» (продолжительность – 13/31 дней). Сумма осадков за период – 0/0,45 мм; средняя температура воздуха – 31,9/22,3; средний уровень ультрафиолета – 27,0/25,2 Вт · м-2.
- Фаза «Колошение» (продолжительность – 9/11 дней). Сумма осадков за период – 3,9/0 мм; средняя температура воздуха – 26,6/25,4; средний уровень ультрафиолета – 19,2/28,9 Вт · м-2.
- Фаза «Цветение» (продолжительность – 15/11 дней). Сумма осадков за период – 0,5/2,1 мм; средняя температура воздуха – 35,1/28,7; средний уровень ультрафиолета – 28,2/26,8 Вт · м-2.
- Фаза «Созревание» (продолжительность – 23/11 дней). Сумма осадков за период – 0,6/5,5 мм; средняя температура воздуха – 27,6/31,3; средний уровень ультрафиолета – 13,7/28,8 Вт · м-2.
В условиях Орловской области оптимальное снабжение яровых зерновых влагой считается, когда в период «посев-всходы» выпадает 30–45 мм осадков. Результаты полевых наблюдений показали, что в период с 7 по 23 мая выпало около 25 мм осадков. Это повлекло некоторое увеличение срока всходов (11–12 дней), притом, что по многолетним данным всходы появляются на восьмой-девятый день.
В период «всходы-цветение» выпало 95 мм осадков, что является достаточным водоснабжением для обеспечения нормального развития посевов яровых и озимых.
В фазы развития «цветение-созревание» зерновых выпадение осадков зафиксировано суммарно всего 0,6 и 5,5 мм – для ярового ячменя и озимой пшеницы соответственно. Тем не менее, показатели запасов почвенной влаги в корнеобитаемом слое составили 67–73 мм, что является оптимальным для обеспечения формирования урожая зерновых культур.
В целом температурный режим и солнечная активность сезона 2021 г. мало отличались от среднемноголетних показателей по Орловской области и не оказали аномального влияния на процесс развития растений и активность бактерий биостимулирующего препарата.
Сезонная изменчивость индекса NDVI по данным многолетних наблюдений демонстрирует следующее. Цикличный тренд значений временных рядов культур аналогичен и в большей степени соответствует закону нормального распределения. Меньшие значения индекса вегетации, не зависимо от культуры, характеризуют зимний период. Конец весны – середина лета обусловлены максимальными значениями NDVI. Наибольшие среднемноголетние значения индекса вегетации отмечены в июне. Для озимой пшеницы максимальные значения NDVI приходятся на июнь и составляют 0,72; для ярового ячменя – 0,56 [23].
С помощью сезонных снимков с БПЛА получены временные ряды значений индекса NDVI для опытных участков озимой пшеницы и ярового ячменя, обработанных новым биопрепаратом, и контроля. Тренд изменения значений вегетационного индекса по фазам развития растений приведен на примере озимой пшеницы урожая 2021 г. (рис. 1). Как видно, значения NDVI возрастают от ранних сроков вегетации культуры к фазам колошения и цветения. Максимальное значение NDVI, равное 0,78, отмечено 18 июня, что соответствует фазе начала колошения, зафиксированной 17 июня. Данный характер динамики изменения вегетационного индекса по отношению к зафиксированным в исследованиях срокам фенофаз растений можно считать типовым для обеих культур.
Установлено, что развитие вегетации растений опытного участка с озимой пшеницей на 10–15 дней опережает контроль в период с 15.04 вплоть до начала фазы колошения (13–15.06). Пиковые значения вегетационного индекса NDVI для опытного и контрольного участков приходятся на 18.06. При этом максимальное значение NDVI для участка, обработанного новым биопрепаратом, составляет 0,83, что превышает пиковое значение контрольного участка (0,78) на 5–6 % (рис. 2). Аналогичные сравнительные данные были получены для посевов ярового ячменя.
Для решения об использовании пиковых значений вегетационного индекса NDVI в качестве управляемой переменной x функции У = f(x) выполнен анализ статистических данных урожайности обеих культур и соответствующих им пиковых значений NDVI в ретроспективе с 2016 по 2020 г. (табл. 1).
Рис. 1. Сезонная динамика вегетационного индекса NDVI по фазам
развития озимой пшеницы (2021 г.) (использованы материалы [8])
Рис. 2. Сравнительная оценка значений и сезонной динамики NDVI озимой пшеницы
по результатам использования нового биопрепарата
Коэффициентами корреляции по данным признакам 0,79 и 0,75 – для пшеницы озимой и ярового ячменя, соответственно, с высокой степенью надежности подтверждается возможность прогноза урожайности культур на основе данных их пиковых значений NDVI.
Используя данные таблицы 1 получены следующие прогностические уравнения урожайности культур:
– для пшеницы озимой
У = 104,09–297,13x + 264,91x2; (1)
– для ячменя ярового
У = 38,9–80,39x + 71,22x2, (2)
где x – максимальное значение NDVI в текущем сезоне; У – ожидаемая урожайность культуры, ц∙га-1.
Таблица 1
Численные значения урожайности и максимумов индекса NDVI
культур с расчетом коэффициента корреляции по данным отдельных
производственных участков в 2016–2020 гг.
|
Год наблюдений |
Производственный участок |
NDVI |
Урожайность, ц·га-1 |
|
Пшеница озимая |
|||
|
2016 |
25 |
0,79 |
27,46 |
|
2017 |
23 |
0,82 |
45,50 |
|
26 |
0,89 |
46,32 |
|
|
31 |
0,85 |
44,91 |
|
|
2018 |
36 |
0,73 |
30,76 |
|
2019 |
22 |
0,66 |
24,00 |
|
33 |
0,77 |
26,71 |
|
|
2020 |
23 |
0,81 |
40,47 |
|
24 |
0,80 |
31,53 |
|
|
26 |
0,80 |
44,97 |
|
|
2021 |
38 |
0,78 |
32,10 |
|
|
k = 0,79 |
||
|
Ячмень яровой |
|||
|
2016 |
27 |
0,84 |
21,22 |
|
30 |
0,81 |
20,50 |
|
|
2017 |
54 |
0,87 |
23,21 |
|
2018 |
37 |
0,54 |
19,73 |
|
38 |
0,50 |
14,02 |
|
|
39 |
0,61 |
15,20 |
|
|
2019 |
27 |
0,57 |
18,16 |
|
34 |
0,57 |
18,01 |
|
|
2020 |
13 |
0,63 |
14,26 |
|
29 |
0,71 |
18,27 |
|
|
|
k = 0,75 |
||
Коэффициент детерминации R2 и MAPE составили соответственно 0,74 и 10,3 – для модели урожайности пшеницы озимой; 0,66 и 8,4 – для модели ячменя ярового.
С целью оценки степени достоверности прогностических моделей урожайности, полученных в ходе настоящего исследования, в них произведена подстановка пиковых значений NDVI сезона 2021 г. Важно отметить, что указанные расчеты выполнялись непосредственно в ходе формирования графика динамики изменения NDVI текущего сезона, т. е. до уборки урожая.
Подстановка в уравнения (1) и (2) значений NDVI, полученных на сравниваемых участках опытных культур, показала следующее. Расчетная продуктивность озимой пшеницы на участке, обработанном препаратом «Нигор++», составила 40,0 ц · га-1. Прогноз урожайности контрольного участка находился на уровне 33,5 ц · га-1, что на 16,2 % меньше. Для ярового ячменя спрогнозирована урожайность 28,9 и 28,4 ц · га-1 – соответственно для участков опыта и контроля.
Позднее были получены значения фактической урожайности культур, принимавших участие в эксперименте. Так, продуктивность опытных участков озимой пшеницы и ярового ячменя находилась на уровне 41,3 и 31,3 ц га-1, тогда как урожайность по контролю составила для этих культур 36,1 и 28,7 ц · га-1 соответственно.
Полученные расчетные и фактические значения урожайности культур для различных вариантов эксперимента представлены в таблице 2. Как видно, погрешность выполненной на основе предложенных регрессионных моделей прогнозной оценки урожайности культур не превышает 6,7 %. Это позволяет рекомендовать предложенные математические модели для прогноза урожайности культур в условиях опытно-производственного хозяйства Орловского ГАУ. Такой прогноз позволит определять возможные объемы реализации растениеводческой продукции и рынки сбыта в «удачном» году, а также планировать превентивные мероприятия для минимизации ущерба, например, при неблагоприятно складывающихся погодных условиях.
Таблица 2
Расчетные и фактические значения урожайности культур
в различных вариантах полевого эксперимента
|
Сравниваемые участки |
Пиковое значение NDVI |
Расчетная урожайность, ц · га-1 |
Фактическая урожайность, ц ∙ га-1 |
Погрешность оценки, % |
|
Озимая пшеница |
||||
|
Опыт |
0,83 |
40,0 |
41,3 |
3,2 |
|
0,78 |
33,5 |
36,1 |
6,2 |
|
|
Яровой ячмень |
||||
|
Опыт |
0,76 |
28,9 |
31,3 |
6,7 |
|
Контроль |
0,74 |
28,4 |
28,7 |
1,04 |
Заключение
1. Сравнением временных рядов NDVI, полученных по среднемноголетним и текущим данным сезона 2021 г., обнаружено опережение динамики вегетативного развития культур в 2021 г. на 7–10 дней и незначительное снижение максимального (пикового) значения вегетационного индекса.
2. Установлено, что, независимо от культуры, развитие вегетации растений, обработанных новым биологическим СЗР защитно-стимулирующего действия «Нигор++», опережает контроль вплоть до начала фазы колошения. Характеризующее наибольший объем вегетативной массы максимальное значение индекса NDVI для культур, обработанных новым биопрепаратом, составляет 0,83, что превышает пиковое значение контрольного участка на 5–6 %.
3. Коэффициенты корреляции между максимальными сезонными значениями индексов NDVI и продуктивностью культур составили 0,79 и 0,75 для пшеницы озимой и ярового ячменя соответственно. Это подтверждает надежность прогноза урожайности культур, выполняемого на основе данных их пиковых сезонных значений нормализованного разностного вегетационного индекса.
4. Получены прогностические модели урожайности культур на основе полиномиальных (второй степени) функций. Коэффициенты детерминации и средние абсолютные ошибки моделей составляют 0,74 и 10,3; 0,6 и 8,4 соответственно для озимой пшеницы и ярового ячменя, что подтверждает достоверность предложенных регрессионных моделей и возможность их применения для планирования производственно-экономических задач в опытно-производственном хозяйстве университета.
5. Сравнение прогнозных и фактических значений урожайности культур в полевом сезоне 2021 г. выявило погрешность модельной оценки, не превышающую 6,7 %.
6. Теоретический подход к использованию показателей сезонной динамики индекса NDVI для сравнительного анализа эффективности новых систем защиты может применяться в научных и производственных организациях, реализующих принципы цифрового сельского хозяйства.
