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

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ МИГРАЦИОННЫХ ПРОЦЕССОВ В ЛАНДШАФТНО-ГЕОХИМИЧЕСКИХ СИСТЕМАХ

ВАК РФ 25.00.23, Физическая география и биогеография, география почв и геохимия ландшафтов

Автореферат диссертации по теме "МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ МИГРАЦИОННЫХ ПРОЦЕССОВ В ЛАНДШАФТНО-ГЕОХИМИЧЕСКИХ СИСТЕМАХ"



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

КОШЕЛЕВА Наталья Евгеньевна

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ МИГРАЦИОННЫХ ПРОЦЕССОВ В ЛАНДШАФТНО-ГЕОХИМИЧЕСКИХ СИСТЕМАХ

25.00.23 — физическая география и биогеография, география почв и геохимия ландшафтов

АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора географических наук

Москва 2003

Рабата выполнена на кафедр« геохимии ландшафтов и географии почв географического факультета Московского государственного университета им. М.В .Ломоносова

член-корреспондент РАН, доктор географических наук Касимов Н.С.

доктор географических наук В исто ров А.С.

доктор физико-математических наук, профессор Крапивин В.Ф.

доктор биологических наук Зейлкгер A.M.

Ведущая организация: Институт географии РАН

Защита состоится «__»_2003 г. в 15 часов на заседании

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

С диссертацией можні <v- <:-усилиться б библиотеке географического факультета МГУ на 21-м этаже.

Автореферат разослан _ _ 2003 г.

Научный консультант;

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

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

Горбунова И. А.

ВВЕДЕНИЕ

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

Количественное изучение условий, форм и интенсивности миграции и трансформации разнообразных химических веществ в ландшафтах привело к попыткам создания математических моделей, в той или иной степени отражающих теоретические представления и наблюдаемые особенности поведения элементов в почвах, поверхностных и подземных водах зоны гипергенеза (Козловский, Корнблюм, 1972; Пачепский, 1976, 1990; Гильманов, 1978; Сысуев, 1986, 2003; Рыжова, 1987, 1992; Козловский, 1991; De Vries et aL, 1991,1997; Parton, 1996; Назаровидр., 1999,2001; Крайнов, 1999 и многие др.).

Методологической базой исследований миграционных процессов является системная концепция Б.Б. Полынова о структурообразующем значении потоков вещества в геохимических ландшафтах. Водным потокам в большинстве природных зон принадлежит ведущая роль: вода является носителем многих веществ и той средой, где протекают разнообразные физико-химические и биологические процессы их трансформации. Поэтому для нынешнего этапа развития математического моделирования миграционных процессов характерны поиски наиболее адекватных и эффективных способов формализации и интеграции геохимических знаний в модели массоперсноса (Schnoor, 1996), Основой для интеграции являются модели гидр о логического цикла, описывающие взаимодействие водных объектов с природными комплексами и влияние различных природных факторов и структур на интенсивность процессов миграции (Linsley, Crawford, 1960; Jonch-Clausen, 1979; Кучмент и др., 1983; Назаров и др., 1988, 1996, Кондратьев, 1992; Ландшафтно-гидро логический анализ..., 1992; Корытный, 2001).

Развитие количественной теории миграционных процессов в ЛГС локального уровня связано с разработкой математических моделей, воспроизводящих основные типы геохимических сопряжений в ландшафтах (Глазовская, 1964). Не менее важна разработка моделей, описывающих геохимическую нагрузку миграционных потоков, а также физико-химические превращения и биологические процессы в различных ландшафтных условиях. Их объединение позволит описать перемещение, накопление, химические взаимодействия и распад различных веществ и таким образом проследить их судьбу s компонентах ландшафта. Этой проблеме и посвящена диссертационная работа.

Обосновать и практически

Цель работы - разработать", теоретически <

использовать систему математических1 Йодё&й, отраж' ■ фонд научной литературу

4v U^j^CJ^xJZLxD-—

иощих процессы движения t

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

1. Провести классификацию моделей миграции и трансформации вещества в ландшафтах с учетом структуры ЛГС, природы ландшафтно-геохнмнческих процессов и свойств участвующих в них веществ,

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

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

4. Построить комплексные модели влаго- и солепереноса в ЛГС степной зоны; путем численных экспериментов с моделями обосновать оптимальный вариант использования водных и земельных ресурсов.

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

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

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

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

соотношения их пространственно-временных масштабов.

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

Объекты исследования и исходные материалы. .Для реализации моделей были выбраны хорошо изученные объекты, по которым накоплен значительный фактический материал. В пределах лесной зоны моделировался водный цикл бассейна р. Медвенки, левого притока р. Москвы, на территории которого расположена Подмосковная воднобалансовая станция. Часть параметров была определена по данным многолетних наблюдений Истринского опорного пункта ВНИИЛМ, а также Проблемной лаборатории охраны геологической среды МГУ и ПГО "Центргеология", в которых автор принимала участие В 1982-1985 гг.

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

Исследования миграции и трансформации тяжелых металлов проводились на территории Московской и Калужской областей. Моделирование подвижности тяжелых металлов в почвах Смоленско-Московской возвышенности основано на материалах, собранных O.A. Самоновой и В.Б. Вагнером на территории Сатин-ского научно-учебного полигона МГУ в 1991-1995 гг. Для реализации балансовой модели круговорота тяжелых металлов в агроландшафтах использованы материалы полевых опытов Центральной опытной станции ВНИИ удобрений и агро-почвоведения, проводимых H.A. и И.Н. Черных, а также литературные данные.

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

1, Разработаны теоретико-методологические основы математического моделирования ландшафтно-геохимическнх процессов в ЛГС локального уровня, включающие

- классификацию математических моделей миграции и трансформации вещества в ЛГС, учитывающую пространственную организацию ЛГС, характер процессов и свойства участвующих в них веществ;

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

трансформации вещества в ЛГС;

- принципы построения комплексных моделей;

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

2. На основе интеграции частных моделей лакдшафтно-геохимнческих процессов созданы комплексные модели биогеохимических циклов;

- водного цикла в каскадной системе 1-го порядка лесной зоны;

- водного и солевого режима оросительной системы в степной зоне;

- формирования качества воды в водохранилище лиманного типа;

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

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

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

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

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

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

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

Результаты моделирования были использованы в проектных разработках Нижне-Днестровской и Явки не кой оросительных систем на юге Украины (проектирующий институт Укрюжгипроводхоз). Компьютерная программа, реализующая динамическую модель формирования поверхностного и подземного стока, использовалась при создании постоянно действующей модели геологической среды зоны влияния Московского градопромышленного комплекса. Результаты исследований отражены в базе данных по состоянию территориально-аквальных комплексов малых и средних рек Московского региона (Per. свид. №7782 в Государственном регистре баз данных «Информрегистр»),

Основные направления исследования были поддержаны Государственной программой "Экологическая безопасность России" (1993-1995), Международным фондом Сороса (1994-1996), Российским фондом фундаментальных исследований (1999-2001), Федеральной программой "Университеты России -фундаментальные исследования".

Материалы диссертации нашли отражение в курсе лекций "Основы математического моделирования" для аспирантов ВНИИ удобрений и атропочвоведения (1989-1992) и используются автором при чтении курсов "Моделирование почвенных и ландшафтно-геохимических процессов" и "Математические методы обработки ландшафтно-геохими ческой информации" на географическом факультете МГУ, Разработанные в ходе исследования теоретические и методические вопросы вошли в учебное пособие "Моделирование почвенных и ландшафтно-геохимических процессов" (1997).

Апробация работы н публикации. Основные результаты работы докладывались, обсуждались и опубликованы в трудах международных, всесоюзных, региональных конференций, совещаний и семинаров: "Математические вопросы динамики атмосферы и гидросферы" при МГУ (МГУ, 1977; 1981); "Повышение эффективности использования орошаемых земель" (Ставрополь,

1980); "Математическое моделирование гидрогеологических процессов" (Москва,

1981); "Системные исследования водных проблем" (Вороново, 1986); заседания гидрогеологической секции МОИП (Москва, 1985, 1987); Всесоюзный научный семинар по использованию подземных вод на орошение (Баку, 1986); "Статистические методы в теории передачи и преобразования информационных сигналов" (Киев, 1988); II Международный симпозиум по методам математической обработки агрохимической информации (Германия, Эберсвальде, 1990); "Применение ЭВМ при гидрогеохимическом моделировании" (Ленинград, 1991); "Инженерная экология-91" и "Проблемы эконнформатики" (Звенигород, 1991, 1992); 'Тяжелые металлы и радионуклиды в агрожосистемах" (Немчнновка, 1992); XX и XXII ассамблеи Европейского Геофизического союза (Германия, Гамбург, 1995; Австрия, Вена, 1997); 3-ая международная конференция по биогеохимии (Фракция, Париж, 1995); "Diffuse Pollution" (Чехия, Прага, 1995); "Тяжелые металлы в окружающей среде" (Пущино, 1996); 6-ая международная конференция В.М. Гольдшмидта (Германия, Гейдельберг, 1996); "Теория и техника передачи, приема и обработки информации" (Харьков-Туапсе, 1996); "Технологии конструирования адаптивных агроэкосистем и агроландшафтоа" (Лобня, 1998); 16-й Всемирный конгресс по почвоведению (Монпелье, Франция, 1998); 5 и 7-ая Международные конференции по биогеохимии микроэлементов (Австрия, Вена, 1998; Швеция, Уппсала, 2003); Ломоносовские чтения (Москва, МГУ, 1999; 2001); "Геохимия биосферы" (Новороссийск, 1994,1999,2001); "Роль почвы в функционировании экосистем" (Польша, Люблин, 1999); "География на рубеже веков: проблемы регионального развития" (Курск, 1999); "Университеты России - фундаментальные исследования: география" (Новороссийск, 2000); III съезд Докучаевского общества почвоведов (Суздаль, 2000); 29-й Международный географический конгресс (Сеул, Корея, 2000); "Функции почв в геосферных-биосферных системах" (Москва, 2001); "Геохимическая экология и биогеохимическое изучение таксонов биосферы" (ГЕОХИ РАН, Москва, 2003); "Современные методы эколого-геохимической оценки состояния и изменений окружающей среды" (Новороссийск, 2003).

