Бесплатный автореферат и диссертация по сельскому хозяйству на тему

Преимущественные потоки почвенной влаги в условиях стационарной фильтрации

ВАК РФ 06.01.03, Агропочвоведение и агрофизика

Автореферат диссертации по теме "Преимущественные потоки почвенной влаги в условиях стационарной фильтрации"

На правах рукописи

Ананьева Ксения Аркадьевна

ПРЕИМУЩЕСТВЕННЫЕ ПОТОКИ ПОЧВЕННОЙ ВЛАГИ В УСЛОВИЯХ СТАЦИОНАРНОЙ ФИЛЬТРАЦИИ

Специальность 06.01.03 - агропочвоведение, агрофизика

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата биологических наук

Москва - 2004

Работа выполнена на кафедре физики и мелиорации почв факультета почвоведения Московского государственного университета им. М.В. Ломоносова Научный руководитель: доктор биологических наук

профессор Е.В. Шеин

Официальные оппоненты: доктор биологических наук

АИ.Поздняков кандидат биологических наук В.Ф.Басевич

Ведущее учреждение: Владимирский научно-исследовательский институт сельского хозяйства

Защита диссертации состоится_в 15.30 в аудитории М-2 на

заседании диссертационного совета К501.001.04 в МГУ им. М.В. Ломоносова.

С диссертацией можно ознакомиться в библиотеке факультета почвоведения МГУ

Автореферат разослан «_»_2004 г.

Приглашаем Вас принять участие в обсуждении диссертации на заседании Диссертационного Совета в МГУ, а отзывы на автореферат в двух экземплярах, заверенные печатью, просим направлять по адресу: 119992, ГСП-2, Москва, Ленинские горы, МГУ, факультет почвоведения, Ученый Совет.

Ученый секретарь диссертационного совета

Богатырев Л.Г.

Актуальность темы. В настоящее время уделяется большое внимание исследованиям так называемых преимущественных путей движения влаги (от англ. - preferential flow, fingers). Образование их связывают с неравномерностью движения (в песчаных почвах), со слоистостью почвенного профиля, с наличием гидрофобных прослоек, макропор различной природы. Значительно меньше изучены вопросы образования преимущественных путей фильтрации влаги, связанных с неоднородностью пространственного распределения физических свойств. Одним из основоположников такого рода работ был профессор МГУ Е.А.Дмитриев. В его работах приведены качественные результаты по преимущественным путям миграции, обусловленным в основном генетическими особенностями почв. Количественная взаимосвязь формирования преимущественных путей фильтрации и пространственного распределения физических свойств в масштабе почвенного индивидуума не была изучена.

Цель работы: изучить закономерности формирования преимущественных путей фильтрации воды в серых лесных почвах Владимирского ополья.

Задачи:

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

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

3. Выявить закономерности фильтрации воды в серых лесных почвах в связи с пространственной неоднородностью плотности.

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

Практическая значимость. Представлены количественные характеристики форми-

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

Апробация. Основные результаты работы были доложены на конференциях студентов и аспирантов «Ломоносов-2000» (Москва, МГУ, 2000 г.), Ш съезде почвоведов (Суздаль, 2000), на 2-й международной конференции по теме "Подпочвенное уплотнение" (Godollo, Венгрия, 2000).

Публикации. По теме диссертации опубликовано 6 работ.

Структура и объем работы. Диссертационная работа изложена на_страницах,

включает_таблиц и рисунков; состоит из введения,_глав, выводов, списка

литературы, включающего_наименований, и приложения.

Глава 1. СОСТОЯНИЕ ПРОБЛЕМЫ

К настоящему времени многочисленными исследователями показано, что вода в почве движется специфическими путями за счет разнообразных механизмов. Все специфические процессы переноса влаги в почве можно весьма условно разделить на 3 характерные группы гидрологических почвенных явлений:

1) формирование отдельных "водных каналов", линий преимущественного потока, "пальчатых структур" (англ. - preferential flow, fingers);

2) быстрый перенос по "проводящим11 зонам порового пространства (макропорам, трещинам) с последующим обменом с "застойными" зонами;

3) неравномерность переноса влаги, связанная с пространственной вариабельностью почвенных свойств.

Формирование "водных каналов" отмечали, как правило, в слоистых почвах. Исследователи подчеркивали, что существование слоев в почве является одной из причин неравномерности фронта увлажнения, возникновения «пальчатых структур» движения влаги (Качинский НА., Дмитриев ЕЛ., Зайдельман Ф.Р., Hillel D. и Baker R.S.,Van.Ommen H.C с соавт. и др.): нестабильность потоков воды возникает в двухслойных (гидрофобный/гидрофильный, насыщенный/сухой, высоко/плохо фильтрующий слои) почвах (Flury М. с соавт., Ritsema C.J. с соавт., Steenhuis T.S. с соавт.). Большинство подобных явлений отмечалось в песчаных почвах, в которых сами слои были гомогенными, но заметно отличались по гидрологическим свойствам.

Быстрый перенос влаги по "проводящим" и последующее медленное ее перераспределение в "застойных" зонах связывают с развитым межагрегатным (трещинным) и агрегатным поровым пространством. Этот раздел специфического почвенного переноса включает такие физические процессы, как набухание/усадку педов (агрегатов), обмен веществом межпедного и внутрипедного пространства. Этот тип переноса приурочен, как правило, к структурным почвам. В большом количестве работ (Дмитриев ЕА., Зайдельман Ф.Р., Шеин Е.В. с соавт. и др.) подчеркивалась специфика быстрого переноса влаги по макропорам в виде: 1) отсутствия адсорбции ионов (Skopp J.); 2) стабильности существования макропор (в основном, биопор) (Шеин Е.В. с со-авт.) и др.

