Влияния внесения биоугля на теплофизические свойства дерново-подзолистой почвы и составляющие энергетического баланса яровой пшеницы в Ленинградской области при различных условиях почвенного увлажнения

Биоуголь является продуктом пиролиза биомассы в бескислородной среде. Он имеет пористую структуру и содержит множество различных функциональных групп (Lehmann, Joseph, 2015). Биоуголь может содержать некоторое количество гумино- вых веществ (Lin et al., 2012). Кроме того, его молекулярная структура имеет высокую степень химической и микробиологической стабильности (Cheng et al., 2008). Физические и химические свойства биоугля сильно зависят от температуры и продолжительности пиролиза, типа исходного сырья (Joseph et al., 2010; Bruun et al., 2011). В качестве исходного материала для производства биоугля используется широкий спектр сырья, включая древесную щепу, органические отходы, растительные остатки и птичий помет (Sohi et al., 2010). Элементный состав биоугля обычно включает углерод, азот, водород и вспомогательные питательные элементы, такие как K, Ca, Na и Mg (Zhang et al., 2015). Биоуголь имеет большую удельную площадь поверхности и содержит ряд полярных и неполярных молекул, обладающих сильным сродством к неорганическим ионам, например, к ионам тяжелых металлов, фосфатам и нитратам (Schmidt et al., 2015; Kammann et al., 2015). Внесение биоугля может улучшить структуру почвы, увеличить пористость, уменьшить плотность и увеличить влагоемкость почвы (Baiamonte et al., 2015). Биоуголь может увеличить электропроводность почвы (Oguntunde et al., 2004) и емкость катионного обмена (Laird et al., 2010), а также снизить кислотность почвы (Oguntunde et al., 2004). Кроме того, внесение биоугля способствует увеличению микробной биомассы (Liang et al., 2010), состава бактериального сообщества почвы (Grossman et al., 2010), базального дыхания (Steiner et al., 2008). Секвестрация углерода за счет применения биоугля приводит к повышению плодородия почв с низким содержанием питательных веществ (Jiang et al., 2012; Liu et al., 2012). Существуют различные механизмы повышения доступности питательных веществ для растений в агроэкосистемах: 1) минерализация лабильной фракции биоугля, которая содержит органически связанные питательные вещества (Sohi et al., 2010; Lehmann et al., 2012); 2) снижение выщелачивания питательных веществ благодаря физико-химическим свойствам биоугля (Liang et al., 2006); 3) снижение потерь азота при денитрификации

(Cayuela et al., 2013); 4) удержание азота, фосфора и серы, связанное с усилением биологической активности почвы или изменением бактериального сообщества (Pietikäinen et al., 2000). Согласно результатам некоторых полевых экспериментов применение био- угля улучшило качество почвы, увеличило урожайность и способствовало росту растений (Lehmann et al., 2006; Major et al., 2010; Zhang et al., 2010; Uzoma et al., 2011). Однако в работе (Asai et al., 2009) приводятся данные о снижении урожайности при применении биоугля. Снижение урожайности может быть связано с высоким содержанием летучих, а также токсичных и вредных веществ в биоугле, которые способствуют снижению поглощения питательных веществ, и, как следствие, подавляют рост растений. Таким образом, повышение урожайности может зависеть от свойств биоугля и свойств почвы.

Теплофизические свойства почвы (теплопроводность, объемная теплоемкость и температуропроводность) влияют на энергетический баланс подстилающей поверхности, аккумуляцию и передачу тепла в почве. Внесение биоугля может снижать теплопроводность и температуропроводность почвы (Zhang et al., 2013; Usowicz et al., 2016; Zhao et al., 2016) и, как следствие, снижать количество тепла, передаваемого через почву (теплопроводность). Терморегулирующий потенциал биоугля выглядит многообещающим, так как внесение биоугля может иметь стратегическое значение для управления климатическими рисками, такими как засухи и волны тепла, которые в последнее время остро ощущаются в Северо-Западном регионе России. Теплофизические свойства почвы в значительной степени коррелируют с другими свойствами почвы, такими как плотность и влажность почвы, содержание органического вещества, концентрация растворенных веществ. Например, теплопроводность почвы может увеличиваться с увеличением плотности и влажности почвы и уменьшаться с увеличением органического вещества почвы и концентрации солей (Abu-Hamdeh, Reeder, 2000; Zhao et al., 2016) . Снижение плотности при применении биоугля может снизить теплопроводность (Zhang et al., 2013) , тогда как увеличение влагоемкости почвы при применении биоугля может увеличить теплопроводность (Usowicz et al., 2016) .

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

Бюллетень Почвенного института им. В.В. Докучаева. 2023. Вып. 116

