Heterogeneity of pore space properties at the pedon scale of the phaeozems humus horizon

Автор: Fomin D.S., Yudina A.V., Romanenko K.A.

Журнал: Бюллетень Почвенного института им. В.В. Докучаева @byulleten-esoil

Рубрика: Статьи

Статья в выпуске: 117, 2023 года.

Бесплатный доступ

The variability of soil parameters depends on the chosen method of measurements, the genesis and type of soil land use, and the level of hierarchy of soil structure organization. Computed tomography of soils is an actively growing method of soil structure study, for which many methodological issues remain relevant. The aim of this work was to examine the variability of the main parameters of the pore space (total porosity, number and average pore size) by the example of the humic horizon of a Phaeozem soil. For this purpose, an excessive number (15 microcores of 2 × 3 cm volume) was sampled from soil profile. Based on statistical evaluation of parameter variation, the objective was to determine the optimal number of replicates allowing full characterization of the soil pore space microstructure at the pedon scale. The smallest difference in heterogeneity between pedon and representative elementary volume REV was observed for total porosity (~12 times), while this ratio is larger (~14 times) for number and average pore size. On average, the threshold level, at which the dispersion of properties stopped decreasing, was 7.3 ± 0.6 monoliths for total porosity, 6.5 ± 0.6 monoliths for pore number, and 7.5 ± 0.4 monoliths for LT. Thus, minimal number of replicates necessary for full characteristic of soil structure is 7 monoliths. Sampling and analysis of microcores in 3 repetitions allows to describe the heterogeneity of the structure of the upper pedon horizon only by 25-30%.

Еще

Soil pore space, soil structure, soil heterogeneity, computed tomography, soil sampling

Короткий адрес: https://sciup.org/143181144

IDR: 143181144   |   DOI: 10.19047/0136-1694-2023-117-208-231

Текст научной статьи Heterogeneity of pore space properties at the pedon scale of the phaeozems humus horizon

Вариабельность почвенных свойств зависит от метода измерений, генезиса и типа землепользования почвы (Дмитриев, 1999; Урусевская и др., 2000) , а также от уровня иерархии структурной организации почвы (Yudina, Kuzyakov, 2019; Yudina et al., 2018) . Вариабельность почвенных свойств как в пространстве ( пространственная ), так и во времени ( временная ) имеет разные закономерности для разных масштабных уровней. Оба вида вариабельности свойств почв необходимо учитывать при оценке качества почв и, соответственно, при разработке методов измерения, критериев оценки (Vogel et al., 2019) .

Оценка вариабельности параметров почвы проводится на основе дескриптивной статистики с помощью статистических характеристик, обозначающих меры местоположения (центрированности), разброса (рассеивания), формы распределения случайной величины (Дмитриев, 1999; Самсонова, 2007) . Наиболее часто используемыми в почвенных исследованиях являются: среднее ( m ), стандартное отклонение ( standard deviation , SD ) и коэффициент вариации ( coefficient of variation , CV , % ) (Самсонова, 2007; Pen nock et al., 2008) . Для оценки степени вариабельности свойств почв было предложено использовать следующие градации CV : от

0 до 15% – низкая вариабельность, 16–35 – средняя, и > 36 – высокая (Wilding, 1985) . Считается, что такие физические свойства как плотность почвы, общая пористость, содержание ила и песка характеризуются низкой–средней, низкой, средней–высокой и низкой–средней вариабельностью соответственно (диапазоны CV равны 3–26, 7–11, 16–53 и 3–37% соответственно) (Mulla, McBrat- ney, 2001) . Однако в последнее время все больше внимания уделяется динамичности физических параметров и разработке на их основе индикаторов функций почв (Rabot et al., 2018; Yudina, Kuzya- kov, 2019) .

