Почвоведение, 2022, № 6, стр. 713-727

ВВЕДЕНИЕ

Изучение гидрологического режима длительно осушаемых почв в агроландшафтах является актуальным направлением, которое имеет фундаментальную и прикладную стороны. Целенаправленное изменение водного режима путем осушительной мелиорации влечет определенные изменения в функционировании почв [13, 14, 26]. Однако режимным исследованиям длительно осушенных почв уделяется недостаточное внимание. В то же время при разработке современных агротехнологий необходимо знание конкретных гидрологических особенностей почв. Для этого требуется проведение мониторинговых исследований такой длительности, которая вскроет закономерности динамики влажности в годы с разным количеством осадков на фоне современных климатических изменений [25]. Результаты научных работ должны иметь, в том числе количественное выражение для последующего построения моделей гидрологического режима с определенным варьированием границ показателей.

Важным итогом мониторинга гидрологического режима почв является возможность диагностики степени гидроморфизма по морфологическим признакам в профиле. С генетической точки зрения – это возможность воплощения неодокучаевской формулы: факторы–процессы–свойства [9] с дополнением: факторы–режимы–процессы–свойства. Ведь признаки гидроморфизма являются результатом процессов, протекающих в определенных условиях водного режима. А режим складывается в результате сочетания факторов почвообразования (как природных, так и антропогенных).

Чаще всего исследования направлены на выявление связей между факторами (климатом, рельефом, породами, растительностью) и водным режимом [5, 7, 13, 16, 19, 20, 22]. Давно ведется поиск связи между морфологическим проявлением гидроморфизма и развитием процесса оглеения [15, 27, 29]. Существует набор аналитических характеристик гидроморфизма [13, 17, 28].

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

Основной причиной оглеения является переувлажнение в анаэробных условиях. Переувлажнение наступает, когда объемная влажность становится больше предельной полевой влагоемкости (ППВ). Однако различия почв по гранулометрическому составу приводят к неточности использования величины ППВ в качестве единой критической границы переувлажнения в песчаных, супесчаных, глинистых и хорошо оструктуренных почвах. Поэтому необходимо найти показатель, который отражает наступление экологического переувлажнения, отрицательно влияющего на урожай сельскохозяйственных культур. Для решения этой задачи Зайдельман [16] использовал границу воздухоносной пористости (ВП) – 8% в пахотном слое и 6% в подпахотном. Шеин [24] в качестве критической предлагает границу ВП 10%.

Актуальным направлением является выявление количественных индикаторов гидроморфизма при изучении режима влажности почв (послойной динамики влажности). Т.А. Романовой на основании многолетних исследований и обобщения массива данных сети гидрометеостанций предложено в качестве критерия гидроморфизма почв Беларуси считать количество дней с влажностью почвы больше наименьшей влагоемкости (НВ) (состояние переувлажнения) в слое 0–20 см [22].

Этот показатель достоверно коррелирует с интенсивностью оглеения и уровнем грунтовых вод. Поэтому Романова [11] полагает, что предложенный критерий может быть основанием для дифференциации почв по степени гидроморфизма. Однако остается открытым вопрос: насколько достоверен показатель при переходе, например, от Белорусской почвенной провинции к Прибалтийской?

В Калининградской области осушительные мелиорации сельскохозяйственных земель начаты в XVII–XVIII вв. [23]. В настоящее время 596 тыс. га (около 80% от общей площади сельскохозяйственных угодий) осушается разными способами. При этом гидрологический режим почв изучен крайне слабо, особенно на фоне климатических изменений [2].

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

ОБЪЕКТЫ И МЕТОДЫ

Исследования проводили на производственном пахотном поле (ключевой участок “Перелески”) в Зеленоградском районе Калининградской области, в пределах Самбийской холмисто-моренной равнины (юго-восточная Прибалтика). Площадь поля 100 га. Участок исследования по ландшафтным условиям является типичным для Самбийской равнины [8]. Рельеф представляет собой чередование асимметричных моренных и камовых холмов и понижений между ними (рис. 1). Почвенный покров поля контрастный. На вершинах холмов и склонах сформировались буроземы разной степени оглеения (рис. 2). Почвообразующими породами для них являются валунные моренные и водно-ледниковые глубоковыщелоченные супеси и суглинки с прослойками глин [12].