Заметно меньшее внимание уделяли взаимосвязи образования отдельных "водных тяжей" и пространственной изменчивости физических свойств почв в масштабе почвенного индивидуума. Здесь следует отметить цикл работ Е.А.Дмитриева, в котором на качественном уровне показана высокая неоднородность путей движения влаги, возникновение преимущественных путей фильтрации. Однако в этих работах не изучались физические закономерности формирования неоднородного водного потока, наличие взаимосвязей путей фильтрации с физическими свойствами почв.

Глава 2. ОБЪЕКТЫ И МЕТОДЫ ИССЛЕДОВАНИЯ

Изучение путей миграции влаги проводились на территории Владимирского НИИСХ во Владимирском ополье. Для данной территории характерна высокая комплексность почвенного покрова при доминировании серых лесных почв. Учитывая различную степень оподзоленности серых лесных почв, мы руководствовались подробной классификацией Е.А.Дмитриева, используемой при исследовании серых лесных почв Владимирского ополья по длинномерным траншеям, в которой выделялись следующие градации по оподзоленности: мелко оподзоленные - мощность оподзо-ленной толщи не более 20 см (сЛ20П1), неглубоко оподзоленные - мощность 20-50 см (сЛ2оп2), глубоко оподзоленные - мощность более 50 см (сЛ2оп3). Выделялись также почвы со вторым гумусовым горизонтом (ВГТ) и остаточно-карбонатные.

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

серой лесной неглубоко оподзолениой почвы со 2-м гумусовым горизонтом (ВГТ) сЛ20п2'ВГГ и 6 серых лесных почв разной степени оподзоленности без BIT.

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

Полевые исследования пространственной картины миграции воды проводились с использованием крахмальной метки (2%-ый раствор водорастворимого крахмала). На поверхности почвы устанавливалась рама и фильтровалась вода в течение 6-8 часов. Затем, сразу после фильтрации воды, "ступенькой" подается меченная растворимым крахмалом вода; слой на поверхности, как правило, составлял 5 см. Через 10-12 часов послойно (через 5-10 см) на горизонтальных площадках отмечались пути фильтрации воды по появлению синих пятен при увлажнении поверхности йодной водой. Морфология и характерное расположение зон фильтрации крахмала зарисовывались с помощью полевого почвенного пантографа в масштабе 1:5. После зарисовки. основных путей фильтрации влаги на площадках размером не менее 1 м2 с шагом 10 см по сетке определялась влажность и плотность почвы (объем бура 4 см3). На рис. 1 представлены, статистики. распределения плотности по глубине для всех объектов. Все распределения различаются. Характерным является проявление так называемой «плужной подошвы» - увеличение медианных значений плотности при изменении параметров варьирования на глубине 20-30 см. Особое место занимаем серая лесная почва с BIT. Для нее характерно отсутствие «плужной подошвы»: очень низкие значения плотности отмечаются на глубине 35 см. Отметим также высокую вариабельность плотности, уменьшающуюся от поверхностных, пахотных, к глубинным слоям.

В лабораторных исследованиях для получения основной гидрофизической характеристики (ОГХ) использовался метод тензиостатов. Для построения "верхней" части ОГХ использовался метод десорбции влаги над насыщенными растворами солей: КС1 (pF=5,3) и K2SO4 (pF=4,4)) (Вадюнина АФ., Корчагина ЗА, 1986).

Рис.1 Статистики распределения плотности серых лесных почв по глубине: а) мелко оподзоленной; б) неглубоко оподзоленной с ВГТ; в) мелко оподзоленной; г) неглубоко оподзоленной; д) серой лесной; е) мелко оподзоленной; ж) подпахотных ' "■'■"*" • горизонтов неглубоко оподзоленной

глубина, см •

Глава 3 РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЯ

На рис. 2 приводится пример формы основных путей миграции по крахмально му окрашиванию для серой лесной мелко оподзоленной почвы. Движение воды происходит практически вертикально. Пятна крахмального окрашивания располагаются друг под другом, уменьшаясь с глубиной по площади. Отчетливо видны "водные тяжи". Растекание по площадке небольшое, большая часть крахмального окрашивания. приходится на область внутри рамы. Однако, в большинстве исследованных почв растекание более значительно, и на глубине 15-20 см основное окрашивание приходится на область вне рамы.

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

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

1. Показатель удельной площади различных категорий плотности (51). Для этого сначала были выделены пять условных категорий: <1,3, 1,3-1,4, 1,4-1,5, 1,5-1,6 и >1,6 г/см3. На каждой площадке с помощью процедуры кригинга (программа Surfer) прорисовывались изоплеты плотности. Удельная площадь каждой категории рассчитывается как процент, занимаемый данным диапазоном ко всей исследованной площади площадки:

2. Показатель вертикального соответствия категорий плотности (52). Рассчитывается как процент площади совпадения соответствующей категории плотности верхней и нижней площадок к площади категории на верхней площадке. Указывает на «вертикальность» пространственного расположения различных зон по плотности: чем он выше, тем лучше совпадение зон одинаковой плотности, расположенных одна под другой.

3. Показатель переноса влаги по соответствующей категории плотности (53). Рассчитывается как процент площади соответствующей категории плотности, приходящейся на отмеченную по крахмальному окрашиванию зону фильтрации влаги. Указывает на долю соответствующей категории плотности в переносе влаги на данной глубине: чем он выше, тем больше участие данной категории в переносе воды.

53 = ^1-100% Бокр