Представление о вариабельности почвенных свойств рассматривается через призму конкретного свойства и методики измерения (Дмитриев, 1999) . Классический отбор средней пробы дает ограниченный набор информации прежде всего потому, что нарушается целостность почвенной структуры. Исследования почвенных свойств в образцах, имеющих такое же строение, что и в естественных условиях, позволяет расширить понимание неоднородности почвы на различных уровнях. Компьютерная томография (CT) почв – активно развивающийся метод изучения структуры почв в ненарушенном состоянии (Абросимов и др., 2021; Mooney et al., 2022) , который в зависимости от используемого разрешения съемки и размера образца почвы позволяет описать неоднородность строения почвы в масштабах агрегатов (мм – см), педона (см – м) и ландшафта (м – км). Характеристика строения агрегатов почв размерами 0.25–10 мм обычно проводится при разрешении в единицы и десятки мкм (например, Papadopoulos et al., 2009; Zhao et al., 2017; Yudina et al., 2022) . Горизонты почв на уровнях педона и ландшафта исследуют с помощью компьютерной макротомографии с разрешением от 10 до 100 мкм и более в монолитах цилиндрической формы объемом от 100 см3 и крупнее (например, Rab et al., 2014; Luo et al., 2010; Lucas et al., 2019) .

Использование компьютерной томографии для изучения почвенной структуры и порового пространства всегда является задачей оптимизации качества изображения и оптимального объема образца. Несмотря на перспективность использования метода СТ, томографическое оборудование мало распространено в почвенных лабораториях по причине довольно высокой стоимости и высоких требований к навыкам оператора и аналитика полученных данных. Обычно качество изображения, т. е. разрешение съемки, определяется исследовательской задачей. Разрешение изображения определяет границы возможного размера образца, и, как правило, для получения изображения с разрешением пикселя 30 мкм и меньше требуется отбор микромонолитов диаметром 2– 3 см. Вопрос о количестве повторностей, необходимых для полноценной характеристики структуры с помощью микромонолитов, остается актуальным, так как обычно исследователи используют минимальное число повторностей, равное 3 или 4.

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

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

Объект исследования – темно-серая почва со вторым гумусовым горизонтом под широколиственным лесом (Каширский район, Московская обл., Россия; 54°46′37.52′′ N; 38°01′55.34′′ E). Название почвы по WRB (2014) – Luvic Retic Greyzemic Phaeozems. Базовые показатели и дополнительная информация об исследованном почвенном разрезе: описание профиля, физические свойства и характеристики порового пространства по генетическим горизонтам, – представлены в работе Скворцовой и др. (2021) ; агрегатный состав и механическая устойчивость гумусового горизонта – Fomin et al. (2021) .

В работе был исследован гумусовый горизонт AU, имеющий мощность 0–40 см. Все микромонолиты почв размерами 2 см в диаметре и 3 см в высоту были отобраны с глубины 10–15 см в 15-кратной повторности в пределах ~3 м2 (по трем стенкам разреза). Микромонолиты сразу после отбора были плотно упакованы и убраны на хранение в холодильную камеру при 4 °С для сохране- ния естественного состояния влажности.

Микротомографическая съемка. Съемка монолитов производилась на рентгеновском компьютерном микротомографе Bruker SkyScan 1172G (США). Образцы ненарушенного строения снимались с разрешением 11 мкм/пиксель, с оборотом образца на 360° вокруг своей оси, угловым шагом 0.4°, параметрами тока на источнике рентгеновского тока 70 кВ и 130 мА, а также с 0.5 мм Al фильтром рентгеновского излучения. Экспозиция каждого кадра составляла 900 мс. Рентгенограмма на каждой угловой позиции является усреднением трех кадров для повышения качества рентгенограммы и снижения шума.

Реконструкция микротомограмм почв. Реконструкция изображений проведена в программе NRecon (США) v 1.7.4.6 от производителя микротомографа. Для каждого набора рентгенограмм была проведена коррекция несоосности ( misalignment compensation ), чтобы скомпенсировать дрейф оси вращения образца относительно оптической оси пучка рентгеновских фотонов. Были использованы алгоритм маскирования дефектных пикселей, алгоритм подавления кольцевых артефактов и алгоритм нивелирования эффекта ослабления пучка, выражающийся в виртуальном увеличении рентгеновской плотности периферической части образца и уменьшении рентгеновской плотности центральной части образца. Для реконструкции изображений был выбран динамический диапазон черного и белого, где граница черного соответствовала значению рентгеновской плотности равной -1 000 единиц Хаунсфилда (HU), а граница белого соответствовала 3 000 HU. То есть все компоненты почвы, например, газовая фаза, имеющая плотность ниже -1 000 HU, на изображении были отображены черным цветом (интенсивность цвета в серых оттенках I = 0); а компоненты почвы, такие как железистые образования или массивные зерна кварца или полевых шпатов, имеющие высокую плотность (3000 HU), отображались белым цветом (интенсивность цвета в серых оттенках I = 255). К данному диапазону рентгеновских плотностей принадлежат все имеющиеся в образцах фазы почвы.

