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

ФОРМИРОВАНИЕ РЕЖИМА И РЕСУРСОВ ПОЧВЕННЫХ ВОД В ЗИМНЕ-ВЕСЕННИЙ ПЕРИОД

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

Автореферат диссертации по теме "ФОРМИРОВАНИЕ РЕЖИМА И РЕСУРСОВ ПОЧВЕННЫХ ВОД В ЗИМНЕ-ВЕСЕННИЙ ПЕРИОД"

А-зобэ?

МОСКОВСКИЙ ОРДШ ЛЕНИНА, ОРДЕНА ОКТЯБРЬСКОЙ РЕВОЛЩШ И ОРДЕНА ТРУДОВОГО КРАСНОГО ¡ЗНАМЕНИ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

ИН.Ы.В.ЛОМОНОСОВА ФАКУЛЬТЕТ ПОЧВОВЕДЕНИЯ

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

ГУСЕВ ЕВГЕНИЙ МИХАЙЛОВИЧ

ФОРМИРОВАНИЕ РЕЖИМА И РЕСУРСОВ ПОЧВЕННЫХ ВОД В ЗИМНЕ-ВЕСЕННИЙ

ПЕРИОД

06.01.03 - агрошчвоведекие и агрофизика

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

Москва - 1992

Работа выполнена в Институте водных проблем Российской академии наук.

доктор биологических наук, профессор Й.И.СТДНИЦУН доктор сельскохозяйственных наук, профессор А.м.ГЛОБУС доктор сельскохозяйственных наук Н.А.МГРОЩВВ

Ведущее учреждение - ИНСТИТУТ ГЕОГРАИЮ РОССИЙСКОЙ АКАДВ-

т ник.

Защита состоится фвб^ал^ 1993 г. в 15-2^- в аул. Ы-2 на заседании Специализированного Совета Д.053.05.31 на Факультете почвоведения МГУ им. Ы.В.Ломоносова.

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

Приглашаем Вас принять участие в обсувдении диссертанте» на заседании Специализированного совета Д.053.05.31 в Мсс*.: . ском университете, а отзывы на автореферат в двух экземпляр заверенные печатью, просим направлять по адресу:

II989&, Москва, Ленинские горы, ИГУ, факультет почвоведе яия, Ученый совет.

Официальные оппонент:

Ученый секретарь Специализированного Совета

Л.А.ЛВЕВДВВ1

ВВЕДЕНИЕ 1 О0ДАЯ ХАРАКТЕРИСТИКА РАБОТЫ

. Актуальность проблемы. Почвенные вот являются одним из важнейших элементов биосферы., обеспечивая существование жизни ва Земле..В последние года теш эволюции биосферы в связи с активизацией человеческой популяции на планете резко усилился, проводя к экологическим-катастрофам в ряде регионов. .Ускорение изменений биосферы требует постоянного совершенствования научно-инженерного инструментария во взаимодействии с природой. Это совершенствование не всегда успевает за биосферной ¿камикой. В частности, негативные последствия многих крупных водохозяйственных мероприятий в нашей стране объясняются не только издержками командно-административной системы, но и слабим развитием фундаментальных исследований в учении о водах суш, в частности о почвенных водах. При этом следует отметить тот факт, что само направление учения о почвенных водах, как части вод суши, слс-жаирсь в последние 10-20 лет. Создание этого учения во многом связано с работами д.г.н. А.И.Будаговского.

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

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

Под почвенной водой в настоящей работе понимается вода, содержащаяся в корнеобитаемом слое почвы. Вода, пребьгващэя в ипостаси почвенной воды,' играет в биосфере исключительно важную роль. [Будаговский, Гусев, 19691. Под режимом же почвенных вод понимается временная динамика характеристик, в той или иной мере отражающих их сое ш инда ^См^ММё. ЛИЦ

ра, величины и напр;

стость, температу-ла). Понятие ре-

М»:ч.се;, . зададим :

им, К. , »„.^.^в и.,г, Мо А-"'

сурсов почвенных вод, принятое в настоящей работе» обсуждается в ее первой главе.

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

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

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

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

Задачи исследований:

- экологическое обоснование концепции ресурсов почвенных вод,

- разработка теоретической основы описания переноса воды. в ненасыщенных средах с позиций статистической неравновесной термодинамики,

- экспериментальное исследование процессов промерзания почвы, миграции воды в промерзшую почву со стороны талой зоны, впиты-^ вания воды в почву в период весеннего снеготаяния, оттаивания и

прогревания почвы после схода снежного покрова,

- разработка комплекса -динамических физико-математических моде-• лей процессов формирования режима почвенных; вод в зимне-

весенний период»

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

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

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

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

Подмосковная воднобалансовая станция (лесная зона). Почва

- дерново-среднеподзолистая, средне суглшщстая.

Волховская воднобалансовая станция (лесостепная зона). Почва - серая лесная, среднесуглинистая.

Нижнеденицкая воднобалансовая станция (лесостепная зона). Почва - типичный чернозем, средне- и тяжелосуглинистый.

Придеснянская воднобалансовая станция (лесостепная зоне). Почвы - крупнопылеватые серые и светлосерые лесные, средний, суглинок и супесь.

Курская биосферная станция (лесостепная зона). Почва типичный чернозем,, средне- и тяжелосуглинистый. Научная новизна работы.

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

На основе статистической теории неравновесных процессов получено общее уравнение переноса »ода в ненасыщенных средах, являющееся теоретическим обоснованием используемы! в настоящее время методов описания движения води в ненасыщенной почве. Теоретически обосновано использование в случае слаОонестационарных процессов переноса вод» в почве уравнения Дарси- Букин-гема.

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

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

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

Осдовнре яяпррцявмме т^юуддоя:

- уравнение нестационарного переноса,вода в ненасыщенных средах;

- комплекс <&гаика-матечатиче ских моделей формирования гидротермического режима почвы в зимне-весенний период (модели промерзания почвы в холодный период, формирования температурного црофиля почвы к началу снеготаяния, впитывания талых вод в почву в период весеннего снеготаяния и зимних оттепелей, оттаивания и весеннего прогревания почвы);

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