4. Показатель «вертикальности» движения влаги (54). Рассчитывается аналогично показателю (52), предложенному для характеристики вертикальности расположения соответствующих зон плотности, однако для зон фильтрации влаги, отмеченных по крахмальной метке. Выражается в процентах площади пересечения зон фильтрации на верхней и нижерасположенной площадках исследования к площади зоны фильтрации вышерасположенной площадки: чем он выше, тем лучше совпадение зон фильтрации, расположенных одна под другой.

5. Показатель переноса влаги по соответствующей категории перепадов плотности (5п). Для этого сначала были выделены 4 условные категории перепада плотности: <0,2, -0,2-0,0-0,2 и >0,2 г/см3. Удельная площадь каждой категории рассчитывается как процент, занимаемый данным диапазоном ко всей исследованной площади площадки. Сам показатель переноса влаги по соответствующей категории перепадов плотности рассчитывается как процент площади соответствующей категории, приходящейся на зону фильтрации влаги, отмеченную по окрашиванию крахмала. Этот показатель описывает взаимосвязь распределения перепада плотности и крахмального окрашивания: чем он выше, тем больше вклад категории перепада плотности в проведении воды в профиле.

6. Коэффициент ветвления (5ветв). Рассчитывается как отношение количества пятен окрашивания на площадке, умноженное на ее площадь, к общей площади окрашивания. Этот показатель оценивает влияния слоистости профиля на формирование отдельных зон фильтрации (зон преимущественного потока): чем он выше, тем больше отдельных «языков» сформировалось в процессе движения воды вниз по профилю.

Рассмотрим распределение всех этих показателей по глубине на примере серой лесной мелко оподзоленной почвы и серой лесной неглубоко оподзоленной почвы с ВГП

Для распределений 51 общим является то, что с глубиной неоднородность по плотности уменьшается, что подтверждается приведенными выше рисунками.

В профиле серой лесной мелко оподзоленной почвы (рис. За) в поверхностном слое (глубина 5 см) присутствуют все категории, но их вклад в общее распределение различен. Наибольшую площадь занимает категория <1,3 г/см3 (36%), одинаковая площадь представлена средними категориями плотности - 1,3-1,4 г/см3 и 1,4-1,5 г/см3 (около 23%), наименьшая площадь приходится на самую плотную категорию - >1,6 г/см3, несмотря на то, что она обладает наибольшим вертикальным соответствием (рис. 36). С глубиной отмечается значительное колебание площадей, занимаемых всеми категориями плотности. Но можно отметить несколько особенностей распреде-

ления плотности. С глубиной закономерно увеличивается вклад относительно плотной категории - 1,5-1,6 г/см1 (начиная с глубины 25 см), а площадь, занимаемая наиболее плотной и рыхлыми категориями, стремится к нулю. Кроме того, на глубине 15 см отмечается максимум для самой плотной категории - >1,6 г/см3.

Распределение «S3 носит похожий характер (рис. Зв). Отмечается значительное колебание ее с увеличением глубины, в глубинных слоях увеличивается доля достаточно плотной категории — 1,5-1,6 г/см3, а доли остальных категорий составляют менее 10%. И "наследуется" эта категория (показатель 52), начиная с глубины 25 см,

Рис. 3. Распределение пространственных показателей фильтрации по глубине серой лесной мелко оподзоленной почвы (пояснения в тексте)

в наибольшей степени (рис. 36).

Что касается вертикальности движения воды - показатель 54 (рис. Зг), то максимум его приходится на глубину 30 см (40,8%). В целом по профилю колебание этого показателя не очень большое (от 16,2-28,9%), что свидетельствует о достаточно быстром формировании "водных тяжей" и сохранении их положения с глубиной.

Для остальных исследованных почв без BIT сохраняются те же закономерности в распределении этих показателей. Остановимся подробнее на закономерностях формирования преимущественных потоков в серой лесной почве с BIT.

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

а б

Рис. 4. Основные пути миграции влаги по крахмальному окрашиванию (а) и распределение диапазонов плотности (б) по глубине серой лесной почвы с ВГГ (обозначения см. рис. 2)

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

В профиле в поверхностных слоях присутствуют все категории (рис. 5а), занимая различные площади, причем наименьшую площадь занимает категория >1,6 г/см1, а наибольшую - 1,4-1,5 г/см3 (показатель 81). Резкая смена ситуации наблюдается на глубине 20 см: с глубиной отмечен быстрый рост площади, занимаемой самой рыхлой категорией, а площади остальных категорий постепенно сводятся к 0. Соответствен-

б.

глубина, см глубина, см

—♦—<1,3 -»-1,3-1,4-*-1,4-1,5 -*-1,5-1,6 -*->1,б

Рис. 5. Распределение пространственных показателей фильтрации по глубине серой лесной неглубоко оподзоленной почвы с ВГГ

но, и распределение по глубине доли этой категории плотности в проведении влаги (показатель 83) носит аналогичный характер: с глубины 20 см доминантное значение приобретает самая «рыхлая» категория (рис. 5в). Слабо «растекается» влага в этой почве: 84 изменяется от 56 до 48% (рис. 5г). Связано это, по-видимому, со строгой горизонтальной ориентацией доминирующей градации плотности, в особенности после глубины 20 см пахотного слоя (рис. 56).

Проведенный анализ количественных показателей распределения плотности и формирования преимущественных потоков позволяет утверждать:

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

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

3. Серая лесная почва со 2-м гумусовым горизонтом занимает особое положение среди исследуемых почв. Для нее характерно наличие "точки резкой смены ситуации" (глубина 20-25 см), что соответствует появлению 2-ого гумусового горизонта и резкого уменьшения плотности почвы.

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

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

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

1) не существует заметно выделяющихся слоев по плотности для данных почв;

2) не наблюдается однообразного распределения межслойных перепадов плотности, даже в почвах единого генезиса;