Dokuchaev Soil Bulletin, 2023, 116 ново-подзолистой супесчаной почвы при различных условиях увлажнения и влагообеспеченности посевов яровой пшеницы в Ленинградской области и на составляющие энергетического баланса. Данная работа является актуальной, так как влияние внесения биоугля на теплофизические свойства почвы изучено меньше, чем влияние биоугля на физико-химические свойства и микробиологические показатели. Изменение теплофизических свойств почвы при внесении биоугля может быть различно для почв, отличающихся гранулометрическим составом. Исследования влияния биоугля на теплофизические свойства различных типов почв при разном уровне увлажнения являются важными для задач моделирования роста и развития растений при применении биоугля (Archontoulis et al., 2016) .

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

Полевые эксперименты проводились во время вегетационного периода 2022 г. на Меньковской опытной станции Агрофизического научно-исследовательского института (Гатчинский район Ленинградской области). При проведении экспериментов изучалось влияние внесения биоугля на составляющие энергетического баланса при различных условиях почвенного увлажнения. Возделываемая сельскохозяйственная культура на опытных участках – яровая мягкая пшеница сорта “Дарья” ( Triticum aestivum L.). Эксперимент включал опытный участок с внесением биоугля в дозе 21.9 т га^-1 в пахотный горизонт почвы и контрольный участок. Биоуголь был изготовлен из отходов деревообрабатывающей промышленности методом медленного пиролиза при температуре 600 °С. Биоуголь внесли на поверхность участка и заделали в почву 05.05.2022 г. перед посевом. Посев был произведен 07.05.2022 г. Площадь опытных участков составляла 16 м². Координаты опытного участка с внесением биоугля – 59°25'28.6" с. ш., 30°02'40.3" в. д.; контрольного участка –  59°25'28.6" с. ш.,

30°02'39.6" в. д. Перед посевом в почву дополнительно внесли азофоску (NPK) марки 15 : 15 : 15 в количестве 50 кг д. в. по каждому элементу. Полевые измерения проводились с 11 до 16 часов по местному времени в следующие даты: 19.05.2022, 31.05.2022, 09.06.2022, 23.06.2022, 07.07.2022, 21.07.2022, 04.08.2022.

Для анализа метеорологической информации за весь вегетационный период яровой пшеницы были использованы срочные данные близлежащей стандартной метеорологической станции “Белогорка” (WMO ID = 26069). Почва на опытных участках – супесчаная дерново-подзолистая. Гранулометрический состав пахотного горизонта почвы на опытных участках (Моисеев и др., 2014) : 10.08% – физическая глина; 18.80% – пыль; 70.12% – песок; 1.02% – гравий.

Определялись следующие теплофизические свойства почвы: теплопроводность (K, Вт м^-1 К^-1), объемная теплоемкость (C, МДж м^-3 К^-1), температуропроводность (D, мм² с^-1). Теплофизические свойства почвы измеряли методом линейного импульсного источника теплоты (He et al., 2018; Brunetti et al., 2022) с помощью сенсора SH-3 прибора TEMPOS. Измерения теплофизических свойств проводились в полевых условиях каждые 15 минут. Период измерений был выбран исходя из времени наступления термодинамического равновесия между сенсором и почвой после изменения расположения сенсора, а также времени измерения прибора. Для опытных участков было рассчитано еще одно комплексное теплофизическое свойство, характеризующее водный и тепловой режим верхнего слоя почвы, – тепловая инерция (ТИ, Дж м^-2 К^-1 с^-0.5), или теплоусвояемость (De Vries, 1963) . Тепловая инерция представляет собой сопротивление почвы изменению температуры ее верхнего слоя (несколько сантиметров). Тепловая инерция широко используется в моделях определения влажности верхнего слоя почвы (Lu et al., 2009) :

ти = V к • с.

Уравнение энергетического баланса подстилающей поверхности формулируется как (Choudhury et al., 1987):

LE = Rn - H - G, где LE – скрытый поток тепла (Вт м-2), Rn – радиационный баланс (Вт м-2), H – турбулентный поток тепла (Вт м-2), G – поток тепла в почву (Вт м-2).

Суммарное испарение – сумма испарения почвы и транспирации. Суммарное испарение является важнейшим показателем энерго-массообмена в системе “почва–растение–приземный слой воздуха”. Суммарное испарение (ET, мм час^-1) определялось через скрытый поток тепла (LE, Вт м^-2). Скрытый поток тепла оценивался как остаточный член в уравнении энергетического баланса с помощью специальных агрометеорологических измерений прибором АМПАК (Efimov et al., 2018) .