По теме диссертации опубликовано более 50 печатных работ.

Структура и объем работы. Диссертация состоит из введения, б глав, выводов, списка литературы и пяти приложений. Работа изложена на 429 страницах, содержит 38 таблиц, иллюстрирована 59 рисунками. Библиография включает 525 наименований.

Автор выражает глубокую признательность научному консультанту работы члену-корреспоиденту РАН Н.С. Касимову. Автор искренне благодарен за ценные советы, поддержку и критические замечания профессору М.А. Глазов-ской. Автор признателен коллегам, во взаимодействии с которыми в разные годы проводились исследования, нашедшие отражение в диссертации: Н-И. Алексе-евскому, В.Б. Вагнеру, М.И, Герасимовой, Д.Л, Голованову, В.А. Жуку, И.Ф.Комиссаровой, В.Ф. Ладонину, С.Н, Лебедеву, Е.М. Мироновой, В.Н. Михайлову,

Н,А. Назарову, А.И. Новикову, И.С. Пашковскому, Т.Ф. Прилепиной, O.A. Само-

новой, ИЛ. Храновнчу, И.Н. Черных, RA. Черных, Е.В. Шеину, В.П. Яковлеву.

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

Глава I. Теоретические основы моделирования процессов в ландшафтно-геохнмнческнх системах

Практическим инструментом исследований ландшафтао-геохимических процессов является системный анализ, который на основе количественного описания связей и отношений в ЛГС позволяет построить их математические модели. Элементы системного анализа стали использоваться еще в 40-50-е годы прошлого века в работах Б.Б. Полынова, в 70-е годы он окончательно утвердился благодаря работам М.А. Глазовской о ЛГС и А.И. Перельмана о видах и формах миграции элементов в ЛГС. В настоящее время системный подход к изучению географической оболочки в целом и ее отдельных компонентов прочно вошел в теорию и практику географических исследований (Роде, 1971; Сочава, 1975; Проблемы регионального географического прогноза, 1982; Фридланд, 1984; Пузаченко, 1985; Мотузова, 1999 и др.). Изучение процессов функционирования геосистем через дифференциацию и миграцию химических элементов позволило выявить закономерности пространственного распределения твердой фазы вещества, особенности динамики жидкой, газообразной и живой фаз, разработать методы ландшафтно-геохимического анализа, синтеза, диагностики состояний, прогноза и оценки геосистем (Сысуев, 1986, 2003; Ландшафтно-геохимические основы,.., 1989; Дьяконов и др„ 1995) и применить их для оптимизации природопользования (Моделирование роста ..., 1986; Солнцева, 1988; Пегов, Хомяков, 1991; Биогеохимические основы..., 1993; Елпатьевский, 1993; Бе Упез ег а1., 1997; Геоэкологическое прогнозирование ..., 2002).

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

В зависимости от природы процессов в ЛГС их математические модели объединены нами в четыре группы с собственным математическим аппаратом (табл.1). Модели отдельных фаз водного цикла описывают движение миграционной фазы (воды, воздуха). Транспортные модели позволяют рассчитать перенос растворенных и взвешенных элементов и их соединений в водных и воздушных потоках, а также механический перенос с головым или гравитационным материалом, солифлкжционными массами. Термодинамические и кинетические модели отражают разнообразные превращения веществ в процессе миграции. На основе термодинамических уравнений рассчитывают равновесные соотношения между исходными реагентами и продуктами химической реакции. Если равновесие недостижимо, с помощью кинетических уравнений моделируется скорость реакции или физико-химического процесса. Для элементов-биофилов важное значение имеет биологическая миграция, которую отражают модели продукционных процессов и разложения органики, а также модели биологического круговорота веществ.

Таблица 1

Математические модели миграции и трансформации вещества вЛГС

лг системы Модели отдельных фаз водного цикла Транспортные модели Термодинамические и кинетические модели Мол ели биологической миграции Комплексные модели

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

Катены • процессов вертикального влагопереноса в ЭЛГС; • формирования склонового стока. • латерального переноса веществ в виде растворов и коллоидов; • водной и ветровой эрозии; • солнфлгокции; • гравитационного переноса в твердой фазе. • химических равновесий; • скорости химических реакций и фюи ко-химических процессов в латеральных потохах. • продукционных процессов; • разложения органического вещества; • биологического круговорота и пищевых цепей в ЭЛГС и их сопряжениях. биогеохимического круговорота

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

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

В классификации моделей три уровня пространственной организации геохимических ландшафтов: элементарные системы (ЭЛГС), катены и каскадные системы (КЛГС), в ряду которых возрастает степень открытости и уменьшается роль биологической миграции (Снакин, 1987). В ЭЛГС доминирует вертикальный влаго- и массоперенос, обусловливающий биогеоценотический круговорот веществ. Его незамкнутость является причиной латеральной миграции в катенах. В моделях катенарного уровня протекающие в ЭЛГС процессы описывают упрощенно, интерес представляет лишь та часть веществ, которая выносится из элементарных ландшафтов и вовлекается в биогеохимический круговорот. Модели КЛГС воспроизводят миграционные процессы в совокупности катен и гидрографической сети, а также в бассейнах подземных вод. Взаимодействие в каскадных ЛГС является причиной геохимического круговорота, основную работу в котором выполняют водные потоки.

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

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

Глава II. Методы моделирования миграции и трансформации вещества в ландшафтно-геохимическнх системах

На основе классификации моделей миграции и трансформации вещества в ЛГС в главе проанализированы различные подходы и обобщен опыт моделирования отдельных видов ландшафтно-геохимических процессов и построения комплексных моделей.

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

Процессы вертикального тепло- и влагопереноса в ЭЛГС описывает семейство^ моделей с общим названием SWAP (Soil-Water-Atmosphere-Plant). Этте модели в динамике воспроизводят компоненты водного баланса на поверхности раздела земли и атмосферы: физичесхое испарение и транспирацию, перехват осадков растительностью, рост растений, поверхностное задержание, впитывание воды почвой и поверхностное стекание, влагоперенос в почве и локальные колебания уровня грунтовых вод (УГВ).

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

Модели КЛГС, как правило, описывают миграцию в катенах и русловое движение. Отдельно моделируется миграция подземных вод. Гидрогеологические модели строятся на основе схематизации, в ходе которой определяются условия питания и разгрузки, а также выявляется структура потока подземных вод. Условия формирования подземного потока представляются в фильтрационной форме с обоснованием исходной гидродинамической модели процесса (Шестаков, 1995), которое включает выяснение режима и мерности течения, задание граничных условий и параметров уравнения геофильтрации. Чаще всего геофильтрационные потоки схематизируются как плоские потоки в вертикальном сечении (одномерное течение) или в плане (двумерное плановое течение).

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

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

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

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

вещества). Конвективно-дисперсионное уравнение образует теоретическую основу математического моделирования процессов миграции многих химических элементов и их соединений в компонентах как элементарных, так и КЛГС: почвенных растворах, грунтовых водах, поверхностном и речном стоке, в водоемах.

Когда миграция веществ в ландшафтах происходит преимущественно в форме взвесей, используют модели водной эрозии. При отсутствии общепринятых физически обоснованных моделей эрозии количественная оценка эродирующей способности водных потоков проводится на основе эмпирических уравнений. Как правило, эти модели связывают количество смьггого материала с морфом етрическнм и характеристиками склона, типом растительности и свойствами почв и слагающих склон грунтов, а также с параметрами дождя (Швебс, 1974, 1981; Ларионов, 1993). Универсальное уравнение эрозии представляет собой произведение основных влияющих на процесс факторов. Модели размыва русла и транспорта взвешенных и донных наносов в каскадных Л ГС учитывают гидравлическую крупность частиц, слагающих русло, и транспортирующую способность потока (Мирцхулава, 1976; Ермаков, Мартынов, 1984; Алексеевский, 1998).