Рис. 1.

Расположение и топографическая карта ключевого участка с площадками почвенно-гидрологического мониторинга в ареалах буроземов (1–8 – номера разрезов).

Рис. 2.

Профили буроземов на площадках мониторинга.

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

Все оглеенные буроземы и дерново-глеевые почвы осушаются системой закрытого гончарного дренажа со сбросом вод в открытые каналы. Последнюю реконструкцию дренажа провели в 1975 г. Глубина заложения дрен 80–110 см, междренное расстояние в среднем 12 м. С конца 90-х гг. ХХ в. прекращен текущий ремонт и обслуживание дренажных систем. Поэтому в настоящее время устьевые оголовки коллекторов в откосах каналов частично засыпаны. Это блокирует около половины объема стока из закрытых осушителей в магистральный канал.

До 2012 г. участок исследований распахивался. С 2012 г. по июль 2016 г. поле временно выведено из севооборота в связи с прокладкой газопровода. В августе 2016 г. почвы снова начали распахивать и по 2020 г. чередовали посев озимой пшеницы, озимого рапса и ярового ячменя.

Строение почв изучали в разрезах глубиной 1.2–2.4 м. Морфологическое описание и диагностику проводили по классификации почв СССР (1977 г.) [17]. Для оглеенных горизонтов использовали систему индексов Зайдельмана [15]. В глееватых горизонтах полностью или частично сохраняется цветовая гамма генетического горизонта или породы. Интенсивность проявления глееобразования диагностируется по площади “холодной” глеевой окраски: в слабоглееватых горизонтах до 20%, в среднеглееватых – 20–50% кутаны сизовато-коричневые, в сильноглееватых – 50–80% и кутаны коричневато-сизые или сизые [14]. Под почвенным гидроморфизмом в настоящей работе понимается комплекс морфологических признаков, предложенных Зайдельманом [15].

Площадки многолетнего мониторинга расположены на четырех вершинах холмов и на четырех склонах. В основе мониторинга лежит методика изучения водного режима, предложенная Роде [21]. Дополнительно исследовали динамику воздухоносной пористости.

Режим влажности изучали буровым методом. Периодичность бурения 2 раза в месяц. Отбор образцов выполняли послойно каждые 10 см до глубины 1 м. Период ежегодных непрерывных наблюдений: с апреля по ноябрь. Мониторинг проводили с 2012 по 2020 гг. Бурение приурочено к междренным участкам. Для установления уровня верховодки применяли бурение до 2 м (1–2 раза в месяц круглогодично в 2012, 2017–2020 гг.).

Данные о среднесуточной температуре и количестве осадков брали с ближайшей к полю метеостанции г. Калининграда.

Свойства почв изучали стандартными и рекомендованными методиками: гранулометрический состав пипет-методом с подготовкой почвы путем обработки 4%-ной Na₄P₂O₇, плотность твердой фазы (ρ_s, г/см³) пикнометрически, плотность сложения (ρ_b, г/см³) методом режущих колец (цилиндров) объемом 100 см³, общую пористость (ε, %) и пористость аэрации (ε_air, %) – расчетными методами, максимальную гигроскопическая влажность и наименьшую влагоемкость (НВ) – по Николаеву, влажность завядания расчетным методом, влажность разрыва капиллярной связи (ВРК) принята для песков и супесей 0.6 НВ, для суглинков 0.7 НВ, полевая влажность почв – термостатно-весовым методом (объемная влажность получали в результате умножения полевой влажности на плотность сложения), полная влагоемкость принята равной величине общей пористости [6, 14]. Воздухоносную пористость (ВП, %) рассчитывали как разницу между общей пористостью и фактической объемной влажностью. Анализы выполняли в 4-кратной повторности. Для статистической и графической обработки данных использовали программу Microsoft Excel.

РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ

Калининградская область относится к Прибалтийской провинции дерново-подзолистых и болотно-подзолистых почв южной тайги [11]. Завалишин и Надеждин доказали, что в западной части региона (Самбийская равнина) преобладают буроземные почвы [12] в условиях мягкого приморского климата. Подробные авторские исследования в начале ХХI в. подтвердили это [2].

Среднегодовая температура воздуха за период 1891–1958 гг. составляла +7.1°С [1]. Для начала XXI в. характерен повышенный температурный фон (около +8°С) [4]. Среднемноголетняя сумма осадков для Калининграда (Кенигсберга) за период 1891–1958 гг. равна 710 мм [1], а по обновленным данным (Г.М. Баринова) – 781 мм [4]. Весь период мониторинга можно разделить на три группы по количеству осадков: сухие годы – 2014, 2015, 2018 (осадков меньше среднемноголетней нормы), средние – 2013, 2019, 2020 (осадков на уровне среднемноголетней нормы), влажные – 2012, 2016, 2017 гг. (осадков больше среднемноголетней нормы) (табл. 1).

Три из четырех буроземов в автономных позициях рельефа имеют четкие признаки оглеения в форме пятен и прожилок сизовато-серого цвета (табл. 2). Даже неоглеенный в летний период бурозем (ареал 1) в сырые периоды поздней осени приобретает серовато-сизоватые прожилки в иллювиальных горизонтах (с глубины 90 см). В других почвах граница появления признаков оглеения варьирует от 29 до 50 см. Иллювиальные горизонты и почвообразующая порода являются глееватыми от слабой до сильной степени (по терминологии Ф.Р. Зайдельмана). Авторский опыт обследований почв моренных равнин и публикации по региону [2, 10] показывают, что такая картина является типичной для агроландшафтов Калининградской области. Завалишин и Надеждин [12] отмечали, что в области проще найти оглеенную почву, чем неоглеенную. Таким образом, в автономных позициях рельефа на Самбийской равнине при норме осадков 781 мм преобладают глееватые почвы.

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

Буроземы образовались на валунных моренных и перемытых суглинках. Как отмечает Зайдельман [15], камни в мелкоземе резко уменьшают объем фильтрующих пор, что приводит к уменьшению коэффициента фильтрации. В результате анализа физических свойств выяснилось, что заиленные горизонты и горизонты с высокой плотностью сложения отличаются низкой общей пористостью. Поэтому уже при влажности, равной НВ, пористость аэрации в них близка к 10% или меньше.

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

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

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

Признаки оподзоливания присутствуют в единичных почвах. Особенностью склоновых буроземов является отсутствие горизонта слабого оглеения. В большинстве почв после горизонтов А1 и В1 следует среднеглееватый пятнистый горизонт, в котором площадь оглеенных участков более 20%. Ортштейны встречаются как в гумусовом, так и подгумусовых горизонтах.

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

Общей чертой всех буроземов на поле является ритмичность динамики влажности. После зимнего и ранневесеннего периода сквозного промачивания в засушливом (как правило) апреле начинается просыхание верхней части профиля. Уже в конце мая возникают зоны биологического иссушения (влажность меньше ВРК), которые в июне распространяются до глубины 40–70 см в зависимости от количества осадков в конкретном году. В октябре (в сырые годы в сентябре) влажность в метровой толще буроземов начинает увеличиваться. В ноябре чаще всего значения влажности близки к ранневесенним величинам. Цикл заканчивается, и снова наступает период сквозного зимнего промачивания.

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