Результаты реконструкции изображений были сохранены в формате 8bit BMP в виде стеков файлов с послойными последовательными изображениями двумерных срезов.

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

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

Рассчитанные томографические параметры. Для анализа порового пространства были рассчитаны следующие параметры – общая пористость, количество пор и средний размер пор.

Общая томографическая пористость ( P t , %) – это процентное соотношение количеств пикселей, относящихся к поровому пространству, ко всему объему образца.

Количество пор ( NOP , шт./мм-3) – это количество отдельных кластеров пикселей, относящихся к поровому пространству, отнесенное к объему образца.

Средний размер пор ( LT , мм) – это средний размер сферы (шара), которую можно вписать в поры. Данный параметр также называется “локальная толщина пор” ( local thickness of pores ).

Все исходные данные приведены в Приложении 1.

Анализ вариограмм. Для выявления количества повторностей, достаточных для описания свойств порового пространства почвы, использовался следующий алгоритм: a) выбиралось конкретное свойство; b) общая выборка разделялась на подвыборки разного размера; c) для каждой подвыборки рассчитывалось среднее значение; d) между средними значениями подвыборок для каждого размера рассчитывалась дисперсия между подвыборками; e) анализировалась взаимосвязь между размерами подвыборки и дисперсией между подвыборками с помощью метода вариограмм; f) определялись значение и доверительный интервал порога ва- риограммы.

Данные для анализа были подготовлены с помощью алгоритма бутстреп ( Bootstrap ). Бутстреп – это метод передискретизации выборки с заменой чисел, который позволяет оценить распределение генеральной совокупности даже по одной выборке. С помощью алгоритма бутстреп были созданы множественные подвыборки с последовательным увеличением размера подвыборки от 2 до 15 повторностей для каждого свойства, т. е. с увеличением числа повторностей отбора образцов. Количество подвыборок для каждого числа повторностей каждого свойства составляло 1 000 шт. Для каждой подвыборки было рассчитано среднее значение, после чего между подвыборками рассчитывалась дисперсия. Для набора дисперсий строились вариограммы, проводились их аппроксимация и анализ. Важно отметить, что по оси абсцисс вместо традиционного для вариографического анализа расстояния между точками, было выбрано количество повторностей отбора образцов.

Аппроксимация вариограмм проводилась экспоненциальной функцией по формуле (Демьянов, Савельева, 2010) :

у(^) = с 0 + (с — с 0 )[1 — exp(1 — 3^ )], (уравенение 1)

где у(^) - дисперсия параметра, h - размер подвыборки, т. е. число повторностей; с 0 – наггет или самородок, т. е. точка пересечения функции с осью ординат; с – плато, т. е. минимальная дисперсия свойства; a – действительный радиус или порог, рассчитанный с поправкой (Chiles, Delfiner, 1999) .

Количество повторностей, достаточное для описания свойств порового пространства почвы, определялось на основе величины и доверительного интервала значения порога. Диапазон повторностей определялся на основе округления доверительного интервала до целых чисел. Количество повторностей, достаточное для описания порового пространства профиля темно-серой почвы с 95%-ной вероятностью, рассчитывалось как максимальное значение для всех исследуемых свойств нижнего предела доверительного интервала.