Комбинированные транспортные модели описывают вынос веществ с поверхности водосбора как в растворенном, так и во взвешенном виде. Чаще всего они находят применение в оценке загрязнения поверхностного стока из рассредоточенных источников (Назаров, 1996),

Моделирование переноса и диффузии примесей в атмосфере предполагает задание трехмерного поля скоростей воздушных потоков, которое можно определить, опираясь либо на данные наблюдений, либо на результаты расчетов по моделям циркуляции. Перенос взвешенных в воздухе частиц описывает уравнение турбулентной диффузии, которое в общем случае содержит члены, характеризующие трансформацию и взаимодействие различных веществ, а также их гравитационное осаждение (Ньистадт, Ван Доп, 1985; Берлянд, 1985; Марчук, Кондратьев, 1992; Белолипешшй, Шохин, 1997; Аргучинцев, Аргучинцева, 2001).

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

Кинетические и термодинамические модели. Физико-химические процессы оказывают заметное влияние на баланс многих веществ в ЛГС, сопоставимое с процессами массопереноса.

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

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

Для'ЛГС характерны параллельные реакции и цепочки последовательно протекающих реакций. Замкнутые цепочки описывают циклы химических элементов, таких как азот, углерод, сера, фосфор, тяжелые металлы (Schnoor, 1996; De Vriesetal., 1997).

Когда характерное время протекания химических реакций (например, комплексообразования и кислотно-щелочных реакций) в несколько раз меньше характерного времени транспортных процессов, с помощью термодинамических моделей рассчитываются равновесные концентрации ионов в растворах и в других термодинамических системах. Система нелинейных алгебраических уравнений относительно активностей ионов имеет в своей основе закон действующих масс. Ее дополняют уравнения стехиометрического баланса масс и баланса электрических зарядов. Результирующая система уравнений, как правило, решается итерационным методом Ньютона-Рафсона. Характеристика компьютерных программ, используемых при термодинамическом моделировании, дана в (Schnoor, 1996; Крайнов, 1999).

Изотерма адсорбции представляет собой простейшую и одновременно широко используемую термодинамическую модель, которая дает интегральную равновесную характеристику способности почвы к адсорбции. Обобщением изотерм на случай переменных физико-химических характеристик почвы являются многофакторные регрессионные модели ПОДВИЖНОСТИ химических элементов (Касимов и др., l995;De Vriesetal., 1997; Кошелева и др., 2002).

Термодинамические модели важнейших физико-химических процессов обеспечены необходимыми для расчетов параметрами — константами равновесия химических реакций и свободными энергиями веществ, участвующих в процессах (Наумов, Рыженко, ХодаковскиЙ, 1971; Батлер, 1973). Для кинетических моделей при усилении техногенной нагрузки на ландшафты требуются экспериментально установленные скорости превращений искусственных соединений в различных природных средах и условиях, некоторые из которых можно измерить в установках MICROCOSM.

Модели биологической миграции обычно реализуются на уровне ЭЛГС, где биологический круговорот веществ наиболее интенсивен. Он объединяет противоположно направленные процессы образования живого вещества и разложения органических остатков.

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

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

тепло- и водообеспеченность, запасы биогенных веществ. Для дикорастущих растений наиболее важны процессы, связанные с биогеохимическими циклами С, N, микроэлементов (De Vries et al., 1997) и формированием зональных запасов органического вещества (Рыжова, 1992), в том числе гумуса в луговых (Гильманов, 1978) и лесных (Сысуев, 1998) экосистемах.

В современных условиях заметное влияние на ростовые процессы и качество продукции оказывают техногенные выбросы, содержащие загрязняющие вещества. Их угнетающее действие моделируется с помощью функций, связывающих микробиологическую и ферментативную активность почв, продуктивность естественных ценозов, потери урожая и ухудшение его качества с концентрацией токсичных веществ, например, тяжелых металлов (Гузев и др., 1980; Кошелева и др., 1995; Черных, Овчаренко, 2002).

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

Частные модели продукционных процессов и разложения органических остатков объединяются в модели биологического круговорота химических элементов и их соединений в ЛГС, В их основе - балансовые уравнения относительно содержания элементов в почве или донных фунтах водных объектов (Богатырев, Рыжова, 1994).

Модели биологического круговорота позволяют исследовать как кратковременную динамику, так и многолетнюю и даже многовековую эволюцию биогеоценозов. Для круговорота элементов-биофилов характерен незначительный вынос вещества из системы, поэтому их балансы являются эффективным средством мониторинга природных и техногенно измененных ЛГС (Гильманов, Базилсвич, 1983; Geoderma, 1997, vol.81). Модели круговорота углерода используются для исследования устойчивости системы "почвенный гумус-растительность", что имеет большое значение при оценке плодородия почв и планировании агротехнических мероприятий (Рыжова, 1992).

Модели пищевых цепей отражают динами^ численности популяций и взаимодействия между различными трофическими уровнями (Schnoor, 1996; Назаров, Леонов, 1999). Результатом исследования таких моделей является определение устойчивых состояний системы (Пегов, Хомяков, 1991; Кошелева, Яковлев, 2002).

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

Комплексные модели. Увеличение геохимической нагрузки на ландшафты выдвигает на первый план комплексные модели, охватывающие полные

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

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

Для региональных территориальных комплексов с линейными размерами от десятков до первых сотен км в ПСА РАН созданы модели миграции и трансформации солей и загрязнителей (Геоэкологическое моделирование..., 2002), которые функционируют в рамках автоматизированной системы регионального экологического прогноза (АСРЭП), В основе моделей АСРЭП лежит концепция дискретной смены состояний биогеосистем, с помощью которой их динамика представляется в виде последовательной смены режимов функционирования. Однако при ориентации на минимальное информационное обеспечение модели АСРЭП не дают пространственно-временного распределения миграционных потоков в ЛГС,

Наименее развиты комплексные модели биогеохимических циклов для элементарных и каскадных ЛГС локального уровня.

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

Если изучаемые процессы образуют незамкнутый цикл, как, например, в моделирующей системе MIKE SHE (http://www.dhisoftware.com), он воспроизводится цепочкой моделей, в которой выход из одной модели является входом для другой. Когда некоторый процесс является результирующей нескольких параллельно протекающих процессов, последние моделируются независимо и их выходные функции служат входными данными для расчета результирующего процесса. Такая структура характерна для балансовых моделей (Кучмент и др., 1983; Кошелева и др., 1989; Пачепскнй, 1990; De Vrieset al„ 1997).

Если в системе проявляются сильные обратные связи, взаимозависимость между процессами может быть описана системой уравнений (Рыжова, 1992; Хомяков и др., 2000; Кошелева, Яковлев, 2002 и др.), В случае многокомпонентной системы прибегают к организации итерационного процесса, в ходе которого происходит согласование полученных в отдельных модулях решений и постепенно уточняются значения всех переменных (Simunek, Suarez, 1994).

Важную роль при создании комплексных моделей играет принцип согласования пространственного и временного масштабов изучаемых процессов, который предопределяет выбор размеров элементов ЛГС и соответствующих шагов по времени. Времени^ иерархия процессов при динамическом моделировании выступает как ограничение: в качестве основного объекта моделирования должны быть выбраны процессы в одном узком диапазоне характерных времен (Арманд, Таргульян, 1974). Более "медленные" переменные вводятся з модель в виде констант, а более "быстрые" - как шум, вызывающий случайные отклонения. Систематизация основных биогеоценогических процессов по величине характерных времен выполнена В.В. Сысуевым (1998).

Степень пространственного усреднения определяется исходя из специфики задачи и пространственной структуры изучаемого объекта. Как правило, последний представляется в виде совокупности относительно однородных участков (Кошелева и др., 1986; 1989; Vaughan, Corwin, 1994; Назаров, 1996). Дня некоторых приложений важны лишь осредненные значения выходных функций и характеристики варьирования. Тогда изучаемый объект может быть представлен в виде статистического ансамбля участков, для которых соответствующие факторы - случайные величины с известными статистическими параметрами (Назаров, Сирин, 1988; Гусев и др., 1998).

Способ усреднения параметров в пределах ландшафтно-однородных участков зависит ОТ характера их пространственных колебаний (Brus, de Gruijter, 1997). При независимости отклонений от среднего в соседних точках участка вычисляется среднее арифметическое (Pachepsky, Acock, 1998; Seyfrîed, 1998), при пространственной скоррелированности показателя для осреднения используется процедура блочного кригинга. Она сглаживает высокочастотное варьирование исходных данных, обусловленное ошибками измерений, с одной стороны, и варьированием на близких расстояниях (менее размера блока) - с другой (Кузякова и др., 2001).

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

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

- соответствия характерных размеров элементов ЛГС характерным временам миграционных процессов, определяющий выбор пространственно-временных масштабов описания процессов;

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

Глава III. Моделирование водного цикла в каскадных ландшафтно-геохнмических системах

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

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

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

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

где - плановые координаты (м); Т - водопроводимость (м2/сут.) основного горизонта (м); №" - переток из слабопроницаемых покровных отложений (м/сут.); кИ,тн - коэффициент фильтрации (м/сут) и мощность толщи, разделяющей основной и нижележащий горизонты (м); Иц - уровень воды в нижележащем горизонте (м); Р1 - водопроводимость, определяющая связь водоносного горизонта с рекой (сут."1); Нр - уровень воды в реке (м).

