Модифицированный метод прогнозирования эрозии почвы и ее последствий

УДК 631.6.02

Модифицированный метод прогнозирования эрозии почвы и ее последствий

Ю.П. СУХАНОВСКИЙ, доктор сельскохозяйственных наук, зав. лабораторией (e-mail: soil-er@kursknet.ru) А.В. ПРУЩИК, кандидат сельскохозяйственных наук, старший научный сотрудник С.И. САНЖАРОВА, кандидат биологических наук, ведущий научный сотрудник Ю.А. СОЛОВЬЕВА, кандидат географических наук, старший научный сотрудник Всероссийский научно-исследовательский институт земледелия и защиты почв от эрозии, ул. Карла Маркса, 70Б, Курск, 305021, Российская Федерация

Цель исследований - модификация разработанного ранее метода прогнозирования эрозии почвы и ее последствий. Метод прогнозирования основан на использовании комплексной модели, которая включает динамику эрозии почвы, трансформацию органического вещества и почвообразование (образование гумусового слоя). Проведен анализ математической модели почвообразования применительно к целинным и пахотным черноземам Центрального Черноземья. В результате был обоснован другой выбор значений параметров, уменьшающий ее погрешность, который был использован в модифицированной комплексной модели. Для проведения расчетов разработана специальная компьютерная программа (на языке программирования Visual Basic Microsoft Excel 2010). С помощью модифицированных моделей были разработаны прогнозы для четырех сценариев использования черноземов Курской области: типичный в XIX-XX вв.; почвозащитный; консервация (переход к многолетним травам); целина (залежь, переходящая к условиям целины). Используя только модифицированную модель почвообразования, были получены следующие результаты: начальная скорость почвообразования будет в интервале 0 - 0,16 мм/год; продолжительность восстановления эродированных почв будет в интервале 5,6 - 8,7 тыс. лет. Модификация комплексной модели практически не повлияла на прогноз трендов мощности гумусового слоя и запасов гумуса. Отклонение от предыдущих прогнозов для техже самых условий не превысило 4,6%. Это меньше погрешности исходной комплексной модели, равной 15%. Модификация метода прогнозирования эрозии почвы и ее последствий практически только уменьшила погрешность модели почвообразования. Оценить саму погрешность невозможно, поскольку отсутствуют необходимые для

этого данные.

Ключевые слова: эрозия почвы, последствия, методы прогнозирования, почвообразование, восстановление почвы, землепользование, чернозем.

Для цитирования: Модифицированный метод прогнозирования эрозии почвы и ее последствий / Ю.П. Сухановский, А.В. Прущик, С.И. Санжарова, Ю.А. Соловьева // Земледелие. 2016. № 2. С. 29-32.

Из всех видов деградации почвы главной причиной сокращения площади пашни служит эрозия [1, 2]. Для поиска решений о рациональном использовании оставшихся почвенных ресурсов требуются методы прогнозирования, в частности, эрозии почвы и ее последствий. К числу главных из них относится проблема восстановления почв. Разработанные долгосрочные прогнозы [3, 4] для разных сценариев использования эродируемых черноземов привели к выводам: на современном уровнеземледелия почвенные ресурсы будут продолжать сокращаться, но можно существенно замедлить этот процесс; для восстановления запасов гумуса в эродированных почвах потребуются столетия, а для восстановления мощности гумусового слоя почвы - тысячелетия.

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

Сейчас в разных странах темп почвообразования для условий целины и пашни оценивается, как правило, в интервале сотых и десятых долей миллиметра в год [6-8]. Оценить погрешность темпа почвообразования невозможно, поскольку отсутствуют необходимые для этого данные. Средняя (за 200 лет после распашки целины) скорость сокращения гумусового слоя эродированных черноземов Курской области составляет 1,0-2,5 мм/год [9], что намного больше интенсивности почвообразования. При таком землепользовании сокращение гумусового слоя было обусловлено практически только эрозионными потерями почвы, то

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