Радиационный баланс над участком с внесением биоугля измерялся с помощью балансомера (Foken et al., 2021) , радиационный баланс для участка без внесения биоугля рассчитывался через изменение уходящего длинноволнового излучения в уравнении радиационного баланса (Brutsaert, 1975) .

Турбулентный поток тепла определялся аэродинамическим методом (Liu et al., 2007) :

„  P a C p (T - T a )

H =-----------,

rah где ρa – плотность воздуха (кг м-3); Cp – удельная теплоемкость воздуха при постоянном давлении (~ 1013 Дж кг-1 К-1) (Allen, 1998); Ta – температура воздуха на высоте 2 м (°C); Ts – температура подстилающей поверхности (°C); rah – аэродинамическое сопротивление (с м-1).

Поток тепла в почву оценивался по данным измерений радиационного баланса и проективному покрытию (Su, 2002) . Проективное покрытие определялось с помощью фотографий растительности (Bennett et al., 2000) . На фотографиях выделялись пиксели, соответствующие растительности, через фильтр преобладания зеленого цвета в матрице RGB. Проективное покрытие рассчитывалось как площадь растительности по отношению к общей площади поверхности:

G = R,-[Г, +(1-/, )-(Fs -Г,)], где Γc – коэффициент для полностью сомкнутого растительного покрова (= 0.05); Γs – коэффициент для оголенной почвы (= 0.315); fc – проективное покрытие (-).

Температура подстилающей поверхности является интегральной величиной температуры излучения почвы и листьев растений. Интегральную температуру подстилающей поверхности на опытных участках измеряли бесконтактным методом с помощью пирометров Optris CT LT (Khanal et al., 2017) . Измерения температуры подстилающей поверхности проводились с интервалом 90 секунд. Различие температуры подстилающей поверхности и температуры воздуха характеризует интенсивность теплообмена между поверхностью Земли и атмосферой. При этом температура листьев растений с учетом состояния атмосферы характеризует водный стресс растений (Jones, Leinonen, 2003) и интенсивность транспирации (Zhao et al., 2022) . Для определения температуры листьев была проведена съемка тепловизионной камерой Testo 868 в каждый экспериментальный день.

Объемная влажность почвы определялась на участке с внесением биоугля с помощью емкостного датчика влажности почвы через измерения диэлектрической проницаемости почвы Decagon 10HS (Visconti et al., 2014) . Измерения объемной влажности почвы проводились с интервалом 90 секунд.

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

В результате анализа информации о погоде за вегетационный период 2022 г. по стандартной метеорологической станции “Белогорка” (рис. 1) было выделено две экспериментальные даты с контрастными метеорологическими условиями:

• 1)    07.07.2022 г. сложились засушливые погодные условия, которые характеризуются предшествующим длительным периодом высоких температур воздуха без выпадения осадков (17 дней без осадков со среднесуточной температурой 19.7 °С);

• 2)    04.08.2022 г. отмечается как дата высокого увлажнения, так как за двое суток (02.08.2022 г.) выпало максимальное количество осадков (71.3 мм) за весь вегетационный период.

В таблице 1 представлены осредненные значения параметров микроклимата сельскохозяйственного поля, измеренные при- бором АМПАК: температура воздуха (Ta), относительная влажность воздуха (RH), атмосферное давление (P), скорость ветра (u); данные объемной влажности почвы (W), измеренной на глубине 10 см сенсором Decagon 10HS; данные фенологических измерений – высота растений (hc) и оценка проективного покрытия (fc) по результатам обработки фотографий. Высота растений (hc) необходима для оценки параметра шероховатости количества движения и потока тепла, которые используются для определения аэродинамического сопротивления (rah). Данные измерений объемной влажности почвы (W) (табл. 1) согласуются с данными количества осадков, измеренного на метеорологической станции “Белогорка”:

• 1)    07.07.2022 г. наблюдаются засушливые условия с отсутствием влаги в почве (W = 0.01 м³ м^-3);

Рис. 1. Среднесуточная температура воздуха (°C) и количество осадков (мм) по данным метеорологической станции “Белогорка”.

Fig. 1. Average daily air temperature (°C) and precipitation (mm). Data from the meteorological station “Belogorka”.

• 2)    04.08.2022 г. наблюдается максимальное за даты проведения экспериментов значение влажности почвы (W = 0.26 м³ м^-3).

19.05.

31.05.

09.06.

23.06.

07.07.

21.07.

04.08.

т

a

11.7

15.7

19.5

22.2

20.0

24.2

25.3

RH

57

70

64

45

57

60

51

P

101.3

101.1

100.7

101.0

100.6

101.1

101.7

u

2.0

1.7

2.2

1.7

1.6