На первом этапе попытались выявить связь между количеством осадков по каждому году мониторинга и индивидуальным количественными показателями режима влажности и ВП на восьми ареалах почв методом корреляционного анализа. Связь оказалась прямолинейной во всех случаях. Устойчивая значимая связь за период 2012–2020 гг. наблюдается с показателями: влажность >НВ в слое 0–20 см и влажность почв >НВ по всему слою 1 м. Однако коэффициент детерминации сильно варьирует от 0.35 до 0.77. Это позволяет сделать вывод о том, что переувлажнение пахотного слоя и метровой толщи буроземов только на 35–77% зависит от количества осадков. К неучтенным факторам, влияющим на влажность буроземов на вершинах холмов, относятся термический режим, положение почвы в рельефе, водопроницаемость, водоудерживающая способность почв. Выяснилось, что коэффициент корреляции сильно варьирует в зависимости от временного отрезка. При рассмотрении четырехлетнего периода (2017–2020 гг.) коэффициенты корреляции в большинстве случаев значительно увеличиваются. Одна из причин заключается в том, что этот отрезок включал резко контрастные по увлажнению годы (сырой 2017 и сухой 2018). В девятилетнем периоде (2012–2020 гг.) связи ослабляются из-за увеличения разнообразия погодных условий.

На втором этапе объединили почвы по элементам рельефа (вершины и склоны), а все годы мониторинга разделили на три группы по степени увлажненности (табл. 4).

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

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

На третьем этапе разделили буроземы на вершинах холмов по степеням оглеения. Между неоглеенной и сильноглееватой почвами различия достоверны по частным показателям в отдельные группы лет: влажность <ВРК в слое 0–20 см в средние годы, влажность <ВРК в слое 0–100 см в сухие годы, влажность >НВ по всему слою 1 м в средние и сырые годы. И только один показатель период ВП <10% в слое 0–100 см статистически значимо отражает различия между буроземами разной степени оглеения вне зависимости от влажности года (табл. 5).

На четвертом этапе сравнили показатели режима влажности и воздухоносной пористости оглеенных буроземов на разных элементах рельефа (вершины и склоны) в годы с разным количеством осадков. Оказалось, что элемент рельефа достоверно влияет на длительность существования верховодки в метровой толще почвы (в глееватых буроземах на склонах она больше). Критерием отличия буроземов по степени оглеения является среднестатистический период с ВП <10%.

Тип водного режима (по А.А. Роде) изученных буроземов промывной дренажный, подтип атмосферного питания. Глееватые почвы выделяются в класс периодического капиллярного насыщения с почвенной верховодкой.

Круглогодичные наблюдения за верховодкой позволили установить, что это явление характерно для всех глееватых почв на склонах и сильноглееватого бурозема на вершине холма (рис. 2, 3). В профиле неоглеенного бурозема верховодки не было с 2012 по 2020 гг. В среднеглееватых буроземах на вершинах холмов верховодка образовывалась только в ареале 2 в экстремально сырые месяцы осени 2017 г. в иллювиальных горизонтах на короткий срок (до 10 дней).

Рис. 3.

Динамика верховодки в сильноглееватом буроземе на вершине холма.

Причина образования верховодки на вершине холма – скопление влаги над уплотненными заиленными слоями, глинистыми прослойками при низких коэффициентах фильтрации пород и затрудненном боковом стоке. На склонах верховодка формируется из натечных вод с повышений. Основная зона локализации верховодки находится в пределах глубин 80–130 см. Толщина слоя верховодки не превышает 30 см.

Верховодка начинает формироваться в глееватых буроземах в конце октября–ноябре в условиях быстрого снижения среднесуточных температур воздуха и как следствие, резкого сокращения испарения и транспирации. Источником формирования верховодки является влага атмосферных осадков и туманов. В зимний период верховодка присутствует постоянно в условиях неглубокого промерзания почв (до 20 см в морозные периоды) или отсутствия такового при положительных температурах (зимы 2014–2015 и 2019–2020 гг.). Наименьшие среднестатистические глубины верховодки характерны для февраля, когда наблюдается пик стока. В конкретных ареалах почв могут встречаться отклонения (с разбросом максимума стока от января до апреля в отдельные годы) в связи с локальными особенностями накопления свободной влаги.