Прогнозирование и поиск рациональных вариантов землепользования основывается на математических моделях, описывающих динамику рассматриваемых процессов. Для черноземов Центрального Черноземья разработана единственная комплексная модель [3], которая включает динамику эрозии почвы, трансформацию органического вещества и образование гумусового слоя. Проверка этой модели показала, что через 200 лет после распашки целины она занизила мощность гумусового слоя в среднем на 13% [10], а запасы гумуса в слое 0-50 см - на 15% [3]. Следовательно, ее можно использовать для долгосрочного прогнозирования эрозии черноземов и ее последствий.

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

Чтобы уменьшить погрешность оценки, мы модифицировали модель динамики мощности гумусового слоя для не подверженных эрозии черноземов Центрального Черноземья. Она входит в состав комплексной модели [3], которая используется для долгосрочного прогнозирования динамики эрозии почвы, трансформации органического вещества и образования мощности гумусового слоя для тех же территориально расположенных черноземов. Оценить погрешность модели почвообразования (образования мощности гумусового слоя) практически невозможно по причине отсутствия необходимых данных. Поэтому обоснование уменьшения погрешности проведено, используя математический анализ уравнения е динамики мощности гумусового Л слоя и сопоставления решений этого д уравнения с данными для черноземов е Курской области [9]. При разработке | прогнозов для разных сценариев землепользования необходимо за- ю давать для неэродированных черно- м земов среднемноголетние значения © поступления в почву растительных 6

остатков. Эти значения получены с использованием данных многофакторного полевого опыта ФГБНУ «ВНИИ земледелия и защиты почв от эрозии» (Медвенский район, Курская область) без внесения удобрений и представлены в работе [3]. Для прежних условий [4] разработаны новые прогнозы с целью получения более точных данных и оценки, насколько модификация модели почвообразования повлияла на прогнозы. Расчеты проведены по разработанной компьютерной программе на языке программирования Visual Basic Microsoft Excel 2010.

Далее показано, как можно уменьшить погрешность модели образования гумусового слоя, которая используется в комплексной модели [3]. Для условий целины тренд мощности гумусового слоя почвы H(t) описывается дифференциальным уравнением [6]

0],

где t - возраст почвы, год; H(t) -мощность гумусового слоя, соответствующая ее возрасту, см; X - постоянная величина, год - 1; Hmax - максимальное (предельное) значение, к которому с увеличением возраста стремится мощность гумусового слоя, см. Из уравнения (1) следует, что с увеличением возраста почвы (с увеличением Н(^) скорость образования гумусового слоя dH(t)/dt уменьшается и стремится к нулю. Образно говоря, чем моложе почва, тем с большей скоростью она образуется. В дальнейшем под почвообразованием понимается образование гумусового слоя почвы.

Решение уравнения (1) записывается

"(0 = Нтах[1-*ехр(-А0],

/с-1-Н0/Нтах,

где exp - экспонента, H0

«О

о

СЧ СМ

О

S

ф

и

ф

s

ф

со

интервале 0 < t < 8 тыс. лет. Важно то, что эти данные не имеют одинаковую календарную дату начала почвообразования.

Если в уравнении (1) под величиной t понимать календарное время и за t = 0 принять начало голоцена (10 тыс. лет назад), то уравнения (1) и (2)будут описывать тренд (изменение во времени среднего значения) мощности гумусового слоя целинной почвы после окончания ледникового периода. Строго говоря, это может быть справедливым только при условии, что для всех использованных данных измерений процесс почвообразования начинался одновременно 10 тыс. лет назад. Как отмечено выше, это не так, а получить необходимые данные измерений практически невозможно.

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

Выбор значений параметров в уравнении (2). Для современных (t = 10 тыс. лет) целинных почв рассчитанные значения Н^) для разных выборок практически одинаковые (100 и 96 см, см. табл. 1). В начале почвообразования (t = 0) для выборки 1 новообразованный гумусовый слой отсутствовал (H0 = 0), что логично для условий после окончания ледникового периода. Для выборки 2 такой слой (H0 = 12 см) уже присутствовал. Это можно объяснить следующим. В выборке 2 (в отличие от выборки 1) добавлено большое количество данных для почв с возрастом десятки и сотни лет. Однако для этих данных начало почвообразования было не 10 тыс. лет назад, то есть были условия более благоприятные для почвообразования, чем после окон-чанияледникового периода. Поэтому при аппроксимации данных измерений зависимостью (2) это привело к появлению H0 > 0 и, как следствие, к уменьшению значения X. Таким образом, для описания эволюции гуму-