2.4

1.3

W

0.19

0.25

0.17

0.18

0.01

0.17

0.26

^hc

-

0.07

0.15

0.30

0.60

0.65

0.65

f c

0

0.15

0.3

0.6

0.7

0.75

0.7

После предобработки исходных данных и удаления выбросов была получена 91 точка измерений теплофизических свойств пахотного горизонта почвы (50 точек для участка с внесением биоугля и 41 точка для контрольного участка). В целом в течение одного дня теплофизические свойства на обоих участках практически не изменялись, но в измерениях присутствовали выбросы. Сильные выбросы могли происходить из-за наличия камня или гравия в радиусе измерения сенсора. Затем было проведено осреднение за время проведения каждого эксперимента. Были осредне-ны значения теплопроводности (K), объемной теплоемкости (C), температуропроводности (D), тепловой инерции (ТИ) для участка с внесением биоугля (табл. 2) и контрольного участка (табл. 3).

Таблица 1. Параметры микроклимата сельскохозяйственного поля: температура воздуха (T a , °C), относительная влажность воздуха (RH, %), атмосферное давление (P, кПа), скорость ветра (u, м с^-1); объемная влажность почвы (W, м³ м^-3); высота растений (h_c, м); проективное покрытие (f_c, -) в даты проведения экспериментов

Table 1. Microclimate parameters of an agricultural field: air temperature (T_a, °C), relative air humidity (RH, %), atmospheric pressure (P, kPa), wind speed (u, m s^-1); volumetric soil moisture (W, m³ m^-3); crop height (h_c, m); projective cover (f c , -) on the experiment’s dates

Таблица 2. Осредненные значения теплофизических свойств пахотного горизонта почвы на опытном участке с внесением биоугля: теплопроводность (K, Вт м^-1 К^-1), объемная теплоемкость (C, МДж м^-3 К^-1), температуропроводность (D, мм² с^-1), тепловая инерция (ТИ, Дж м^-2 К^-1 с^-0.5)

Table 2. Average values of thermal soil properties on the experimental plot with the biochar application: thermal conductivity (K, W m^-1 K^-1), volumetric heat capacity (C, MJ m^-3 K^-1), diffusivity (D, mm² s^-1), thermal inertia (ТИ, J m^-2 K^-1 s^-0.5)

Дата	K биоуголь	биоуголь	D биоуголь	I ll „ „ биоуголь
19.05.2022	0.358	1.072	0.334	619.6
31.05.2022	0.435	1.295	0.336	750.5
09.06.2022	0.245	1.014	0.241	498.2
23.06.2022	0.389	1.189	0.327	679.9
07.07.2022	0.158	0.965	0.164	390.2
21.07.2022	0.505	1.210	0.417	781.6
04.08.2022	1.155	1.840	0.628	1457.8

Таблица 3. Осредненные значения теплофизических свойств пахотного горизонта почвы на контрольном участке: теплопроводность (K, Вт м^-1 К^-1), объемная теплоемкость (C, МДж м^-3 К^-1), температуропроводность (D, мм² с^-1), тепловая инерция (ТИ, Дж м^-2 К^-1 с^-0.5)

Table 3. Average values of thermal soil properties on the control site: thermal conductivity (K, W m^-1 K^-1), volumetric heat capacity (C, MJ m^-3 K^-1), diffusivity (D, mm² s^-1), thermal inertia (ТИ, J m^-2 K^-1 s^-0.5)

Дата	K контроль	г контроль	D контроль	ТИ контроль
19.05.2022	0.290	1.057	0.274	553.9
31.05.2022	0.493	1.385	0.356	826.3
09.06.2022	0.394	1.288	0.306	712.5
23.06.2022	0.513	1.355	0.379	833.9
07.07.2022	0.225	1.185	0.190	515.8
21.07.2022	0.618	1.643	0.376	1008.0
04.08.2022	1.06	1.889	0.559	1412.1

Минимальные значения всех теплофизических свойств пахотного горизонта почвы приходятся на засушливый период с минимальной влажностью почвы – 07.07.2022 г. Максимальные значения теплофизических свойств почвы приходятся на 04.08.2022 г., когда наблюдалась максимальная влажность почвы.

Статистический анализ значимости различий теплофизических свойств почвы выполнен с помощью непараметрического U критерия Манна–Уитни–Уилкоксона, так как данные не прошли проверку на нормальность распределения по результатам теста Шапиро–Уилка, а также количество наблюдений в каждой выборке невелико (n < 30).

