Некоторые характерные физические свойства почв различного гранулометрического состава (наиболее вероятный диапазон – в скобках)

* ВРК определены по характерным основным гидрофизическим характеристикам (ОГХ) на основании метода А.Д.Воронина

** Природные пески почти всегда слоисты. Вследствие этого приведенные данные весьма ориентировочны

В процессе полевых исследований проводятся полевые измерения физических свойств почв и отбор образцов для лабораторных исследований, рассчитываются агрофизические показатели и создаются сводные таблицы изученных свойств, включающие усредненные (или медианные) значения для каждой точки исследований по всем почвенным разностям. Полевые исследования рекомендуется проводить с помощью экспресс-методов, которые позволяют быстро и эффективно оценить физические свойства почв, а соответственно агрофизическую ситуацию почвенного покрова в ландшафте (17, 160, 216). Обязательными, прежде всего, являются послойные (через каждые 10 см до глубины 40-50 см) определения плотности почвы буровым методом с объемом бура не менее 100 см³, водопроницаемости методом трубок и сопротивления пенетрации (твердости). Такой набор быстро определяемых агрофизических свойств почвы дает возможность впоследствии выделить участки полей с повышенной уплотненностью (как с поверхности, так и в подпахотном слое - так называемое «подпочвенное уплотнение»), зоны возможного застоя и быстрого (инфлюкционного) движения влаги. Указанные параметры, в том числе влажность почвы – это минимальный список необходимых агрофизических свойств. Этот набор агрофизических свойств должен дополняться лабораторными определениями, прежде всего, гранулометрического состава, порозности агрегатов, структурного состава и, в случае необходимости, других агрофизических характеристик.

Полученные результаты заносятся в таблицы (например, в формате xls), которые служат основой для создания, импортируемой в ГИС, базы данных.

На основе полученной базы данных послойно по каждому параметру агрофизического состояния с помощью той или иной предусмотренной в ГИС процедуры интерполяции (кригинг, сплайн и т.п.) строятся изоплеты исследованных свойств почв в рамках заданного участка исследований. Послойное отображение свойств с одной стороны дает полную картину изменчивости физических параметров по площади исследуемого поля, а с другой стороны, с помощью изоплет позволяет выделить участки (зоны) с соответствующими градациями агрофизических свойств. В результате прорисовки изоплет в ГИС’е получается карта исследованного свойства, на которой выделяются градации с неблагоприятными и оптимальными свойствами.

При агрофизических пространственных исследованиях часто бывает весьма затруднительно в достаточном количестве для прорисовки изоплет определить те или иные агрофизические свойства для всего массива, например, диапазон доступной влаги. Это важнейшее агрофизическое свойство предполагает независимое экспериментальное определение наименьшей влагоемкости (НВ) и влажности завядания (ВЗ). Определение указанных агрофизических свойств – дело трудоемкое, а в случае пространственных площадных определений с многочисленными повторениями – практически невозможное. Поэтому предложен подход пространственного восстановления подобных трудноопределяемых свойств. Метод на основе так называемых «педотрансферных функций» (ПТФ), которые позволяют с помощью статистических (регрессионных) зависимостей находить искомое свойство по ряду других с ним взаимосвязанных. Физики почв любят говорить: «Педотрансферные функции превращают данные, которые у нас есть, в данные, которые нам нужны!».

Итак, задача формулируется следующим образом: создан банк данных для территории по определяемым в полевых условиях свойствам почв (плотность, сопротивление пенетрации, водопроницаемость и др.). Необходимо с помощью ГИС построить карту неизвестного, трудноопределяемого параметра например, запасов влаги при диапазоне доступной влаги в корнеобитаемой толще.

Решение этой задачи можно разбить на несколько этапов:

1. Необходимо выбрать 8-12 точек на исследуемом ландшафте по возможности с наиболее контрастными физическими свойствами. Если на участке есть почвы, заметно различающиеся по гранулометрическому составу (на несколько градаций: глина – суглинок - песок), то лучше эти группы почв выделить в самостоятельные выборки. Основой для составления такого набора точек должна быть почвенная карта. В зависимости от площади исследования и его сложности (комплексности почвенного покрова) на участке располагается 8-10 ключевых разрезов, в которых исследуются основные физические свойства почв.

2. В этих точках необходимо определить как традиционные для почвенно-физического обследования свойства (указанные выше плотность, сопротивление пенетрации, водопроницаемость и возможные иные экспрессно определяемые свойства), а также и искомые свойства (НВ и ВЗ). В результате получается массив из 8-12 пространственно сопряженных свойств: НВ, ВЗ, плотности, сопротивления пенетрации, водопроницаемости соответствующих слоев почвы. Этот массив помещается в табличном виде в среду EXCEL или STATISTICA.

3. С помощью процедуры множественного линейного регрессионного анализа находятся уравнения связи между традиционными свойствами (плотность, сопротивление пенетрации, водопроницаемость) и трудноопределяемыми необходимыми агрофизическими свойствами (в примере НВ, ВЗ). В результате получаются уравнения связи, в которых в виде функции выступает трудноопределяемое свойство, а в виде аргументов – традиционные агрофизические свойства. С помощью этих уравнений можно рассчитать, или, как говорят «восстановить» значение искомого свойства по известным. Это уравнение и носит название педотрансферной функции (ПТФ).