сового слоя целиннои почвы в уравнениях (1) и (2) лучше использовать значения параметров, полученные по выборке 1. В дальнейших расчетах использовали эти значения.

В комплексной модели [3] использовали только уравнение (1) при X =

0,00027 год - 1, а в качестве H при' ' max ^

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

Замедление почвообразования на пашне. Принято, что перед распашкой целины состояние почвы было близко к равновесному состоянию. После распашки целины изменился растительный покров. Уравнение (1) справедливо и для пашни, если вместо H использовать H - мак-

max max, a

симальную (предельную) мощность гумусового слоя для пашни. Для почв в равновесном состоянии справедливо соотношение [6]

-]0,37

H = H

maх,а max

мощность гумусового слоя в начальный момент времени, то есть H0 = H(t =

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

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

1). Выборка 2 [7] включает в себя выборку 1 [6], дополненную новыми данными. В обеих выборках использованы данные для возраста почвы в

1. Значения параметров почвы для лесостепных черноземов

Величины am и a можно интерпретировать как среднее ежегодное поступление растительных остатков в почву, соответственно, перед и сразу после распашки целины (эти величины не зависят от времени t). После распашки уменьшилось поступление в почву растительных остатков, что привело к снижению содержания гумуса. Поэтому из уравнения (3) следует, что для пашни выполняется

неравенство Hmax, a < Hmax. С учеТОм

этого, из уравнения (1) следует, что при одинаковой мощности гумусового слоя H(t) скорость почвообразования на пашне была меньше, чем на целине.

В комплексной модели [3] величина a зависела от времени t и рассчитывалась как количество растительных остатков, поступающих в почву в момент времени t (это занижало значение H ).

max, a

Изменение мощности гумусового слоя после распашки целины. Как отмечено выше, для оценки значений параметров H , k и X использовали

max

данные для возраста почвы в интервале 0 < t < 8 тыс. лет, то есть не использовали данные о мощности гумусового слоя современных целинных почв (при t = 10 тыс. лет). Далее показано, что мощность гумусового слоя перед распашкой целины можно приближенно оценивать по мощно-

Выборка X, k H , Рассчитанные по уравнению (2) значения H(t), см

год - 1 см t = 0 (H0) t = 10 тыс. лет

Выборка 1 0,00027 1,0 107* 0 100

Выборка 2 0,00015 0,902 120 12 96

* - рассчитано по графику на рис. 1.14-2 [6]

сти современных неэродированных пахотных почв. Изменения черноземов Курской области произошли ориентировочно за 200-250 лет после распашки целины [9]. Используя уравнение (2) со значениями параметров для выборки 1 (см. табл. 1), были проведены расчеты изменения мощности гумусового слоя целины за 300 лет. В результате получено, что приращение гумусового слоя равняется 6 мм. Поскольку, как показано выше, скорость почвообразования на пашне уменьшается, то и приращение гумусового слоя было еще меньше.

Следовательно, с погрешностью не более 6 мм можно принять, что мощность гумусового слоя целинной почвы с возрастом 9,7-10 тыс. лет равняется мощности современной неэродированной пахотной почвы с возрастом 10 тыс. лет.

Уменьшение погрешности уравнения (1) для конкретной почвы. Значения параметров Н , X и к интер-

^ ^ тах' ^

претируются как средние значения для использованной выборки данных измерений. Это служит источником погрешности при рассмотрении конкретной почвы. По данным [11] для Центрального Черноземья мощность неэродированных черноземов изменялась в интервале 25-115 см (объем выборки 68), а для Курской области - в интервале 60-100 см (объем выборки 195, использованный в работе [9]) со средним значением 74 см, соответствующим возрасту почвы 10 тыс. лет.