Теплопроводность почвы 07.07.2022 г. на 29.7% меньше на участке с внесением биоугля (Кбиоуголь), чем на контрольном участке (Kконтроль). Объемная теплоемкость почвы 07.07.2022 г. на 18.5% меньше на участке с внесением биоугля (Cбиоуголь), чем на контрольном участке (Cконтроль). Температуропроводность почвы 07.07.2022 г. на 13.7% меньше на участке с внесением биоугля (Dбиоуголь), чем на контрольном участке (Dконтроль). Тепловая инерция почвы 07.07.2022 г. на 24.3% меньше на участке с внесением биоугля (ТИбиоуголь), чем на контрольном участке (ТИконтроль). Различия между средними значениями теплофизических свойств почвы на участке с внесением биоугля и контрольном участке 07.07.2022 г. статистически значимые (p-value < 0.05) согласно тесту Манна–Уитни–Уилкоксона. Снижение всех теплофизических свойств почвы при внесении биоугля происходит из-за того, что значения теплофизических свойств биоугля меньше, чем соответствующие значения минеральной части твердой фазы почвы (Liu et al., 2018). При этом внесение биоугля может увеличивать влагоемкость почвы (Basso et al., 2013; Suliman et al., 2017), особенно на легких почвах, вследствие чего значения некоторых теплофизических свойств почвы при внесении биоугля могут быть и выше, чем без внесения биоугля, что подтверждается результатами измерений теплофизических свойств почвы 04.08.2022 г. после выпадения большого количества осадков. Теплопроводность почвы 04.08.2022 г. на 9.4% больше на участке с внесением биоугля (Кбиоуголь), чем на контрольном участке (Kконтроль). Объемная теплоемкость почвы 04.08.2022 г. на 2.6% меньше на участке с внесени- ем биоугля (Cбиоуголь), чем на контрольном участке (Cконтроль). Температуропроводность почвы 07.07.2022 г. на 12.3% больше на участке с внесением биоугля (Dбиоуголь), чем на контрольном участке (Dконтроль). Тепловая инерция почвы 04.08.2022 г. на 3.2% больше на участке с внесением биоугля (ТИбиоуголь), чем на контрольном участке (ТИконтроль). Различия между средними значениями теплофизических свойств почвы на участке с внесением биоугля и контрольном участке 04.08.2022 г. статистически значимые (p-value < 0.05) согласно тесту Манна–Уитни–Уилкоксона.

Биоуголь, изменяя теплофизические свойства почвы, может влиять на составляющие энергетического баланса и суммарное испарение через уменьшение отраженной радиации (Genesio, 2012) , увеличение собственного излучения Земли (Genesio, 2012) , увеличение температуры подстилающей поверхности (Feng et al., 2021) , а также снижение испарения с поверхности почвы (Wang et al., 2018) . В таблицах 4 и 5 приведены оценки влияния биоугля на составляющие энергетического баланса и суммарное испарение.

Таблица 4. Составляющие энергетического баланса (радиационный баланс (R n , Вт м^-2), турбулентный поток тепла (H, Вт м^-2), поток тепла в почву (G, Вт м^-2), скрытый поток тепла (LE, Вт м^-2) и суммарное испарение (ET, мм час^-1) на опытном участке с внесением биоугля

Table 4. Energy balance components (radiation balance (R_n, W m^-2), turbulent heat flux (H, W m^-2), soil heat flux (G, W m^-2), latent heat flux (LE, W m^-2)) and evapotranspiration (ET, mm h -1 ) on the experimental plot with the biochar application

Дата	R n биоуголь	H биоуголь	G биоуголь	LE биоуголь	ET биоуголь
19.05.2022	204.8	163.1	64.5	13.7	0.02
31.05.2022	336.3	76.1	92.6	168.1	0.25
09.06.2022	324.6	151.0	76.5	106.8	0.16
23.06.2022	404.3	144.7	63.1	197.0	0.29
07.07.2022	301.13	164.6	39.0	114.8	0.17
21.07.2022	446.0	82.7	51.9	311.5	0.46
04.08.2022	483.7	25.5	62.6	395.6	0.58

Анализ результатов определения составляющих энергетического баланса над участком с внесением биоугля и контрольным участком показал:

• 1)    Радиационный баланс на участке с внесением биоугля (R n биоуголь ), по сравнению с радиационным балансом на контрольном участке (R n контроль ), в среднем уменьшился на 2.5%. Среднее значение для R n биоуголь = 360.0 Вт м^-2, для R n контроль = 369.2 Вт м^-2. Значение критерия Манна–Уитни–Уилкоксона для оценки различий средних значений (p-value > 0.1) свидетельствует об отсутствии достоверных различий средних значений R_{n биоуголь} и R n контроль .