3) во всех случаях отмечается наличие экстремума (максимума или минимума) значений на глубине 20-30 см (в зависимости от почвы), что, по-видимому, связано с формированием "плужной подошвы" и, как следствие, смены условий.

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

На рис. 7 представлены графики распределения этого показателя по глубине на примере серой лесной мелко оподзоленной почвы и серой лесной неглубоко оподзо-ленной почвы с ВГГ.

В серой лесной мелко оподзоленной почве (рис. 7а) в поверхностном слое доминируют две категории перепадов плотности - <-0,2 и -0,2-0. С глубиной основной в проведении влаги становится категория 0-0,2, а затем - -0,2-0. Значение же категории <-0,2 в проведении воды с глубиной резко падает и не превышает в дальнейшем 1%. Категория >0,2 практически не участвует в этом процессе (значения Sn по профилю для нее не превышают 0,7%).

В серой лесной почве неглубоко оподзоленной почве с ВГГ (рис. 76) распреде-а б

5 15 25 35 45 5 15 25 35

глубина, см глубина, см

Рис. 7. Распределение перепадов плотности (5п) и коэффициента ветвления (5ветв) для серых лесных почв: а) мелко оподзоленной; б) неглубоко оподзоленной с ВГГ

ление Sn иное. В верхнем слое доминирующей категорией является -0,2-0, значительно меньше вклад в проведении влаги категории 0-0,2. С глубиной основная проводящая категория постоянно меняется: на 15 см примерно одинаково значение 8п для всех категорий, кроме <-0,2, глубже 20 см остаются только две категории (>0,2 и 0-0,2), которые сменяют друг друга в лидерстве проведения влаги.

Следует отметить, что значения 8п редко превышают 10-15% (максимальное значение - около 35% для ¡сЛ2вгг на глубине 10 см). Можно сделать вывод, что перепад плотности в исследуемых условиях не оказывает значительного влияния на передвижение воды в толще почвы в виде преимущественных потоков влаги. Но в той части, где градиент плотности совпадает с крахмальным окрашиванием, доминирующими являются средние категории - -0,2-0 и 0-0,2 г/см3, т. е. фильтрация осуществляется по тем участкам, которые имеют либо одинаковую по вертикали плотность, либо плотность должна отличаться в соседних горизонтах не более чем на 0,2 г/см3.

Для оценки влияния слоистости профиля на формирование отдельных зон фильтрации (зон преимущественного потока) нами был рассчитан коэффициент ветвления (5ветв). Полученные зависимости распределения 5ветв по глубине представлены на рис. 7. Эти графики для большинства объектов не совпадают, но общим является то, что почти все они имеют один максимум на глубине 20-30 см в зависимости от почвы.

Максимальные значения различаются достаточно сильно. В серой лесной мелко оподзоленной почве (рис. 7а) отмечается равномерное возрастание коэффициента ветвления с глубиной до глубины 25 см (максимальное значение - 107,2%), а затем его постепенное уменьшение. Внизу профиля 5ветв постоянен.

Наименьший "пик" обнаруживается в серой лесной почве неглубоко оподзо-ленной почве с ВГТ (рис. 66) - около 15% (кроме того, в этой почве большое значение 5ветв отмечается с поверхности - 22%). Для этой почвы распределение имеет слабо выраженный синусоидальный характер.

При попытке сопоставить распределение коэффициента ветвления с распределением перепадов плотности можно сделать вывод об отсутствии совпадении, качественно лучше соответствие наблюдается с распределением плотности. Для всех исследуемых почв можно отметить следующие особенности:

- графики распределения коэффициента ветвления более всего соответствуют графикам распределения проводящей категории плотности;

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

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

Следующей задачей была попытка нахождения соответствия путей фильтрации с категориями порозности - с распространением мезо- и микропор. Для этого в лабораторных условиях на насыпных образцах были получены ОГХ - основные гидрофизические характеристики (образцы отбирались по сетке 20x20 см). По ОГХ возможно определить объемы микро-, мезо- и макропор. Объем микропор вычислялся по объему влаги, содержащейся при давлении воды рр>4,18 (верхняя часть ОГХ). Объем ме-зопор вычислялся как разница между нижней границей макропор, которая рассчитывалась как влажность при рр 1,5, и объемом микропор, т.е. по влажности в диапазоне рр 1,5-4,18. Объем макропор не вычислялся, так как в насыпных образцах не приходится говорить о реальном их количестве.

Далее полученные значения были разбиты на три класса. Поскольку не было возможности создать единую шкалу, для каждого слоя были свои границы классов. Эти градации были условно названы как "мало", "средне" и "много" мезо- и микро-пор.

Аналогично показателю переноса влаги по соответствующей категории плотности (53) был рассчитан показатель переноса влаги по соответствующей категории мезо- и микропор. Рассчитывается как процент площади соответствующей категории пор, приходящейся на зону фильтрации влаги.

Указывает на долю соответствующей категории пор в переносе влаги на данной глубине. Результаты расчетов представлены на рис. 8.

Из приведенных рисунков видно, что ни для распределения микропор, ни для распределения мезопор нет определенной закономерности в проведении воды. Во

всех слоях представлены все три категории, но в разных соотношениях. а б

100

100

ОТНОСИТЕЛЬНЫЙ ВКЛАД КАТЕГОРИИ ПОР ■♦-"мало" Я- "средне" А- "много"

Рис 8. Изменения показателя переноса влаги БЗ по различным категориям пор: а) мезопор; б) микропор с глубиной ■

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