Метаанализ. Для сравнения вариабельности свойств порового пространства с другими почвами и образцами был проведен метаанализ. Основным требованием к отбираемым статьям – это наличие нескольких пространственных повторностей в рамках педона. Из исследовательских работ (Zhao et al., 2017; Liang et al., 2019; Zhao et al., 2020; Menon et al., 2020; Guo et al., 2020) были собраны и проанализированы данные о количестве повторностей, среднем значении и стандартном отклонении для общей пористости, количества пор и среднго размера пор. На основе собранной информации рассчитывался коэффициент вариации ( CV ) и соотносился с результатами данной работы.

Статистическая обработка. Статистическая обработка и визуализация выполнялись в среде R 4.2.1 . Расчет базовых статистик выполнялся с использованием пакета broom (Robinson, 2014) . Бутстрепинг выполнялся с использованием пакета mosaic (Prim et al., 2017) . Визуализация – с использованием пакета ggplot2 (Wickham et al., 2016) . Тест Шапиро–Уилка выполнялся с уровнем значимости a = 0.1 (Royston, 1995) с использованием пакета stats . Доверительный интервал рассчитывался как двусторонний с уровнем значимости a = 0.05.

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

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

Общая пористость ( P t ) имеет самый низкий коэффициент вариации (17.1%) среди изученных свойств, который позволяет отнести общую пористость к свойствам со средней вариабельностью (Wilding, 1985) . Средне значение P t (27.65%) отличается от медианы (26.6%) менее чем на 4%. Минимум и максимум общей пористости различаются в 1.8 раза (табл. 1).

Количество пор ( NOP ) на основе величины коэффициента вариации относится к свойствам почв с высокой вариабельностью (36.2%). Среднее значение NOP (49.59) отличается от медианы (48.0) на 3.2%. Разброс количества пор между образцами внутри педона составляет 57 пор на мм3.

Средний размер пор ( LT ) – единственное свойство, не имеющее нормального распределения (табл. 1). Кроме того, LT характеризуется самым высоким коэффициентом вариации (46.2%). Среднее значение и медиана LT различаются на 25% (0.39 и 0.31 мм соответственно). В то же время отношение максимума к минимуму LT лишь немногим выше, чем у NOP , и составляет 4.65 раза.

Таблица 1. Вариабельность параметров порового пространства исследованной серой лесной почвы (n = 15)

Table 1. Pore space parameters variability of studied Phaeozem soil (n = 15)

Статистический показатель

Параметр

Pt (%)

NOP (мм-3)

LT (мм)

Среднее ± СО*

27.65 ± 4.72

49.59 ± 17.94

0.39 ± 0.18

Коэффициент вариации

17.1%

36.2%

46.2%

Медиана

26.6

48.0

0.31

Минимум

20.19

16.49

0.19

Максимум

36.20

73.50

0.89

Нижний доверительный интервал (Q2.5)

20.83

21.46

0.20

Верхний доверительный интервал (Q97.5)

34.90

73.21

0.76

Тест Шапиро–Уилка (W/p-value)

0.944/0.44

0.939/0.37

0.847/0.02

Примечание. *СО – стандартное отклонение.

Note. *СО – standard deviation.

Определение оптимального количества повторностей при описании порового пространства почвы. Один из основных методов анализа пространственной неоднородности почвенных свойств – это анализ вариограмм. Анализ основан на описании дисперсии свойства с помощью модели вариограммы, которая включает параметры: самородок (c0), плато (c) и порог (a). В контексте данной работы самородок равен полной дисперсии свойства в горизонте на уровне педона и равняется квадрату стандартного отклонения свойства из таблицы 1. Плато – это дисперсии репрезентативного элементарного объема (representative elementary volume, REV), внутри которого данное свойство в горизонте стационарно. Идеального значения плато можно достичь при использовании неограниченно большого количества повторностей. Однако количество повторностей, которое достаточно для описания дисперсии свойства в горизонте на уровне REV, определяется параметром a, т. е. порогом.

Самородок модели аппроксимации вариограммы P t показывает значение, описывающее значение дисперсии свойства горизонта на уровне педона (рис. 1A). Неоднородность пространственного распределения общей пористости в горизонте на уровне педона (параметр c 0 ) более чем в 12 раз выше, чем на уровне репрезентативного элементарного объема (параметр c ). Оптимальное количество повторностей, которое требуется для достоверного описания неоднородности (p = 0.95) общей пористости внутри горизонта на уровне педона, находится в диапазоне доверительного интервала от 6.11 до 8.80 и в среднем составляет 7.31 шт. При округлении до целого числа минимальное количество повторностей равно 6 шт., а максимальное – 9 шт.