Таблица 5. Составляющие энергетического баланса (радиационный баланс (R_n, Вт м^-2), турбулентный поток тепла (H, Вт м^-2), поток тепла в почву (G, Вт м^-2), скрытый поток тепла (LE, Вт м^-2)) и суммарное испарение (ET, мм час^-1) на контрольном участке

Table 5. Energy balance components (radiation balance (R_n, W m^-2), turbulent heat flux (H, W m^-2), soil heat flux (G, W m^-2), latent heat flux (LE, W m^-2)) and evapotranspiration (ET, mm h^-1) on the control plot

Дата	R n контроль	H контроль	контроль	EE.„ контроль	контроль
19.05.2022	215.7	125.1	67.9	31.3	0.05
31.05.2022	343.7	53.3	94.6	195.8	0.29
09.06.2022	335.6	106.3	79.0	151.1	0.22
23.06.2022	426.0	81.3	66.5	278.2	0.41
07.07.2022	310.9	135.5	40.3	140.3	0.21
21.07.2022	446.6	80.5	51.9	314.2	0.46
04.08.2022	485.5	21.5	62.9	401.1	0.59

• 2)    Турбулентный поток тепла над участком с биоуглем (H_{биоуголь}), по сравнению с турбулентным потоком тепла над контрольным участком (H контроль ), в среднем увеличился на 35.5%. Среднее значение для H биоуголь = 113.5 Вт м^-2, для H_{контроль} = 83.8 Вт м^-2. Значение критерия Манна–Уитни– Уилкоксона для оценки различий средних значений

(p-value < 0.05) свидетельствует о наличии статистически значимых различий средних значений H биоуголь и H контроль .

• 3)    Поток тепла в почву на участке с биоуглем (G биоуголь ), по сравнению с потоком тепла в почву на контрольном участке (G контроль ), в среднем уменьшился на 2.8%. Среднее значение для G биоуголь = 66.1 Вт м^-2, для G контроль = 68.0 Вт м^-2. Значение критерия Манна–Уитни–Уилкоксона для оценки различий средних значений (p-value > 0.1) свидетельствует об отсутствии достоверных различий средних значений G контроль .

• 4)    Скрытый поток тепла над участком с биоуглем (LE_{био-уголь}), по сравнению со скрытым потоком тепла над контрольным участком (LE контроль ), в среднем уменьшился на 17.0%. Среднее значение для LE биоуголь = 180.5 Вт м^-2, для LE контроль = 217.5 Вт м^-2. Значение критерия Манна–Уитни–Уилкоксона для оценки различий средних значений (p-value < 0.05) свидетельствует о наличии статистически значимых различий средних значений LE биоуголь и T F

контроль .

• 5)    Суммарное испарение над участком с биоуглем (ET_био- уголь ), по сравнению с суммарным испарением над контрольным участком (ET_{контроль}), в среднем уменьшилось на 13.9%. Среднее значение для ET_{биоуголь} = 0.28 мм час -1 , для ET контроль = 0.32 мм час -1 . Значение критерия Манна–Уитни–Уилкоксона для оценки различий средних значений (p-value < 0.05) свидетельствует о наличии статистически значимых различий средних значений ET_{биоуголь} и FT ™ контроль .

Увеличение турбулентного потока тепла над участком с внесением биоугля (Hбиоуголь) происходит из-за увеличения температуры подстилающей поверхности на участке с внесением биоугля (рис. 2). Особенно ярко различия температуры подстилающей поверхности проявляются в начале и середине периода вегетации. Различия в турбулентном потоке тепла обусловлены более высокой температурой подстилающей поверхности с внесением биоугля. Средняя температура подстилающей поверхности на участке с биоуглем составила 25.9 °C, на контрольном участке – 24.4 °C (увеличение температуры подстилающей поверхности на участке с биоуглем на 6.4% по сравнению с контролем). Увеличение турбулентного потока тепла может также происходить из-за умень- шения теплопроводности почвы, так как радиационный баланс будет в большей мере расходоваться на передачу энергии от подстилающей поверхности к атмосфере, чем на передачу энергии вглубь почвы.

Биоуголь обладает меньшей отражательной способностью (альбедо) по сравнению с почвой (Genesio, 2012) , но при этом биоуголь увеличивает собственное излучение вследствие увеличения температуры подстилающей поверхности. В результате два противоположных потока радиационного баланса компенсируют друг друга и значимых различий средних значений радиационного баланса на участке с внесением биоугля (R_{n биоуголь}) и контрольном участке (R n контроль ) обнаружено не было.

Рис. 2. Температура подстилающей поверхности (°C) на участке с внесением биоугля и на контрольном участке.

Fig. 2. Crop surface temperature (°C) at the plot with biochar application and at the control plot.

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