Рис. 4.

Динамика верховодки в сильноглееватых буроземах на склонах: а – ареал 5, в – ареал 6, с – ареал 7, d – ареал 8.

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

Вначале вегетационного периода (апрель) верховодка нередко присутствует в слое 1 м или находится на глубине 1.0–1.2 м. Испарение и транспирация способствуют расходованию влаги, и уже к июню верховодка исчезает или перетекает в прослойки рыхлых супесей глубже 130 см. В летний период в большинстве лет верховодка отсутствует в профиле глееватых буроземов (до 130 см). При этом тонкие слои свободной гравитационной влаги могут обнаруживаться в глубоких горизонтах (140–180 см) почвообразующих пород вплоть до августа. На основании этого заключили, что капиллярное насыщение профиля почвы происходит практически круглогодично. Это является причиной критических значений воздухоносной пористости в нижней части профиля буроземов.

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

Результаты корреляционного анализа показали, что граница оглеения в почвах совпадает с верхней границей зоны с ВП <10%. Коэффициент корреляции для буроземов на вершинах холмов 0.97, а на склонах 0.99. Также тесная связь выявлена между границей оглеения и длительностью периода с ВП <10% (коэффициенты корреляции соответственно –0.84 и –0.92).

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

ЗАКЛЮЧЕНИЕ

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

Предлагается характеризовать гидрологический режим почв в вегетационный период с апреля по октябрь конкретными количественными показателями: период с влажностью >НВ и <ВРК в слоях 0–20 и 0–100 см, период сквозного промачивания (влажность >НВ по всему слою 1 м), длительность периода и глубина зоны с воздухоносной пористостью <10% в слое 0–100 см, вероятность формирования и длительность существования верховодки в слое 0–100 см.

В автоморфных легкосуглинистых и супесчаных буроземах на вершинах холмов признаки оглеения не наблюдаются (или выражены слабо в сырые сезоны поздней осени и зимы) при отсутствии в годичном цикле верховодки и возникновении зон с воздухоносной пористостью меньше 10% на глубине 70–100 см не более 20 дней в период с апреля по октябрь во влажные по осадкам годы. Для среднеглееватых осушенных супесчаных и легкосуглинистых буроземов на вершинах холмов характерно развитие зон с критической воздухоносной пористостью с глубины 50–63 см продолжительностью 50–65 дней в течение вегетационного периода. Верховодка возникает эпизодически при выпадении экстремального количества осадков.

Сильноглееватые осушенные буроземы на вершинах холмов и склонах характеризуются постоянным присутствием верховодки в профиле в осенний, зимний и ранневесенний периоды. Вероятность поднятия верховодки в слой 0–100 см с апреля по октябрь в почвах на склонах составляет от 33% в сухие по осадкам годы до 75% во влажные. Длительность нахождения верховодки в метровой толще от 5 до 23 дней. Зоны с критической воздухоносной пористостью глубже 40–60 см существуют в почвах в среднем 127–175 дней в течение вегетационного периода.

На основании многолетнего мониторинга (2012–2020 гг.) установлено, что показатель период с влажностью >НВ в слое 0–20 см в условиях гумидного климата не имеет устойчивой достоверной связи со степенью оглеения в осушаемых буроземах пахотного агроландшафта. Критерием различия буроземов разных степеней оглеения (от неоглеенных до сильноглееватых) служит показатель длительность периода с воздухоносной пористостью <10% в слое 0–100 см.

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

ФИНАНСИРОВАНИЕ РАБОТЫ

Исследования проводили в рамках инициативно-поисковой госбюджетной научно-исследовательской темы кафедры агропочвоведения и агроэкологии ФГБОУ ВО “КГТУ” “Почвенные ресурсы Калининградской области: оценка, использование, продуктивность, управление”. Работа выполнена при финансовой поддержке гранта РФФИ 19-29-05277 мк “Цифровой структурно-функциональной анализ ландшафта в системе адаптивно-ландшафтного земледелия”.