Для общего количества пор пространственная неоднородность в горизонте для педона (параметр c 0 ) почти в 14 раз больше, чем для REV (параметр c ), что еще больше, чем для общей пористости (рис. 1B). Однако для достижения дисперсии уровня REV потребуется меньше повторностей. Доверительный интервал параметра a составляет 5.11–8.32 шт., а среднее значение равно 6.54. При округлении до целого числа минимальное количество повторностей для достоверного определения количества пор равняется 5, а максимальное – 8.

Масштаб различий пространственной неоднородности среднего размера пор в горизонте для педона и REV практически равняется количеству пор (рис. 1C).

2 4 6 8 10 12 14

Число повторностей

Рис. 1. Модели вариограмм свойств порового пространства гумусового горизонта серой почвы. По оси X – количество повторностей отбора образцов (шт.). По оси Y – дисперсии свойств: ( A ) – общая пористость (%), ( B ) – количество пор (мм-3), ( C ) – средний размер пор (мм). Сплошная линия – аппроксимация данных дисперсии на основе экспоненциальной функции (уравнение 1). Вертикальная пунктирная линия обозначает порог (коэффициент a ). Горизонтальная пунктирная линия обозначает уровень плато (коэффициент c ).

Fig. 1. Variogram models of the pore space properties of the Greyic Phaeo-zems humus horizon. X-axis is the number of sampling repetitions (pcs). Y-axis: ( A ) – total porosity (%), ( B ) – number of pores (mm-3), ( C ) – average pore size (mm). The solid line is an approximation of variance data based on an exponential function (Equation 1). The vertical dashed line denotes the coefficient a . The horizontal dashed line denotes the plateau level (the coefficient c ).

Параметр c 0 почти в 14 раз больше, чем параметр c . Однако нижняя граница доверительного интервала порога (параметра a ) для LT на 22% выше, чем для NOP , и на 8% выше, чем для P t , и составляет 6.66 шт. Верхняя граница доверительного интервала равна 8.6 шт., а среднее значение равно 7.54. При округлении до целого числа минимальное количество повторностей для достоверного определения количества пор равняется 7, а максимальное – 9.

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

В масштабе почвенного покрова (весь ареал серых лесных суглинистых и глинистых почв на территории бывшего СССР) неоднородность и среднее значение содержания органического углерода гумусового горизонта серых лесных почв сильно различается между лесными и пахотными вариантами (Урусевская и др., 2000) . Для серых лесных почв под лесом содержание органического углерода характеризуется разбросом от 1.53 до 13.91% (среднее значение и стандартное отклонение равно 5.55 ± 1.10%) и коэффициентом вариации 19%. Анализ углерода во всех горизонтах был оценен для 980 образцов, из которых гумусовый горизонт был представлен 146 образцами.

В масштабе водораздела химические и физико-химические свойства пахотного горизонта серой лесной почвы Владимирского ополья характеризуются различной степенью варьирования (Жарова и др., 2002) . Актуальная кислотность (pH водный) и удельная поверхность (УП) почвы характеризуются крайне низким коэффициентом вариации (КВ) – 3.22 и 8.9% (среднее значение и стандартное отклонение равны 6.82 ± 0.22 и 85.5 ± 7.6 м2/г для pH и УП соответственно). Наибольшее значение КВ характерно для содержания P 2 O 5 и равно 59.2% (9.8 ± 5.8 мг/100 г почвы), в то же время КВ содержания K 2 O более чем в 2 раза ниже и равен

22.9% (8.3 ± 1.9 мг/100 г почвы). При этом количество проб для pH составляло 403 образца, для удельной поверхности – 202 образца, для содержания P 2 O 5 и K 2 O – 279 образцов.

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

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