- концепция экологической роли почвенных вод, заключающаяся в том, что в процессе своего движения в системе почве-

- s -

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

- методы оценки блеяния растительной мульчи sa гидротермический режим почвы в зимне-весенний период.

- методы определения коэффициентов эффективной теплопроводности г эффективное диффузии пара в мульче.

Практическая значимость работы.

Полученные в работе результаты были использованы:

- в плановых работах Института водных проблем АН СССР и представлены в отчетах по плановой тематике Академии наук СССР (тема "Водные ресурсы и водный режим страны в отдаленной перспективе", 1985 г., гос. регистрац. ЖЯ827018477; тема "Изменение ресурсов и режима вод суши под влиянием природных колебаний климата и антропогенной деятельности", 1990 г., гос. регистрац. ÜOI.9.00018345);

- при реализации Программы биосферных и экологических исследований (Блок "Экология сельскохозяйственного производства"), 1989г.;

- при выполнении теми "Прогноз природопользования и экологии до 2015 г.*, разрабатывавшейся по заданию ГКНТ СССР и Госкомприрода СССР, 1991 г.;

- при разработке международных проектов, проводившихся по задания Международного комитета по планетарным геофизическим исследованиям (КАЛГ) (субпроект 17.I "Физика почвенных вод", 1980-1985 г.г.; субпроект Ш-3.1 "Изучение переноса воды в системе атмосфера - растительный покров - ненасыщенная зона -грунтовые воды в целях разработки научных основ управления почвенными водами", 1986-1990 г.г.).

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

Результаты работа также могут Сыть использованы при оценке эффективности различных мелиоративных приемов, а также в практике гидрологических расчетов и прогнозов.

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

- б -

Московский городской гидрологический семинар, (1977 г.); рабочие■ встречи участников субпроекта 17.I Международного комитета по планетарным геофизическим исследованиям (КАПП "Физика почвенных вод" (1978 г., СССР, Курск; 1982 г., ЧССР, Требишов; 1983 г., НРБ, София); Всесоюзная школа-семинар "Моделирование гидрологически! процессов", (Звенигород, 1982 г.î; Чехословацко-польский симпозиум "Физика почвенных вод* (1380 г., ПНР. Закопане); 11-ая Всесоюзная конференция по применению математичесюп: методов и ЭВМ в почвоведении (1963 г., Пуйино); Доклада в Отделе гидродинамики пористых сред Института гидрологии и гидравлики САН (1985, 1X6 г.г., ЧССР, Братислава); рабочая встреча участников проекта III-3 КАПГ "Геофизические основы управления водами суши" (1986 г., ЧССР, Стара Лесна); XX Совещание КАПГ (1989 г, ВНР, Шопрон); Международный семинар "Движение воды в системе атмосфера-ненасыщенная зона-грунтовые воды и его моделирование" (1990, ЧОФР, Стара Лесна); Всесоюзное совещание "Определение и использование параметров влаго- и солепереноса в почвах" (1991 г., Пущино), Заседание кафедры физики и мелиорации почв МГУ (1993 г., Москва); Секции Ученого совета Ш1 АН СССР (1985 - 1992 г.г.).

Структура и объем работы.

Диссертация состоит из введения, пяти глав, основных выводов и приложения. Она включает 250 стр. машинописного текста, 32 рисунка, 7 таблиц и список литературы из 175 наименования (4.9- на иностранных языках).

Ще раз подчеркнем, что создание общего учения о почвенных водах в основном связано с именем А.И.Вудаговского. В то -же время большое значение для развития наших работ имели предшествующие исследования различных частных аспектов этого учения, которые получили свое отражение в трудах очень многих исследователей, таких как Нерпин C.B., Глобус A.M.. Пачепский Я.А.. Калюжный И.Л., Кучмент Л.С., Мотовилов Ю.Г., Менжулин Г.В., Еудыко М.К., Росс Ю.К., Тооминг Х.Г., Еихеле З.Н., Мол-дау Х.А., Сиротенко Ю.Д., Бойко А.П., Чураев Н.В.. Воронин А.Д., Судницын И.lt., Лыков A.B., Мичурин Б.Н., Ничипорович A.A., Галямин Е.п., Вазилевич Н.И., Родин Л.Е., Коронкевач H.H., Кренке А.Н., Кузьмин п.П., Монин A.C., попов Е.Г., Ершов

Э.Д., Кулак В.Я., Павлов ¿.В., Роде I.A., Свирепев D.U., Виноградов Ю.Б., Денисов D.M., Чудаовский А.Ф., Полуэктов P.A., Слейчер Р., Кари Дж.в., Лархер. В., Адлер ы., Христов И.Д.. Phylip J.R., Witt С.!., Taylor S.A., Childs B.C., Youngs E.G., Rose D.A.» Anderson D.H., Rubin J., Gardner W.R., Nualem Y., Da^n G.t Montelgji J.L., Freeze R.A., Sutor J., Noväk V., Spoalto G., Eagleaon P.S., Yosida Z, и многих других,

СОДЕРЖАНИЕ РАБОТЫ

ГЛАВА I. ПОЧВЕННЫЕ ВОДЫ И ЭКОЛОГИЯ НАЗМШ РАСТИТЕЛЬНЫХ СООБЩЕСТВ.

Воды суш принято подразделять на поверхностные, подземные, почвенные воды, ледники и снежники (Великанов, 1943). Понятия о ресурсах поверхностных и подземных вод употребляются в гидрологии, гидрогеологии и водном хозяйстве уже давно и общеприняты. Что же касается понятия ресурсов почвенных вод, то оно было введено сравнительно недавно (Будаговский, 1985} и связано с тем, что почвенные воды - важнейший фактор существования и развития растительного покрова суш. в качестве мери ресурсов почвенных вод А.И.Будаговским на основе анализа уравнения водного баланса водосбора предложено использовать величину суммарного испарения с поверхности суши за безморозный период Ес, представляющего собой сушу испарения воды почвой Ед и поверхностью растительного покрова - транспкрации Ет.

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

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

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

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

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

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

Циркуляция вода в структурах гидрологического цикла в конечном итоге поддерживает непрерывную циркуляцию биоэлементов суши. Швтому собственно их и можно рассматривать в качестве биосферных ресурсов. В большей степени это относится к почвенным водам, поскольку они приводят в движение самое мощное звено в циркуляции биоэлементов наземных экосистем - звено зеленых растений суши. Причем, в настоящее время 80-90» так называешь "безвозвратных потерь" стока, т.е. изымаемых человеком ресурсов поверхностных а подземных вод, используются в

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

Жизненные процессы на Земле имеют четко выраженный сезонный ход. При этой величина суммарного испарения за период вегетации определяется не только условиями теплого периода. В таких странах, как СНГ, СОЦ., Канада, в странах Северной Европы и др. годовая величина суммарного испарения в значительной мере обусловлена величиной снегонакопления на поверхности почвы и последуодим разделением выпавших зимних осадков между впитыванием их в почву и поверхностным стоком в' период весеннего снеготаяния.

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

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

зование аналитических методов неизбежно связано с выделением основных черт явления и пренебрежением второстепенными. Одна-

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

Удобным при исследовании оказалось * выделение следувдих характерных периодов.

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

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

3. Весеннее снеготаяние и впитывание талых вод в почву.

4. Оттаивание почвы после схода снежного покрова.

5. Зимние оттепели.

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

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

ГЛАВА 2. ПЕРЕНОС ВОДЫ В НШАСЩПИП ПОРИСТЫХ СРЕДАХ.

В совокупности исследований, посвященных построению ОГО1-

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

»- я — , (I)

вг

которое связывает поток воды д в конкретной точке системы г и в данный момент времени t с градиентом потенциала почвенной воды ф в этой же точке и в тот же момент времени. Параметр К в (I) называется коэффициентом влагопроводности. Эмпирический характер приведенного уравнения и появляющиеся время от времени сведения о его несоответствии экспериментальным данным заставили обратиться как к поиску теоретического обоснования этого уравнения, так и к.определению области его применимости.

В качестве пути решения задачи был выбран третий из указанных подходов, основанный на применении статистической тории неравновесных процессов. В результате было получено уравнение нестационарного переноса воды в ненасыщенной среде [Гусев, 1979а, 1979в, 1981а): г «

где Щ - влажность среды, г1 - переменные интегрирования.

Полученное уравнение переноса в отличии от уравнения Дарси-Букингема показывает, что поток воды д в данный момент времени г и в конкретной точке среды г определяется градиентом водного потенциала во все предыдущие моменты времени tA (эффект запаздывания) и-в некоторой окрестности рассматриваемой точки ¡2-г,I (эффект нелокальности}.

Отмеченные эффекты описываются'на основе релаксационной функции £, являющейся в данном уравнении основной гидрофизической характеристикой среды.

В случае пренебрежения эффектами нелокальности (они имеют

- а -

место лишь при наличии очень резко меняющихся в пространстве градиентов потенциала)уравнение (2) упрощается да

(3)

*1

Поскольку получение теоретического описания структуры функции £ в настоящее время вряд ли возможно, был разработан метод ее экспериментального определения (Гусев 19636}. Метод основан на получении опытных данных о динамике влажности среды и потоков воды в ней, сведений о зависимости от влажности среды ее кашллярно-сорбшюнного потенциала, и последующей теоретической обработке этой информации. При этом задача сводится х решению интегрального уравнения, близкого по структуре к интегральным уравнениям Фредгольма 1-го рода. При атом в качестве решения берется функция впг^-г), которая доставляет " минимум следутему функционалу:

»гв; = - я)г<ЗР + а Ыи.) . (4)