Используя значения указанных параметров для выборки 1 (см. табл. 1), рассчитанная по уравнению (2) средняя мощность современных (при t = 10 тыс. лет) целинных (и неэроди-рованных пахотных) черноземов приближенно равняется 100 см, что для курских черноземов больше реальной на 26 см. Используя то же уравнение с теми же значениями параметров, получено, что мощности 74 см соответствует возраст почвы t = 4,4 тыс. лет, что на 5,6 тыс. лет меньше реального

возраста. Из уравнения (1) следует, что с уменьшением возраста почвы (с уменьшением H(t)) увеличивается скорость почвообразования dH(t)/dt. Следовательно, занижение возраста почвы ведет к завышению скорости почвообразования.

Другим источником погрешности служит отсутствие в использованных выборках данных измерений мощности гумусового слоя современных (при t = 10 тыс. лет) целинных и неэро-дированных пахотныхчерноземов. Из всех величин, входящих в уравнение

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

Чтобы уменьшить погрешность, предложено при рассмотрении конкретной почвы значение Hmax для целины рассчитывать по уравнению (2) при известном измеренном значении H (t = 10 тыс. лет) современной целинной или неэродированной пахотной почвы. Например, для курских черноземов (при H = 74 см) получено, что Hmax = 79 см. В таком случае средние скорости за 10 тыс. лет, рассчитанные по уравнению (1) и по измеренной мощности гумусового слоя с возрастом 10 тыс. лет, будут одинаковыми. А изменение скорости во времени будет описываться уравнением (1). Заметим, что средняя мощность гумусового слоя современных неэродиро-ванных черноземов Курской области (74 см) меньше максимальной (H

* ' у max

= 79 см) на 6%, то есть эти целинные почвы через 10 тыс. лет приблизились к их равновесному состоянию, что и было принято выше.

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

[3] приведены значения ежегодного поступления растительных остатков в неэродированную почву при современном земледелии: сценарий

типичный - 2,9; почвозащитный -5,6; консервация - 7,4 т/га в год. Урожайность сельскохозяйственных культур (при прочих равных условиях) прямо пропорциональна запасам гумуса в слое 0-50 см [9], а за период земледелия после распашки целины содержание гумуса в пахотном слое почвы уменьшилось примерно в два раза [12, 13]. Это означает, что урожайность этих культур, следовательно, и поступление растительных остатков в почву при ее равновесном состоянии (сразу после распашки целины) должно быть примерно в два раза больше. Поступление растительных остатков для сценария целина остается без изменения, то есть аш = 12,3 т/ га в год [3]. Для неэродированной почвы и сценария «типичный» в XIX-XX вв. скорость почвообразования принята равной нулю. Объясняется это тем, что имеющаяся мощность гумусового слоя (74 см) была уже больше максимально возможной Нтаха = 60 см. Следовательно, можно приближенно принять, что после распашки целины мощность гумусового слоя неэродированной почвы практически не изменилась. Это означает, что мощность гумусового слоя почвы перед распашкой целины можно принять равной мощности современной неэродированной почвы.

Для эродированных почв приведены средние скорости сокращения гумусового слоя за 200 лет (см. табл. 2). Эти скорости в 25-100 раз больше скорости почвообразования для сценария «типичный» в XIX-XX вв. Следовательно, сокращение гумусового слоя в прошлом происходило практически только за счет эрозионных потерь почвы.

Восстановление мощности гумусового слоя. Существует весьма актуальная задача: как можно восстановить эродированные почвы и что для этого потребуется? Сейчас единственный известный способ, гарантирующий восстановление -

2. Параметры почвообразования, характеризующие черноземы Курской области

Параметр Неэродированные Степень эродированности

слабая средняя I сильная

Гумусовый слой, см 74 55 35 24

Средняя скорость сокращения гумусового слоя за 200 лет, мм/год 0 1,0 2,0 2,5

Сценарий: Начальная скорость почвообразования, мм/год

типичный в XIX-XX вв. (зерно-паро-пропашной севооборот, Нтаха = 60 см) 0 0,01 0,07 0,10

почвозащитный (зерно-травяной севооборот, Нтаха = 76 см) 0,01 0,06 0,11 0,14