Суммируя все полученные данные, можно предложить следующую схему переноса влаги при низконапорной фильтрации (рис. 9). Как видно из приведенной схемы, полагаем, что основные зоны фильтрации и, соответственно, наиболее глубоко проникающие водные тяжи, приурочены к зонам почвы с постепенным переходом плотности. А отдельные водные внутрипочвенные потоки обусловлены еще и формированием их в зоне заметного перепада плотности - в данном случае, в зоне "плужной подошвы". Таким образом, в исследованных почвах в условиях стационарной низконапорной фильтрации горизонтальная и вертикальная изменчивость плотности будет являться доминантным фактором формирования преимущественных путей миграции влаги.

0 ^ Ш ■ ш

рыхлая слаба уплотненная плотная

уплотненная

1 поток влаги (величина стрелки соответствует « интенсивности потока)

Рис. 9. Схема формирования преимущественных путей переноса влаги при фильтрации

ВЫВОДЫ

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

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

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

3. Коэффициент ветвления (5ветв), показывающий "рассечение" общего потока с поверхности на отдельные водные тяжи в процессе фильтрации, увеличивается на глубине 25-45 см, что связано с образованием плужной подошвы. Ниже ее Sветв резко падает, что обусловлено формированием и стабилизацией постоянных потоков.

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

СПИСОК ОПУБЛИКОВАННЫХ РАБОТ:

1. В.М. Гончаров, МА. Мазиров, М.А. Бутылкина, Е.В. Фаустова, А.В. Кириченко, КА. Марченко. Основные характеристики уплотнения почв/Яезисы докладов Международной научно-практической конференции "Проблемы воздействия движителей на почву и эффективные направления ее решения", Москва, 1998, с. 66-67

2. Марченко КА. Пространственное распределение плотности почвы и путей движения влаги по профилю серых лесных почв/ЛГезисы докладов VII международной конференции студентов и аспирантов по фундаментальным наукам "Ломоносов-2000", 2000, секция "Почвоведение", с. 49.

3. Shein E.V., AB.Umarova, A.V.Dembovetsky, X.A. Marchenko. Subsoil density distribution and matter movement in soils//Proceedings of 2nd Workshop and International Conference on SUBSOIL COMPACTION, Godollo, Hungary, 2000, p. 105-111.

4. Умарова А.Б., Милановский Е.Ю., Марченко K.A., Яковлева Т.Ю., Калмыкова О.В. Исследование движения влаги и веществ в серых лесных почвах методами меток и встроенных лизиметров//Материалы III съезда почвоведов, Суздаль, 2000.

5. Шеин Е.В., Марченко К.А. Взаимосвязь путей движения влаги и пространственного распределения плотности почвы Владимирского ополья//Почвоведение, №7, 2001, с. 823-833

6. Шеин Е.В., Марченко К.А. Преимущественные пути миграции влаги//Вест-ник Московского университета, сер. 17, Почвоведение, 2002, № 1, с. 45-49.

Типография ордена «Знак почета» издательства МГУ 117234, Москва, Ленинские горы Заказ № 1558 Тираж 100 экз.

,»--50*

РНБ Русский фонд

2004-6 1893

Содержание диссертации, кандидата биологических наук, Ананьева, Ксения Аркадьевна

ВВЕДЕНИЕ.

ГЛАВА 1. СОВРЕМЕННЫЕ ПРЕДСТАВЛЕНИЯ О ДВИЖЕНИИ ВОДЫ В ПОЧВЕ И УСЛОВИЯХ, ЕГО ОПРЕДЕЛЯЮЩИХ.

1.1. Движение воды в условиях стационарной фильтрации.

1.2. Преимущественные пути фильтрации влаги.

1.2.1. Почвенная структура.

1.2.2. Преимущественный поток растворенных веществ.

1.2.3. Изменение фильтрации в связи с пространственной вариабельностью физических свойств почв.

• 1.3. О Владимирском ополье.

ГЛАВА 2. ОБЪЕКТЫ И МЕТОДЫ ИССЛЕДОВАНИЯ.

2.1. Характеристика объекта исследования.

2.2. Основные методы исследования.

2.2.1. Полевые методы.

2.2.2. Лабораторные методы.

ГЛАВА 3. КОЛИЧЕСТВЕННАЯ ОЦЕНКА КАЧЕСТВЕННЫХ ПРЕДСТАВЛЕНИЙ О ПРЕИМУЩЕСТВЕННЫХ ПУТЯХ ФИЛЬТРАЦИИ ВОДЫ

ВЫВОДЫ.

СПИСОК ОПУБЛИКОВАННЫХ РАБОТ.

Введение Диссертация по сельскому хозяйству, на тему "Преимущественные потоки почвенной влаги в условиях стационарной фильтрации"

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

Опыт показывает, что адекватные модели в почвоведении появляются, когда исследуемые явления, процессы или свойства глубоко поняты и изучены экспериментально [Дмитриев и др., 1986].

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

В последнее время уделяется большое внимание формированию так называемых преимущественных путей движения влаги. Образование их связывают с неравномерностью движения (в песчаных почвах), со слоистостью. Значительно меньше изучены вопросы образования преимущественных путей фильтрации влаги, связанных с неоднородностью пространственного распределения физических свойств. Основоположником такого рода работ был Е.А.Дмитриев. Однако, и в этих работах приведены качественные результаты по преимущественным путям миграции, обусловленным в основном генетическими особенностями почв. Количественная взаимосвязь формирования преимущественных путей фильтрации и пространственного распределения физических свойств в масштабе почвенного индивидуума не была изучена.

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

Этому и посвящена данная работа, которая проводилась на кафедре физики и мелиорации почв МГУ под руководством профессора, доктора биологических наук Евгения Викторовича Шеина.

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

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