4. По полученным педотрансферным функциям рассчитывают трудноопределяемые свойства для всего массива по известным традиционным свойствам. Получается пространственно распределенный массив данных по трудноопределяемому свойству, с которым поступают по описанной выше процедуре: помещают в ГИС, рисуют карты изоплет и выделяют на картах участки с неблагоприятными/оптимальными диапазонами рассматриваемых свойств.

В качестве примера получения и анализа агрофизической информации в ГИС-технологиях рассмотрим построение карты запасов продуктивной влаги на участке комплекса серых лесных почв Владимирского ополья (опытное поле Владимирского НИИСХ). Для этого на исследуемом поле по равномерной сетке с шагом 21 м было заложено 62 разреза глубиной до 60 см. В каждом разрезе на глубинах 0, 10, 20, 30, 40 и 60 см определялись плотность почвы, сопротивление пенетрации, водопроницаемость, а также в лаборатории - плотность твердой фазы, содержание органического углерода. Кроме того, на контрастных по морфологическому описанию почвенных разностях, но относящихся к суглинистым почвам (от легкого до тяжелого суглинков) было выделено 10 точек, в которых определялась также величина НВ методом залива малых площадей. По этим 10 наблюдательным разрезам на основе процедуры множественного регрессионного анализа были получены уравнения, связывающие величину НВ с другими свойствами почв (плотностью почвы, водопроницаемостью, плотностью твердой фазы, содержанием углерода). Оказалось, что уравнение связи величины НВ (в долях от абсолютно сухого веса почвы) для исследуемой территории выглядит следующим образом (по данным В.Г. Тымбаева):

, (2.37)

где ρ_b – плотность почвы, ρ_s – плотность твердой фазы (г/см³), c - содержание углерода (%).

Это уравнение имело множественный коэффициент корреляции 0.58, и высокую его достоверность по критерию Фишера (с уровнем значимости <0.05). Кроме статистического исследования достоверности коэффициентов и значимости это уравнение необходимо исследовать и экспертно, на основании известных фактов и взаимосвязей в физике почв. Действительно, и увеличение органического вещества в почве, и плотность твердой фазы увеличивают влагосодержание при НВ, а вот увеличение плотности эту величину снижает. Это отражено и знаках приведенного уравнения. Следовательно, оно удовлетворяет и известным взаимосвязям физических свойств почв и данное уравнение использовано как ПТФ для определения величины НВ и в оставшихся 42 точках исследуемого массива. Например, если в какой-либо точке получены значения плотности почвы 1.4 г/см³, содержания углерода 1.2% и 2.65 г/см³ – плотности твердой фазы, то в этой точке величина НВ составит 0.073-0.046*1.4 + 0.025*1.2 + 0.059*2.65 = 0.195 или 19.5 % к весу почвы. Полученный массив данных для исследуемой площади позволяет привлечь аппарат ГИС для прорисовки изоплет НВ, а при соответствующем определении ВЗ и диапазона доступной влаги для всех исследуемых слоев. На рис. 2.10 приведен пример изоплет диапазона доступной влаги для слоев 10-15 и 30-35 см рассмотренного сельскохозяйственного опытного поля. На нем четко выделяются районы с высоким диапазоном доступной влаги, а также с низкими его значениями.

Рис. 2.10. Изоплеты пространственного изменения варьирования запасов влаги, соответствующих диапазону доступной влаги (мм водн.слоя):

а) в слое 10-15 см, б) в слое 30-35 см (участок ВНИИСХ). Области с изоплетой >25 мм – оптимальный диапазон, <21 мм – неудовлетвортельный. (По данным В.Г.Тымбаева)

На этой карте пространственного распределения диапазона доступной влаги выделяются зоны наиболее вероятного иссушения почв, т.е. зоны с низкими запасами влаги при диапазоне доступной влаги (<21 мм водного слоя), а также зоны с наиболее высокими запасами (>25 мм).

Рассмотрим пример исследования пространственного распределения плотности почвы на участке светло-каштанового солонцового комплекса (Волгоградская область, Ергенинская возвышенность,). На территории в 200 га (2х1 км) было заложено 12 ключей на основных элементах рельефа: на плакорном участке со светло-каштановой почвой (3 разреза), на склонах со светлокаштановой солонцеватой почвой (4 разрезов), солонцами на склонах балки (3 разреза) и солончаками по днищу балки (2 разреза). В этих точках были изучены послойно плотность почвы, водопроницаемость, сопротивленияе пенетрации в естественном состоянии и после определения водопроницаемости при влажности равной НВ. Были получены соответствующие ПТФ, позволяющие по сопротивлению пенетрации восстанавливать плотность почвы. Затем по сетке с шагом в 200 м были определены сопротивления пенетрации и рассчитана плотность почвы.

Рис. 2.11. Картограмма медианных значений плотности светло-каштановой почвы в классификации Бондарева (1985)

На картограмме (рис. 2.11) отражены медианные значения плотности почвы. На ней отчетливо выделяются зоны сильного переуплотнения (>1.5 г/см³). Эти зоны в основном охватывают часть склона вблизи реки и вблизи дорог, что связано с антропогенным воздействием и более сильно выраженными процессами осолонцевания. Причем переуплотнение идет по всему профилю практически с поверхности. Также выделяются зоны, связанные с эрозионными процессами – смытые почвы в нижней части балки имеют высокую плотность, а в более пологой ложбине отмечается намытость почв и соответственно более низкая плотность. Для остальной территории плотность светло-каштановых почв составляет 1.3-1.4 г/см3.