?

Здесь Р - пространственно-временная область определения штоков, т.е. интервал времени, в течение которого измерялись потоки, а также пространственная область их фиксации, А -

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

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

Методика определения реелаксационной функции была реали-

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

На ркс.1 показан общий вид релаксационной функции &(1Г,Т'}> где Г^.-Г.

Eii« ¿„к1 (Hirt

Рис.1. Общий вид зависимости-релаксационной функции кварцевого песка от влажности V и времени запаздывания Т

b ъ

5

4

Рис.2. Зависимость коэффициента влагопроводности кварцевого песка от влажности, рассчитанная по данным о релаксационной функции (I) и определенная методом - "квазистационарной инфильтрации" (2)

в.Й ар. о,м>

цои1"^3

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

циональная зависимость %(Т) была бы близка к дельта-функции Ь(Т), в величина релаксационной функции значимо отличалась- бы от нуля только при Г « 0.

Представляет интерес проверить согласованность полученных в опыте с капиллярным опусканием данных по Л(И,Т) с какими-то другими экспериментальными материалами по влагопереносу в исследуемой системе, полученными при экспериментальном анализе других процессов в этой же системе. В частности, можно показать ГГусев, 196363, что коэффициент влагопрср^днооти, определенный, например, методом "квазистационарной инфильтрации", должен Сыть связан с релаксационной функцией следушим соотношением

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

Дальнейший анализ полученных уравнений переноса воды в ненасыщенных средах (2), (3) показал, что возможны различные способы их упрощений для разных частных случаев движения воды в ненасыщенных средах IГусев, 1979а, 197Эв1, вплоть до уравнения Дарси-Букингема. Таким образом, уравнения (2), (3) могут в некоторой мере служить обоснованием этого уравнения. Проведенный анализ также показал, что значительных, отклонений фактических данных от рассчитанных результатов, основанных на использовании уравнения Дарси-Букингема, можно ожидать лишь в случае сильно нестационарных процессов, которые отличаются резкими изменениями влажности во времени. При миграции воды в талой зоне к фронту промерзания и в смоченной части мерзлой зоны при впитывании воды с поверхности почвы быстрые изменения влажности практически не имени места в связи с относительно малыми величинами интенсивности влагопереноса. Поэтому основой для описания переноса вода в почве в этих случаях в настоящей работе послужило уравнение Дарси-Букингема.

о

(5)

- п* -

ГЛАВ! 3. ФОНДОВШИЕ РЕХША ПОЧВЕННЫХ ВОД В ЭНШИИ ПЕРИОД.

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

Вопросу промерзания почвы посвящено очень много работ Однако работ, основанных на применении аналитических методов при совместном решении задач промерзания и миграции воды со стороны талой зоны в мерзлую почву, нет. Поэтому в настоящей работе, были разработаны два метода расчета глубины промерзания почва [Гусев, 1965, 1988; Гусев, Ясинский. 1986], отличительной чертой которых является одновременный учет фактического хода метеоэлементов и фактической динамики снежного покрова в течение зимнего сезона, а также учет миграции влаги со стороны талой зоны. Первый способ основан на рекуррентном пошаговом во времени вычислении положения фронта промерзания, второй - на получении аналитического .выражения искомой характеристики на любой момент времени.

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

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

При решении задачи промерзания в мерзлой зоне почвы , а также в вышерэсположенном слое снега принималась линейная аппроксимация температуры,' в талой зоне - параболическая (рис.3) [Гусев, 1985, 1988]. Аппроксимировалось также отклоне- • ние б текущего значена влажности и от ее начального перед промерзанием значения В соответствии с законом дарск-Букингема полагалось, что при отклонении ЬЬ/Ь% от нуля поток

- 1С -

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

(6>

где паршетр Т>щ близок по смыслу ж значении коэффициенту диффузии почвенной воды при

а)

О V

__

< І® —г

г і

А

"7

5)

Рис.3. Схематизация распределения температури (а) и влажности (б) при промерзании почвы. I - снежный покров, 2 - мерзлая почва, 3 - Талая почва.

Выделение отклонения 0 удобно тем, что для него возможен простой вид аппроксимации в пределах талой зоны до некоторой глубины г=ц (рис.30):

« «Ш [ ££ ] . <7>

Ответим, что аппроксимация (7) справедлива до тех пор, пока граница депрессии влажности ц не достигла уровня грунтовых -вод. После этого необходим другой вид апроксимации (Гусев, 19861. Здесь же будем анализировать случай достаточно глубоко-

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

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

с, е т сг^

- £ "2 5

PJ.ni ~ и)--г -

" ■ 2 Е+^Л,

гв 2 Г*

<31

(8)

^ 1 / -Г ?+12аз(и*0> - М

где - плотность воды, Ъ - теплота фазового перехода вода-лед, Го - температура поверхности снега, , \г - теплопроводности соответственно снега и мерзлой почвы, о3 - температуропроводность талой почвы, сг- эффект Аная теплоемкость мерзлой зоны, Г* - температура почвы на глубине затухания сезонных колебаний, tQ - параметр, характеризуют«! влияние предзимья СГусев, 19051, Л - высота снежного покрова, количество неза-

-е. Л

мерзшей воды в почве, п -~ргг<3г - среднее по мерзлой зоне

^о

общее влагосодержание.

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

для [Гусев, 19881

-С ~ £ -5 —| *

о

1«

где й2 = - 0£(32- среднее приращение влажности в мерзлой ^ о

"б * £

зоне, а 0 (I) = - среднее за период промерзания

^ о

отклонение влажности у фронта промерзания от ее начального

эначения, 3(1) ~ поток води у фронта промерзания со стороны талой зоны.

Таким образом, чтобы определить приращения влагозапасов в промерзшей зоне, надо знать /Ш . а также изменение влажности у фронта промерзания. Согласно многочисленным аксперлменталь-ным данным (например, (Ершов, 1979; Глобус, 1963)) величину потока воды можно определить, исходя из следующих соотношений:

й|>

/(5) = - КСи)— дг

дф от

= -К(и)-

г=£ М О* ^

где Ыи) - коэффициент термовлагопроводности этого выражения и закона Дарси-Букингема для итоге дает (Гусев, 1988)

1

дт

-Ми)— дг

(Ю)

Использование талой почвы в

■

1 +

ЗЙ

<2£_ «а*

г . цря Я0.(£) +

=

при И1 Е

Ш) г/3 + ?