Полевые материалы были получены в результате экспедиций кафедры физики и мелиорации почв 1998-2000 годов на серых лесных почвах Владимирского ополья на территории Владимирского НИИСХ.

Лабораторные исследования проводились в лабораториях кафедры физики и мелиорации почв МГУ по окончании полевого периода.

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

Заключение Диссертация по теме "Агропочвоведение и агрофизика", Ананьева, Ксения Аркадьевна

выводы

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

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

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

3. Коэффициент ветвления (,9ветв), показывающий "рассечение" общего потока с поверхности на отдельные водные тяжи в процессе фильтрации, увеличивается на глубине 25-45 см, что связано с образованием плужной подошвы. Ниже ее £ветв резко падает, что обусловлено формированием и стабилизацией постоянных потоков.

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

СПИСОК ОПУБЛИКОВАННЫХ РАБОТ:

1. В.М. Гончаров, М.А. Мазиров, М.А. Бутылкина, Е.В. Фаустова, А.В. Кириченко, К.А. Марченко. Основные характеристики уплотнения почв//Тезисы докладов Международной научно-практической конференции "Проблемы воздействия движителей на почву и эффективные направления ее решения", Москва, 23-25 сентября 1998 г., с. 66-67

2. Марченко К.А. Пространственное распределение плотности почвы и путей движения влаги по профилю серых лесных почв//Тезисы докладов VII международной конференции студентов и аспирантов по фундаментальным наукам "Ломоносов-2000", 12-15 апреля 2000 г., секция "Почвоведение", с. 49.

3. Shein E.V., A.B.Umarova, A.V.Dembovetsky, Х.А. Marchenko. Subsoil density distribution and matter movement in soils//Proceedings of 2nd Workshop and International Conference on SUBSOIL COMPACTION, Godollo, Hungary, 29-31 May 2000, p. 105-111.

4. Умарова А.Б., Милановский Е.Ю., Марченко K.A., Яковлева Т.Ю., Калмыкова О.В. Исследование движения влаги и веществ в серых лесных почвах методами меток и встроенных лизиметров//Материалы III съезда почвоведов, Суздаль, 2000.

5. Шеин Е.В., Марченко К.А. Взаимосвязь путей движения влаги и пространственного распределения плотности почвы Владимирского опо-лья/Шочвоведение, №7, 2001, с.823-833

6. Шеин Е.В., Марченко К.А. Преимущественные пути миграции влаги// Вестник Московского университета, сер. 17, Почвоведение, 2002, № 1, с. 45-49.

Библиография Диссертация по сельскому хозяйству, кандидата биологических наук, Ананьева, Ксения Аркадьевна, Москва

1. Bachmat Y., Elrick D.E. Hydrodynamic instability of miscible fluid in a vertical porous column. Water Resour., 1970, №6/

2. Corey J.C. Evaluation of dyes for tracing water movement in acid soils. Soil.Sci., 1968, v. 106, №5

3. Flury, M., H. Fluhler, W.A.Juiy, and J.Leuenberger. Susceptibility of soils to preferential flow of water: a field study. Water Resour. Res., 1994, 30:1945-54

4. Glass R.J., Steenhuis T.S., Parlange J.Y. Wetting front instability as a rapid and farreaching hydrological process in the vadose zone. J. Contam. Hydrol. №3. 1988.

5. Glass R.J., Steenhuis T.S., Parlange J.Y. Wetting front instability. Water Resours. Res. 25,1989.

6. Hill E.D., Parlange J.Y. Wetting front instability in layered soils. Soil Sci. Am. Proc., 1972, №36.

7. Hillel, D., R.S.Baker. A descriptive theory of fingering during infiltration into layered soils. Soil.Sci., 1988, 146:51-56

8. Kueper B.H., Flind E.O. An overview of immiscible fingering in porous medium. J. Contam. Hydrol, 1988, №2

9. Lacey W.J., Laguna W. De. Method of preparing radioactive cations for tracing-ground water. Science, 1956, v. 124, №3218

10. O.Lin C.C. "The Theory of Hydrodynamic stability". Cambridge Univ. Press, London and New York, 1955.1 l.Peck A.J. Moisture profile development and air compression during water uptake by bounded porous bodies: 3 vertical columns. Soil Sci., 1965, 100.

11. Philip S.R. Stability analysis of infiltration. Soil Sci. Am. Proc., 1975, №39.

12. Raats P.A.C. Unstable wetting fronts in uniform and nonuniform soils. Soil Sci. Am. Proc., 1973, №37.

13. H.Ritsema, C.J., L.W.Dekker, J.M.H.Hendrickx, and W.Hamminga. Preferential flow mechanism in water repellent sandy soil. Water Resour. Res., 1993, 29:21832193

14. Saffman P.G., Taylor G.I. The penetration of a fluid into a porous medium or Hele-Show cell containing a more viscous liquid. Proc. R. Soc., 1958, London, ser. A 245.

15. Smith R.E., Diekkruger. Field-Scale Soil Water Flow in Heterogeneous Soil I. Modeling Statistical Soil Variation and Large-Scale Constituent Reletions

16. Steenhuis, T.S., C.J.Ritsema, L.W.Dekker, and J.Y.Parlange. Fast and early appearance of solutes in groundwater by rapid and far-reaching flows. In: Transactions of the 15th World Congress of Soil Science, Acapulco, Mexico, July 1994, 184-203

17. Wooding R.A. Growth of fingers at an unstable diffusing interface in a porous medium or Hele-Show cell. J. Fluid Mech., 1969, №39.

18. Александровский A.JI., Гольева А.А. История развития почв со вторым гумусовым горизонтом. Тезисы докладов 2-го съезда общества почвоведов, книга 2-я, С-Пб, 1996, с. 3.