Зона аэрации так же представлена в виде двухслойной толщи, в которой положение УГВ Н,р определяется водообменом с нижележащим горизонтом IV и притоком к свободной поверхности У„:

<1НГ„ М-^-Уп-*.

где м " коэффициент гравитационной емкости. У„ вычисляется в зависимости от влажности почвы, коэффициента фильтрации и мощности зоны аэрации.

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

ЩНп {-х<У> ~ £(!'Х'У) ~ К'

где тп,{1п - мощность (м) и коэффициент гравитационной емкости активного почвенного слоя; & - (в-вд)/- капиллярная насыщенность активного слоя; во - максимальная молекулярная влажность. Для зимнего периода учитывается накопление осадков в снежном покрове, Е(:,х,у) оценивается по схеме с сосредоточенной емкостью. Меняющиеся от года к году погодные условия описываются вероятностной моделью климата, которая генерирует ряды средних месячных температур воздуха и осадков.

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

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

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

Путем численных экспериментов нами установлены относительные ошибки определения ¡V в зависимости от погрешностей определения напоров в скважинах. Получена формула для оценки математического ожидания этой ошибки, которая использована для определения числа скважин т в пределах каждого района. Установлено, что для исследованного случая, когда ошибка определения IV не превышает 20%, т~5. Эта методика дает возможность спроектировать режимную сеть скважин, обеспечивающую заданную точность оценки и рационально разместить их на территории.

Реализация модели в бассейне р. Медвенки. С помощью модели была получена детальная картина пространственно-временного распределения потоков поверхностных и подземных вод в пределах бассейна р. Медвенки, левого притока р. Москвы, и его изменения под воздействием природных и антропогенных факторов. Параметры и входные функции модели задавались до данным наблюдений на Истринском опорном пункте ВНИИЛМ и на Подмосковной воднобалансовой станции. На территории бассейна р, Медвенки выделено 92 ячейки размером 483x483 м (рис,1) с тремя видами землепользования: лесным, полевым (сельскохозяйственные угодья) и пойменным, различающимися по условиям миграции поверхностных и подземных вод.

граница бассейна

□-1 сп-2 а-з

ОТ-4

-5

-7

Рис.1. Результаты моделирования водного цикла в бассейне р. Медвенки: а) распределение по площади бассейна среднего многолетнего стока весеннего половодья; и б) модуля питания подземных вод (мм): 1 - 0-20, 2 - 20-40, 3 - 40-60. 4 - 60-80, 5 - 80-100, 6 -¡00-120, 7-Цг< О

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

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

Для описания водного режима элементарных ландшафтов и интенсивности латеральных потоков в КЛГС проведены численные эксперименты с моделью и определена миграционная структура КЛГС (рис.1) водосбора р. Медвенки. На вход модели подавались случайные реализации осадков и температур воздуха длиной в 25 лет. Установлено, что модуль питания грунтовых вод для всей территории водосбора практически не меняется во времени и составляет 61,2 мм/год. Средний модуль подземного стока составляет 54 мм/год при о=10 мм/год.

Распределение поверхностного склонового стока по территории водосбора неравномерно - большие модули стока весеннего половодья характерны для открытых участков, а максимальные модули дождевого стока приурочены к областям с близким залеганием УТВ. Полученное распределение позволило выявить места концентрации стока, приуроченные к тальвегам логов и балок. Изменчивость годовых объемов склонового и речного стока характеризуют эмпирические кривые обеспеченности, близкие к лог-нормальному закону распределения.

Разработанная модель была использована для вероятностных прогнозов изменения водного цикла КЛГС под влиянием хозяйственной деятельности. Был получен прогноз возможных изменений водного режима и структуры водного баланса бассейна р. Медвенки при сведении леса и распашке территории, при восстановлении леса и при работе водозабора подземных вод с расходом 2000 м3/сут. Полученные результаты свидетельствуют о существенном влиянии ландшафтных условий и водоотбора на все компоненты баланса, которое выразилось в изменении соотношения поверхностного и подземного стока в реку (рис.2), максимальных расходов воды и объемов речного стока, питания грунтовых вод и глубин их залегания. При восстановлении леса уменьшалась величина весеннего поверхностного стока и увеличивалось питание грунтовых вод. При сведении леса и распашке территории в речном стоке возрастала доля поверхностной составляющей при сокращении испарения и питания грунтовых вод. Эксплуатация водозабора приводила к уменьшению подземного стока в реку и сокращению общего объема речного стока.

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

Рис,2. Рассчитанные кривые обеспеченности среднемесячного модуля питания грунтовых вод (1) и среднемесячного модуля подземного стока в р.Медвенку (2) при изменении ландшафтных условий: I - ненарушенные условия, 2 — при облесении бассейна, 3 - при вырубке леса. 4 - при работе подземного водозабора

Глава IV. Моделирование водной миграции солей в ландшафтно-геохнмических системах с преобладанием химического стока

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

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

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

из почвы, либо к накоплению в профиле почв солей, содержащихся в поливной воде н в лессовой толще (Ковда, 1946-1947; Ковда, Егоров, 1968; Зимовец, 1974; Новикова, 1975; Моргун, 1976; Николаенко, 1991*, Куликов, Магнатаев, 1997).

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

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

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

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

Концентрация ионов С определяется с шагом 1 сутки по уравнению конвективной диффузии (Лукнер, Шестаков, 1986);

где О - коэффициент конвективной диффузии; V - скорость фильтрации; г -вертикальная координата; I - время. Затем находят объемы влаги, перетекающие из одного слоя в другой за декаду, и новые значения концентрации ионов в каждом слое. Последние уточняются в зависимости от процессов ионного обмена и растворения-осаждения. Кроме того, моделируется разложение почвенного органического вещества и выделение СОг, а также его миграция в почвенном воздухе. Интенсивность выделения СО2 линейно связана с содержанием органики и температурой почвы.

3(С0) _ д

а ~ &

d(Cv) &

Ионные равновесия в системе "раствор (Са2+, М§2+, Ма+, вО^", СГ, НСОз", СОз2") - газ (СОг) - твердая фаза (гипс, кальцит) - обменный комплекс" рассчитываются по модели, представляющей собой модификацию модели Я.А. Пачеп-ского (1976, 1990). Реакции растворения-осаждения гипса и кальцита, ионного обмена, образования и диссоциации угольной кислоты при растворении СО; моделируются на основе термодинамических уравнений относительно равновесных концентраций компонентов, содержащихся в различных фазах. Активности ионов в каждой фазе вычисляются в предположении соответствия катионного обмена между и линейной изотерме сорбции, для разновалентных ионов справедлив закон действующих масс. Равновесные концентрации СОг в растворе и в составе почвенного воздуха определяются на основе закона Генри,

Компьютерная программа при суточном расчетном шаге вычисляет средние месячные значения УГВ, а также все составляющие годового баланса: поступление воды с осадками, требуемое количество оросительной воды, поверхностный сток, изменение вяагозапасов в зоне аэрации, переток через водоупор. При включении дренажа определяются средние месячные и годовые значения модуля дренажного стока и одновременно вычисляются потоки влаги и объемы перетоков за декаду между слоями почвы и зоны аэрации. Раз в декаду рассчитываются изменения концентраций Саг+, 304г", СГ, а также количество гипса и кальцита. Такие же расчеты проводятся перед каждым поливом.

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

На территории НДОС в основном распространены черноземы южные слабозасоленные с сульфатно-гндрокарбонатным типом засоления. В ненарушенных условиях содержание солей в водной вытяжке изменялось от 3 мг-экв/100г на глубине 1 м до 10 на глубине 10 м при содержании кальцита на этих глубинах 1,7% и 9% от веса сухой породы. Параметры, характеризующие свойства почв и пород зоны аэрации, определялись по литературным данным (Зборищук и др., 1980; Баер, Лютаев, 1980), а также на основе материалов полевых изысканий. Калибровка модели проводилась по данным наблюдений на 1-ой очереди НДОС.

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

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

В ходе моделирования определена динамика всех составляющих водного баланса агроландшафгов в многолетнем разрезе, установлены темпы изменений качества поровых растворов и засоления (рис.3). Рассмотрено два варианта задачи. В первом варианте моделировались автономные ландшафты, где подъем УГВ вызывался только фильтрационными потерями при орошении. Во втором варианте дополнительно учитывали латеральную составляющую, обусловленную потерями,из распределительной сети. Результаты моделирования за 50-летний период показали, что в автономных ландшафтах УГВ поднимался медленно, со скоростью 0,3-0,5 м/год. При этом в почвенном слое на протяжении всего времени, включая период, когда грунтовые воды приближались к поверхности земли, сохранялся промывной режим. Эта особенность водного режима приводила к постепенному опреснению почвенной толщи, пока уровни находились на глубине более 3,5 м. В отдельные маловодные годы происходило временное накопление солей и повышение содержания ионов Na в ППК до 1,1 мг-экв/100 г.