При исследовании Mollisols под вторичным лесом Китая (Liang et al., 2019) в 4-кратной повторности коэффициент вариации общей пористости и количества пор находится в диапазоне среднее ± стандартное отклонение CV серой почвы (рис. 2A, 2B). Вероятно, неоднородность порового пространства, сформированного с участием лесной растительности, имеет схожий порядок вариабельности у разных почв.

Почвы Gleysols, Ultisols и Calcicorthic Aridisols на пастбищах при 3-кратной повторности (Zhao et al., 2020) были исследованы с использованием тех же параметров порового пространства, что и в нашей работе: Pt, NOP, LT. Коэффициент вариации общей пористости выше, чем в среднем у серой почвы, и находится в диапазоне между средним значением и максимумом. CV количества пор практически равен серой почве. Для среднего размера пор CV существенно ниже, по сравнению со средним значением для серой почвы. В нашей работе вариабельность параметров возрастала в ряду Pt – NOP – LT, в то время как у Zhao et al. (2020) в таком же ряду вариабельность снижалась (рис. 2A, 2B, 2C). Поскольку в работах рассматривались по умолчанию разные почвы, то, вероятно, различия обусловлены типом землепользования. Поэтому можно предположить, что для сельскохозяйственных объектов тенденции вариабельности свойств порового пространства могут в противоположную сторону отличаться от почв под естественными экосистемами.

Гумусовый горизонт почв Mollic Fluvisols под пашней, лугом и лесом при 7–8-кратной повторности показал более низкие значения варьирования общей пористости, чем для серой почвы (Menon et al., 2020) . Все CV общей пористости всех исследуемых в работе объектов находятся в диапазоне между средним значением и минимумом CV серых почв. Данный факт может быть обусловлен тем, что Fluvisols являются молодыми почвами на аллювиальных отложениях с коротким профилем. Вероятно, слабая зрелость аллювиальных почв и близость их состояния к материнской породе обуславливает их однородность и низкую пространственную вариабельность.

В работе Guo et al. (2022) изучались почвы Typic Hapludoll Китая при традиционной вспашке и при использовании технологии no-till в 3-кратной повторности. Общая пористость пахотных объектов крайне однородна, в то время как CV при no-till больше в 28 раз, что существенно превышает средние значения CV серых почв (рис. 1A). Вариабельность общей пористости и количества пор пахотных почв Calcaric Cambisol Китая при выборке 5 повторностей (Zhao et al., 2017) показали существенно более низкие значения CV, чем для серой почвы (рис. 1A, 1B). В обеих работах были исследованы пахотные почвы, в которых естественная вариабельность порового пространства нарушена и обуславливается только воздействием сельскохозяйственной техники. Регулярная вспашка приводит к гомогенизации гумусового горизонта (Bronick, Lal, 2005). При технологии no-till вспашка почвы не про- изводится, и корневая система сельскохозяйственных растений за несколько лет может сформировать поровое пространство с высокой пространственной неоднородностью.

Рис. 2. Зависисмость коэффициента вариации ( CV ) свойств порового пространства для выборок разного размера. ( A ) – общая пористость (%), ( B ) – количество пор (мм-3), ( C ) – средний размер пор (мм). Ось X – число повторностей (шт.). Ось Y – коэффициент вариации (%). Черные точки – среднее значение коэффициента вариации. Черные усы – стандартное отклонение коэффициента вариации. Красные усы – минимум и максимум коэффициента вариации. Цветные точки – данные метаанализа.

Fig. 2. Relationship between the coefficient of variation ( CV ) of pore space properties for different number of replicates of the sample. ( A ) – total porosity, ( B ) – number of pores, ( C ) – average pore size. X-axis – number of replicates (pcs). Y-axis – coefficient of variation (%). Black points – mean value of the CV. Black whiskers – standard deviation of the CV . Red whiskers – minimum and maximum of the CV . Colored dots – meta-analysis data.