19. Алифанов В.М. Изменение серых лесных почв при сельскохозяйственном использовании. Почвоведение, 1979, №1

20. Алифанов В.М. Палеокриогенез и современное почвообразование. Пущино, 1995, с. 159-163.

21. Алифанов В.М., Гугалинская Л.А. Основные палеокриогенные свойства и их роль в почвообразовании. Тезисы докладов 2-го съезда общества почвоведов, книга 2-я, С-Пб, 1996, с.4.

22. Андроников В.Л. Об изменении свойств серых лесных почв при окультуривании. В сб.: Изменение почв при окультуривании, их классификация и диагностика. Колос, 1965.

23. Березин П.Н., Воронин А.Д., Шеин Е.В. Физические основы слитогенеза. Почвоведение, 1989, №1, с.31-38.

24. Бондарев А.Г. Теоретические основы и практика оптимизации физических условий плодородия почв. Почвоведение, 1994, с. 10-15.

25. Бондарев А.Г., Медведев В.В. Изучение физических свойств почв сельскохозяйственного использования Франции. Почвоведение, 1988, №8, с.36-43.

26. Бондарев А.Г., Силаков С.Н. Изменение физических свойств серых лесных почв при окультуривании. Почвоведение, 1993, № 7, с. 107-112.

27. Вадюнина А.Ф., Корчагина З.А. Методы исследования физических свойств почв. М.: Агропромиздат, 1986.

28. Васильев И.С., Роде А.А. Опыт «мечения» почвенной влаги ионом С1 с целью изучения закономерностей ее передвижения в полевых условиях. Почвоведение, 1960, №4

29. Величко А.А. Вопросы геохронологии лессов Европы. Известия АН СССР, 1965, география, №4

30. Величко А.А., Морозова Т.Д., Нечаев В.П., Порожнякова О.М. Позднеплейстоценовый криогенез и современное почвообразование в зоне южной тайги (на примере Владимирского ополья). Почвоведение, 1996, №9, с.1056-1063.

31. Величко А.А., Морозова Т.Д., Нечаев В.П., Порожнякова О.М. Палеокрио-генез, почвенный покров и земледелие. М. «Наука», 1996, с. 92-121. Всего с. 152

32. Воронин А.Д. Основы физики почв. Изд. МГУ, 1986

33. Воронин А.Д. Структурно-функциональная гидрофизика почв. Изд. МГУ, 1984.

34. Воронин А. Д., Губер А.К., Шеин Е.В. Использование почвенно-гидрологических констант для расчета параметров гидрофизических характеристик. Почвоведение, 1996, №5, с. 630-634.

35. Глобус A.M. О применении радиоактивных комплексонатов в качестве метки воды для почвенно-гидрологических исследований. Почвоведение, 1961, №9

36. Грин A.M. Зависимость инфильтрационной способности от физических и химических свойств почв. Почвоведение, 1972, №895

37. Грин A.M., Назаров Г.В. Сравнительная характеристика впитывающей способности почв лесостепной зоны европейской части СССР. Почвоведение, 1965, №3

38. Губер А.К., Шеин Е.В. Адаптация и идентификация математических моделей переноса влаги в почвах. Почвоведение, 1997, № 9, с. 1107-1119.

39. Губер А.К., Шеин Е.В., Ван Ицюань, Умарова А.Б. Экспериментальное обеспечение математических моделей переноса воды в почвах, оценка адекватности и надежности прогноза. Почвоведение, 1998, №9, с. 1127-1138.

40. Гумматов Н.Г., Пачепский Я.А., Щербаков Р.А. Пространственно-временная изменчивость водоудерживания серой лесной почвы. Почвоведение, 1991, № 9, с. 169-175.

41. Гуссак В.Б., Зухуров Н., Ахметжанов И., Раджапов А. Опыт исследования варьирования влажности сероземов. Почвоведение, 1967, №4.

42. Димитров И., Борисова М. Влияние уплотнения сельскохозяйственными машинами на структуру пахотного горизонта почвы и урожайность культур. Почво-знание, агрохимия и экология, 1997, 32, № 3, с. 36-42.

43. Дмитриев Е.А. Закономерности пространственной неоднородности состава и свойств почв. Диссертация на соискание уч. ст. д. б. н. в форме научного доклада. Издательство Московского университета.

44. Дмитриев Е.А. К методике полевого изучения путей передвижения в почве жидкой влаги. Научные доклады высшей школы, 1971, №5, биол. науки

45. Дмитриев Е.А. Некоторые особенности формирования поверхностей смачивания при дождевании (модельные опыты). Почвоведение, 1976, № 10

46. Дмитриев Е.А. Полевой почвенный пантограф. Почвоведение, 1977, №9

47. Дмитриев Е.А., Карпачевский JI.O. и др. О происхождении неоднородности почвенного покрова в лесных биогеоценозах. В кн.: Проблемы почвоведения, 1978

48. Дмитриев Е.А., Макаров И.Б. О понятии «равновесная плотность почв». Почвоведение, 1993, № 8, с. 94-98.

49. Дмитриев Е.А., Сибуль Р.А. Объемная плотность верхнего слоя дерново-подзолистой почвы и ее связь со структурой. Почвоведение, 1980, №4, с. 8391.

50. Дмитриев Е.А., Хохрина Т.К. О путях передвижения впитывающейся в почву влаги. В сб.: Проблемы сельскохозяйственной науки в МГУ, Изд-во МГУ, 1975

51. Дмитриев Е.А., Щеглов В.Н., Басевич В.Ф. Морфология движения впитывающейся во влажную почву влаги и определяющие ее факторы. Почвоведение, 1985, №1