консервация (переход к многолетним травам, Нтаха = 85 см) 0,03 0,08 0,13 0,16

целина (залежь, переходящая к условиям целины, Н = 79 см) 0,01 0,06 0,12 0,15

Восстановление эродированных черноземов по сценарию целина

Возраст мощности гумусового слоя для условий целины, тыс. лет 10 4,4 2,2 1,3

Время на восстановление гумусового слоя для условий целины, тыс. лет 0 5,6 7,8 8,7

Ы (D 3 ь

(D д

(D Ь 5

(D

М 2

О ^

О)

это перевод пашни в залежь, которая со временем перейдет к условиям целины. Для сценария «целина» в уравнении (2) были приняты следующие значения параметров: Нтах = 79 см, X = 0,00027 год - 1, к = 1. Сначала для этих условий по уравнению (2) был рассчитан возраст мощности гумусового слоя неэродированной и эродированных почв. Время на восстановление рассчитано как разность между этими возрастами.

Ранее [4], используя значения параметров для выборки 2 (см. табл. 1), рассчитано, что для восстановления эродированных почв потребуется значительно меньше времени (2,34,9 тыс. лет). Объясняется это тем, что рассчитанный возраст неэродированной почвы (Н = 74 см) равнялся 5,7 тыс. лет, что меньше реального на 4,3 тыс. лет. Как было показано выше, занижение возраста почвы ведет к завышению скорости почвообразования, следовательно, к занижению времени на восстановление эродированных почв.

Таким образом, для уменьшения погрешности уравнений (1) и (2) для лесостепных черноземов необходимо использовать X = 0,00027 год - 1, к = 1, а значение Нтах рассчитывать по уравнению (2) для измеренной мощности гумусового слоя современной целинной или неэродированной пахотной почвы.

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

Проверка исходной модели [3] показала, что через 200 лет после распашки целины она занизила мощность гумусового слоя в среднем на 13% [10], а запасы гумуса в слое 050 см - на 15% [3]. Используя исходную комплексную модель, ранее [4] были разработаны прогнозы трендов мощности гумусового слоя и запасов гумуса в слое 0-50 см для склона длиной 500 м с углом наклона 3° и для четырех сценариев землепользования. Для тех же условий были разработаны такие же прогнозы с использованием модифицированной комплексной модели. Резуль-«о таты расчетов показали, что через о 200 лет внизу склона модифицированная модель (по сравнению с исходной) ^ завысила мощность гумусового слоя о» на 1,8-4,6%, а запасы гумуса в слое | 0-50 см - на 0,9-2,4%. То есть различия результатов были в пределах погрешай ности исходной модели. 5 Таким образом, модификация мо-$ дели почвообразования уменьшила ее

погрешность и позволила для разных сценариев землепользования оценить начальную скорость почвообразования в интервале 0-0,16 мм/год, а время на восстановление эродированных почв - в интервале 5,6-8,7 тыс. лет. Использование этой модели практически (в пределах 4,6%) не повлияло на прогнозы трендов мощности гумусового слоя и запасов гумуса. Для разных сценариев землепользования модифицированные модели могут быть использованы для оценки скорости почвообразования, продолжительности восстановления эродированных почв и для прогнозирования трендов мощности гумусового слоя и запасов гумуса.

Литература.

1. Добровольский Г.В. Педосфера как оболочка высокой концентрации и разнообразия жизни на планете // Почвы в биосфере и жизни человека: сб. научн. тр. М.: МГУЛ, 2012. С. 20-34.

2. Кирюхина З.П., Пацукевич З.В. Эрозионная деградация почвенного покрова России // Почвоведение. 2004. № 6. С.752-758.

3. Сухановский Ю.П., Санжарова С.И., Прущик А.В. Модель динамики содержания гумуса в эродированном черноземе Центрального Черноземья // Агрохимия.

2011. № 12. С. 45-52.

4. Оценка трендов эродируемых черноземов пахотных земель Центрального Черноземья / Ю.П. Сухановский,

A.В. Прущик, Ю.А. Соловьева, С.И. Санжарова // Земледелие. 2015. № 6. С.19-22.

5. Проблема обоснования допустимых эрозионных потерь почвы и подход к ее решению / Ю.П. Сухановский, А.В. Прущик, Ю.А. Соловьева, С.И. Санжарова // Бюллетень Почвенного института им.

B.В. Докучаева. 2015. Вып. 78. С. 3-19.

6. Лисецкий Ф.Н. Пространственно-временная организация агроландшафтов. Белгород: Белгор. гос. ун-т, 2000. 304 с.

7. Голеусов П.В., Лисецкий Ф.Н. Воспроизводство почв в антропогенно нарушенных ландшафтах лесостепи. М.: ГЕОС, 2009. 210 с.

8. Лисецкий Ф.Н., Светличный А.А., Черный С.Г. Современные проблемы эрозиоведения. Белгород: Константа,

2012. 456 с.

9. Сухановский Ю.П., Прущик А.В., Санжарова С.И. Оценка и динамика ресурсов черноземов и серых лесных почв сельскохозяйственных угодий Курской области // Вестник Курской государственной сельскохозяйственной академии. 2013. № 6. С. 53-59.

10. Сухановский Ю.П. Вероятностный подход к расчету эрозионных потерь почвы // Почвоведение. 2013. № 4. С. 474-481.

11 .Справочник агрогидрологических свойств почв Центрально-Черноземных областей. М.: Гидрометеорологическое изд-во,1964. 114 с.

12. Адерихин П.Г. Изменение черноземных почв ЦЧО при их использовании в сельском хозяйстве // Черноземы ЦЧО и их плодородие. М.: Наука, 1964. С.61-69.

13. Афанасьева Е.А. Черноземы СреднеРусской возвышенности. М.: Наука, 1966. 224 с.

Modified Forecasting Method for Soil Erosion and Its Consequences

Yu.P. Sukhanovskii, A.V. Prushchik, S.I. Sanzharova, Yu.A. Solov'eva

All-Russian Research Institute of Soil Management and Soil Erosion Control, ul. Karla Marksa, 70b, Kursk, 305021, Russian Federation

Summary. The aim of the research was to modify the previously developed method of forecasting of soil erosion and its consequences. The forecasting method is based on a complex model, which includes models of the dynamics of soil erosion, organic matter transformation and soil formation (humus horizon formation). It was carried out an analysis of the soil formation model applied to virgin soil and arable chernozem of the Central Chernozem Region. As a result, it was proposed another parameter choice for the model, which reduces its error. The new choice of parameter value was used in the modified complex model. It was developed a special computer program (in the programming language Visual Basic Microsoft Excel 2010) to carry out the calculations. Using the modified models the forecast was developed for four land use scenarios in Kursk region: typical one, representative of XIX-XX centuries; soil protective; conservation at transition to perennial grasses; virgin land (restoration of long-continued neglected fields). Using only the modified model of soil formation, the following results were obtained: the initial rate of soil formation will be in the range of 0-0.16 mm/year; the duration of the eroded soil restoration will be in the range of 5.6-8.7 thousand years. The modification of the complex model has almost no impact on the forecast trends of the humus horizon and humus reserves. The deviation from the previous forecasts at the same conditions did not exceed 4.6%. This is less than the error of the original complex model (15%). The modification of the forecasting method of soil erosion and its consequences in practice only reduced the error of soil formation model. Error estimation of the soil formation model is impossible because there is no necessary data.

Keywords: soil erosion, consequences, forecasting methods, soil formation, soil restoration, land utilization, chernozem.

Author Details: Yu.P. Sukhanovskii, D. Sc. (Agr.), head of laboratory (e-mail: soil-er@kursknet.ru); A.V. Prushchik, Cand. Sc. (Agr.), senior research fellow; S.I. Sanzharova, Cand. Sc. (Biol.), leading research fellow; Yu.A. Solov'eva, Cand. Sc. (Geogr.), senior research fellow.

For citation: Sukhanovskii Yu.P., PrushchikA.V., Sanzharova S.I., Solov'eva Yu.A. Modified Forecasting Method for Soil Erosion and Its Consequences. Zemledelie. 2016. No 2. Pp. 29-32(in Russ.).