Scientific journal
Advances in current natural sciences
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,775

MODELING OF MOISTURE INTERPROFILE MOVEMENT IN CHERNOZEMS E OF AGROLANDSCAPES OF UIMON INTERMOUNTAINE DEPRESSION (KATUN BASIN, CENTRAL ALTAI)

Baboshkina S.V. 1 Puzanov A.V. 1 Elchininova O.A. 1 Rozhdestvenskaya T.A. 1 Kuznetsova O.V. 1 Balykin S.N. 1 Troshkova I.A. 1
1 Institute for Water and Environmental Problems SB RAS
Water retention curves (WRC) and hydraulic conductivity functions of ordinary chernozems of Uimon intermountain depression (r. Katun basin, Central Altai) were obtained, using semi-empirical calculation methods of pedotransfer functions. Approximation parameters of WRC by Van Genuhten`s equation were determined in RETC program, on the basis of experimental data: particle size distribution of the soil, its density, wilting point and field capacity. The boundaries of soils grain-size fractions determined by the method of Kachynskii, for use in the RETC, were made equal to the International FAO classification by the linear interpolation method. Using a computer program HYDRUS-1D (authors: J. Simunek, M.Th. Van-Genuchten, M. Seina) the modeling of moisture vertical movement in agrochernozems was executed and volumes of moisture reserves reduce due to its infiltration into the lower layers during the dry 10-days period with root water uptake factor were calculated. It is shown that in all soil layers of arable land chernozem wet faster and the general moisture loss from its upper horizon after initial atmospheric wetting up to 10 days without rainfall are more significant (than in the soil of hayfields), that determine better conditions for leaching, during good wetting of A gorizont, which represent by pH change to acidic condition. Vertical outflow of moisture to the deeper layers of the haying soil prevent by high density of B and C horizons with alkaline reaction and less coefficient of hydraulic conductivity. The simulation of moisture movement with the root plants uptake consumption showed that transpiration factor corrects the dynamics of moisture in chernozem of haymaking largely, than in chernozem of plowed field.
chernozem
Uinon intermontane depression
soil moisture
water retention curve
RETC
HYDRUS-1D
hydraulic conductivity
moisture vertical distribution
Van Genuhten function
infiltration
transpiration

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

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

Процесс поступления в почву воды и ее последующее распределение является важной составляющей режима любой почвы. В почвах сухих степей, где практически единственным источником поступления воды является атмосферная влага, особый интерес представляет процесс впитывания воды и ее дальнейшее перераспределение в ненасыщенной влагой сухой почве (инфильтрация) [8]. Количественной характеристикой водоудерживающей способности почв является основная гидрофизическая характеристика (ОГХ), как зависимость между капиллярно-сорбционным давлением влаги (Рк-с) и влажностью. Вертикальное движение влаги в ненасыщенных влагой степных почвах характеризуется функцией влагопроводности – зависимостью между коэффициентом влагопроводности (Квл) и Рк-с [8].

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

Материалы и методы исследования

Объектами данного исследования были выбраны черноземы обыкновенные Уймонской котловины – одного из межгорных понижений Центрального Алтая (бассейн р. Чендек, левый приток р. Катунь). В настоящее время почвенный покров котловины представлен в основном черноземами обыкновенными распаханными. Рассматриваемые нами агроландшафты занимают выровненный и топографически однородный участок с уклоном поверхности около 0 °. Первый вариант – пашня, посевы овса, стадия – выход в трубку. Второй вариант – залежь, разнотравно-злаковая ассоциация. На этих двух полях были заложены почвенные разрезы на расстоянии не более 100 м друг от друга. Отбор проб проводился 8 июля 2014 года. По архивным данным сайта www.rp5.ru 5 июля 2014 года шли ливневые дожди и выпало в сумме 14 мм осадков, а с 6 по 8 июля выпало еще 3,5 мм осадков. На момент наблюдений верхний горизонт изучаемых нами черноземов был отчетливо более влажный, чем нижние слои.

Полевые и лабораторные исследования были проведены общепринятыми стандартизированными методами. Полевая влажность определялась термостатно-весовым методом, плотность почвы – методом режущего кольца, плотность твердой фазы – пикнометрически. Гранулометрический состав почв определяли пипеточным методом по Качинскому. Общую пористость определяли через соотношение плотности сложения почвы и плотности ее твердой фазы [1]. Определение наименьшей влагоемкости в поле проводили через 24 ч после окончания заливочных экспериментов на водопроницаемость [1, 9]. Содержание гумуса в почве определяли методом мокрого сжигания по Тюрину, рН – потенциометрически [2].

Для аппроксимации функции ОГХ (основной гидрофизической характеристики) агрочерноземов использовали наиболее широко известное в настоящее время в мировом почвоведении уравнение Ван Генухтена:

bab01.wmf,

где m = 1 – 1/n, θ – равновесная влажность (влажность, соответствующая определенному давлению почвенной влаги P, θs – влажность, близкая к влажности насыщения, θr – остаточная влажность, α – величина, обратная величине капиллярно-сорбционного давления, приближающегося к давлению входа воздуха, n – крутизна кривой [8, 9]. Восстановление ОГХ проводилось полуэмпирическим расчетным методом с помощью педотрансферных функций (ПТФ) – зависимостей, позволяющих восстановить ОГХ по традиционно определяемым базовым свойствам почв (как правило, это гранулометрический состав, а также плотность, порозность, содержание органического вещества и др.) [8]. В данной работе использовались ПТФ Rosetta Lite 1.1, заложенные в программе RETC [10]. Для получения параметров аппроксимации ОГХ функцией Ван-Генухтена использовались экспериментально полученные данные гранулометрического состава, плотности, порозности, наименьшей влагоемкости и влажности завядания. Границы фракций гранулометрического состава были приведены к границам международной классификации FAO методом линейной графической интерполяции [7].

Процесс вертикального перемещения влаги из черноземов обыкновенных Уймонской степи различного сельскохозяйственного использования в ситуации первоначального насыщения только верхних горизонтов влагой описан с помощью системы компьютерного моделирования влагопереноса в почвах HYDRUS-1D (США, авторы: J. Simunek, M.Th. Van-Genuchten, M. Seina, адрес: www.hydrus2d.com). Расчеты проводили для почв, каждый горизонт которых обладал своими водно-физическими свойствами. Гидрологический режим черноземов моделировали для периода в 10 (июльских) дней без последующих атмосферных осадков. Верхнее граничное условие задавалось как «Постоянный поток», а на нижней границе принималось «Постоянное давление». Моделирование отрицательного (нисходящего) передвижения влаги проводилось с учетом поглощение воды корнями растений. В этом случае выбрали подмодель этого процесса в Main process, а в команде The Root Water Uptake Parameters – параметры корневого потребления для определенного вида культуры из встроенной базы данных: в случае пашни под посевами овса – «Corn – vegetative period», в случае сенокоса – «Grass». Распределение плотности корневых систем строилось в графическом редакторе (Soil Profile – Graphical Editor). Далее, в команде Water Flow (Constant Boundary) на верхней границе задавали интенсивность потенциальной транспирации [7] – максимальное испарение, которое может происходить при данных погодных условиях над конкретным растительным покровом сельскохозяйственной культуры, находящейся в оптимальных агрономических условиях (без стрессов). Величину Root water uptake 1,4 мм/сут. мы принимали, исходя из сравнений с величинами потенциальной транспирации (а также эвапотранспирации, испарения) и расчетов через величины урожайности и коэффициенты водопотребления, приводимыми в работах других авторов [3–6].

Результаты исследования и их обсуждение

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

Таблица 1

Физико-механические и физико-химические показатели агрочерноземов Уймонской котловины

Объект

Гори-зонт

Гранулометрический состав ( %) и классификация по Качинскому

Гранулометрический состав ( %) и классификация FAO

Плот-ность, г/см3

Пороз-ность,

%

Содер-жание гуму-са, %

рН

Ил, < 0,001

Физич. глина, < 0,01

Принадлежность

Clay < 0,002

Silt 0,002–0,05

Sand > 2

Принадлежность

Чернозем (сенокос)

Ад

7,1

29,6

л-сугл.

12,0

51,1

36,9

Silt loam

0,88

62,9

4,7

7,1

В

12,2

40,8

ср.-сугл.

19,8

42,6

37,6

Loam

1,26

51,1

2,1

8,0

С

16,2

38,1

ср.-сугл.

21,0

40,2

38,8

Loam

1,36

48,0

1,1

7,9

Чернозем (пашня)

Ап

7,2

34,5

ср.-сугл.

15,3

45,0

39,7

Loam

1,07

58,7

9,5

6,8

В

7,7

27,0

л-сугл.

13,0

41,9

45,1

Loam

1,14

55,4

5,4

6,6

С

4,8

16,9

супесь

10,0

33,4

56,6

Sandy loam

1,17

57,8

3,6

7,7

Для расчетов и моделирования процесса движения влаги в HYDRUS 1D было решено применить параметры аппроксимации ОГХ функцией Ван Генухтена и коэффициенты фильтрации, полученные по ПТФ Rosetta программы RETC с использованием данных гранулометрического состава (приведенные к международной классификации FAO методом линейной интерполяции) (SSC), плотности (BD), наименьшей влагоемкости (TH33), влажности завядания (TH1500) (табл. 2–3).

Таблица 2

Гидрофизические свойства чернозема обыкновенного под залежью. Параметры аппроксимации ОГХ функцией Ван Генухтена и коэффициенты фильтрации

Горизонт (глубина, см)

θs, см3/см3

θr, см3/см3

n

α

Кф, см/сут

Ад (2–16)

0,500

0,0255

1,3196

0,0245

161,2

В (16–26)

0,4248

0,0365

1,4145

0,0112

31,4

С (26–40)

0,4112

0,0556

1,4228

0,0165

23,72

Таблица 3

Гидрофизические свойства чернозема обыкновенного под пашней. Параметры аппроксимации ОГХ функцией Ван Генухтена и коэффициенты фильтрации

Горизонт (глубина, см)

θ, см3/см3

θr, см3/см3

n

α

Кф, см/сут

Ап (0–26)

0,459

0,0354

1,5362

0,0054

88,26

В (26–40)

0,4427

0,0283

1,3922

0,013

76,95

С (40–55)

0,4241

0,0193

1,4281

0,010

84,67

Кривые водоудерживания изученных нами суглинистых черноземов Уймонской степи, как правило, имеют выраженную S-образность и наиболее дифференцированы в области капиллярной и гравитационной форм влаги (рис. 1).

baboh1.tif

Рис. 1. Кривые ОГХ различных горизонтов агрочерноземов Уймонской степи различного режима использования – под пашней (посевы овса) и под залежью (сенокос, разнотравье)

В состоянии, близком к насыщению (Р < – 50), наибольшей влагоудерживающей способностью отличается верхний гумусово-аккумулятивный горизонт чернозема под залежью – наиболее легкий по гранулометрическому составу и наиболее рыхлый. В области капиллярной (pF = 1,7–3) и пленочно-капиллярной малоподвижной влаги (pF = 3–4,5) наибольшей влагоудерживающей способностью отличается среднесуглинистый, наиболее гумусированный и довольно пористый верхний Ап горизонт пашни. Наилучшими влагоудерживающими свойствами в области пленочной и адсорбированной форм влаги (pF > 4,5) обладает один из наиболее тяжелых по гранулометрическому составу и наиболее плотный Ск горизонт под сенокосом, что объясняется увеличенным содержанием в нем мелких пор, которые удерживают влагу при низких значениях давления [8].

В области полного насыщения почвы влагой наибольшей влагопроводностью (и наибольшим коэффициентом фильтрации) отличается наиболее легкий по гранулометрическому составу и наиболее рыхлый Ад горизонт почвы под сенокосом (рис. 2). При этом значительно более плотные В и С горизонты этого варианта чернозема отличаются самыми низкими влагопроводящими свойствами, что, очевидно, будет затруднять вертикальное нисходящее движение влаги в этих почвах. В области от капиллярного механизма удержания влаги (P > – 10 см водн. ст.) лучшими влагопроводящими свойствами отличается чернозем, используемый под посевы овса (пашня). Почва под сенокосом в условиях иссушения будет проводить влагу заметно хуже.

Для расчета движения влаги в почве в качестве начальных условий взяли природную полевую ситуацию, когда, на момент отбора проб (на 8 июля 2014 года) влажность верхних горизонтов была порядка 30 %, а подпахотных В горизонтов – 16–18 %. Судя по ОГХ (график (рис. 1) и матрица из программы RETC), в рассматриваемых почвах эти значения влажности соответствовали значениям давления порядка – 270 – – 290 см водн. ст. в верхних горизонтах и – 640 – – 660 см водн. ст. в подпахотных В горизонтах. Однако, при практически одинаковых величинах капиллярно-сорбционного давления более суглинистый, гумусированный и более плотный Ап горизонт почвы под посевами овса изначально будет содержать больше объемной влаги (0,325 см3/см3), т.к. отличается большей влагоемкостью (рис. 1, 4).

Характер и динамика распределения влаги в рассматриваемых условиях в агрочерноземах различного сельскохозяйственного использования Уймонской котловины заметно различаются (рис. 3, а–б). Можно отметить, что в почве пашни за первые двое суток за счет большего Квл нижележащих горизонтов В и С происходит довольно быстрое промачивание всех слоев почвы (рис. 3, а) – нисходящее движение влаги в почве пашни по макропорам и трещинам, очевидно, происходит более интенсивно, чем в почве под залежью.

baboh3a.tif baboh3b.tif

а) чернозем под пашней б) чернозем под залежью

Рис. 3. Динамика профильного распределения влажности (см3/см3) в черноземах обыкновенных различных агроландшафтов Уймонской котловины (результат моделирования за 10-дневный период без осадков, с учетом корневого потребления растениями)

В последующие дни (в отсутствие дополнительного атмосферного увлажнения) нижележащие горизонты В и С почвы под пашней полученную из верхнего горизонта влагу постепенно теряют, происходит постепенное иссушение всех слоев почвы, причем до значений ниже исходной величины. Что касается чернозема под залежью, то, несмотря на более легкий гранулометрический состав верхнего Ад горизонта, в первый день влага не проникает на всю глубину почвенной толщи, а только до глубины 32 см (рис. 3, б). Очевидно, вертикальному внутрипочвенному стоку препятствуют существенно более низкие Квл в менее пористых и более плотных В и С горизонтах (рис. 2).

baboh2.tif

Рис. 2. Функции влагопроводности для А, В и С горизонтов черноземов агроландшафтов Уймонской степи

В целом за 10 суток, последующих после дождей, с учетом корневого потребления влаги растениями, почва под пашней теряет из верхнего Ап горизонта глубиной 26 см 0,114 см3/см3 влаги или 29,5 мм водного слоя что на 5,5 мм больше в случае моделирования простого вертикального внутрипочвенного оттока, без учета транспирации. Объемная влажность изменяется в Ап горизонте с учетом корневого потребления от 0,325 до 0,211 см3/см3 или (на весовую влажность) от 30,4 % до 22,5 %. Потери влаги в почве под сенокосом, особенно в Ад горизонте – более рыхлом и легком по гранулометрическому составу и с более плотным распределением корней, оказались более существенным, по сравнению с простым вертикальными оттоком, без учета транспирации, что особенно заметно по Ад и В горизонтам. За 10 суток влажность Ад горизонта с учетом корневого потребления изменится от 0,279 см3/см3 до 0,183 см3/см3 (или с 31,8 % до 16,1 %), влагопотери составят 0,096 см3/см3 или 15,4 мм водного слоя (для 16 см почвенной толщи), что на 6,8 мм (почти в 2 раза) больше, если корневое потребление при моделировании не учитывать.

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

Заключение

Моделирование вертикального нисходящего движения влаги в черноземах обыкновенных Уймонской котловины в ситуации начального атмосферного увлажнения верхних горизонтов и корневого потребления влаги растениями с применением пакета HYDRUS-1D показало, что из двух сравниваемых участков (под залежью и под посевами овса) в черноземе пашни промачивание всех слоев почвы более быстрое и происходит уже в первые сутки. Полученная во время атмосферных осадков влага сравнительно быстро двигается вниз по почвенным порам пахотного варианта чернозема. В результате в почве пашни формируются более благоприятные условия для выщелачивания, что находит свое отражение, например, в снижении уровня рН до нейтральных величин, по сравнению с почвой под залежью. В почве под разнотравьем, в которой рН существенно выше и наблюдается щелочная реакция среды всех горизонтов, вертикальному оттоку влаги в более глубокие слои препятствуют уплотненные В и С горизонты, кроме того, значительная доля поступившей в нее атмосферной влаги поглощается корнями растений.