Скрытый поток тепла (LE) и суммарное испарение (ET) уменьшались при внесении биоугля, вероятно, из-за снижения испарения почвы (Wang et al., 2018) , а не из-за снижения транспирации. Для проверки данного предположения был проведен анализ температуры листьев растений с помощью тепловизионных изображений, так как температура листьев сильно связана с транспирацией растений (Gates, 1968) . В таблице 6 представлены результаты обработки съемок тепловизионной камеры Testo 868, в таблице 6 представлена температура листьев на участке с внесением биоугля (Т_{c биоуголь}) и на контрольном участке (Т_{c контроль}). Согласно критерию Манна–Уитни–Уилкоксона значимых различий между температурой листьев на участке с внесением биоугля и на контрольном участке обнаружено не было (p-value > 0.1). На рисунке 3 представлен пример съемки тепловизионной камеры 23.06.2022 г. для участка с внесением биоугля и на контрольном участке.

Таблица 6. Температура листьев растений (T c , °C) на участке с внесением биоугля и контрольном участке в даты проведения экспериментов

Table 6. Leaf temperature (T_c, °C) on the plot with the biochar application and on the control plot on the experiment’s dates

" "——_	31.05.	09.06.	23.06.	07.07.	21.07.	04.08.
Время съемки	11:57	12:03	11:54	11:29	13:14	12:34
Т c биоуголь	16.7	26.3	27.1	23.1	31.7	29.9
т c контроль	17.0	26.4	27.4	23.2	31.4	30.3

Рис. 3. Тепловое изображение участка с внесением биоугля (слева) и контрольного участка (справа) 23.06.2022 г.

Fig. 3. Thermal image of the plot with the biochar application (left) and the control plot (right) on June 23, 2022.

Интегральная температура подстилающей поверхности по съемке тепловизионной камерой 23.06.2022 г. в 11:54 равна 31.9 °C для участка с внесением биоугля и 30.4 °C для контрольного участка. Снижение собственного испарения почвы при внесении биоугля может быть вызвано изменением структуры почвы и изменением ее гидрофизических и теплофизических свойств. Биоуголь увеличивает количество микропор, которые имеют плохую влагопроводность. В результате чего снижается влагопе-ренос из нижележащих слоев к поверхности почвы (Wang et al., 2018) .

В настоящее время активно развивается направление определения влажности почвы и влагообеспеченности растений по результатам оценки тепловой инерции с помощью данных дистанционного зондирования Земли и тепловизионных съемок с беспилотных летательных аппаратов (Lu et al., 2009; Soliman et al., 2013; Qin et al., 2013). На рисунке 4 представлены следующие регрессионные зависимости: зависимость между тепловой инерцией на участке с биоуглем (ТИбиоуголь) и объемной влажностью почвы (Wбиоуголь); между тепловой инерцией на участке с биоуглем (ТИбиоуголь) и суммарным испарением (ETбиоуголь); между тепловой инерцией на контрольном участке (ТИконтроль) и суммарным испарением на контрольном участке (ETконтроль). Ниже представлены значения параметров моделей взаимосвязи теплофизических свойств пахотного горизонта почвы (теплопроводность (K), объемная теплоемкость (C) и температуропроводность(D)) и тепловой инерции (ТИ) с объемной влажностью почвы (W) и суммарным испарением (ET), а также их коэффициенты детерминации (R2) с p-уровнем значимости (p-value):

• •    Взаимосвязь объемной влажности почвы и теплопроводности для участка с внесением биоугля (W биоуголь ~ K биоуголь ):

W биоуголь = 0.16599 × K биоуголь + 0.09877 (R²= 0.32, p-value = 0.107)

• •    Взаимосвязь суммарного испарения и теплопроводности для участка с внесением биоугля (ET биоуголь ~ K биоуголь ): ET биоуголь = 0.47346 × K биоуголь + 0.05623 (R²= 0.59, p-value = 0.027)

• •    Взаимосвязь суммарного испарения и теплопроводности для контрольного участка (ET контроль ~ K контроль ):

ET контроль = 0.58317 × K контроль + 0.01924 (R²=0.75, p-value = 0.007)

• •    Взаимосвязь объемной влажности почвы и объемной теплоемкости для участка с внесением биоугля (W биоуголь ~ C биоуголь ):

W биоуголь = 0.19047 × C биоуголь – 0.05788 (R²= 0.36, p-value = 0.09)

• •    Взаимосвязь суммарного испарения и объемной теплоемкости для участка с внесением биоугля (ET биоуголь ~ C биоуголь ):

ET биоуголь = 0.5250 × C биоуголь – 0.3682 (R²= 0.59, p-value = 0.026)