— Е

> и-..

<

(II)

9

Т

е2 + 12 V -г £

Здесь А0 -

/ /«ЛИ. о

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

Модель (7) - (II) была проверена на данных лабораторных опытов (Зигог, 1971), а также на материалах наблюдений на ряде воднобалансовах и Курской биосферной станций. Для проведения сопоставления с данными лабораторных оштов был использован аналитический вариант реализации модели, при сравнении с данными полевых наблюдений использовался численный вариант с шагом по времени, равным одной декада.

При моделировании тепловлагопереноса для конкретных полевых объектов возникает проблема нахождения таких па^метров как коэффициент термовлагопроводности й0 и коэффициент диффузии воды в почве С . Поскольку данных о зависимостях К(и), ф(и),

и(Т), позволяющих определить указанные параметрі для конкретных почв, еще недостаточно, возникает вопрос об их оценке путем решения обратных задач. Методика определения подобным образом параметров к0 и была разработана в настоящем исследовании [Гусев, 1968]. Необходимую информация іфи этом представляют данные о температуре поверхности снега (примерно равной температуре воздуха), динамике высот» снежного покрова, глубины промерзания и профиля влажности почвы, имеющиеся, например, в материалах наблюдений воднобалансовых станций. Вообще говоря, необходимы также данные о количестве воды, поступившей в почву в период значительных оттепелей, но можно выбрать зимние сезоны без оттепелей (что и было сделано в данной работе).

Определение 1>я и й0 по разработанной методике было осуществлено на материалах наблюдений, проведенных на полевых участках Подмосковной воднобалансовой станции (почва - дерно-во-подзо- листая, средний суглинок) и Волховской воднобалансовой станции (почва - серая, лесная, средний суглинок). Параметры и й0 определялись в зимние сезоны 1576 -1978 гг. и оказались равными соответственно для Подмосковной станция 2,5'ІСГ4 с^/с и 1,5'Ю"5 см2/с-град, для Волховской 2,0-10~* см2/с и 0,5'ІСҐ5 см2/с-град.