Понимание функционирования почвенного покрова большой территории и на уровне педона основано на предположении о стационарности (однородности) свойств внутри выбранной единицы масштаба (Козловский, 2003; Gerke, Karsanina, 2020). Описание неоднородности почвы в целом в определенном масштабе является крайне сложной задачей, и для этого используют наборы интегральных свойств или индикаторов (Козловский, 2003; Pachepsky, Hill, 2017). Индикаторы должны обладать важным свойством – чувствительностью к разнообразию условий в рамках заданного масштаба (Priorr, 2003). Общая пористость, количество пор и средний размер пор являются интегральными параметрами порового пространства и почвенной структуры в целом (Reynolds et al., 2009; Rabot et al., 2018).

Чувствительность индикаторов порового пространства к пространственной неоднородности не означает чувствительность к динамике почвенной структуры во времени и к воздействию внешних факторов на нее (Vogel et al., 2019) . Ранее нами было показано для 3 генетических горизонтов дерново-подзолистой почвы, что общая пористость является плохим индикатором динамики порового пространства при изменении влажности почвы (Fomin et al., 2023) . Как раз одной из причин слабой пригодности общей пористости к оценке динамики структуры является ее высокая пространственная неоднородность.

Важной частью оценки пространственной неоднородности горизонта на уровне педона является отбор достаточного количества образцов для анализа (Pennock et al., 2008) . По результатам проделанной работы на примере гумусового горизонта серой лесной почвы мы можем дать рекомендации по количеству образцов для анализа порового пространства. Минимальный размер выборки определялся как максимальное округленное до целого значение нижнего доверительного интервала порога, т. е. коэффициента a . Для общей пористости минимальное число повторностей равняется 6, для количества пор – 5, для среднего размера пор – 7. Итого, минимальное количество образцов, которое стоит отбирать из профиля почвы для оценки порового пространства, составляет 7 шт. Более того, классический отбор проб для анализа порового пространства почвы (и других свойств в том числе) в 3-кратной повторности даст возможность описать пространственную неоднородность только на 25–30% (для P t – 29%, для NOP – 25%, для LT – 30%).

ЗАКЛЮЧЕНИЕ

В данной работе мы изучили неоднородность распределения на уровне педона общей пористости, количества и среднего размера пор гумусового горизонта серой почвы под лесом. Сравнение результатов при разном размере выборки показало снижение дисперсии при увеличении размера выборки. Анализ вариограмм поз- волил оценить различия вариабельности свойств для уровня педона и репрезентативного элементарного объема. Наименьшая разница неоднородности между педоном и репрезентативным элементарным объемом REV наблюдается для общей пористости (~12 раз), в то время как для количества и среднего размера пор это соотношение больше (~14 раз). Анализ порога позволил выявить различия в количестве повторностей, необходимых для полного описания неоднородности на уровне педона для каждого свойства. В среднем уровень порога, при котором дисперсия свойств прекращала снижаться, для Pt равняется 7.31 ± 0.56 шт., для NOP – 6.54 ± 0.61 шт., для LT – 7.54 ± 0.41 шт.

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

Результаты работы показали, что общая пористость, количество пор и средний размер пор являются подходящими индикаторами порового пространства и почвенной структуры, которые имеют высокую чувствительность к пространственной неоднородности. По итогам данной работы мы можем дать рекомендации по количеству повторностей микромонолитов для анализа порового пространства на уровне педона, на примере гумусового горизонта серой почвы. Минимальное количество повторностей составляет 7 шт. Отбор и анализ микромонолитов в 3-кратной повторности позволяет описать неоднородность педона только на 25–30%.

ПРИЛОЖЕНИЕ

Приложение 1. Исходные данные

Повторность

NOP

LT

Pt

1

35.34732

0.461487

31.83155

2

16.49299

0.18891

32.1104

3

70.23075

0.311226

28.96115

4

48.0025

0.888957

36.19992

5

62.38746

0.515515

31.97288

6

73.49665

0.224321

26.60423

7

37.39783

0.298239

22.01071

8

42.79712

0.214597

25.39324

9

72.67089

0.29925

23.46264

10

30.69543

0.399024

30.84838

11

71.79816

0.280292

20.1869

12

58.00355

0.460414

24.69886

13

31.7377

0.273979

32.49419

14

51.6884

0.49043

24.57977

15

41.07253

0.506368

23.39107

Статья научная