мгокв

л

....^.í.r^r.'.^AAr.«^) а ...■•yjlVL*'»«MgJ -

,*■" •* г

с. \ ¡

1-1-1-1-1-1 ) ' ■ I Г--1

г ч ь & «з и <* к п го

Рис.3. Прогноз водного и солевого режима автономных (I) и подчиненных (II)

агроландшафтов при орошении: а) динамика содержания ионов в паровом растворе почвы (слой OSO см); 6') динамика УГВ и дренажного стока в подчиненном ландшафте (3)

МО-3020-<0-

В подчиненных ландшафтах происходил быстрый подъем уровня, вызванный боковым притоком, и постепенное засоление почвы. Содержание ионов Са1*, М§2+ и 5С>42" увеличивалось в 1,2-2 раза, Ыат и СГ - в 4,5 раза. Минерализация норового раствора достигала в конце расчетного срока 7 г/л при преобладании иона СГ, недопустимо возрастало содержание иона Ка+ в ППК - до 4 мг-эквЛООг, или 9% от емкости поглощения. В почвенной толще появлялся гипс (0,33%), увеличивалось содержание кальцита (до 2%). Отмечен]гые изменения проявлялись несмотря на наличие дренажа, который поддерживал УГВ в вегетационный период на глубине около 2 м. Расчеты для подчиненного ландшафта обнаружили крайне неблагоприятные условия и тем самым подтвердили необходимость применения наиболее современных методов полива и водораспределительной сети.

Таким образом, водный режим почв является наиболее важным фактором, определяющим направленность процессов засоления-рассоления в ландшафтах юга черноземной зоны. Создание оптимального водно-солевого режима ОС предполагает забор значительных объемов воды из водных источников. С помощью имитационной модели ОС получена интегральная оценка водного баланса рассматриваемой территории за многолетний период. Имитация функционирования НДОС при поэтапном вводе III и 1У-ОЙ очередей позволила построить кривые обеспеченности объемов водных дефицитов и недобора урожая. Из их анализа следует, что проектируемая ОС полностью удовлетворяет нормативным требованиям, при высоком к.п.д. водораспределительной сети прогноз водного и солевого режима орошаемых массивов благоприятен.

Глава V. Моделирование водной миграции солей в аквальных ландшафтах

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

Формирование качества воды в замкнутых водоемах рассмотрено на примере вдхр. лиманного типа. Лиман Сасык площадью 210 км2 и глубиной 2-3 м расположен на побережье Черного моря. Как вдхр. сезонного регулирования Сасык эксплуатируется с 1979 г., когда морской залив перегородили дамбой и соединили с Дунаем. По каналу Дунай-Сасык в вдхр. поступает ежегодно 400800 млн.м3 пресной воды, которая, смешиваясь с водами лимана и опресняя их, расходуется на орошение и промывки и частично сбрасывается в море. С севера в вдхр. впадают реки Когильник и Сарата, с минерализацией воды от 3 до 6 г/л. Под вдхр. и вдоль его берегов в подстилающих породах мощностью более 50 м

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

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

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

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

т ¿5 & & а

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

= (2)

где к - коэффициент фильтрации донных отложений, g - ускорение силы тяжести, р - давление. Система уравнений (1) и (2) решалась итерационным методом Гаусса-Зейделя.

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

Объем и качество воды в водохранилище характеризовались уравнениями водного и солевого балансов северного и южного плесов за декадные интервалы времени. Качество воды из плесов зависит от привноса основных ионов (Саг+, М^1*, 3042", СГ, НСО/) поверхностными и подземными водами, водозабора, испарения и водообмена между плесами (рис.4). При сложной гидрохимической обстановке рассмотрен баланс каждого иона в отдельности. Входные функции в

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

Рис.4, Составляющие водного и солевого баланса лимана Сасык: ш----- регулируемые; -- нерегулируемые

Прогноз качества воды в вдхр. Сасык. С помощью модели проанализированы различные проектные варианты подачи и отъема воды из водоема на орошение. Прогноз солевого режима выполнен для двух основных вариантов. Первый вариант соответствовал реальным условиям подачи воды по каналу (720 млн. м /год), во втором варианте рассматривалась возможность увеличения водоподачн вдвое.

Водный и солевой режимы водохранилища моделировались с момента начала его эксплуатации. Первый вариант использовался для калибровки модели. В приходных и расходных статьях баланса учитывались осадки, речной сток и испарение. Согласно полученным результатам (рис.5), в первые годы эксплуатации водохранилища происходило опреснение воды, затем режим стал квазистационарным, при котором колебания минерализации и содержание Na+, SAR имели сезонный характер и зависели от поступления воды из источников.

г. мес.

Рис.5. Рассчитанные по модели колебания уровня воды (1), содержания иона Ка+ (2) и БАЯ (3) в южном плесе в начальный период эксплуатации водохранилища (первый вариант)

В первом варианте средний за 15-летниЙ период вынос солей составил 262 тыс. т/год при разгрузке 1,2 млн. м3/год, во втором - 216 тыс. т/год при разгрузке 0,83 млн. м3/год. Доля донного привноса солей в солевом балансе вдхр, равнялась соответственно 30% и 22%. Рассчитанная величина минерализации воды южного плеса была близка к 0,8 г/л при Иа/Са = 2,3-2,5 в первом варианте и к 0,6-0,7 г/л при Ыа/Са =1,4 во втором. Таким образом, в обоих вариантах значения Ыа/Са не удовлетворяли требованиям, предъявляемым к оросительной воде. Уменьшение этого соотношения возможно только при существенном увеличении водообмена. Последнее вряд ли реально, поскольку для этого требуется практически шестикратное увеличение водоподачи по каналу - 4000-5000 млн. м3/год.

Результаты моделирования показали, что качество воды в вдхр, Сасык в значительной степени определяется донной разгрузкой минерализованных вод. Запасы растворенных солей, накопившихся под дном водохранилища, таковы, что при современной интенсивности разгрузки для полного замещения рассолов водой с меньшей минерализацией потребуются столетия. Поэтому использование лимана Сасык в качестве вдхр, нецелесообразно, орошение водой из него приведет к деградации почв {Куликов, Магнатаев, 1997).

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

Глава VI. Моделирование трансформации н миграции тяжелых металлов в природных и агроландшафтах

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

Моделирование трансформации тяжелых металлов в почвах. Модели процессов трансформации "1М в почвах относятся либо к кинетическим, либо к термодинамическим моделям, по детальности описания процессов их можно разделить на механистические и эмпирические. В механистических моделях описывают равновесное распределение ТМ между твердой фазой, почвенным раствором, органическим веществом и неорганическими комплексами, опираясь на физико-химические свойства почвы и данные химических анализов различных форм ТМ й а1., 1996; Бе Упеэ е1 аЦ 1997). Так как формы ТМ,

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

Эмпирический подход состоит в построении и исследовании обобщенных

статистических зависимостей между формами нахождения ТМ в почвенном профиле и свойствами и режимами почв различных геохимических провинций (van der Zee et al., 1988; Third Codata Conference, 1993; Мотузова, 1999). Этот подход, по нашему мнению, позволяет преодолеть принципиальные трудности информационного обеспечения механистических моделей при условии репрезентативности экспериментальных данных.

Связи между содержанием подвижных форм ТМ и комплексом ландшафт-но-геохшшческих факторов и условий описываются методами многомерного анализа данных, С помощью регрессионных моделей можно с достаточной точностью предсказать содержание подвижных форм ТМ в том диапазоне значений факторов, который отмечен в натурных наблюдениях (Мотузова, 1988, 1999; Puis, Powell, 1991; Пачепский, 1992; Малкина-Пых, 1998).

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

Ниже приводятся результаты исследования подвижности ТМ в почвах смешанных лесов, для которых установлена фоновая геохимическая структура ТМ, а также выявлен комплекс факторов, определяющих поведение ТМ в различных ландшафтных условиях (Васильевская, Якущевская, 1972; Тойкка, 1973; Каури-чев и др., 1977; Микроэлементы в почвах СССР, 1981; Ландшафтно-геохимиче-ские основы, 1989; Мотузова, 1992; Обухов, Попова, 1992; Цаплина, 1994 и др.).

Подвижные формы ТМ характеризовались их концентрациями в вытяжке 1н HCI вследствие простоты методики их определения, высокой информативности и возможности сопоставления с результатами других исследований (Цюрупа, 1961; Соловьев, 1989). Валовые содержания металлов определялись количественным спектральным методом в лаборатории ГГЭ г. Бронницы, подвижные формы - методом атомной абсорбции на спектрофотометре фирмы "Хитачи 170/80" в лаборатории географического факультета МГУ (аналитики Е.А. Шахпекдерян, Т.Г. Сухова). Количество гумуса (по Тюрину) и рН водный измерялись общепринятыми методами.

Для моделирования использовались данные по Смоленско-Московской возвышенности с пределах Сатинского полигона МГУ. ГИС содержит координаты и описания 137 почвенных разрезов, а также результаты химических анализов: содержание гумуса, рН, концентрации валовых и подвижных форм Мп, 2n, Си, Pb, Со, Ni, Сг и подвижных форм Fe в генетических горизонтах почв, распространенных на изучаемой территории. К ним относятся дерново-подзолистые почвы на покровных, делювиальных и моренных суглинках и флювиогляциальных отложениях, дерновые на делювиальных суглинках и балочном аллювии, пойменные дерновые на отложениях различного гранулометрического состава и болотные.

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

Для большинства металлов построенные регрессионные модели имеют удовлетворительную точность. Их валидация проведена на основе опубликованных ранее результатов экспериментальных исследований механизмов мобилизации ТМ в почвах.

Интерпретация моделей позволила выявить провинциальные особенности поведения ТМ в почвах и оценить роль факторов, обусловливающих дифференциацию подвижных форм ТМ (рис,б). Установлено, что подвижность Мп и Си в изученных почвах определяют в основном их органо-минерапьные соединения. Связь подвижных форм РЬ с содержанием гумуса проявляется только в восстановительной среде. Соединения Сг с органическим веществом легкоподвижны и поэтому выносятся из почвенного профиля. Зависимость между содержанием гумуса и подвижных форм N1 нелинейна, у Со и 2п ее не обнаружено.

--Мп = 200;----Мп = Ю00мг/кг почт

Рис.6. Концентрация подвижных форм (Сп) марганца (А) и хрома (Б) в зависимости от их валового содержания (Св) и основных почвенных свойств: А - содержание гумуса, %, рИ - реакция почвенного раствора, Ре и Мп -содержание подвижных железа и марганца, мг/кг почвы, £=/ - наличие оглеения, $=0 — его отсутствие

Щелочные условия способствуют уменьшению содержания подвижных форм 2п, РЬ, N1, Со и увеличению - Мп, У большинства ТМ (Мп, Си, РЬ, N1) в восстановительной пгеевой среде концентрация подвижных форм снижается, у 2п - несколько увеличивается. Ассоциированность Си, РЬ, Сг, № с

соединениями Fe, как правило, ведет к накоплению подвижных форм этих элементов с ростом содержания подвижного Fe. Значимый эффект на поведение ТМ оказывают подвижные соединения Мп, с увеличением концентрации которых растет количество подвижных форм Со, Zn( Pb, Ni и уменьшается - Cu и Сг. Подвижные формы соединений Мл, Сг, Ni имеют прямую положительную связь с их валовыми концентрациями. Для Си и Со характерна обратная зависимость между этими показателями, а у РЬ и Zn она отсутствует.

Улучшение точности регрессионных моделей, вероятно, может дать учет пространственного положения точек опробования с характеристикой рельефа через локальные уклоны и кривизны дневной поверхности (Pachepskyet al., 2001). Такое исследование может обеспечить точная цифровая модель рельефа, в частности, Сатинекого полигона. Другой способ учета пространственного положения - определение геохимических ландшафтов через качественную переменную с последующим применением алгоритма регрессионных деревьев (Ralws, Pachepsky, 2002).

Построенные модели могут быть эффективно использованы I) при планировании более широких по территориальному охвату зональных исследований; 2) при проведении эколого-геохимической экспертизы устойчивости почв; 3) При картографировании условий миграции ТМ в ЛГС.

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