• •    Взаимосвязь суммарного испарения и объемной теплоемкости для контрольного участка (ET контроль ~ C контроль ): ET контроль = 0.60799 × C контроль – 0.53279 (R²= 0.88, p-value = 0.001)

• •    Взаимосвязь объемной влажности почвы и температуропроводности для участка с внесением биоугля (W биоуголь ~ D биоуголь ):

W биоуголь = 0.41111 × D биоуголь + 0.03200 (R²= 0.45, p-value = 0.059)

Рис. 4. Регрессионные линии взаимосвязи тепловой инерции (ТИ, Дж м^-2 К^-1 с^-0.5) с суммарным испарением (ET, мм час^-1) и с объемной влажностью почвы (W, м³ м^-3) для участка с внесением биоугля (слева) и для контрольного участка (справа).

Fig. 4. Regression lines between thermal inertia (ТИ, J m -2 K -1 s^-0.5) and evapotranspiration (ET, mm h -1 ), and soil volumetric moisture (W, m³ m^-3) for the plot with the biochar application (left) and for the control plot (right).

• •    Взаимосвязь суммарного испарения и температуропроводности для участка с внесением биоугля (ETбиоуголь ~ Dбиоуголь):

ETбиоуголь = 1.03178 × Dбиоуголь – 0.08497 (R2= 0.56, p-value = 0.032)

• •    Взаимосвязь суммарного испарения и температуропроводности для контрольного участка (ETконтроль ~ Dконтроль):

ETконтроль = 1.3371 × Dконтроль – 0.1475 (R2= 0.66, p-value = 0.017)

• •    Взаимосвязь объемной влажности почвы и тепловой инерции для участка с внесением биоугля (Wбиоуголь ~ ТИбиоуголь):

Wбиоуголь =  0.0001654  × ТИбиоуголь +  0.0533856 (R2 =  0.38, p-value = 0.08)

• •    Взаимосвязь суммарного испарения и тепловой инерции для участка с внесением биоугля (ET биоуголь ~ ТИ биоуголь ):

ETбиоуголь = 0.0004479  × ТИбиоуголь – 0.0555895  (R2 = 0.59, p-value = 0.026)

• •    Взаимосвязь суммарного испарения и тепловой инерции для контрольного участка (ET контроль ~ ТИ контроль ):

ETконтроль = 0.0005430  ×  ТИконтроль – 0.1361763  (R2 = 0.81, p-value = 0.003).

Данные регрессионные модели могут найти применение при определении влажности пахотного горизонта почвы на начало вегетации и суммарного испарения, данные величины являются крайне важными при моделировании роста и развития растительности (Steduto et al., 2009) . Однако для построения более полных и точных моделей необходимо провести дополнительные полевые и лабораторные измерения с ненарушенными почвенными образцами во всем диапазоне изменения влажности исследуемой почвы, так как взаимосвязь теплофизических свойств с влажностью почвы и суммарным испарением не является линейной, если рассматривать весь диапазон изменений влажности почвы (Qin et al., 2013) .

ВЫВОДЫ

В целом анализ показал:

• 1.    Внесение биоугля в условиях влажности почвы в абсолютно сухом состоянии (объемная влажность почвы W = 0.01 м³ м³) уменьшило на 29.7% (p < 0.05) теплопроводность, уменьшило на 18.5% (p < 0.05) объемную теплоемкость, уменьшило на 13.7% (p < 0.05) температуропроводность, уменьшило на 24.3% (p < 0.05) тепловую инерцию пахотного горизонта супесчаной дерново-подзолистой почвы.

• 2.    Внесение биоугля при значениях влажности почвы, соответствующей наименьшей влагоемкости (объемная влажность почвы W = 0.26 м³ м³) увеличило на 9.4% (p < 0.05) теплопровод-

• ность, уменьшило на 2.6% (p < 0.05) объемную теплоемкость, увеличило на 12.3% (p < 0.05) температуропроводность, увеличило на 3.2% (p < 0.05) тепловую инерцию пахотного горизонта супесчаной дерново-подзолистой почвы.

• 3.    Внесение биоугля незначимо (p > 0.1) уменьшило радиационный баланс и поток тепла в почву, значимо (p < 0.05) увеличило на 35.5% турбулентный поток тепла, значимо (p < 0.05) уменьшило на 17.0% скрытый поток тепла, значимо (p < 0.05) уменьшило на 13.9% суммарное испарение.

• 4.    Внесение биоугля значимо (p < 0.05) на 6.4% увеличило температуру подстилающей поверхности.

• 5.    Внесение биоугля не оказало значимых (p > 0.1) изменений на температуру листьев растений.

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