Сопоставление результатов расчета глубин' промерзания почвы и приращения влагозапасов в мерзлой зоне с опытными данными, полученными в результате проведения экспериментальных исследований водно-термического режима почвы на Курской биосферной станции ИГ АН СССР (1979-1980, 1983-1984 г.г.), данными наблюдений Подмосковной (I972-1981 гг), Волховской (19761984 г.г.) И Нижнедевицкой (1969-1972, 1976-1980 г.г.) водно-балансовых станций показало их удовлетворительное согласие (часть результатов приведена на рис.4).

ГЛАВА 4. ПОПОЛНЕНИЕ ЗАПАСОВ ПОЧВЕННЫХ ВОД В ПЕРИОД СНЕГОТАЯНИЯ

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

0,0!< 0,06

V

0,0$

<кач -

л- б], «Уем»

о о,оН цов Э'

Рис.4. Сопоставление рассчитанных и измеренных данных по динамике промерзания и изменении влагозапасов в промерзшей зоне почвы для Подмосковной (а) и Волховской (б) воднооа-лансовнх станций, {р и - соответственно рассчитанные и

измеренные глубины промерзания; Ср, рассчитанные и измеренные приращения влажности в мерзлой зоне почвы.

чета этого профиля [Гусев, Ясинский, 1986). В основу решения был также положен метод интегрального теплового баланса.

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

А, сг1г- суПг ог1гП * с:Пг£к}/\г

= + с2 л2 - ---г--; (12)

У У •

где

е - 2

У -т- 5 * Ь.

Ч

- 2 i -

r° азс Г p(t-t) 3(t -тЛ

- h\ ^r [VP ° j

arc f p(to-i) q(ta-i)]

о

Bf, - функции X,. Xg, h, ot, o,, p, <j 1Гусев,

Ясинский, 1966], а p и q - корни следущего квадраптого уравнения:

г2 + гЕ + 5 » 0 , (13)

В - 3

г 1 5 1 f h i

[нз--] ft + 3 —I.

-+ а1-:—

£г

36

a

Здесь a, ii oj - коэффициента темературопроводкости соответственно снега и мерзлой почвы, tQ - продолжительность предвесеннего периода (Гусев, Ясинский, 1986).

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

Звание температурного профиля почвы перед снеготаянием дает возможность подойти к расчету интенсивности впитывание талых вод. Особенностью данной работы в отличии от уже имеющихся по данной задаче (например Мотовилов, 1977; Зарецкий, Лавров, 1986)является то, что предлагаемая методика расчета впитывания талых вод основана на моделировании этого процесса с помощью аналитических методов [Гусев, 1989]. Численный характер решения задачи проявляется лишь в отношении изменения искомых величин во времени.

Используемое при атом дифференциальное уравнение переноса

воды в почве имеет известный вид в форме Дарси-Вукингема:

б ф„

q-* - к-8 + к , (14)

ъ ö г

где qB - шток воды в почве, капиллярно-сорбционная составляющая потенциала почвенной воды. Гидрофизические характе-

«О . $

Рис.5, распределение температури в почвеннш профиле. а)Нижнедевицкая водаобалвнсовая станция, б) Курская биосферная станция. I. Данные по термометрам сопротивления, 2.Данные по вытяжным термометрам, 3. Рассчитанные данные.

ристики X и в данном случае являются функциями содержания жидкой воды в почве и и ее льдистости Л, связанными с общей влажностью почвы ІГ следующим соотношением

Рд

1Г = —Л + и , (15)

Р*

где рд , рж - плотности соответственно льда и воды.

Интегрируя (14) от г*0 (поверхность почвы) до г=г0 (поло-юние фронта впитывания) и учитывая при этом, что в пределах колонны впитывания Я^ " Я0> где д0 - интенсивность впитывания, получим

?в - Ж О + в/г0), (16)

г0

К = | Кйг - средний коэффициент влагопроводности смо-о

Ф(0)

- эффективный капиллярный по-

ченной зоны, а В

Ф(г)

тенцаал почвы. В первом приближении можно считать, что интегральные параметры Ї и Й зависят от среднего содержания жидкой

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

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

*г Г®. гГн"о'1 Г 120,11

Ч [ 5- ^ ] ч 1?г <17>

а * -2г0Г1; + /9Ыоа)\г * 12 а^ , (18)

где ¿А - приращение льдястости в смоченной зоне, га - температура почвы перед снеготаянием.

Присоединяя к (17) уравнение движения воды в смоченной зове (16), уравнение динамики глубины проникания о (18), а также условие водного баланса смоченной зоны

| (я0- Кн)& - г0 [ ~ й + и - ия | , (19)

• ■ « о

где и^ - начальное содержание незамерзшей воды в мерзлой зоне, £д = Л0.) - коэффициент влагопроводности почвы перед фронтом впитывания, получаем систему четырех уравнений относительно пяти неизвестных ÍЛ(t)t U(t)l Яоа), 20(г), и оа). Для замыкания этой системы необходимо привлечь еще одно уравнение.

№< служит в случае выполнения условия - _ *

Л + и < Р , (20)

где Р* - эффективная пористость почвы ( с учетом доли пор, занятых защемленным воздухом), уравнение

Я0'<Г0> <21>

где .Тф - интенсивность водоотдачи снежного покрова. Если ке условие (20) нарушается, замыкающим уравнением является

Л + 5 - Р , (22)

В этом случае q0 < <Га (поверхность почвы затоплена и имеет место водообразование).

Решение уравнений (1б)-(19) и (21) или (22) для заданных моментов времени дает возможность рассчитать динамику интенсивности впитывания и общее количество впитавшейся за период снеготаяния воды.

Автором была разработана методика информационного обеспечения предлагаемой модели впитывания талых вод в мерзлую почву. Прежде всего это касается ^задания основных тепло- и гидрофизических параметров почвы ин, а^ Р*, Я(и,Л), Я(и,Л),