Интенсивность миграции и аккумуляции ТМ в агроландшафтах целесообразно оценивать с помощью комплексной модели биогеоценотического круговорота ТМ. Уравнение баланса ТМ в кор необитаем ом слое почвы за многолетний период записано как:

dS = Q-I-P di м-р

где р - плотность почвы (кг/м3); т ~ мощность пахотного слоя (м). Переменная состояния S(t) характеризует среднее содержание (накопление) металла (мг/кг) в пахотном слое в í-ый интервал времени. Его динамика определяется несколькими процессами, поэтому в правой части уравнения задаются Q ~ поступление металла с атмосферными осадками, пылью и удобрениями (мг/мг'-год); / -выщелачивание металла из пахотного слоя; Р - вынос с растительной продукцией (рис.7). Реализация модели предполагает непосредственные измерения S(0) всех интересующих металлов в почвах и оценку правой части уравнения (3) путем моделирования составляющих баланса.

Источниками ТМ в агроландшафтах являются атмосферные выпадения техногенного происхождения и поступления с удобрениями и мелиорантами. Интенсивность Q(t) привноса ТМ определяется суммой

І-1 -1 І-И -2 С3> -3 -► -4

Рис,7. Структура балансовой модели ТМ в агроландшафте: І - переменные состояния; 2 - процессы; 3 - внешние факторы; 4 - потоки ТМ

Q(i) = h ■ C} + (Ts • C2 + 2 Z>i ■ С Д > • ;o-'vA

где Ti и A - объем твердых выпадений (кг/га-год) и слой атмосферных осадков (м/год); С/ и С: - содержание в них ТМ (мг/кг), D, и С/, - доза (кг/га-год) /-го вносимого удобрения или мелиоранта и концентрация в нем ТМ (г/м3), к -количество видов удобрений и мелиорантов в используемой технологии возделывания сельскохозяйственных культур.

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

I(t) = (h-evp)Cp.

1(0 уточняется при калибровке всей балансовой модели. Эвапотранспирация evp определяется по модифицированной формуле Девита (DeWit, 1958).

Вынос металлов с растительной продукцией P(t) (мг/га) рассчитывается в зависимости от урожайности возделываемой культуры 7(кг/га) и концентрации R в ней элемента (мг/кг сухого вещества): P(t) = Y'R. R является функцией поглощения R » /fCV, индивидуального для каждой сел ьс кохо~ йственной культуры, В качестве аргумента С„ (мг/кг) используется концентрация ТМ в ацетатно-аммонийной вытяжке (рН=4,8) из почвы, которая по своим свойствам максимально приближается к корневым выделениям растений. Для ее определения использованы зависимости, аналогичные рассмотренным выше. Величина R оценивается с помощью банка прикладных моделей, обобщающих результаты полевых, михрополевых и вегетационных опытов с разными ' сельскохозяйственными культурами (Кошелева и др., 1995; Черных, 1995). j

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

Из-за пространственных изменений в интенсивности потоков ТМ, формирующих результирующий баланс ТМ, численное моделирование различных сценариев антропогенного воздействия проводится по отдельным участкам исследуемой области. Эта задача решалась нами на примере Московской области, региональный фон ТМ которой формируется под воздействием атмотехногенного привноса. Средства химизации сельского хозяйства усиливают загрязнение почв, поскольку обогащены Hg, Ag, Cr, Zn, Cu, Ni, Co, V, Pb. Результирующие нагрузки на агроландшафты области отражает карта-схема (Kosheleva, Garanushkin, 1999).

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

области, на междуречье p.p. Москвы и Овей, где ведутся полевые и микрополевые опыты с различными сельскохозяйственными культурами на дерново-подзолистых почвах суглинистого состава. Район испытывает значительную техногенную нагрузку, связанную с промышленностью а транспортом. К приоритетным загрязнителям относятся Pb, Cd, а также Мл, Zn, Си.

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

С помощью модели было определено время, за которое уровень содержания ТМ достигнет ПДК. Установлено, что с максимальной скоростью аккумулируется Zn и, если в ближайшие десятилетия не произойдет заметного снижения темпов атмотехногеиного поступления ТМ в почву, то уже через 30-70 лет содержание Zn достигнет ПДК, Удвоение концентрации Zn произойдет через 75-190 лет. РЬ, Си н Ni накапливаются не столь интенсивно и поэтому существующую нагрузку на агроландшафты можно считать допустимой.

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

Защищаемые положения в выводы

На защиту выносится совокупность теоретических и практических результатов, обеспечивающих количественное описание и прогноз процессов миграции и трансформации вещества в ЛГС.

1. Разработаны теоретико-методологические основы математического моделирования ландшафтно-геохимических процессов в элементарных и каскадных ЛГС локального уровня.

=5 Проведена классификация математических моделей миграции и трансформации вещества в ЛГС. В зависимости от природы ландшафтно-геохимических процессов их математические модели объединены в четыре группы: модели отдельных фаз водного цикла; транспортные модели; термодинамические и кинетические модели; модели биологической миграции. Соотношение между процессами зависит от структуры и масштабов ЛГС, а также от свойств контролируемых элементов и их соединений. Пространственная организация геохимических ландшафтов представлена ЭЛГС, катенами и

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

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

=> Сформулированы принципы построения комплексных моделей биогеохимических циклов: (1) модульный, предполагающий разработку и последующее "сшивание" моделей, описывающих процессы разной природы; (2) соответствия характерных размеров элементов ЛГС характерным временам миграционных процессов; (3) агрегирования пространственных переменных в зависимости от характера их неоднородности.

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

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

- водный цикл в каскадной системе 1-го порядка лесной зоны;

- водиый и солевой режим оросительной системы в степной зоне;

- формирование качества воды в водохранилище лиманного типа;

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

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

=> Верификация комплексной модели свидетельствует о том, тго она корректно отражает закономерности формирования и движения радиальных и латеральных потоков поверхностных и подземных вод в лесной зоне. Анализ

миграционной структуры КЛГС 1-го порядка в бассейне р. Москвы показал

- влияние погодных условий на разгрузку подземных вод в пойме реки;

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

- увеличение питания горизонта подземных вод в лесных ландшафтах по сравнению с полевыми.

=> Предложенная модель учитывает ландшафтную структуру каскадной системы, поэтому численные эксперименты с моделью позволили оценить различные природные ситуации и хозяйственные мероприятия, еще не наблюдавшиеся или не реализованные в данном бассейне. Получен прогноз возможных изменений миграционной структуры КГЛС в бассейне р. Москвы под воздействием антропогенных факторов, согласно которому;

- при восстановлении лесной растительности на пашне произойдет уменьшение на 40% величины весеннего поверхностного стока и увеличение на 20% питания грунтовых вод;

- при сведении леса и распашке территории увеличится на 25-35% доля поверхностной составляющей в речном стоке и уменьшится на 5-10% испарение и питание грунтовых вод;

- эксплуатация водозабора подземных вод приведет к сокращению на 20% общего объема речного стока и на 40-50% подземного стока в реку.

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

=> Блок вертикального влагопереноса в ЭЛГС с удовлетворительной точностью описывает все пять типов геохимических сопряжений с преобладанием химического стока, так как допускает произвольные соотношения между скоростями поступления и инфильтрации воды в почву при любом УГВ, С его помощью рассчитан режим увлажнения почв, дренажного стока и грунтовых вод на оросительных системах юга черноземной зоны в начальный и эксплуатационный периоды.

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

=> Апробация модели на Нижне-Днестровской оросительной системе выявила качественные и количественные различия в солевом режиме автономных и подчиненных ландшафтов: постепенное рассоление почвенной толщи первых, при котором В ППК снижалось содержание Иа и возрастало содержание М§, и постепенное засоление почвы вторых под влиянием инфиль-трационных потерь, которое сопровождалось увеличением содержания Ыа и С1вППК.

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

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

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

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

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

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

Список работ, опубликованных по теме диссертации

1. Кошелева Н.Е., Хранович ИЛ. Модель выбора оптимальных параметров ВХС, обеспечивающих жесткий график водопотребления// Водные ресурсы, 1983, №4, с.61-74.

2. Кошелева Н.Е., Михайлов В.Н., Назаров H.A. Ожидаемые изменения стока р.Терек под влиянием роста безвозвратного водопотребления-'/ Вестник Моск. ун-та, Сер. География, 1984, N3, с.37-45.

3. Кошелева Н.Е. Оценка надежности и корректировка гидрологических данных при обосновании крупных водохозяйственных мероприятий// Географический анализ при планировании народнохоз. программ. Мл юд-во МГУ, 1986, с.87-93.

4. Кошелева Н.Е., Пашковский И.С. Определение иифильтрационного питания по данным режимных наблюдений с использованием функций влияния// Рац. использование водных ресурсов, вып.7, М.: Наука, 1986, с.46-51.

5. Кошелева Н.Е,, Пашковский И.С., Прилепина Т.Ф. Геофильтрациопная модель прогноза подземного и поверхностного стока речного бассейна// Охрана и гидрогеол. нсслед. подзем, вод вЦешр. районах РСФСР, М., 1986, с.98-107.

6. Кошелева Н.Е., Пашковский U.C. Численное моделирование влагопереноса в связи с прогнозом водного режима зоны аэрации при орошении// Рациональное использование водных ресурсов, вып.8, М,: Наука, 1987, с.56-63.

7. Кошелева Н.Е., Пашковский U.C., Прилепина Т.Ф. Особенности водного режима зоны аэрации при сезонном промерзании 1рунтов// Оценка и прогноз изменения гидрогеологических условий на застраиваемых территориях. М., 1987, с.69-78.

8. Кошелева Н,Е„ Пашковский И.С. Оценка потребности воды на орошение по расчетам влагопереноса в зоне аэрации// Состояние и перспективы использования подземных вод для орошения. М.: Наука, 1988, с.44-48.

9. Кошелева U.E., Пашковский И.С., Прилепина Т.Ф. Динамико-стохастическая модель формирования поверхностного и подземного стока на водосборе'/ Результаты исследований в обл. компл. инж. изысканий на орош. землях. Сб. науч. тр. В/О "Союзводпроекг". М., 1989, с.9-24.

10. Arutjunova L.V., Kirikoj J.T., Koscheleva N.E. Informationssystem des agrookologichen Monitorings auf der Grundlage des Geographischen Netzes von Feldversuchen. Feldversuchswesen, Koordiiüerungsstelle für Feldversuchswesen Bad Lauchstadt, nl5. Berlin 8(1991) l,pp.3-20.

11.Кошелева H.E., Арутюнова Л.В, Система математических моделей для экологической экспертизы технологий возделывания сельскохозяйственных культур (по данным ахроэкологического мониторинга)// Материалы международного симпозиума "Инженерная экология-91", Звенигород, 1991. М., 1991, с, 114-117,

12. Арутюнова Л,В., Кошелева Н.Е. О разработке экологического блока информационно-вычислительной системы агроэкологического мониторинга// "Автоматизация научных исследований в области растениеводства. Проблемы фитстронии", Научно-тех. бюлл. ВИР, 1991, N215, С.5-15.

13.Кошелева Н.Е., Пашковский И.С. Оценка водообеспеченности и прогноз засоления на оросительной системе по расчетам водного и солевого режимов зоны аэрации// Оценка и прогноз изменений инж.-геол. условий осваиваемых терр. М.: ЦБНТИ "Водстрой", 1991, вып.1, 4.1, с.51-82.

14. Исследования в полевых и производственных опытах по изучению эффективности систем питания растений в комплексе с др. средствами химизации. Методика и

программа// Ладоннн В.Ф., Цимбалист H.H., Самойлов Л.Н. и др. М.: ВИУА, 1991, ч.1.186с.

15.Научные основы повышения эффективности комплексного применения средств химизации (удобрений, пестицидов, регуляторов роста) в интенсивных технологиях возделывания зерновых культур, обеспечивающих их высолю окупаемость и экологическую безопасность// Милащенхо Н.З., Ладоннн В.Ф., Цимбалист Н.И, и др. М.:ВИУА, ¡991,4.2. 143с.

16. Кошелева U.E. Оценка загрязнения окружающей среды средствами химизации// Химизация сельского хозяйства, 1991, N 11,с.91-97,

17.Кошелева И.Е., Комиссарова 14.Ф., Орешина ЕЛ., Арутюнова Л,В. Методические основы экологической паспортизации сельхозпредприятий с использованием базы данных// Вестник Российской академии сельхознаук, 1992, N6, с.56-62.

18.Кошелеэа Н.Е., Яковлев В.П. Расчет спектра импульсного случайного процесса с независимыми шпервалами// Радиотехника и электроника, 1993, вып.9, с. 1597-1604.

19. Кошелева Н.Е, Научно-производственная компьютерная система для опенки загрязнения почв и сельскохозяйственных культур тяжелыми металлами и обоснования способов снижения их токсичности// Тяжелые металлы н радионуклиды в агроэкосисгемах. M.: Агроэколас, 1994, с.270-277,

20.Кошелева Н.Е., Комиссарова И.Ф., Черных НА., Черных И.Н„ Большакова Л.С. Прикладная модель накопления тяжелых металлов в некоторых культурах на дерново-подзолистых почвах// Тяжелые металлы и радионуклиды в агроэкосистсм&х. 1994, М.; Агроэколас, с.237-244.

21.Касимов U.C., Кошелева Н.Е., Самонова O.A. Подвижные формы тяжелых металлов в почвах лесостепи Среднего Поволжья (опыт многофакториого регрессионного анализа)// Почвоведение, 1995, N6, с.705-713.

22. Кошелева Н.Е,, Комиссарова И,Ф„ Лебедев С.Н. Компьютерная система оценки почв, загрязненных тяжелыми металлами, для выращивания сельскохозяйственных культур// Агрохимия, 1995, N4, с.65-73.

23.Кошелева Н.Е., Пашковский И.С, Моделирование формирования качества воды в водохранилище Сасык//Водные ресурсы, 1996, том 23, Nl,c.66-72.

24. Kasimov N.S., Kosh eleva N,E„ Samonova O.A. Mobile forms of heavy metals in soils of Middle Volga Forest-Steppe (Experience of multivariate regressional analysis}// Eurasian Soil Sei. 28/9 (Sept. 1996), p.47-61.

25. Кошелева Н.Е. Моделирование почвенных и ландшафтно-геохимических процессов. Учебное пособие. М.: МГУ, 1997,109 с.

26.Самонова OA., Касимов Н.С., Кошелева U.E. Ассоциации тяжелых металлов в профиле дерново-подзоластых почв южной тайги// Вестник Моск. ун-та. Серия 17. Почвоведение, 1998, N2, с.14-19,

27.Kasimov N.S., Kosheleva N.E., Samonova O.A. Mobile forms of heavy metals in soddy-podzolic soits of southern taiga. Proceedings of 16-th World Congress of Soil Science, Montpellier, France, 1998. Symp.6,6 pp.

28. Самонова O.A., Касимов H.C., Кошелева Н.Е, Об оценке результатов мониторинга тяжелых металлов в дерново-подзолисгых почвах южно-таежных ландшафтов// География на рубеже веков; проблемы регионального развития. Материалы междуиар. науч. конф. 22-25 сент, 1999, Курск, т.З, с.219-225.

29. Кошелева Н.Е. Моделирование миграции тяжелых металлов в агроландшафтах// География на рубеже веков: проблемы регионального развития. Материалы

междунар. науч. конф, 22-25 сект. 1999, Курск, т. 2, с. 88-93.

30.Samonova O.A., Kosheleva N.E., Kasimov N.S. Variations of heavy metal concentrations in soddy-podzolic soils under forest and meadow plant communities// Polish journal of soil science, Ossolineum, Wroclaw, 2000, vol.33, no,2,12 pp.

31.Кошелева H.E., Касимов H.C., Самонова O.A. Миграционная способность тяжелых металлов в почвах смешанных лесов// География и окружающая среда, М.: ГЕОС, 2000, с.415-428.

32.Самонова O.A., Касимов Н.С., Кошелева Н.Е, Пространственно-временное варьирование содержаний тяжелых металлов в дерново-подзолистых почвах южной тайги// Всстник Моск. ун-та, сер.17, почвоведение. 2000, №2, с.20-26.

33.Кошелева Н.Е. Моделирование биогеохимических циклов тяжелых металлов в агроландшафтах на основе балансового подхода// Геохимия ландшафтов и география почв. Смоленск: Ойкумена, 2002, С.389-405.

34. Kosheleva N.E., Kasimov N.S., Samonova O.A. Methods of multivariate statistical analysis in landscape and geochemical research// Proceedings of 17-th World Congress of Soil Science, Bangkok, Thailand, 2002, Symp.2,10 pp.

35.Кошелева H.E., Касимов H.C., Самоиова O.A. Регрессионные модели поведения тяжелых металлов в почвах Смоленско-Московской возвышенности// Почвоведение, 2002, №8, с.954-966.

36. Кошелева Н.Е., Яковлев В.П. Нелинейное моделирование в экологии// Успехи современной радиоэлектроники, 2002, №11, с.35-43.

37. Состояние тсрригориально-аквальных комплексов малых и средних рек Московского региона// Алексеевский Н.И., Калинкин И.В., Кошелева Н.Е. и др. Государственный регистр баз данных, ФГУП Научно-технич. центр «Информрсгистр». Рсг.свнд. Лг7782, дата per. 23.04.2002, perJfc 0220208320.

38. Гаранюшкин M.H., Кошелева Н.Е. Моделирование распределения микроэлементов в почвах и донных отложениях западного побережья Каспийского моря// Геоэкология Прикаспия,М.:изд-воМГУ, 2003. Вып.6,15 с. (в печати).

Депонированные рукописи и тезисы конференций

39.Кошелева Н.Е. О применении системного анализа для решения проблем рационального использования и охраны водных ресурсов// Материалы VII конф.асп. и мол. ученых. Секция "Охрана геол. срсды", Деп.ВИНИТИ N5521 or 26.12.1980, с.21-30.

40. Кошелева Н.Е. Метода моделирования в проектировании водохозяйственных систем// Материалы VIH конф.асп. и мол. ученых. Секция "Охрана геол. среды", Деп.ВИНИТИ N 5854-81ДЕП, с.34-50,

41.Кошелева Н.Е., Фролова H.JI. Возможности применения системного подхода при анализе и моделировании гидрологических систем// Материалы IX конф. асп, и мол, ученых. Серия "Охрана геол. среды", М„ 15.3.1982. Деп ВИНИТИ N 782- 83ДЕП от 14.2.1983, с.33-49.

42.Кошелева Н.Е., Назаров IIA. Об условиях применения уравнения Сен-Венана без инерционных членов и кинематической волны для описания водных потоков// Материалы конф. асп. и мол, ученых. Серия "Охрана геол. среды и водных ресурсов", М„ Деп. ВИНИТИ Ы6482-83Деп. от 2.12.1983, с. 188-199.

43. Кошелева U.E., Миронова Е.М, Асимптотический метод решения уравнений движения водного потока// Материалы XII конф.асп. и мол. ученых геол. ф-та МГУ. Серия "Охрана геол.среды и водных ресурсов". Деп.ВИНИТИ N3715 от 22.5.1986. С.74-84.

44.Бурлин MJO., Кошелева Н.Е., Пашковский И.О. Прогнозное моделирование влажностного и солевого режима зоны аэрации массива орошения на персональном компьютере. Тез, докл. Всес. семинара "Применение ЭВМ при гидрогеохимическом моделировании". Л.: изд-во ЛГУ, 1991, с.27-29.

45.Кошелева Н.Е., Комиссарова И.Ф., Черных Н.А., Черных И,Н. Экологическая экспертиза сельскохозяйственных территорий, загрязиетшых тяжелыми металлами. Матер. Междунар. симпозиума "Проблемы экоинформатики". Звенигород, 1992, с.138-139.

46.Кошелева Н.Е., Комиссарова И.Ф„ Шеин Е.В., Губер А.К. Математическое моделирование миграции и транслокации в растения тяжелых металлов для прогноза загрязнения почв// Материалы Междунар, симпозиума "Проблемы экоинформатики". Звенигород, 1992, с.140-141.

47. Вагнер В.Б., Кошелева Н.Е., Герасимова М.И., Исаченкова Л.Б. Разработка локальной почвенно-геохимической геоинформационной системы многоцелевого назначения. Тез. кон ф. "Геохимия биосферы", Новороссийск, 1994, с.68-69.

48. Кошелева Н.Е. Исследование радиальной геохимической дифференциации почв методами многомерного статистического анализа. Тез. конф. "Геохимия биосферы", Новороссийск, 1994, с.72-73.

49. Kosheleva N.E., Komissarova I.F., Chcmikh N.A, Heavy metal balance in arable soils: field experiments and computer modelling. Annates Geophysicae, part II, Supplement II, vol.13,1995. C298.

50. Kosheleva N,E,, Komissarova i.F., Chernikh N.A., Lebedev S.N. Response of plant species to concentration levels of Zn,Pb and Cd as influenced by soil conditions. 3nd Internal, conf, on Biogeochcmistry of Trace Elements, Paris, 15-19 May 1995.

51. Kosheleva N.E., Komissarova 1.Г., Chernikh N.A. The assessment of long-term effects of heavy metal loads on soils and crops: balance approach, Proc. of the 2d Inlernat. IAWQ Spec. Conf.and Symp. on Diffuse Pollution, Bmo&Prague, 13-18 aug, 1995, part II, 593594.

52. Кошелева Н.Е. Об использовании количественных показателей для оценки радиальной дифференциации почв// "Теория и техника передачи, приема и обработки информации". Тез. докл. 2-оЙ междунар. конф. Харьков-Туапсе, ХТУРЭ, 1996, с.121-122.

53.Кошелева Н.Е. Прогнозная модель миграции тяжелых металлов в агролапдшафтах// "Тяжелые металлы в окружающей среде". Сб. тез., Пушино, 1996, с. 146-147, 272.

54. Kosheleva N.E., Garanushkin M.N. Heavy metal contamination of agricultural landscapes in Moscow region. Proc .of ext. abstracts of 5-th Intemat. Conference on the Bicgeochemistry of Trace Elements, July 11-15,1999. Vienna, Austria, Vol.11, p.968-969.

55.Kasimov N.S., Samonova O.A., Kosheleva N.E. Behaviour of trace elements in soils of Southern Taiga Ecosystems// 29th International geographical congress. 14-18 august 2000. Seoul, Korea. Abstracts, p.227-228.

56.Кошелева Н.Е. Применение метода регрессионных деревьев при обработке результатов ландшафтно-геохимических исследований// Тезисы Междунар. школы "Современные методы эколого-гесхимической оценки состояния и изменений окружающей среды", Новороссийск, НИИ геохимии биосферы, 2003, с.67-68.

Усл. ист, д. 2,1

Зак. 482

Тираж 100 экз.

АНО «Издательство МСХА» 127550, Москва, ул. Тимирязевская, 44