52. Дмитриев Е.А., Щеглов В.Н., Басевич В.Ф. Характер миграции воды во влажных почвах. Почвоведение, 1985, №8

53. Дубенский Н.Я. О почвах Владимирской губернии. Ж. Моск. Общества с.х., 1855, №4-5

54. Емельянов В.А., Маслов В.П. О варьировании плотности почв на орошаемом поле. 1985

55. Иванов JI.A. Ботанические и почвенные исследования в Юрьевском и Суздальском уездах Владимирской губернии. М., 1897

56. Иванов JI.A. Еще о Владимирском черноземе. Почвоведение, 1899, №2

57. Корсунская Л.П., Гумматов Н.Г., Пачепский Я.А. Сезонная динамика содержания корневой биомассы, углеводов и показателей структуры серой лесной почвы. Почвоведение, 1994, № 10, с. 75-78.

58. Костяков А.И. О динамике коэффициента просачивания воды в почвогрунты и необходимости динамического подхода к его изучению в мелиоративных целях. Почвоведение, 1932, №3

59. Кошелева Н.Е. Моделирование почвенных и ландшафтно-геохимических процессов. Изд-во Московского ун-та, 1997, 110 с.

60. Красюк А.А. Почвенные районы Иваново-Вознесенской, Костромской и Владимирской губерний. 1925

61. Кузнецова И.В. Изменение дифференциальной пористости и водно-физических свойств почв при уплотнении. Плодородие почв и его изменение при уплотнении и разуплотнении. М., 1984, с. 18-24.97

62. Кузнецова И.В., Данилова В.И. О разуплотнении почв под влиянием процессов набухания-усадки. Почвоведение, 1988, №6, с. 59-70.

63. Макеев А.О., Дубровина И.В. География, генезис и эволюция почв Владимирского ополья. Почвоведение, 1990, №7, с. 5-25.

64. Марков К.К., Величко А.А. Четвертичный период. Недра, 1967, т. 3

65. Мешалкина Ю.Л., Самсонова В.П. Использование методов геостатистики для выявления пространственных неоднородностей сельскохозяйственных угодий. Тезисы докладов 2-го съезда общества почвоведов, книга 2-я, С-Пб, 1996, с. 382.

66. Мильков Ф.Н. О природе ополий на Русской равнине. В сб.: Вопросы регионального ландшафтоведения и геоморфологии. Изд-во Львовск. Ун-та, вып.8, 1964.

67. Михович А.И. Метод контроля инфлюкции. Почвоведение, 1961, №8

68. Полевые и лабораторные методы исследования физических свойств и режимов почв (под ред. Шеина Е.В.). Изд-во Московского университета, 2001.

69. Рубилин Е.В., Долотов В.А. Генезис и география лесных почв Европейской лесостепи СССР. Изд. ЛГУ, 1970

70. Рубцова Л.П. О генезисе почв Владимирского ополья. Почвоведение, 1974, №6, с. 17-27.

71. Рубцова Л.П. О неоднородности почвенного покрова Владимирского Ополья. Природа и сельское хозяйство Калужской области, 1970

72. Рубцова Л.П. О неоднородности почвенного покрова Владимирской опытной станции. Труды Калужской государственной областной сельскохозяйственной опытной станции, 1970, т. 7

73. Самсонова В.П., Дмитриев Е.А., Витязев В.Г. Пространственная организация морфологических профилей и удельной поверхности в пахотной дерново-подзолистой почве. Почвоведение, 1997, № 10, с. 1200-1206.98

74. Танфильев Г.И. О Владимирском черноземе. Почвоведение, 1899, №1

75. Троицкий А.И. Серые лесные почвы северной части Среднерусской возвышенности (путеводитель 1 Делегатского съезда почвоведов). Изд. АН СССР, 1958

76. Троицкий А.И., Воробьева Э.С., Калнина В.А. Серые лесные почвы северной лесостепи и их агропроизводственная характеристика. В сб.: Опыт крупномасштабной картографии и агропроизводственной характеристики почв Нечерноземной полосы, Колос, 1970

77. Тюрюканов А.Н., Быстрицкая T.JI. Ополья центральной России и их почвы. М., Наука, 1971

78. Хантулев А.А., Гагарина Э.И. и др. Особенности строения и причины формирования почв со сложным гумусовым профилем в зоне бореальных лесов Евразии. В сб.: География и генезис почв, Изд. ЛГУ, 1969

79. Хантулев А.А., Гагарина Э.И. Материалы о генезисе и географии почв со сложным гумусовым горизонтом. Почвоведение, 1972, №2

80. Шеин Е.В. и др. Лабораторные методы исследования физических свойств почв. ГЕОС, 2000

81. Шеин Е.В., Березин П.Н., Капинос В.А. Задачник по физике почв. Изд. МГУ, 1988

82. Шеин Е.В., Гудима И.И., Бандина Л.Ю. Исследование воздухоносной пористости в почвах при орошении. Почвоведение, 1989, №12, с. 60-68.

83. Шеин Е.В., Капинос В.А. Сборник задач по физике почв. Изд. МГУ, 1994

84. Шеин Е.В., Салимгареева О.А. Пространственная вариабельность водного режима чернозема типичного. Тезисы докладов 2-го съезда общества почвоведов, книга 2-я, С-Пб, 1996, с. 388.

85. Шеин Е.В., Салимгареева О.А. Пространственная вариабельность физических свойств и водного режима чернозема типичного. Почвоведение, 1997, № 4, с. 484-492.

86. Якушевская И.В. Почвы Владимирского Ополья. М., 1956, автореф. дис.

87. Ярошенко П.Д., Юрова Э.Л. Об остепненных лугах Владимирского ополья. Бюлл. МОИП, отд. биол., 1970, т.57, вып.499