Рассмотренная методика расчета впитывания и была апробирована на материалах наблюдений Подмосковной водно-балансовой станции и данных экспериментальных исследований, проведенных нами совместно с сотрудниками ИГ АН СССР на Курской биосферной станции. Сопоставлялись рассчитанные и наблюденные гидрографы склонового стока (в сипу малости размеров стоковых площадок полагалось равенство интенсивности водообразования и склонового поверхностного стока): на двух полевых площадках подмосковной воднобалансовой станции с 1971 по 1983 г.г. и четырех -Курской биооферной станции с 1979 по 1985 г.г.Кроме того, для Курской биосферной станции было проведено сопоставление рассчитанных и измеренных с помощью электротермометров температурных профилей почвы в период снеготаяния.

Временной расчетный шаг был равен суткам. Величины тепло- ■ физических параметров почвы приняты на основе' результатов исследований, проведенных Голицыной др.(1972), Павловой, Голицыной (1974), Калюкным, Павловой (1981), с учетом имеющейся на станциях информации б влажности и плотности почвы (средних по слою 0-30 см).

Такие гидрофизические параметры почвы как ее коэффициент фильтрации Я0 и эффективный .капиллярный потенциал Як для Подмосковной водно-балансовой станции оценивались на основе результатов специальных исследований впитывания воды в талую почву, проведенных на станции в 1961 г., и были приняты Яо=0,3 мм/мин, Вк=10 см.

