На сегодняшний день одной из наиболее актуальных задач оценки экологического состояния природных ландшафтов является необходимость разработки методов экологической оценки продуктивности ландшафтов, которые должны включать частные оценки продуктивности его составляющих, лежащие в основе комплексной или интегральной характеристики климатических, почвенных и других факторов, оказывающих влияние на продуктивность природных систем. Целью научного исследования было создание компьютерной модели оценки влияния эрозионных процессов на продуктивность почвенного покрова. При этом экологическая оценка продуктивности ландшафтов должна, как правило, иметь региональную привязанность [5], т.е. опираться на закономерности, характеризующие функционирование определенных техноприродных систем.
Математическое моделирование экологической оценки продуктивности ландшафтов. Продуктивность растений может быть определена по методике, предложенной И.С.Шатиловым и А.Ф.Чудновским [8], в основе которой лежит оценка природных ресурсов почвенно-климатической зоны по ряду факторов, а именно: обеспеченности теплом и влагой, фотосинтетической активности солнечной радиации, уровня плодородия почв и др.
Продуктивность биотической составляющей ландшафта (ПУ) характеризуется уровнем энергетических ресурсов природной системы (R), коэффициентом использования свободной энергии (hэн) и коэффициентом влагообеспеченности территории (hв). Значение экологической продуктивности биотической составляющей ландшафта с учетом его влагообеспеченности определяется из соотношения [8]:
где 
– радиационный «индекс сухости», С – калорийность единицы урожая органического вещества; hэн – коэффициент использования свободной энергии
,
где 
– коэффициент использования растениями активной фотосинтетической радиации.
Энергия, затрачиваемая на почвообразование, определяется из соотношения [2]
,
где Qi – энергия, затрачиваемая на почвообразование, кДж/см2; A – показатель, учитывающий состояние поверхности почвы (коэффициент эродированности почвы), R - уровень энергетических ресурсов природной системы, 
– радиационный «индекс сухости».
Для расчета показателя эродированности почвы с учетом рельефа в модели была проведена классификация исследуемой территории на однородные участки по форме и площади водосбора. Такая классификация позволяет учитывать долю антропогенной нагрузки на ландшафт, по-разному отражающейся на эрозионных процессах и зависящая от различных факторов, протекающих на исследуемой территории. Так, для горных территорий характерны ландшафты с преобладанием пастбищных угодий, в меньшей степени затронутых человеческой деятельностью. Для предгорных территорий, с большими территориями позволяющими их использовать под агроландшафты, доля антропогенной нагрузки на ландшафт сильно завышена, и факторы, влияющие на эрозионно–аккумулятивные процессы ранжированы иначе. В связи с этим в модели предлагается использовать различные модели.
Для агроландшафтов предгорной территории республики используется методика картометрических вычислений по крупномасштабным картам.
Схема картографического моделирования следующая:
1) составляются факторные карты по всем показателям и коэффициентам уравнения;
2) путем совмещения различных тематических карт выполняется районирование территории (выделение ареалов по однородности эрозионных процессов);
3) по каждому ареалу рассчитывается смыв и процент земель различной категории эрозионной опасности;
4) окончательные расчеты представляются в виде карты районирования территории по степени эрозионной опасности.
При моделировании процесса эрозии почвы для расчета показателя ее эродированности (А) нами была построена функция, зависящая от трех параметров:
где A – показатель эродированности почвы; 
– фактор влияния климата условной территории на проявление эрозии; 
– фактор влияния рельефа на проявление эрозии; 
– фактор влияния растительного покрова на предотвращение эрозии.
Из климатических факторов прямое воздействие на интенсивность эрозионных процессов оказывают осадки, вызывающие сток. Для оценки эрозионной опасности обильных дождевых осадков используется эрозионный индекс осадков Уишмейера-Смита:
где 
– максимальная интенсивность дождя за период 30 мин.; 
– кинетическая энергия дождя. Годовой эрозионный индекс осадков вычисляется суммированием эрозионных индексов отдельных дождей.
Фактор рельефа рассчитывается на основе уравнения Уишмейера —Смита:
LS = Lq·Wk,
где L - длина склона. W - крутизна склона, q и к -эмпирические коэффициенты, принимающие соответственно значения 0.3-9.5 и 1.4-1.5 Эрозионные процессы более интенсивны на крутых и протяженных склонах из-за возрастания потенциальной энергии текущей воды.
Для определения опасности развития эрозии в зависимости от состава возделываемых культур, вычисляется средневзвешенное значение проективного покрытия почвы культурами (Рср.взв., %):
где 
– проективное покрытие почвы данной культурой в эрозионноопасные периоды, в %; 
– площадь, занимаемая данной культурой, в % от общей площади севооборота.
По полученному значению параметра А можно определить степень эродированности почвы рассматриваемого района в соответствии со шкалой, представленной в таблице:
|
Значение А |
Степень эродированности почвы |
|
0..0,1 |
Не подвержена |
|
0,11..0,3 |
Слабо эродированная |
|
0,31..0,6 |
Средне эродированная |
|
0,61..0,9 |
Сильно эродированная |
|
0,91..1 |
Деградированная |
Для определения параметра, характеризующего экологическую продуктивность почвы, необходимо найти величину потенциально возможной энергии, затраченной на почвообразование - QП [2]
.
Таким образом, экологическая оценка продуктивности ландшафтов (КЭ) определяется соотношением таких осредненных индикаторных величин, как коэффициент продуктивности растений (КР) и почвы (КП)
где 
– коэффициент, характеризующий экологическую оценку продуктивности растительного сообщества
,
– коэффициент, характеризующий экологическую продуктивность почвы
Целью нашего вычислительного эксперимента было определение оценки влияния эрозионных процессов на продуктивность почвенного покрова. Оценивались различные уровни эродированности почвенного покрова.
В качестве испытательного полигона рассматривалась территория Баксанского района, в частности сельхозугодия, расположенные вблизи с. Исламей. Разработан программный продукт, реализующий предложенный алгоритм оценки влияния эрозионных процессов на продуктивность почвенного покрова агроландшафтов рассматриваемой территории. На рисунках 1 и 2 приведены результаты работы программы.
Рис.1. График зависимости продуктивности почвы от состояния почвы
Рис.2. Графики зависимости энергии, затраченной на почвообразование Qi и QП от состояния почвы
Разработанная компьютерная модель позволяет оценить влияние эрозионных процессов на продуктивность почвенного покрова и может быть использована в исследовательских целях.
Рецензенты:
Анахаев К.Н., д.т.н., профессор, Заслуженный деятель науки КБР, зам. директора по селевой проблематике ФГБУ «Высокогорный геофизический институт», г. Нальчик;
Ламердонов З.Г., д.т.н., профессор, Кабардино-Балкарский аграрный университет, г. Нальчик.