Значения К0 для площадок Курской биосферной станции с различной подстилаицей поверхностью принимались на основе результатов экспериментальных определений этого параметра (Ясинский, 1982), а Ик для черноземной почвы - на основе собственных исследований автора [Гусев, 198161.

На рис.6 приведены рассчитанные и измеренные данные по

динамике температуры в почвенном профиле.

1° 1-5 -Ц- '3 -<? -1 О Рис.6. Динамика температуры почвы (Курская биосферная станция, 1980). 1-4 - экспериментальные данные соответственно на 25 и 31 марта (начало снеготаяния). 2 и 4 апреля; & -рассчитанные данные на те же сроки; 6 - рассчитанный профиль температуры на I апреля (начало водоотдачи снежного покрова).

На рис.7 приведены промеры рассчитанных и измеренных

гидрографов поверхностного склонового стока. В качестве пока-

■ *

Г&ШО

лг «иг , ¡лж хж

• I

ш

МЛ *1ж

ша иж ма

аж

Рис.7. Рассчитанные (I) и измеренные (2) гидрографы поверхностного стсАа с различных площадок Подмосковной водно-балансовой (а) и Курской биосферной (б) станций.

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

- 2 6 "

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

станции гк=0,82 , а/*Фв=0,60 , для Курской биосферной - соответственно 0,85 и 0,52. Коэффициент корреляции между рассчитанными и измеренными слоямЬ стока половодья равен 0,95, сред-неквадратическое отклонение - 13 мм.

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

Суть упрощений заключалась в применении средних по времени оттепели характеристик гкдротермического режима почвы. Сопоставление рассчитанных и измеренных объемов впитавшейся воды было проведено на материалах наблюдений Подмосковной (1956-1974 Г.Г.), Нижнедевицкой (1955-1974 Г.Г.) И ПрИДвСВЯЕ-ской (1952-1974 г.г.)водаобалансовых станций, разумеется, для тех лет, когда оттепели имели место.

Среднее квадратичное отклонение суммарных за период оттепели рассчитанных и измеренных слоев стока для всех трех указанных станций при общей за оттепель водоотдаче снежного покрова , лежащей в диапазоне от нескольких миллиметров до нескольких десятков миллиметров, составляет 7-8 ММ. С ТОЧГС1 зрения оценки баланса почвенных вод за зимне-весенний период такая точность приемлема.

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

ГЛАВА 5. НЕКОТОРЫЕ ПРИВШ ЦЕЛЕВОГО АНТРОПОГЕННОГО ВОЗДЕЙСТВИЯ НА ШМ ПОЧВЕННЫХ ВОД.

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

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

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

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

Основным методологическим прием«! выявления эффекта от применения кулисования в данной работе являлось сравнение пополнения запасов почвенных вод, рассчитанных для объекта с

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

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

В связи с этим были разработаны подходы к возможному уче ту пространственной стохастичности некоторых из указанных выше параметров СГусев, 19780, 1902а,б, Гусев," Ясинский, 19871. Один из таких подходов основан на методе статистических ансамблей. Он представляет собой теоретическую разработку учета пространственной стохастичности параметров почвы, в основу которой положено уравнение эволюции не конкретной гидрофизической характеристики режима почвенных вод, например, влажности почвы, а функции ее распределения [Гусев 1982а,б).

Однако, такой вариант учета пространственной стохастичности параметров природных объектов оказался более удобным*при описании процессов вертикального переноса воды: Для учета горизонтальной пространственной неоднородности гидрологических Факторов, что и требуется при оценке эффективности кулисова-ния, нами был разработан и использован другой подход, основанный на применении метода Монте-Карло [Гусев, 19786; Гусев, Ясинский, 1987). Метод базируется на обработке серии решений, полученных на основе детерминированной модели, в которой значения закладываемых в модель параметров определяются с помощью

связи с коэффициентом вариации этого параметра.

В случае наличия кулис склон раебивается на участки чере-~ дующихся полос с кулисами и межкулиснын пространством. Имитация расположения снегозапасов по склону осуществлялась на основе гипотезы о распределении снегозапасов для каждого участка как авторегрессисяных последовательностей: одной - для участков с кулисами; другой - для межкудисаых. Необходимые для розыгрыша статистические характеристики последовательностей находились на основе данных о распределении снежного покрова на полях с кулисами, приведенных в (Калганый, Суаков, 1966).

Остальные детали алгоритма расчета весеннего пополнения влагозапасов почвы (Яш таким же. как и в случае отсутствия кулис. Необходимые для расчетов данные о температуре воздуха и осадках были взяты на станции Петринка, находящейся вблизи Курской биосферной станции Ж1 АН СССР.

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

Одним из звеньев технологии минимальной обработки почвы является мульчирование ее поверхности растительными остаткам;! (оставление стерни, покрытие соломенной мульчей, торфом я т.д.). С позиций использования ресурсов почвенных вод мульчирование является очень эффективным способом уменьшения непро7 дативного испарения воды почвой и повышения водообеспеченнос-тц растительности. 1

Предварительные и очень осторожные оценки, проведенные для степной и лесостепной зон ЁТС (Будаговский, Григорьева, 1991), показывают, что северо-западнее линии, проходящей примерно '-от шной границы Молдавии через ' Полтаву, Ульяновск, • Бугульму, а также для значительной части Северного Кавказа, лежащей несколько менее Маныча, применение мульчирования позволяет для 8-9 лет из . 10 достигнуть оптимальной влагообеспе-ченности посевов зерновых культур.'-Следовательно, мульчирование является очень дешевой и экологически благоприятной альтернативой их орошения.

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

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

В основу метода оценки влияния мульчирования на температурный режим почвы в зимний период [Гусев, Ясинский. 1990] положено математическое моделирование теплопереноса в системе снежный покров - мерзлая почва - талая почва как при наличии между снежным покровом И мерзлой почвой теплового сопротивления мульчслоя, так и при его отсутствии. Основные положения моделирования промерзания почвы и формирования ее температурного профиля в мерзлой зоне в-отсутствии мульчирующего покрытия были изложены выше. Учет мульчи осуществлялся на основе введения в модель дополнительного термического сопротивления между почвой и снежным покровом.

Созданные модели позволили осуществить математическое моделирование динамики промерзания почвы и формирования в ней термического профиля как при отсутствии, так и при наличии на ее поверхности мульчирукщего покрытия. Результаты моделирования сопоставлялись с данными экспериментальных исследования термического режима почвы, которые проводились ЗИМОЙ 1968-1989 г.г. в районе Курской биосферной станции ИГ АН СССР на трех участках: паровом, паровом, покрытом соломенной мульчей, и участке с естественным травяным покровом степной зоны.

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

- (2,4 t ОД; • 1СГ* кал/fc см °С).

Наличие данных о высотах и коэффициентах теплопроводности слоев соломенной мульчи и войлока отмершей естественной травяной растительности позволило провести расчеты динамики промерзания почвы и формирования в ней температурного профиля. Результаты сопоставления рассчитанных и экспериментально определенных характеристик термического режима почвы приведены на рис.8 и 9.

Для оценки влияния соломенной мульчи на изменение

- зг -

f ^бо'1 I-I-1-г

го

«0 «

° ч» , _<

Рис.8. Динамика промерзания почвы на открытом пару (а), пару, покрытом соломенной мульчей (0), участке с отмершим естественным травяным покровом (в). I - данные по мерзлотомзру Данилина; 2 - данные по электротермометрам ТЭТ-2; з - рассчитанные данные.

Рис.9. Профили температуры в мерзлой зоне почвы на открытом пару (I, 4), пару, покрытом соломенной мульчей (2, 5), естественном степном участке (3, 6). I, 2, 3 - экспериментальные данные; 4, 5, 6 - рассчитанные данные.

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

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

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

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

Построенная в результате модель передвижения верхней границы оттаивания С имеет следующий вид

«Г> - -г Ы

а *

(23)

где С - высота слс* мульчи, - аффективная теплопроводность мульчи, ^ - теплопроводность талой почвы, Рв - температура воздуха, Гп - температура поверхности почвы, с3 - теплоемкость талой почвы.

Решение задачи о нахождении профиля температуры почви после ее оттаивания выглядит следующим образом

г - а

+ Т(0)?

—. ^

2 " Св

+ А(г-а)г приОсжо пряо< г<1

(24)

А е

К V

о"*«;

12 а^

7Гх

от

<га •) « 1

Л1

о «

£

V ' С„ • хэ/сэ"

где гз - продолжительность зимнего периода, Ск - глубина промерзания почвы в конце зимнего периода, а момент 1=0 соответствует моменту полного оттаивания почвы.

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

- -

схода снежного покрова. В связи с этим были проведены специальные полевые исследования для его независимого определения.

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

Оценка с^ для соломенной мульчи проводилась на основе экспериментов по определению испарения с оголенной смоченной почвы Е0 и с почвы, покрытой слоем мульчи £ , с использованием следующего соотношения:

«О I *о _ ^

I Ф* в и--

Ср * 1 ф* I Е .

(25)

где 0 - коэффициент конвективного массообмена между подстилающей поверхностью и стандартной высотой наблюдений метесэле-

А

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

Указанные эксперименты были проведены в районе Курской биосферной станции ИГ АН СССР в течение двух летних сезонов (1966, 1989 г.г.). Всего было поставлено около 30 серий опытов при различной высоте слоя соломенной мульчи и при разных значениях средней за период наблюдения скорости ветра (от 0,5 до 2,5 м/с). Проведенные измерения показали, что указанные выше параметры связаны с коэффициентом турбулентного обмена между подстилающей поверхностью и стандартной высотой измерения метеоэлементов гс Г2м) 0 следующими соотношениями

где Д « 0,8 см, а « й'О, где к'0,24 кал/смг°С.

- 2

-----4

О «

Т/С

Рис.10. Динамика среднесуточной температуры почвы в степи (а), на участке, покрытом соломенной мульчей (б), и на пару (в). I - температура воздуха, 2,3 - рассчитанная температура почвы на глубинах соответственно 10 и 40 см, 4,5 -измеренная температура почвы на тех же глубинах.

В результате оказалось, что при скоростях ветра характерных дня марта, апреля в районе Курской биосферной в 2 • Ю'л кал /(см с °С), что и было принято при проведении модельных расчетов динамики температуры почвы под слоем мульчи с целью сопоставления получаемых результатов с натурными данными.

Найденное' значение коэффициента эффективной теплопроводности мульчи дало возможность провести. расчеты термического режима почва после схода снежного покрова на основе рассмотренных вше моделей для конкретных объектов и провести сопоставление расчетных данных с данными непосредственных наблюдений. В качестве таких объектов были взята три площадки в районе Курской биосферной станции: с оголенной почвой, с почвой, покрытой слоем соломенной мульчи с указанными выше характеристиками, а также на степной площадке, покрытой естественным войлоком отмершей травяной растительности с плотностью и 1 кг / м2 и эффективной высотой - 7 см . Коэффициент для отмершей растительности естественного травяного покрова принимался таким же, как и для соломы.

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

Основные выводы

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

нять и увеличивать свою упорядоченность.

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

3. Общее уравнение переносе воды в ненасыщенных средах, полученное на основе статистической теории неравновесных процессов, является теоретическим обоснованием используемых в настоящее время методов описания движения воды в ненасыщенной почве. В случае слабонестационарных процессов переноса воды в почве для их описания возможно использование уравнения Дарси-Букингема, являвдегося частным случаем'полученного уравнения.

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

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

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

ности протекания указанных процессов и могут быть использованы для оценки влияния мульчи на формирование' гидротермического • режима почвы в зимне-весенний период. > "

7. Экспериментально определенные тетьлогидрофизические параметры'соломенной мульчи служат информационной основой как при оценке эффективности мульчи в качестве депрессора непродуктивного испарения «так и при анализе ее влияния на уменьшение амплитуд колебаний термического режима почвы. Установлено, что применение слоя соломенной мульчи толщиной 5-10 см создает гидротермический режим почвы, близкий к режиму естественной' степной экосистемы.

Основные результаты работа отражены в следгаяих публикация!.

1. Гусев Е.М. Зависимость интенсивности впитывания.веды в почву от ее начальной влажности// Водные ресурсы. 1978. №. С.45-52, (0,8 п.л.)

2. Гусев Е.М. Влияние горизонтальной неоднородности коэффициента фильтрации почвы на интенсивность впитывания// Метеорология и гидрология. 1976. * 7. С. 66-73. (0,9 п.л.)

3. Гусев Е.М. Об уравнении нестационарного переноса влаги в -почвогрунтах // Доклады ВАСХНИЛ. 1979. №3. С.39-41. (0,3 п.л.)

4. Гусев Е.М. Исследование процесса впитывания дождевых осадков. Автореф. канд. техн. наук. Москва, Агрофизический ШИ, 1979. 21 С. (1,5 П.Л.)

5. Гусев Е.М. Об одном уравнении изотермического движения воды в почве // Почвоведение. 1979. ЗЙ7. С.61-68. (0,7 п.л.)

6. Гусев Е.М. Методика расчета впитывания и поверхностного стока на единичном склоне с учетом пространственной изменчивости фильтрационных свойств почвы // География и природные ресурсы. 1960. Я. С.161-163. (0,3 п.л.)

7. Гусев Е.М. Проблемы теории переноса жидкости в ненсыщен-ных пористых средах // Физика почвенных вод. М.: Наука, 1981. (1,7 п.л.)

6. Гусев Е.М. Экспериментальное исследование напорного впитывания// Физика почвенных вод. Ы.: Наука, 1961. 0. 195-2С5. (1.5 п.л.)

9. Гусев Е.М. Вариант динамико-стоха'стического моделирования гидрологических процессов' // Метеорология и. гидрология. 1962. Ж. (0,8 п.л.) . .

10. Гусев Е.М. Стохастическое описание водного обмена на сельскохозяйственном поле//Доклада ВАСХНИЛ. 1962. N10,(0.5 п.л.)

11. Гусев Е.Ы., Будаговский А.И. определение штоков воды в ненасыщенных средах // Вод&е ресурсы. 1962. *2. (1,0 п.л.)

12. Гусев Е.М. Определение релаксационной функции ненасыщенной почв« // II-ая Всесоюзная конференция по применению математических методов и ЭВМ в почвоведении. Пущино. 1983.(0,2 п.л.)

13. Гусев Е.М. Определение релаксационной гидрофизической функции ненасыщенной почвы // Почвоведение. 1983. С.51-56. (0,5 п.л.) •

14. Гусев Б.М. Приближенный численный расчет глубины промерзания почвы // Метеорология и гидрология. 1965. # 6. С. 94102. (0,9 п.л.)

15. Гусев Е.М., Ясинский C.B. Методика расчета динамики термического режима почвы в зимний период// Материалы метеорологических исследований. М.: Междуведомственный геофизический КОМИТЕТ АН СССР, 1966. * 10. С. 87-100. (1,3 П.Л.)

16. Гусев Е.М., Ясинский C.B. Моделирование склонового стока талых вод с учетом пространственной вариации снегозапасов // Изучение и оптимизация водных ресурсов Центральной лесостепи. Московский филиал Географического общества СССР. Курск, 1987. (0,6 п.л.)

17. Гусев Б.М. Уйрощенный расчет глубины промерзания почвы с учетом миграции воды к границе промерзания // Водные ресурсы. 1968. * I. С. 33-42. (1,0 п.л.)

16. Гусев Е.М. Впитывание воды в почву в период снеготаяния// Водные ресурсы. 1989. M 2. С. 19-31. (1,3 п.л.)

19. Гусев Е.М., Будаговский А. И. Почвенные воды // Водные ресурсы. 1989. N5. С.16-27. (1,2 п.л.) - •

20. Гусев Е.М. Экологическая роль почвенных вод и их ресурсы // Водные ресурсы. 1990. Jt 5.'C.II0-I2I. (1,3 п.л.>

21. Гусев Е.М., Ясинский C.B. Влияние мульчи из растительных остатке^ на формирование'термического режима почвы в зимний период // Почвоведение. 1990. *12. С.46-54. (0,9 п.л.)

22. Gusev Е.М. Modelling or water transfer In soils Юг the winter-spring period // IAHS Publication N 204. Hydrologlcal Interaction Between Atmosphere, Soil and Vegetation. 1991. p.235-244. (0,7 п.л. )

- 4 1 -

23. Cusev Ye.M., Busarova O.Ye., Yaslnsky S.V. Impact of mulch layer composed by organic remains on the soil hydro thermic regime in the water-spring period // Annals Geophysical. Eart II. Oceans, atmosphere, hydrology & Nonlinear geophysics. Supplement II to Volume 10. Springer 'International. 1992, P.273. (0,1 n.ji.)