WWW.KN.LIB-I.RU
БЕСПЛАТНАЯ  ИНТЕРНЕТ  БИБЛИОТЕКА - Различные ресурсы
 


Pages:     | 1 |   ...   | 3 | 4 || 6 | 7 |

«УЧРЕЖДЕНИЕ РОССИЙСКОЙ АКАДЕМИИ НАУК УЧРЕЖДЕНИЕ РОССИЙСКОЙ АКАДЕМИИ ИНСТИТУТ ФИЗИКО-ХИМИЧЕСКИХ И НАУК ИНСТИТУТ МАТЕМАТИЧЕСКИХ БИОЛОГИЧЕСКИХ ПРОБЛЕМ ПОЧВОВЕДЕНИЯ РАН ПРОБЛЕМ БИОЛОГИИ РАН ...»

-- [ Страница 5 ] --
Ольчев А.В. 1, Новенко Е.Ю.2 Учреждение Российской Академии Наук Институт Проблем Экологии и Эволюции имени А. Н. Северцова РАН, Москва, Россия aoltche@gmail.com Учреждение Российской Академии Наук Институт Географии РАН, Москва, Россия lenanov@mail.ru Аннотация: уравнения для расчета годовых сумм фактического и потенциального испарения, а также транспирации бореальных лесов по значениям среднегодовой температуры и суммарного количества осадков были получены по результатам численных экспериментов с моделью Н2О и СО2 обмена Mixfor-SVAT. Полученные уравнения были использованы для расчета испарения еловых лесов Центрально-Европейской территории России в Голоцене.

Введение Основной целью данного исследования являлась разработка упрощенных выражений для расчета годовых значений фактического (Е) и потенциального (PE) испарения, а также транспирации (TR) еловых лесов Европейской территории России в голоцене по палеоклиматическим данным. Параметризации были получены на основе анализа значений TR, E и PE, рассчитанных для модельных древостоев с разным участием ели (от 0 до 100%) при разных температурных условиях и режиме увлажнения с помощью Mixfor-SVAT модели. PE рассчитывалось с помощью уравнения Пристли-Тейлора.

Модель Mixfor-SVAT Одномерная процесс-ориентированная математическая модель MixFor-SVAT была разработана для расчета переноса радиации, потоков турбулентного тепла, Н2О и СО2 в разновозрастных смешанных древостоях (Ольчев и др., 2008, Olchev et al. 2008, 2009).

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

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

В общем виде для Е уравнение может быть записано как:

E = a 0 ( a1 6 T + a 2 ) ( a3 N + a 4 ) ( a5 700 p + a 6 ) (1) где T - среднегодовая температура, p - годовое количество осадков в мм год-1, N - доля ели в древостое (0 - 1), a0-a6 - эмпирические коэффициенты. Значения коэффициентов a0-a6 в уравнении 1 были получены методом Левенберга-Марквардта на основе анализа результатов расчетов E, PE и TR по модели MixFor-SVAT (таблица 1).

Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия

–  –  –

Литература Ольчев А.В., Курбатова Ю.А., Варлагин А.В., Выгодская Н.Н. Модельный подход для описания переноса СО2 между лесными экосистемами и атмосферой // Лесоведение.- 2008.- № 3.- С. 3-13.

Olchev A, Ibrom A., Ross T., Falk U., Rakkibu G., Radler K., Grote S., Kreilein H., Gravenhorst G. A modelling approach for simulation of water and carbon dioxide exchange between multi-species tropical rain forest and the atmosphere. // J. Ecological Modelling, 2008.– 212.- P. 122–130.





Olchev A, Novenko E, Desherevskaya O, Krasnorutskaya K, Kurbatova J Effects of climatic changes on carbon dioxide and water vapor fluxes in boreal forest ecosystems of European part of Russia // Environ. Res. Lett.-. 4.- 2009.pp) Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия

МОДЕЛЬ ГЕНЕТИКО-ПОПУЛЯЦИОННОЙ СТРУКТУРЫ

ПАРТЕНОГЕНЕТИЧЕСКОГО ВИДА КАВКАЗСКОЙ СКАЛЬНОЙ ЯЩЕРИЦЫ

DAREVSKIA UNISEXUALIS НА ОСНОВЕ МНОГОМЕРНОГО АНАЛИЗА

Омельченко А.В. 1,2, Вергун А.А. 2,3, Петросян В.Г. 1, Корчагин В.И. 3, Токарская О.Н. 3 1УРАН Институт проблем экологии и эволюции им. А.Н. Северцова РАН, Москва, Россия omi@bk.ru 2ГОУ ВПО Московский педагогический государственный университет, Москва, Россия vermand@mail.ru 3УРАН Институт биологии гена РАН, Москва, Россия tokarskaya@list.ru Аннотация: в настоящей работе представлена математическая модель генетико-популяционной структуры партеновидов, созданная на основе данных о межпопуляционном полиморфизме микросателлитных маркеров генома партеногенетической ящерицы Darevskia unisexualis с использованием математического аппарата многомерного анализа, позволяющая количественно оценить генетический вклад популяций в генофонд вида.

Введение Облигатно-партеногенетические рептилии и, в частности, партеногенетические виды Кавказских скальных ящериц рода Darevskia, система размножения которых исключает какое-либо взаимодействие с мужским геномом, популяции которых представлены исключительно самками, а размножение осуществляется по принципу естественного биологического клонирования, являются уникальными природными объектами для анализа геномов, мониторинга генетической изменчивости и связанной с ней динамики генофондов (Рысков, 1999). Однополые виды животных, и в том числе, виды рода Darevskia, значимо отличаются от двуполых отсутствием генетической рекомбинации при размножении, половым составом популяций, репродуктивным поведением и т.п. (Darevsky et al. 1985), что делает не возможным применение моделей для двуполых видов к видам партеногенетическим. Поэтому, изучение генетических процессов на однополых видах и разработка на основе этих данных математических моделей популяционных структур партеновидов может внести значительный вклад в понимание процессов дифференцировки популяций однополых видов, привести к разработке инструментов количественной оценки генетического вклада каждой популяции в генофонд вида и оценки эффективной численности партеновидов.

Объекты и методы В качестве входных данных использовалась матрица частот встречаемости 6 аллельных вариантов полиморфного локуса Du214 генома D.unisexualis. Общее число особей составляло 65 (из пяти изолированных популяций). Программная реализация модели предполагает использование генератора случайных чисел, при помощи которого были сформированы 5000 выборок из всех популяций партеновида по 10 особей в каждой.

Гетерозиготность определялась по формуле Нея (Хедрик, 2003).

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

Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия (mi1 )2 + (mi 2 )2 +... + (mik )2 di = l l l d s = mi1 + mi 2 +... + mik i =1 i =1 i =1 d s i = i, где: mij - абсолютная частота аллеля j (j = 1, 2, 3…k) в популяции i (i = 1, 2, ds 3,…,l); k – число аллелей в геномах особей исследуемых популяций; l – число популяций, составляющих исследуемый партеновид; - множитель (длина k-мерного вектора), характеризующий распространенность k аллелей в популяции i; - множитель (длина kмерного вектора), характеризующий распространенность k аллелей по виду в целом; - доля всех аллелей популяции j в ДНК всего вида (генетический вклад популяции j в генофонд вида).

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

Рисунок 1. Генетический вклад популяций в генофонд вида D.unisexualis

На основании анализа результатов можно сказать, что наибольший вклад в генофонд вносят популяции наименее изменчивые и наиболее генетически стабильные – популяции Лчап, Загалу и Норатус (см. рисунок 2). Дисперсия полученных значений генетического вклада и гетерозиготности равны нулю, т.к. все особи данных популяций содержат только два аллеля исследуемого локуса. В свою очередь, популяции Кутчак и Такярлу, населяющие высокогорные районы демонстрируют высокий уровень генетического разнообразия и низкий уровень генетического вклада (19,24(±0,4)% и 20,13(±0,2)% соответственно), вероятнее всего обусловленный небольшой численностью особей этих популяций.

–  –  –

Исследование поддержано грантом Президента РФ МK-2697.2011.4.

Литература Рысков А.П. Мультилокусный ДНК-фингерпринтинг в генетико-популяционных исследованиях биоразнообразия // Молекуляр. биология, 1999. Т. 33. № 6. С. 997-1011.

Хедрик Ф. Генетика популяций. М.: Техносфера, 2003.- 592с.

Darevsky I.S., Kupriyanova L.A., Uzzel T. Parthenogenesis in Reptiels // Boilogy of Reptilia. 1985.V.15. P. 412-526.

Schцn I., Martens К., van Dijk P. Lost Sex: The Evolutionary Biology of Parthenogenesis.- Springer, 2009.– P.615.

РАЗЛИЧНЫЕ ПОДХОДЫ К ОЦЕНКЕ ВЛИЯНИЯ ЭКОЛОГИЧЕСКИХ ФАКТОРОВ

ГОРНЫХ ТЕРРИТОРИЙ НА ИЗМЕНЧИВОСТЬ ПРИЗНАКОВ

ИНТРОДУЦИРУЕМЫХ СОРТОВ КАРТОФЕЛЯ

Оплеухин А.А.1, Стрельцова Т.А.2, Федюнина М.В.3 Горно-Алтайский государственный университет, Горно-Алтайск, Россия plymbym@rambler.ru Горно-Алтайский государственный университет, Горно-Алтайск, Россия tomagorny@yandex.ru Горно-Алтайский государственный университет, Горно-Алтайск, Россия marmed08@mail.ru Аннотация: пример комплексного использования нескольких методов статистической обработки для анализа данных экологического сортоиспытания картофеля в течении трех лет в зонах различной высотной поясности. В процессе анализа использованы методы позволяющие оценить взаимодействие факторов и всех признаков, с последующей оценкой вариации отдельных признаков в зависимости от факторов среды.

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

Объекты и методы В 1999-2001 годах было проведено экологическое сортоиспытание и генетический мониторинг 14 сортов картофеля в трех пунктах Горного Алтая разных по высотной поясности — в Майме (350м над уровнем моря), Чемале (630м) и Усть-Кане (1100м).

Испытания проводились согласно "Методическим указаниям по экологическому сортоиспытанию картофеля" (1982), с использованием и других современных методик.

Полевые опыты размещались синхронно по вертикальной зональности (площадь делянок 7,5м, повторность 4-х - кратная, рендомизированно).

Испытуемая коллекция включала генотипы 3-х групп спелости - ранние (Приекульский ранний, Корине, Новосибирский, Алмаатинский и Уральский сувенир);

среднеранние (Невский, Огонек, Адретта, Свитанок киевский, Эскорт и Гибрид 89/18) и среднеспелые – Луговской, Ласунок и Символ. Фиксировались показатели: продуктивность (грамм/куст), число клубней с одного куста, крупность (средняя масса одного клубня), высота стебля, число стеблей. Статистическая обработка была проведена при помощи модулей ППП STATISTICA: факторный анализ, главные компоненты и классификационный анализ, дисперсионный анализ; для построения графика распределения вариантов эксперимента в пространстве главных компонент использовался Microsoft Excel.

-1

-2

-3

-2,5 -2 -1,5 -1 -0,5 0 0,5 1 1,5 2 Рисунок 1. Центроиды, обозначающие положение показателей сортов, в пространстве главных компонент, в зависимости от года и места выращивания (треугольник – 1999 г, ромб

– 2000 г, квадрат – 2001 г: цвет белый – Майма, серый – Чемал, черный – Усть-Кан) Для рассмотрения распределения вариантов эксперимента мы использовали полученные координаты вариантов в пространстве главных компонент для построения графика (рисунок 1), из графика видно, что по оси второй главной компоненты (ось ординат) наиболее отличны варианты за 1999 год в Майме и Усть-Кане. По координатной оси первой главной компоненты (ось абсцисс) распределение наблюдается в основном по годам вегетации: наименьшие значения имеет 2000 г, около нуля 1999, и наибольшие, в целом, 2001 г, однако, с довольно сильными различиями по пунктам выращивания. Обратимся к данным по средним арифметическим признаков, для выяснения причин смещения в Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия пространстве главных компонент указанных вариантов. Продуктивность и высота стеблей были наименьшими в 2000 г (480 г/куст, 59см) особенно в Усть-Кане (399 г/куст. 49см) и наибольшими в 2001 (713 г/куст, разница по высоте стеблей с 2000 г не столь существенна).

В 1999 г в Майме наблюдалось большое число клубней с куста (16) и низкая крупность (39 г), в то время как в Усть-Кане малое количество клубней с куста (5,6), но высокая крупность (103 г), что, однако, не отразилось на продуктивности, в результате обратной корреляции указанные величины скомпенсировали продуктивность.

Для объяснения отклонений крупности и количества клубней в 1999 г в Усть-Коксе и Майме обратимся к метеорологическим данным. По-видимому, критическим было количество осадков в июле: в Майме на 30% ниже нормы, в Усть-Кане на 16% выше нормы.

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

АГРОКЛИМАТИЧЕСКИЙ ПОТЕНЦИАЛ СРЕДНЕЙ СИБИРИ ПРИ ИЗМЕНЕНИИ

КЛИМАТА В XXI ВЕКЕ

Парфенова Е.И. 1, Чебакова Н.М.1, Лысанова Г.И.2 Институт леса СО РАН, Красноярск, Россия lyeti@ksc.krasn.ru Институт географии СО РАН, Иркутск, Россия lysanova@irigs.irk.ru Аннотация: рассматриваются изменения агроклиматических ресурсов территорий, используемых для выращивания сельскохозяйственных культур на юге Средней Сибири. Данные метеонаблюдений с середины XX века и прогнозные сценарии изменения климата для XXI века показывают направленный тренд в увеличении теплоресурсов территории, что может повлечь изменение ассортимента и урожайности возделываемых культур.

Многочисленные исследования подтверждают тенденции в изменении климата с середины прошлого века и особенно в его последней четверти (IPCC 1990, 1995, 2001, 2007).

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

Цель нашего исследования:

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

В Сибири прогнозируется, что при изменении климата в течение XXI века 40% площади может быть занято лесостепью и степью за счет отступания лесов на север (Tchebakova et al., 2011). Именно эти растительные зоны являются пригодными для сельского хозяйства, и потепление климата может способствовать его продвижению на север. Общее агроклиматическое районирование подразделяет территорию на пояса, зоны, подзоны в соответствии с тепло- и влагоресурсами, а специальное агроклиматическое районирование дает основу для оценки возможности возделывания конкретных сельскохозяйственных культур (Синицина и др. 1973). Мы провели общее агроклиматическое районирование Средней Сибири для текущего и будущего климата с учетом трендов изменения климата.

Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия Для 100 метеостанций в регионе исследований были собраны данные (температуры января и июля и годовое количество осадков) за период 1960-2010 гг. (рис. 1а), по которым оценено изменение климата за период 1990-2010 гг. по отношению к базовому периоду 1960-1990 гг. Анализ данных показал, что за последние два десятилетия (1990-2010 гг.) зима и лето потеплели на 1-2°С (рис. 1 б, в) в лесостепной и степной зонах. Осадки в основном увеличились на 10%, хотя в межгорных сухих котловинах юга отмечается уменьшение осадков на 10-20% (рис. 1 г). В агроклиматическом районировании Красноярского края и Республики Хакасия (Атлас…,1994) наивысшей градацией по теплу указана сумма температур выше 10оС, равная 2000-2200, а в течение периода 1990-2010 гг. появляются агрорайоны с суммами температур 2200-2600.

–  –  –

В соответствии со сценариями изменения климата HadCM3А2 и HadCM3В1 Центра Хадли аномалии метеоэлементов на территории Средней Сибири достигнут к 2080 г. 4-5оС для июльской температуры, 5-8оС для январской температуры и 10-20% увеличения годовых осадков (www.ipcc-data.org). При таком изменении к 2080 г. увеличатся площади теплых и очень теплых агрорегионов на 8-10%, которые будут пригодны для использования в сельском хозяйстве при достаточном увлажнении или поливе. Также появятся жаркие и экстремально жаркие агрорегионы (до 10% по сценарию А2), которые окажутся непригодными для использования из-за недостатка увлажнения. Такие изменения в климатических ресурсах приведут к изменению как ассортимента культивируемых видов и их размещения, так и оптимумов их урожайности.

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

, 1973), построены ареалы основных культур в пределах исследуемой территории в современном климате и в 2080 году. Также, на основе собранных статистических данных по урожайности возделываемых агрокультур в течение 1966-2008 гг. в районах с различными климатическими условиями на юге Красноярского края, Хакасии и Тывы были найдены Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия корреляционные зависимости урожайности этих культур от климатических факторов – суммы градусо-дней выше 5оС и годового индекса увлажнения (рис. 2).

Климатические факторы объясняют 20-70% (R2 = 0.2-0.7) вариации урожайности основных культур. При потеплении в течение XXI века урожайность зерновых может увеличиться с 20 до 35 ц га-1, кукурузы на силос – с 300 до 600 ц га-1, картофеля с 100 до 200 ц га-1, ягодников – с 10 до 60 ц га-1. При этом климатический потенциал агрорегионов ограничивается почвенным потенциалом, за который можно принять почвенные условия степи, лесостепи, подтайги и южной тайги. Традиционные культуры (зерновые, картофель) можно постепенно продвигать на север, примерно 50 км за десятилетие, а также интродуцировать на юге новые культуры (кукурузу на зерно, некоторые сорта винограда, абрикос), которые будут определяться условиями зимы, или могут потребовать дополнительного полива. Таким образом, Агроклиматическое районирование Средней Сибири должно быть скорректировано в соответствии с наблюдаемыми и прогнозируемыми климатическими трендами.

Работа выполнена при поддержке РФФИ, грант 10-05000941 Литература Атлас Красноярского края и Республики Хакасия. Ред. В.И.Иванов. Роскартография.- Новосибирск. 1994.– 84 с Синицына Н.И, Гольцберг И.А., Струнников Э.А. Агроклиматология.– Л.: Гидрометеоиздат, 1973.– 344 с.

Tchebakova N.M, Parfenova E.I., and Soja A.J. Climate change and climate-induced hot spots in forest shifts in central Siberia from observed data // Regional Environmental Change.- 2011.

ОПРЕДЕЛЕНИЕ ПЛОТНОСТИ ПОПУЛЯЦИЙ ОРГАНИЗМОВ (НЕ ПОДЛЕЖАЩИХ

ПРЯМОМУ ПЕРЕСЧЕТУ) МЕТОДОМ «НУЛЕВЫХ ПРОБ» НА ПРИМЕРЕ

ФИТОПАРАЗИТИЧЕСКИХ НЕМАТОД

Перевертин К.А.

Центр паразитологии Института проблем экологии и эволюции им. А.Н.Северцова, Москва, Россия perevertink@mai.ru Аннотация: эффективное управление популяциями фитопаразитических нематод (nematode control) в агробиоценозах даже при развитом аппарате математического моделирования при производственном внедрении не может быть реализовано без достаточно простого технологически, но корректного статистически метода оценки главной входной переменной модели – плотности популяции нематод, Pi,(особей/100 см3 почвы). Апробированный нами метод нулевых проб (zero cores) позволяет не пересчитывать микроорганизмы в каждой пробе, а лишь констатировать их наличие/отсутствие. Доля (%) нулевых проб позволяет адекватно оценить Pi.

На долю фитопаразитических нематод приходится до 1/4 потерь мирового урожая от вредных организмов (Maggenty, 1981, Метлицкий, 1981). Хозяйственная значимость потерь сельхозпродукции от фитогельминтозов, неизмеримо возросшая в период интенсификации земледелия, предопределила повышенный научный интерес к этой группе фитопаразитов и, в частности, интерес к математическому моделированию. Успешности последнего способствуют биоэкологические особенности объекта моделирования/управления. Нематода не бегает, не прыгает, не летает, локализуется в слое прикорневой почвы растения-хозяина и имеет достаточно низкий потенциал размножения – сотни--тысячи, альтернативой миллионным порядкам у грибов и миллиардным у вирусов. Однако эйфория работы с «идеальным объектом», выразившаяся в построении множества красивейших матмоделей (имеющих прочный фундамент «от» основных модельеров Ферхюльста и Лотки с Вольтерра (Seinhorst, 1965, 1984, Ferris, 1976, 1984, etc.,)), сменилась определенным унынием, при попытках производственного внедрения, натолкнувшись на необходимость Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия полевой оценки главной входной величины модели – допосевной плотности популяции нематод – Pi. (Для удобства цитирования отметим Pi = ). На прямой вопрос практиков: «Так сколько же проб (cores) и какого объема, необходимо брать на гектар?» приходилось юлить, предлагая известные рекуррентные (секвенциальные) алгоритмы с необходимостью тестирования почвенных проб на каждом шаге еще до объединения их в смешанный образец (sample) по учетной единице.

Занятно, что текущий (2011) год можно назвать юбилейным после выхода фундаментальной статьи Taylor L.R. в журнале «Nature», 1961. За полувековой период многочисленные исследователи, подтверждая наличие степенной связи между средней и дисперсией S2=a расходятся в главном. Правило… или Закон? (Taylor's power low или Taylor’s power rule?) Статистически-приемлемое описание для множества частных случаев, или все же фундаментальное положение справедливое везде и всегда? А ведь Биология очень нуждается в ЗАКОНАХ. Очень заманчиво выглядела Тэйлоровская декларация о видоспецифичности параметра b, пик проверки которой в нематологии пришелся на 90-е (Brown, 1993, Boag &Topham, 1984, Молоканова, 1993 и др.), но она также не нашла убедительного подтверждения.

По нашему мнению, истина лежит, как всегда, посредине. Стаи птиц, стада сайгаков, пятна планктона и даже неравномерное (очаговое) распределение нематод в «равномерных»

агробиоценозах свидетельствуют в пользу гипотезы Тэйлора. Явный минус последней в корреляционном (не функциональном) характере зависимости, но это в огромной степени связано с несовершенством наших учетов, и уж ни в коей мере нельзя предварительно декларировать некорректность Правила для использования в прикладных разработках. Здесь уместно вспомнить о «до-Тэйлоровских» статистических законах распределения, которые многократно употреблялись (и употребляются) в биометрии и, в частности, применительно к нематодным популяциям: и логнормальное, и биномиальное, и Неймана типа А и т.д. и т.п.

Нельзя не упомянуть работы нашего многолетнего соавтора С.В. Васильева (ВИЗР), касающиеся, впрочем, учета не микроскопических вредителей.

Суть метода сводится к внешне парадоксальной рекомендации – учету доли «пустых»

почвенных проб, в которых нематод данного вида не обнаружено. Технологически достаточно простым методом для седентарных нематод является предложенный нами рулонный биотест (roulean bioassay, Перевертин, 1997, 2008). Математическое обоснование предполагает совместное использование правила Тэйлора и одного из законов распределений, наилучшим из которых признано отрицательное биномиальное (Noe, Barker,1985, Метлицкий, 2005).

Если р – вероятность обнаружения нематод исследуемого вида в почвенной пробе, то

q – обратная вероятность выражается следующим образом:

Параметр K отрицательного биномиального распределения (p-q)-K, отражающая степень агрегированности популяций может быть представлен так:

Долю нулевых проб (в %) обозначим через P0:

–  –  –

Таблица 4. Прогнозируемая плотность нематод M.

incognita Полученная аналитическим путем функциональная зависимость при практическом внедрении далеко не всегда оказывается адекватной. Рассчитанные для M.incognita теоретические зависимости «отказывались» работать, например, в Воронежской области.

Настройка факторозависимыми параметрами (К и а) зачастую оказывалась недостаточной.

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

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

ПРИМЕНЕНИЕ РЕГРЕССИОННО-ЛИНЕЙНОЙ МОДЕЛИ ДЛЯ АНАЛИЗА

ДОЛГОВРЕМЕННЫХ РЯДОВ КЛИМАТИЧЕСКИХ ДАННЫХ (на примере метеостанции г. Коломна, Южное Подмосковье) Петров А.С., Курганова И.Н.

Учреждение Российской Академии наук Институт физико-химических и биологических проблем почвоведения РАН, Пущино Московская обл., Россия PetrovAS-69@yandex.ru Аннотация: в результате анализа долговременных рядов климатических данных метеостанции г.

Коломна был выявлен отчетливый тренд увеличения среднегодовой температуры воздуха (с 1913 по 2008 г.) и показано отсутствие ярко выраженного тренда изменения количества осадков (с 1966 по 2008 г.), что свидетельствует об увеличении засушливости климата Южного Подмосковья.

Наиболее заметный рост среднегодовой температуры воздуха наблюдался после 1980 г и составлял более 90% от общего повышения температуры за анализируемый временной период.

Введение Одной из актуальных экологических проблем современности является изменение климата, выражающееся в увеличении средней глобальной температуры воздуха, составляющее у поверхности Земли 0.74оС за последнее столетие (Гулев и др., 2008).

Наиболее ярко изменение климата наблюдается в последние десятилетия и проявляется не только в более выраженном увеличении глобальной температуры, но и в усилении засушливости климата на значительной территории России (Золотокрылин и др., 2007).

Увеличение частоты засух наблюдается как в регионах с прогнозируемым снижением количества осадков, так и в областях, где количество осадков увеличивается в последствие изменения климата (Алифанов и др., 2006). Цель настоящего исследования состояла в

–  –  –

анализе долговременных рядов климатических данных на примере метеостанции г. Коломна (Южное Подмосковье) с применением регрессионно-линейных моделей.

Объекты и методы исследования Долгосрочные ряды наблюдений за температурой воздуха (1913-2008 гг.) и количеством осадков (1966-2008 гг.), полученные на метеостанции г. Коломна и опубликованные на сайте Росгидромета (http://aisori.meteo.ru/ClimateR), были проанализированы методом линейно-регрессионного анализа на месячном, сезонном и годовом уровнях осреднения. Для полученных временных трендов методом наименьших квадратов была оценена их достоверность (R2) и вероятность (Р, %). Интенсивность и знак наблюдаемых трендов характеризовали с помощью коэффициентов линейной регрессии в полученных статистических моделях.

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

–  –  –

В результате проведенного моделирования был выявлен значимый (P99%) тренд увеличения среднегодовой температуры воздуха (СГТВ) на 0.02оС·год-1 (Рисунок 1А) при отсутствии отчетливого тренда изменения количества осадков, что может косвенно свидетельствовать об увеличении засушливости климата Южного Подмосковья. Общий рост СГТВ составил 1.56 оС за период с 1913 по 2008 год и был обусловлен главным образом увеличением средне-зимних и средне-весенних температур на 0.03оС·год-1 (Р99%) в Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия среднем для каждого сезона (Рисунок 1Б). Достоверных трендов изменения температуры в летний и осенний периоды не выявлено. Наиболее отчетливый тренд повышения СГТВ наблюдался, начиная с 1980 года (Рисунок 1А), и составлял более 90% от общего повышения температуры за анализируемый временной период.

В зимний период наиболее значимое повышение температуры воздуха наблюдалось в декабре (0.03оС·год-1; Р97%) и январе (0.04оС·год-1; Р97%). В весенний период значимый рост температуры (Р99%) был характерен для всех трех месяцев и составлял 0.05, 0.03 и

0.01 оС·год-1 в марте, апреле и мае, соответственно (Таблица 1).

–  –  –

Литература Алифанов В.М., Гугалинская Л.А., Иванникова Л.А. Оценка и прогноз гидротермических условий почвообразования серых почв // Почвенные процессы и пространственно-временная организация почв.М.: Наука, 2006.- С. 471–494.

Гулев С.К., Катцов В.М., Соломина О.Н. Глобальное потепление продолжается // Вестник РАН, 2008.- Т. 78.- № 1.- С. 20-27.

Золотокрылин А.Н., Виноградова В.В., Черенкова Е.А. Динамика засух в европейской России в ситуации глобального потепления // Проблемы экологического мониторинга и моделирования экосистем.

СПб.:

Гидрометеоиздат, 2007.- Т. XXI.- С. 160–182.

МОДЕЛИ ДИНАМИКИ ЧИСЛЕННОСТИ И ПРОСТРАНСТВЕННОГО

РАСПРЕДЕЛЕНИЯ ВАЖНЕЙЩИХ РЕСУРСНЫХ ВИДОВ ДИКИХ КОПЫТНЫХ

(ЛОСЯ, КОСУЛИ, КАБАНА) РОССИИ НА ОСНОВЕ ДАННЫХ МНОГОЛЕТНЕГО

МОНИТОРИНГА

Петросян В.Г., Дергунова Н.Н., Бессонов С.А., Омельченко А.В.

Институт проблем экологии и эволюции им. А.Н. Северцова РАН, Москва, Россия petrosyan@sevin.ru Аннотация: представлены модели декомпозиции трендов, периодических и случайных составляющих динамики численности важнейших охотничьи – промысловых видов лося, европейской/сибирской косули и кабана. Показано, что длина интервалов цикличности колебания для лося составляет 10 лет, а для косули и кабана - 11 лет. Приведены модельные оценки репродуктивного потенциала, годового воспроизводства и максимально возможная численность популяций изучаемых видов.

Введение В настоящее время существуют две основные гипотезы, объясняющие причины и механизмы динамики низкой численности важнейших ресурсных парнокопытных животных в России. Сторонники первой гипотезы полагают, что многолетние колебания численности животных вызваны с глобальными изменениями климата и сукцессиями растительности и проявляются с интервалами от 5 до 120 лет[Данилкин, 2007]. Согласно второй гипотезе, основной фактор, определяющий численность копытных в России за последние 50 лет, имеет “антропогенно – хищнический” характер, т.е. утверждается, что динамики численности Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия может быть описана в терминах обобщенной модели “хищник – жертва”. По нашему мнению, динамические изменения численности животных имеют многофакторный характер и разработка походов для анализа и выявления закономерностей этих изменений во времени и пространстве необходимы для понимания и управления ресурсными видами диких копытных России.

Основная цель работы – анализ динамики численности и характера пространственного распределения важнейших ресурсных видов лося (A. alces L.), европейской / сибирской косули (C. capreolus L., C. pygargus Pall.) и кабана (S. scrofa L.) на основе данных многолетнего мониторинга в 1981-2007 гг.

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

Блок Модели МА(p,q):

оценивает тренд методом скользящего среднего (МА) с параметрами шага усреднения p и полиномами степени q, включает сервисные возможности для подбора параметров p и q.

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

Блок АРМА: определяет сезонную и несезонную модель авторегрессии - скользящего среднего (ARMA), используя выборочные коэффициенты АКФ и ЧАКФ. Блок Состояние популяций: включает оценку репродуктивного потенциала популяций, экологической емкости среды обитания, величины годового воспроизводства и интенсивности миграционных процессов на основе модифицированных дискретных моделей БивертонаХолта, Рикера и Мальтуса.

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

В целом существующие оценки по плотности лося наглядно показывают, что количество областей, в которых плотность принимает значение менее 1 особей/1000 равно 37, что составляет 91 % от всей площади ареала. Относительная площадь 32 регионов, на территории которых плотность лося варьирует в диапазоне 1-4 особей на 1000 лесной площади (ЛП), составляет 7% площади ареала. В ограниченном количестве регионов 7 (2 % площади ареала) плотность населения больше 4 особей на 1000 га. Эти оценки позволяют утверждать, что взвешенная оценка средней плотности лося по России составляет 0.71 особей на 1000 га ЛП.

Суммарная площадь областей, на территории которых плотность косули больше 30 особей на 1000 га ЛП, составляет менее 1 % площади ареала. Относительная площадь Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия территории, на которых плотность составляет 0.01-0.091; 1.01 – 4.92; 5.04-9.78;10.1-19.59 особей на 1000 га ЛП, соответственно равны 65, 22, 8, и 4 %%. Эти оценки указывают, что взвешенная оценка средней плотности косули по России имеет значение 2.33 особей на 1000 га ЛП.

Доля площади арела кабана, на которой плотность больше 13 особей, находится в интевалах 8-12, 3-8, 1-3 и меньше 1 особей на 1000 га ЛП составляют 0.3, 0.8, 2.1, 6.9, 89.9 %% соответственно. Средневзвещанная плотность кабана по России равняется 0.84 особей на 1000 га ЛП.

Модельные оценки численности животных в России показывают, что в природных зонах животные распределены неравномерно в широтном направлении. Наиболее высокие плотности животных лося, косули и кабана достигаются в зоне хвойно-широколиственных лесов: 6.7 – 13.6, 27.1-54.3 и 21-42 особей на 1000 га ЛП. В таежной зоне искомые величины для лося, косули и кабана изеняются в диапазонах 4.94-9.88, 19.7-39.5 и 15.3-30.6 особей на 1000 га ЛП, соотвестственно. В лесостепной зоне средние плотности составляют 3.26 – 6.57, 13.1-26.3, 10.1-20.3 особей для лося, косули и кабана, соответственно. В зоне притундровых лесов и редкостойной тайги искомые показатели составляют 0.2 -0.57, 4.9- 9.8, и 3.8-7.68 особей на 1000га ЛП, соответственно.

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

Исследование поддержано программой ОБН РАН «Биологические ресурсы России:

Фундаментальные основы рационального использования» и госконтрактом МинОбрнауки № 02.740.11.0867 Литература Данилкин А.А. К проблеме цикличности в динамике численности диких копытных животных. Вестник охотоведения.- 2007.- Т. 4.- № 1.- С. 89-99.

Петросян В.Г., Дергунова Н.Н., Бессонов С.А., Назарова К.А., Омельченко А. В. Оценка динамики численности и половозрастной структуры лося (A. Alces L.) в России с использованием имитационной модели и данных многолетнего мониторинга. Ж. Известия Самарского научного центра РАН, № 11, 1(7).- 2009.С. 1566-1574.

Петросян В.Г., Дергунова Н.Н., Бессонов С.А. Комплексная имитационная модель оценки динамики численности, половозрастной структуры и оптимальной плотности лося (Alсes alсes L.) на основе данных наземных и дистанционных измерений // Успехи современной биологии.- 2010, т. 130.- № 1.- С.

88-99.

К МЕТОДОЛОГИИ МОДЕЛИРОВАНИЯ ЭКОСИСТЕМ: ПЕРСПЕКТИВЫ СИНТЕЗА

МЕРОЛОГИЧЕСКОГО И ХОЛИСТИЧЕСКОГО ПОДХОДОВ

Печуркин Н.С.

ИБФ СО РАН, Академгородок, Красноярск-36, Россия, nsla@akadem.ru Аннотация: в работе обсуждаются методологические особенности двух основных подходов в экологическом моделировании: мерологического и холистического. Предлагается объединение двух описанных подходов на основе учета циклов лимитирующих веществ в экосистемах.

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

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

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

Метод черного ящика не способен улавливать и дифференцировать циклы, ибо они могут быть не связаны с входами и выходами системы, замыкаясь внутри нее. И здесь опять возникает проблема неоднозначности и множественности траекторий движения и развития так же, как и при имитационном моделировании. Действительно, разные типы систем в определенных условиях могут иметь сходные реакции и в таком случае мы их соотнесем к одинаковому типу. Но при смене условий замаскированные в них разные циклы и подциклы поведут себя совсем по-другому, и вся наша классификация рассыпется в прах. По образному выражению профессора Р. Улановича, «как известно, природа не терпит пустоты, но следует также сказать, что логика не терпит циклов». Таким образом, циклы, при "черноящичном" подходе ускользают из-под изучения.

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

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

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

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

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

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

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

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

ЭНЕРГЕТИЧЕСКИЙ КРИТЕРИЙ ДЛЯ ПРОГНОЗИРОВАНИЯ ДИНАМИКИ

ПОПУЛЯЦИЙ ЧЕЛОВЕКА И ВЫСШИХ ХИЩНИКОВ: ЛОГИСТИЧЕСКАЯ

МОДЕЛЬ

Печуркин Н.С., Сомова Л.А.

ИБФ СО РАН, Академгородок, Красноярск-36, Россия, nsla@akadem.ru Аннотация: в работе обсуждается энергетический критерий для количественного описания динамики популяций человека и высших хищников на основе коэффициента потребления энергии на душу населения. Предлагается простая логистическая модель зависимости показателей прироста популяции от энергообеспеченности.

Описания популяционной динамики, главным образом, основаны на концепции “демографического перехода”, согласно которой численность человеческой популяции, Материалы Второй конференции «Математическое моделирование в экологии»

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

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

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

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

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

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

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

Можно записать уравнение для относительной (удельной) скорости роста популяции человека в виде :

dn / dt = (dN/dt) / N = a( h жук – h) / h жук, Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия где а – максимальная удельная скорость роста, генетически присущая данной популяции, N - численность анализируемой популяции, h жук - величина энергообеспеченности, соответствующей Желаемому Уровню Комфорта, h - величина энергообеспеченности, достигнутая к моменту - t описываемого времени.

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

Таким образом, данная модель описывает все четыре последовательных стадии демографического развития человечества: от самой ранней, с высокими скоростями рождаемости до поздней, с нулевым приростом. Внешне формула дает вполне похожее на логистическую кривую ограничение роста численности. Но «наполнение» ее совершенно иное. Не тривиальный Голод, а Энергетика и Информация правят миром! Конечно, сама величина Желаемого Уровня Комфорта (ЖУК) является многопараметрической (многовекторной, многолучевой), но в данном случае полагается, что именно ее энергетическая составляющая наиболее информативна.

Перечислим основные преимущества введенной меры:

1. Она относится к основным показателям жизнеобеспечения человека, так как энергетика имеет прямое отношение к «жизнесберегающим технологиям».

2. Она имеет строгое количественное выражение.

3. Она понятна всему миру, каждому из нас, независимо от расы и вероисповедания (буквально: киловатт – он и в Африке киловатт).

4. Она не подвержена «инфляционным подвижкам», как это неизбежно происходит с экономическими показателями. Например, при сравнении величин ВВП в денежном выражении всегда приходится прибегать к трудно измеримым паритетам покупательной способности [ППС].

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

ПРИМЕНЕНИЕ ТЕОРИИ БУФЕРНОЙ ЕМКОСТИ ДЛЯ ОЦЕНКИ ПРИРОДНОГО

КОНТРОЛЯ СОДЕРЖАНИЯ ИОНОВ ТЯЖЕЛЫХ МЕТАЛЛОВ В ЗАГРЯЗНЕННЫХ

ВОДАХ Повар И.Г.

Институт химии Академии Наук, Кишинев, Республика Молдова ipovar@yahoo.ca Аннотация: недостаток информации о буферных свойствах водных экосистем по отношению к ионам тяжелых металлов (ТМ) объясняется отсутствием количественной теории, способной предсказать их химическую устойчивость по отношению к изменению равновесного состава.

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

Хотя аналитические методы являются необходимым инструментом для исследования Материалы Второй конференции «Математическое моделирование в экологии»

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

Экосистемы обладают буферным свойством поглощать существенное количество ионных загрязнителей, однако если последние находятся либо возникают в больших количествах, то они способны кумулятивно перегрузить систему и как следствие не смогут быть поглощены. Для оценки риска загрязненных экосистем или естественного подавления загрязнителей необходимо обязательно знать как и до какой степени они распространятся либо их концентрации убавятся в будущем. Таким образом, проблема контроля концентрациями ионов ТМ в естественных системах может быть сформулирована следующим образом: каков будет ответ системы при изменении параметров, которые управляют pМ = -log aM (отрицательный логарифм активности иона металла)? Другими словами, какова величина буферной емкости pМ системы? В отсутствие буферной емкости концентрации загрязнителя водной экосистемы изменялись бы очень быстро. При повышенных концентрациях, когда приток ионных загрязнителей превышает их буферную емкость, имеют место токсичные эффекты. Любая буферная система имеет определенный потенциальный (буферный) запас, который управляет равновесием и поддерживает один из равновесных параметров (например, равновесные концентрации) постоянным.

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

Главной целью нашего сообщения является количественная оценка способов перемещения и преобразования ионов тяжелых металлов (ТМ) в окружающей среде и их распространение в экосистемах при помощи разработанной нами теории буферной емкости (ТБЕ). Ожидается, что ТБЕ поможет предсказывать долгосрочные эффекты в естественном подавлении ТМ в загрязненных водах. Полученные результаты предоставят исследователям инструмент для установления надежных пределов уровней концентраций ионов ТМ в водных объектах окружающей среды. Применение ТБЕ углубит знание и понимание процессов, которые управляют концентрациями компонентов, и сможет служить надежным математическим инструментом для оценки и предсказания долгосрочных эффектов в естественном подавлении ТМ в загрязненных экологических системах. Теория химических равновесий, диаграммный способ распределения химических форм, механизм действия буферной емкости и оценка буферных свойств двухфазных систем были развиты нами в (Fishtik, 2006, Povar, 1996, Povar, 2000, Povar, 2000, Povar, 2003, Повар, 2011, Повар, 2011). В этих работах выведены теоретические соотношения для основных буферных характеристик.

Буферная емкость двухфазной (гетерогенной) системы значительна больше чем аналогичная величина для монофазной (гомогенной) системы в области рМ ее максимального значения.

В работе (Повар, 2011) ТБЕ применена для оценки ион-молекулярной буферной емкости, X, в кислых природных водах в равновесии с минералом гиббсит (твердая фаза) относительно pX, где pX =-log [X] (pX =pH для кислотно-щелочных буферов и pAl для буфера по иону алюминия). Гиббсит Al(OH)3(s) имеет основной контроль над растворимыми алюминиевыми формами при данном значении рН, однако с увеличением степени комплексообразования Al уменьшается его свободная концентрация в природных водах. Эти воды имеют тенденцию поддерживать довольно постоянное значение pAl даже когда ионы Al3+ могут быть введены из внешних источников. Понятие буферной емкости очень полезно Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия в оценке кислотно-щелочного статуса кислых природных вод наряду с оценкой Al буферного эффекта. Развитый подход учитывает оценку буферности H и Al в природных водах, принимая во внимание доминирование Al буферности. Изучено влияние различных факторов (температуры, pH и концентраций главных компонентов в природных водах) на буферные характеристики этих буферных двухфазных систем. Наш подход был проверен на реальных водных качественных данных и был применен к практическому анализу природных водных образцов. Установлено, что природные воды в равновесии с гиббситом обладают определенным потенциальным резервом, посредством которого регулируется равновесие и поддерживается постоянным один из его параметров. Присутствие лигандов, образующих устойчивые комплексы с ионами твердой фазы, способствует усилению буферного действия гетерогенных систем.

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

Литература Fishtik I., Povar I. Buffer Capacity in Multiple Chemical Reaction Systems Involving Solid Phases // Can. J. Chem.– 2006.– V. 84.– P. 1036-1044 Povar I. Buffer properties of heterogeneous systems containing sparingly soluble acids and metal hydroxides // Russ. J.

Inorg. Chem.– 1996.– V. 41.– P. 1167-1172 Povar I. On homogeneous ion buffer systems // Russ. J. Inorg. Chem.– 2000.– V. 45.– P. 1628-1631 Povar I. Buffer properties of heterogeneous water-salt systems in relation to the components of a low-solubility precipitate // Russ. J. Inorg. Chem.– 2000.– V. 45.– P. 1632-1636 Povar I., Luca C. Considerations regarding ionic-molecular buffer solutions // Rev. Chim.– 2003.– V. 54(4).– P.312Повар И., Русу В. Новый тип диаграмм распределения растворимых и нерастворимых форм алюминия в натуральных гетерогенных водных системах // Вода: Химия и экология.– 2011.- №1.– С. 44 – 49.

Повар И., Русу В. Теоретические аспекты ион-молекулярной буферности в природных водах в равновесии с минеральной фазой гиббсит // Вода: Химия и экология.– 2011 (в печати).

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И ГРАФИЧЕСКОЕ ИЗОБРАЖЕНИЕ

РАВНОВЕСНОГО СОСТАВА МНОГОКОМПОНЕНТНОЙ ГЕТЕРОГЕННОЙ

СИСТЕМЫ "НАСЫЩЕННЫЙ ВОДНЫЙ РАСТВОР - МИНЕРАЛ"

Повар И., Русу В.

Институт химии Академии Наук, Кишинев, Республика Молдова ipovar@yahoo.ca Аннотация: моделирование химических равновесий в природной гетерогенной системе сведено к решению системы нелинейных уравнений материального баланса, учитывающих в явном виде распределение химических форм иона алюминия в жидкой и твердой фазах. На основе рассчитываемого равновесного состава многокомпонентной системы разработан новый тип диаграмм гетерогенных равновесий, количественно описывающих распределение химических форм в обеих фазах.

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

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

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

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

Для полной количественной характеристики изучаемых систем нами введены дополнительные переменные и дополнительные уравнения МБ, учитывающие вклад и стехиометрический состав твердых фаз. Дополнительные переменные представляют собой количества осажденных i ионов-компонентов минерала в 1 л раствора, C i, моль/л либо количество молей твердой фазы Р в 1 л раствора. Показано, что эти величины взаимосвязаны с растворимостью твердой фазы S. Предложенная запись учета МБ позволяет свести систему N нелинейных уравнений МБ для известного исходного состава гетерогенной смеси (т.е.

известных исходных (общих, аналитических) концентраций компонентов Ci0 ) с N+1 переменными к системе N нелинейных уравнений с N неизвестными (равновесными концентрациями), которое легко решается численными методами. Остальные неизвестные равновесные концентрации компонентов гетерогенной смеси, входящие в систему нелинейных уравнений, легко выражаются через эту переменную. Одна из особенностей предложенной записи МБ состоит в том, что полученные на их основе результаты расчета равновесного состава могут быть непосредственно быть использованы для построения диаграмм гетерогенных химических равновесий (ДГХР), которые в явном виде отображают распределение химических форм данного компонента в водной и твердой фазах.

В гетерогенных равновесиях мольные доли i зависят также от исходного состава смеси и, следовательно, являются функциями дополнительных переменных: i = f(Ci0, Cj0, …, pH). Такие диаграммы удобно строить в координатах pH - i, или Ci0 - i при фиксированных значениях других переменных. Процедура построения диаграмм гетерогенных химических равновесий включает в себя следующие этапы: (а) термодинамический расчет области устойчивости твердой фазы (природного минерала). Твердая фаза устойчива при Gtot 0.

Условие Gtot = 0 соответствует началу его растворения (осаждения); (б) расчет мольных долей i всех химических форм содержащих данный компонент в твердой фазе; (в) для полноты описания картины вне области устойчивости твердой фазы, т.е. для гомогенного Материалы Второй конференции «Математическое моделирование в экологии»

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

–  –  –

Рисунок 1. Диаграмма гетерогенных химических равновесий в зависимости от рН в системе насыщенный водный раствор– гиббсит.

Концентрации (моль/дм3): COrg, C PO4,

–  –  –

На рисунке 1 приведена ДГХР для заданного состава гетерогенной системы насыщенный водный раствор – гиббсит Al (OH ) 3( S ). Такая диаграмма наглядно отражает распределение различных химических форм алюминия как в твердой фазе, так и в насыщенном растворе. Отметим, что при значении рН растворения-образования твердой фазы мольные доли растворимых и нерастворимых форм алюминия резко изменяются.

Для C Al =·10-5 М в кислой и нейтральной средах равновесные концентрации моноядерных и полиядерных гидроксокомплексов алюминия малы. Для CAl01·10-3 М в кислых и слабокислых средах из растворимых форм преобладают комплексные частицы AlOrg и AlHOrg+. При концентрации CAl0 = 1·10-5 М ощутимый вклад вносят фторидные комплексы алюминия AlF2+ и AlF2+. На основе разработанного подхода возможно изучить, кроме исходного состава смеси, влияние других факторов влияющих на распределение и концентрации растворимых и нерастворимых частиц алюминия в водных системах, как например температура и корректный выбор используемых констант равновесия.

Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия

ИССЛЕДОВАНИЕ СВОЙСТВ ПРОСТРАНСТВЕННО ОДНОРОДНОЙ

МАТЕМАТИЧЕСКОЙ МОДЕЛИ ЧЕТЫРЕХКОМПОНЕНТНОЙ (N-P-Z-D)

ПЛАНКТОННОЙ ПИЩЕВОЙ ЦЕПИ

Подгорный К.А.

Учреждение Российской Академии наук ИБВВ им. И.Д. Папанина РАН, пос. Борок, Россия kap@ibiw.yaroslavl.ru Аннотация: рассматриваются наиболее важные свойства четырехкомпонентной N-P-Z-D-модели планктонной пищевой цепи, в которой используется нелинейная зависимость при описании процесса минерализации детрита, а также учтено отрицательное влияние метаболитов фитопланктона на развитие зоопланктона. В модели учтено также, что существуют определенные пищевые предпочтения в питании зоопланктона: он может питаться не только фитопланктоном, но и детритом.

Цель данной работы состояла в аналитическом и численном исследовании динамического поведения пространственно однородной математической модели четырехкомпонентной планктонной системы. Модель включает в себя биогенные элементы (N), фитопланктон (P), зоопланктон (Z) и планктонный детрит (D). Наше исследование, прежде всего, будет касаться изучения наиболее важных свойств модели в зависимости от изменений порядка ферментативной реакции, которая описывает разложение детрита, а также скорости выделения метаболитов клетками фитопланктона. Ранее в работе (Подгорный, 2009) исследовался более простой вариант модели, в котором предполагалось, что зоопланктон может потреблять только фитопланктон. В данной работе будет рассматриваться более реалистичный вариант модели, который учитывает то, что зоопланктон на самом деле может питаться не только фитопланктоном, но и детритом. В природных экосистемах у зоопланктона существуют определенные пищевые предпочтения в питании, которые могут изменяться с течением времени в зависимости от тех или иных условий. Поэтому важно оценить влияние и этого фактора на характер динамики модельной системы.

Модель представляет собой систему четырех обыкновенных дифференциальных уравнений:

(P + D ) dN N aP dt = e + N b + cP + + P + D Z + qZ + D + k (N 0 N ),

–  –  –

экскреции пищи зоопланктоном; - коэффициент регенерации пищи, не усвоенной организмами высших трофических уровней; - скорость выделения метаболитов (токсических компонентов) клетками фитопланктона. Все переменные модели неотрицательны и выражены в граммах углерода на м3. Параметры модели неотрицательны, причем + 1, 1 и r + s + k a b. Последнее условие означает, что фитопланктон не вымирает полностью при очень низкой концентрации биогенных элементов. Пищевые предпочтения зоопланктона учитываются в новом варианте модели с помощью параметра = p 2 p1, где p1 и p 2 - коэффициенты предпочтения зоопланктоном фитопланктона и детрита соответственно. Включение в модель дополнительного источника питания для зоопланктона добавляет в модель всего один новый параметр - коэффициент предпочтения детрита зоопланктоном. Если = 0, то это означает, что p 2 = 0, p1 0, то есть зоопланктон не питается детритом, а источником пищи для него является только фитопланктон. В этом случае мы получаем модель, рассмотренную ранее в работе (Подгорный, 2009). В природных экосистемах значение параметра, по всей видимости, близко к 0,5. В модельных исследованиях используют значения [0;2] (Popova et al., 1997).

Численные эксперименты с моделью (1) проводились методом Рунге-Кутта-4.

Значения параметров модели выбраны на основе литературных данных (Долгоносов, 2009;

Edwards, 2001; Popova et al., 1997) Нами были приняты следующие значения параметров:

e = 0,03 гС/м3; a = 0,2 (м сут)-1; b = 0,2 м-1; c = 0,4 м2/гС; r = 0,15 сут-1; = 0,33 ; = 0,6 сут-1; = 0,035 гС/м3; = 0,5 ; q = 0,1 сут-1; k = 0,05 сут-1; N 0 = 0,4 гС/м3; s = 0,04 сут-1;

= 0,25 ; [0,0;0,3] м3/(гС сут); = 0,1 сут-1; [0,5;1,2] ; = 0,08 сут-1, [0;2].

Для системы уравнений (1), возможны следующие стационарные состояния:

E 0 = (N 0,0,0,0), E 1 = (N 1, P1,0, D1 ), E 2 = (N 2, P2, Z 2, D 2 ). Аналитическое исследование модели показало, что каждое них является единственными. Для обеспечения их устойчивости необходимо соблюдение определенных условий.

Путем численных экспериментов было показано, что все возможные типы динамического поведения модели (1) происходят в определенных пределах изменения параметра - от нулевого его значения до некоторой критической величины * = f (, ).

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

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

Литература Долгоносов Б.М. Нелинейная динамика экологических и гидрологических процессов / Отв. редактор М.Г.

Хубларян. Предисловие Г.Г. Малинецкого. – М.: Книжный дом «ЛИБРОКОМ», 2009. – 440 с.

Подгорный К.А. Четырехкомпонентная (N-P-Z-D) модель планктонной пищевой цепи // Материалы Национальной конференции с международным участием «Математическое моделирование в экологии»

(1 – 5 июня 2009 г., г. Пущино). – Пущино, ИФХ и БПП РАН, 2009, с. 209 – 210.

Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия Edwards A.M. Adding detritus to a nutrient - phytoplankton - zooplankton model: a dynamical-systems approach // J.

Plankton Res., 2001, V. 23, N. 4, P. 389 - 413.

Popova E.E., Fasham M.J.R., Osipov A.V., Ryabchenko V.A. Chaotic behaviour of an ocean ecosystem model under seasonal external forcing // J. Plankton Res., 1997, v.19, p. 1495 – 1515.

ИССЛЕДОВАНИЕ СРЕДНЕМНОГОЛЕТНЕЙ ДИНАМИКИ БАЛАНСОВ

БИОГЕННЫХ ВЕЩЕСТВ В ЭКОСИСТЕМЕ НЕВСКОЙ ГУБЫ ФИНСКОГО ЗАЛИВА

(ПО ДАННЫМ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ)

Подгорный К.А.

Учреждение Российской Академии наук ИБВВ им. И.Д. Папанина РАН, пос. Борок, Россия kap@ibiw.yaroslavl.ru Аннотация: с помощью имитационной математической модели проведено исследование среднемноголетней динамики балансов БВ, определены наиболее характерные особенности формирования запасов форм азота и фосфора в Невской губе Финского залива за счет естественного круговорота, а также за счет антропогенной нагрузки.

Цель данной работы состояла в том, чтобы путем проведения численных экспериментов исследовать наиболее важные закономерности формирования запасов и среднемноголетней динамики форм азота и фосфора в Невской губе Финского залива за счет естественного круговорота, а также за счет антропогенной нагрузки. Биохимический блок имитационной модели включает в себя описание совместной трансформации форм N и Р (важнейших для экосистемы Невской губы биогенных элементов), поступающих в эту акваторию со сточными водами и речным стоком (Подгорный, 2003, 2009; Подгорный, Леонов, 2011). Модель включает в себя весь известный спектр БВ указанных элементов, а также основных представителей сообщества микроорганизмов, обитающих в водной среде и участвующих в трансформации форм N и Р в естественных условиях. В модели рассматривается 17 переменных состояния водной среды: концентрации растворенных фракций органического N (DON) и P (DOP); минеральных форм P (DIP) и N – аммония ( ) ( ) ( ) NH +, нитритов NO и нитратов NO3 ; взвешенных форм N и P – в составе детрита (ND и PD); биомасс гидробионтов – гетеротрофного бактериопланктона (BN и BP), фитопланктона (FN и FP), простейших (PRN и PRP) и зоопланктона (ZN и ZP), а также концентрация растворенного в воде O 2. Структура взаимодействия этих переменных показана схематично на рисунке 1. Оно осуществляется в естественных условиях под совместным влиянием как природных, так и антропогенных факторов.

В результате обработки данных моделирования было показано, что среднемноголетний внешний приток BN с мая по октябрь изменялся в пределах 220-245 т/мес, PhN – 59-294 т/мес, PrN – 13,6-73,2 т/мес, ZoN – 2,4-12,7 т/мес, ND – 920-1035 т/мес, DON – 1181-2539 т/мес, NH + – 441-627 т/мес, NO 3 – 972-1510 т/мес, NO – 22-37 т/мес.

Таким образом, установлено, что во внешнем притоке в Невскую губу основной вклад в структуре баланса азота вносят DON, ND и NO 3. На долю этих трех составляющих баланса приходится в среднем 80 % внешнего притока всех азотсодержащих биогенных веществ.

Среднемноголетний вынос BN из экосистемы Невской губы в восточную часть Финского залива с мая по октябрь изменялся в пределах 135-198 т/мес, PhN – 83-429 т/мес, PrN – 7,2-94,4 т/мес, ZoN – 4,5-21 т/мес, ND – 1131-1358 т/мес, DON – 1221-2510 т/мес, NH +

–  –  –

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

Показано, что среднемноголетний внешний приток BP с мая по октябрь изменялся в пределах 22,3-24,3 т/мес, PhP – 5,9-29,4 т/мес, PrP – 1,4-7,3 т/мес, ZoP – 0,24-1,27 т/мес, PD – 92-103,5 т/мес, DOP – 22,8-36,0 т/мес, DIP – 12,1-35,1 т/мес. Таким образом, установлено, что основной вклад в структуре баланса фосфора во внешнем притоке в Невскую губу вносят PD, DOP и DIP. На долю этих трех составляющих фосфорного баланса приходится в среднем от 74 до 81 % внешнего притока всех фосфорсодержащих биогенных веществ.

Среднемноголетний вынос из экосистемы Невской губы в восточную часть Финского залива BP с мая по октябрь изменялся в пределах 14,4-27,5 т/мес, PhP – 8,3-40,5 т/мес, PrP – 0,69-9,4 т/мес, ZoP – 0,42-2,07 т/мес, PD – 104,8-126,1 т/мес, DOP – 17,8-42,3 т/мес, DIP – 28,4-53,1 т/мес, Таким образом, основной вклад в структуре баланса фосфора во внешнем притоке в восточную часть Финского залива по-прежнему вносят PD, DOP и DIP. На долю этих трех составляющих фосфорного баланса приходится в среднем от 74 до 86 % общего выноса в Финский залив фосфорсодержащих биогенных веществ.

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

Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия Литература Подгорный К.А. Математическое моделирование пресноводных экосистем нестратифицированных водоемов (алгоритмы и численные методы).- Рыбинск: Изд-во ОАО “Рыбинский Дом печати”, 2003.- 328 с.

Подгорный К.А. Пространственно- неоднородная имитационная модель экосистемы Невской губы Финского залива // Известия Самарского научного центра РАН.- 2009.- Т.- 11.- № 1 (7).- С. 1554–1558.

Подгорный К.А., Леонов А.В. Использование пространственно-неоднородной имитационной модели для изучения процессов трансформации соединений азота, фосфора и динамики кислорода в экосистеме Невской губы Финского залива. Водные ресурсы. 2011 (в печати).

МАТЕМАТИЧЕСКАЯ БИОЛОГИЯ: А.А.ЛЯПУНОВ И СОВРЕМЕННОСТЬ3

Полуэктов Р.А.

Агрофизический Институт, Санкт-Петербург, Россия r_poluektov@mail.ru Аннотация: современный этап развития математической биологии насчитывает более, чем сорокалетнюю историю, хотя ее корни уходят к 1202 году. Понятие математической биологии принято соотносить к тем ее разделам, который занимается развитием динамической теории надорганизменных структур – популяций и биоценозов. В докладе выделяются основные этапы формирования теории, начиная с основополагающих работ А.А.Ляпунова до наших дней.

В книге Liber abacci, написанной в 1202 г., итальянский математик Фибоначчи, сформулировал задачу нахождения потомства от пары кроликов. Это первая известная нам динамическая модель популяции, обладающей возрастной структурой, в которой особи никогда не умирают и ежемесячно приносят постоянное потомство, начиная со второго месяца после своего рождения. Остается только удивляться тому, что следующий шаг в направлении развития модели возрастной структуры популяции – учет смертности особей – был сделан лишь в 1945 г.(Lesli P.H). И все же до второй половины 20-го века имелись лишь разрозненные попытки постановки и решения отдельных задач, относящихся к области моделирования популяций и биоценозов (Malthus, 1820; Ferhulst, 1838; Fisher, 1930; Volterra, 1931, Kostizin 1937). В конце 60-х годов выходят первые отечественные работы по математической теории популяций и биоценозов (Свирежев; Корзухин, Жаботинский;

Базыкин, Гинзбург, Ратнер, Гильдерман, Кудрина, Полетаев). Их появление во многом стимулировалось распространением идей и методов новой тогда науки – кибернетики, а толчком послужила публикация монографии И.А.Полетаева «Сигнал». Однако главным событием, которое стимулировало взрыв исследований по математической биологии, явились летние молодежные школы, проводившиеся в течение ряда лет в Подмосковье под руководством знаменитого биолога Николая Владимировича Тимофеева-Ресовского и чл.корр. АН СССР кибернетика и математика Алексея Андреевич Ляпунова. Первые две школы (1967 и 1968 гг.) проводились в Можайском и Аксакове (под Москвой), а затем несколько раз в Звенигороде и далее они плавно переехали в Пущино-на-Оке. Влияние этих выдающихся ученых на развитие математической биологии невозможно переоценить. Николай Владимирович Тимофеев-Ресовский, будучи ученым-энциклопедистом, так формулировал биологическую проблему, что она оказывалась совершенно готовой к математической алгоритмизации. А Алексей Андреевич Ляпунов в момент становления кибернетики как науки внес решающий вклад в такие ее разделы как теория программирования, математическая лингвистика и математическая биология. Издававшийся А.А.Ляпуновым журнал «Проблемы кибернетики» был единственным изданием, в котором тогда публиковались работы по этим направлениям.

Работа выполнена при финансовой поддержке РФФИ, проект № 04-05-00415 Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия Определенную роль в формировании математической биологии как науки сыграла и специальная школа, получившая название «Математические методы в теории популяций и ценозов» (Р.А.Полуэктов, Ю.М.Свирежев, В.А.Ратнер). Заседания школы приходили поочередно в Ленинграде и его окрестностях, и в Новосибирском Академгородке. В результате оперативного обмена информацией и обсуждения новых результатов была сформулирована концепция и, в последующие годы, развита современная динамическая теория биологических популяций и сообществ, которая к концу 80-х гг. оформилась в самостоятельную дисциплину. Отдельно следует сказать о популяционно-генетических моделях. В периодических встречах автора с А.А.Ляпуновым последний обращал внимание на то, что все разрабатываемые модели носят детерминированный характер. Вместе с тем гаметогенез, скрещивание и образование зигот, выживаемость и тех и других, несомненно, носят случайный характер. Это ставило под сомнение все полученные результаты. Замечания А.А.Ляпунова стимулировали постановку специального исследования, посвященного изучению влияния случайных событий на все этапы генетического оператора. Блестящую работу в этом направлении выполнил сотрудник Агрофизического института И.А.Вол (Вол, 1974). Это исследование явилось предметом его кандидатской диссертации и, к сожалению, до сих пор является недоступным большинству исследователей.

Помимо теории ставились и решались прикладные задачи – управление искусственным культивированием клеточных и микробных популяций, демографические проблемы роста народонаселения, распространение эпидемий, динамика вредителей сельского и лесного хозяйства, динамика лесных насаждений и сельскохозяйственных посевов. Следует сделать несколько замечаний о задачах управления популяциями и сообществами. Проблема заключается в том, что управления в данном случае (так же как и фазовые переменные) носят односторонний характер, что совершенно не свойственно классическим постановкам задач управления. Действительно, собирая урожай (укосы), мы можем только изымать биомассу и особей из популяции, но не добавлять их. Внося удобрения или поливая растения, мы только добавляем агент управления в систему, но не можем извлечь его обратно. Сами переменные (численность особей, их биомасса) также неотрицательны. Это вносит определенные трудности в решение задачи. Даже в линейном случае классические условия управляемости и наблюдаемости Р.Калмана являются только необходимыми, но недостаточными. Не с этим ли связаны все неудачи в создании эффективного метода борьбы с раковыми заболеваниями? Ряд необходимых и достаточных условий управляемости и наблюдаемости при культивировании популяций в хемостатах получил Б.Г.Заславский (Заславский и др., 1988). Однако, они не получили своего дальнейшего развития.

На западе формирование математической биологии во многом повторяло тот путь, который она прошла в нашей стране, иногда с опережением, а иногда с отставанием. В частности, результаты, опубликованные в нашей монографии «Динамическая теория биологических популяций» в 1976 году до сих пор известны только выходцам из нашей страны. Развитие исследований пошло несколько по иному пути. Была сформулирована концепция т.н. «метапопуляций». И хотя она появилась еще в 1970 году как “population of populations” (R. Levins), ее фактическое становление стало возможным лишь недавно.

Последнее обусловлено появлением возможности оперировать с пространственно распределенными данными, используя ГИС-технологи (L.Ginsburg). И еще одно направление современного этапа математической биологии - обоснование фундаментальных законов экологии на основе ее параллелей с классической ньютоновской механикой. Это направление, начатое Л.Р. Гинзбургом еще в годы его работы в АФИ в 70-е годы 20-го века, нашло в последнее время свою законченную формулировку (Ginzburg, 2004).

Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия Литература

Вол И.А. Стохастическая устойчивость и надежность биологических популяций: Дис. канд. техн. наук.– Л.:

АФИ, 1974.- 131с.

Заславский Б.Г., Полуэктов Р.А. Управление экологическими системами.– М.: Наука, 1988.– 295 с.

Ginzburg L., Colvan M. Ecological Orbits: How Planets Move and Populations Grow.– Oxford Univ. Press, 2004.- 166 p.

ПРОБЛЕМЫ СЕМАНТИЧЕСКОЙ МОДЕЛИ "УСТОЙЧИВОГО РАЗВИТИЯ" И

ПУТИ ИХ РЕШЕНИЯ

Пузаченко Ю.Г.

Институт проблем экологии и эволюции им. А.Н.Северцова РАН, Москва puzak@orc.ru, jpuzak@mail.ru Аннотация: формирование концепции «устойчивого развития» рассматривается с позиции максимума производства энтропии. Постулируется необходимость формирования нового метаязыка, обобщающего представления современной термодинамики и синергетики.

Само по себе возникновение, становление и развитие концепции «устойчивого развития» (Sustainable Development) можно рассматривать как пример сложного адаптивного поведения социума, ставящего перед собой естественную цель самосохранения на неопределенно большом интервале времени. Исторически ее возникновение было определено осознанием через слушания в разных странах комиссией Brundtland, начавшей работать в рамках мировой стратегии сохранения среды 1972-1980 г, неразрывности социума, экономики и среды. Широкая поддержка этой идеи международными институтами определила мощный поток как общих работ, направленных на разработку оснований концепции, так и частных, ориентированных на разработку методов достижения желаемой цели. Уровень шума в этом потоке информации, как и во всех сферах международной активности, очень высок и в подавляющем большинстве случаев поиск решений носит стохастический характер. Из этого потока можно выделить фактически одно конструктивное представление о капитале и его емкости и/или пропускной способности (carrying capacity).

Это понятие было распространено из финансовой сферы на промышленность, инфраструктуру, человека, природу. Устойчивое развитие должно строиться на основе приращения того, что необходимо для роста благосостояния, без деградации основного капитала. Понятие удобно, так как на основе простой логистической модели обеспечивает принципиальную измеримость емкости капитала и его приращения как показателя роста благосостояния. С другой стороны «емкость» должна подчиняться закону пропускной способности канала связи Шеннона и закону необходимого разнообразия Эшби. Если первый подход фрагментарно использовался в некоторых моделях новой экономики, то второй, информационный, в явном виде вообще не нашел какого-либо применения. При этом Cobb–Douglas функция, определяющая в частности мультипликативную форму вклада капиталов в конечную продукцию производственной деятельности, как и любые подобные аллометрические соотношения, есть прямой аналог пропускной способности канала связи.

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

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

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

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

Создается впечатление, что необходимо преобразовать аксиоматическую систему и перейти на более высокий уровень метаязыка, способного охватить физические, биологические, социально-экономические и познавательные формы действий, переходы глобальной системы и ее отдельных частей в различные фазовые состояния на общем тренде саморазвития. Очевидно, что это направление широко представлено в современной литературе. Общей его чертой является введение для макроскопического описания в термодинамическую модель информации, в ее соотношении с энтропией. Наиболее общую систему представлений энтропии-информации, по-видимому, дал А.М. Хазен (1998-2000). В его модели принцип максимума производства энтропии перестает быть ограничением для эволюции системы и, напротив, становится ее необходимым условием. Исходные положения можно свести к утверждению «все что наблюдаемо – то устойчиво в смысле Ляпунова на некотором интервале длительности». Естественной «целью», а точнее естественным поведением будет максимизация длительности – локальной стационарности с минимум производства энтропии.

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

Сама концепция «устойчивого развития» эволюционирует в соответствие с моделью:

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

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

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

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

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

ПРОСТРАНСТВЕННО-ВРЕМЕННАЯ ДИНАМИКА ТЕРМОДИНАМИЧЕСКИХ

ПЕРЕМЕННЫХ БИОСФЕРЫ ПО ИЗМЕРЕНИЯМ СО СПУТНИКА MODIS

Пузаченко Ю.Г.1, Сандлерский Р.Б.2 Институт проблем экологии и эволюции им. А.Н. Северцова РАН, Москва, Россия jpuzak@mail.ru, 2srobert_landy@mail.ru Аннотация: демонстрируются оценки термодинамических переменных, включая эксергию и приращение информации на основе спутника Modis Terra. Показана ведущая роль лесной растительности в функционировании биосферы.

Рассматриваются результаты анализа мультиспектральных измерений сезонного варьирования (два измерения за месяц) отраженной солнечной радиации спутником Modis Terra 2002 г. с разрешением 0.5о. Оценивается поглощенная радиация (R), информация Кульбака (I), энтропия (S), эксергия (Ex), по схеме предложенной Ергенсоном и Свирежевым (2004) и тепловой поток (T) от земной поверхности по данным проекта Land Surfase Temperature/Emissivity Modis Terra за тот же год. Внутренняя энергия оценивается через замыкание уравнения баланса.

Каналы спутника захватывают 16.58% энергии солнечной постоянной. В среднем для биосферы эффективность поглощения энергии R = 0.654, эксергии Ex = 0.437, S = 1.74 nit, I = 0.176 nit.

Результаты измерений позволяют выделить четыре фазовых состояния биосферы:

три равновесных состояния с информацией Кульбака меньше 0.14 nit и одно неравновесное.

Равновесные состояния: поверхность, покрытая снегом: I = 0.041 nit, минимум S = 1.07 nit, минимальные Ex = 6.10 вт/м2 и T = -9.77оC; аридные земли: I = 0.131 nit, S = 2.05 nit, Ex на среднем уровне – 90.52 вт/м2, максимальная T = 38.66оC; леса умеренных широт зимой: I =

0.038 nit, S = 2.10 nit, Ex = 18.57 вт/м2, T = -2.0330C. Тропические леса в течение всего года и леса умеренных широт в летний период неравновесны и имеют максимальные I = 0.498 nit и Ex = 137.11 вт/м2, низкую S = 1.54 nit и Т = 31.986оC. Области вечных снегов, пустыни и тропические леса в течении всего года находятся в одном фазовом состоянии. Сезонные изменения термодинамических переменных здесь незначительны и их изменения во времени равновесны. В первой области они связаны с перекристаллизацией снега, во второй с сезонным изменением увлажнения: чем больше увлажнение, тем больше эксергия и ниже энтропия. В лесах умеренных широт фазовый переход резкий. Таким образом, основную полезную работу в биосфере, выражающуюся в интенсификации круговорота влаги и биогенных элементов, осуществляют леса. Область их наибольшей активности совпадает с максимумом приращения информации – рисунок 1.

Материалы Второй конференции «Математическое моделирование в экологии»

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

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

Соответствующее приращение информации для всего массива данных составляет I = 0.083 nit и эффективности эксергии для пространственно-временного транспорта воздушных масс Ex= 0.315. Приращение информации для декабря I = 0.104 nit при Ex = 0.427 и для июля I =

0.439 nit при Ex= 0.631. Таким образом, в целом за год система почти равновесна и существенно неравновесна летом северного полушария, что определяется диспропорцией в площади суши и большим вкладом в полезную работу лесных экосистем умеренных широт.

Рисунок. 1. Приращение информации (светлые тона – максимум приращения информации) Эксергия нелинейно связана с солнечной постоянной (Eint): Ex/dEint =

0.0318exp(0.0036Eint) и в тропических широтах приращение эксергии на 1 вт/м2 солнечной постоянной составляет 4.358 вт/м2. Можно полагать, что увеличение солнечной постоянной будет увеличивать поток скрытого тепла из южного в северное полушарие и приводить к более теплым зимам с более ранней весной и общим потеплением климата. При этом тепловой поток в интенсивно испаряющих дождевых лесах может даже снижаться.

б.

а.

Рисунок 2. Соотношение альбедо и энтропии (а) – по Modis, (б) – по модели на основе данных ERBE (Stephens, O'Brien, 1993) Сопоставление результатов оценки потока энтропии по модели «тепловой машины» с оценками, полученными по спутнику, показывает их полное несоответствие.

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

Более сходные результаты дают измерения длинноволновой и коротковолновой радиации Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия осуществленные со спутников по проекту Earth Radiation Budget Experiment (ERBE) 1986 – 1990 (Stephens, O'Brien, 1993). На рисунке 2 приведено соотношение теоретической модели потока энтропии (б), полученной на основе данных ERBE, с результатами измерения по Modis (а). Если исключить влияние работающей растительности с резким увеличением информации, то соотношение между альбедо и энтропией в этих оценках подобно. Полезная работа лесной растительности, увеличивая неравновесность, снижает производство энтропии. Формально это указывает на то, что лесная растительность неравновесна и стационарна. В этом отношении особый интерес представляет лесная растительность умеренных широт, испытывающая резкий фазовый переход от равновесия с очень большим производством энтропии к стационарному состоянию с минимум энтропии. В приведенные результаты демонстрируют большие возможности мультиспектральной дистанционной информации для анализа термодинамики поверхности земли.

Работа выполнена при поддержке РФФИ № 09-05-00292-а.

Литература Jorgensen S.E., Svirezhev Y.M. Towards a Thermodynamic Theory for Ecological Systems.- Langford Lane Kidlington. Oxford. Elsevier, 2004.- 369 p.

Stephens G.L., O'Brien D.M. Entropy and climate. I: ERBE observations of the entropy production of the earth // Quarterly Journal of the Royal Meteorological Society.-Vol. 119, Issue 509, 1993.- P. 121–152.

МОДЕЛЬНЫЙ АНАЛИЗ ВЛИЯНИЯ ПОЛОВОЙ СТРУКТУРЫ НА ДИНАМИКУ

ЧИСЛЕННОСТИ ПОПУЛЯЦИИ

Ревуцкая О.Л., Фрисман Е.Я.

ИКАРП ДВО РАН, Биробиджан, Россия oksana-rev@mail.ru, frisman@mail.ru Аннотация: исследуется многокомпонентная модель динамики численности популяции, состоящей из новорожденных особей, неполовозрелых особей обоего пола, также самок и самцов, участвующих в размножении. В модели учитываются плотностно зависимые эффекты, действующие на выживаемость младшего возрастного класса. Проведено исследование структуры популяции в равновесном состоянии. Изучаются сценарии перехода к нелинейным режимам динамики.

–  –  –

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

Определены равновесные решения модели (2) и условия их существования.

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

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

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

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

Исследования проведены при частичной финансовой поддержке РФФИ (проекты № 09-04-00146-а, № 11-01-98512-р_восток_а) и ДВО РАН (проект № 11-III-В-01М-002).

КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ АГРОЛЕСОЛАНДШАФТОВ В

ГЕОИНФОРМАЦИОННОЙ СРЕДЕ

Рулев А.С., Юферев В.Г., Юферев М.В.

ГНУ ВНИАЛМИ Россельхозакадемии, г. Волгоград, Россия vnialmi_olp@vlpost.ru Аннотация: в статье рассмотрены вопросы моделирования в среде агролесомелиоративной геоинформационной системы (ГИС), которое обеспечивает решение научных и прикладных задач, связанных с инвентаризацией, анализом, прогнозированием деградационных процессов, а также территориальную организацию и управление ландшафтами.

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

Использование в исследованиях таких программных продуктов, как ENVI, MapInfo, Surfer, "Талка", OzExplorer и др. позволяет решать задачи, связанные с обработкой и анализом космоснимков, построением цифровых моделей рельефа (ЦМР) и ландшафта (ЦМЛ).

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

Большое внимание при моделировании в среде ГИС отводится геокодированию – привязке Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия объектов к цифровой карте территории, расположение которых в пространстве задается сведениями из таблиц баз данных. Эта информация может быть представлена координатами объектов, которые получены приемниками GPS или ГЛОНАСС, цифровыми адресами объектов и др.

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

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

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

Создание же современных ЦМР базируется на использовании матриц высот, которые доступны в геоинформационной среде в виде набора данных (GeoTIFF) радарного обследования поверхности земли зондoм "ASTER, установленным на спутнике TERRA", и зондами "SIR-C" и "X-SAR", установленными на борту "Space Shuttle Endeavour".

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

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

Рисунок 1. Цифровая 3D модель эрозионного ландшафта ключевого участка Материалы Второй конференции «Математическое моделирование в экологии»

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

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

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

Работа выполнена при поддержке РФФИ (проект № 10-05-97000) Литература Рулев А. С. Ландшафтно-географический подход в агролесомелиорации.– Волгоград: ВНИАЛМИ, 2007.– 160 с.

Кулик К.Н., Юферев В.Г. Компьютерное математико-картографическое моделирование агролесоландшафтов на основе аэрокосмической информации. // Доклады Российской академии сельскохозяйственных наук.– 2010.- №1.– С.

52-54.

Рулев А. С., Юферев В. Г. Картографо-геоинформационное моделирование в агролесомелиорации / А.С.

Рулев, В..Г. Юферев / Материалы Всероссийской конференции (с международным участием) «Математические модели и информационные технологии в сельскохозяйственной биологии: итоги и перспективы». 14–15 октября 2010 г., СанктПетербург. – СПб.: АФИ, 2010.- С. 68-71 (288) с.

ЭКОЛОГИЧЕСКАЯ МОДЕЛЬ ПРОДУКТИВНОСТИ СЕЛЬСКОХОЗЯЙСТВЕННЫХ

КУЛЬТУР

Рулева О.В.

ГОУ ВПО Волгоградский филиал РГТЭУ, Волгоград, РФ bifu@mail.ru Дается оценка биологической продуктивности агрофитоценозов как интегрального показателя деятельности полезащитных лесных насаждений. Предлагается методологический подход на основе, которого сформулирован закон о роли лесной полосы как экологического фактора, определяющий развитие агроценозов по экспоненциальной зависимости. Практическая значимость работы — прогноз динамики биопродуктивности на экосистемном уровне.

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

Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия На основе анализа собственных 20-летних исследований по данной проблеме, а также обобщении данных других исследователей в орошаемых агролесоландшафтах был предложен новый метод оценки динамики биопродуктивности агроценозов. (Рулева, 2006).

К растительному сообществу и агроценозам, как частному случаю фитоценозов при всей своей сложности, динамичности развития самого сообщества и его составляющих могут быть «приложимы» подходы и методы оптики мутных сред (Рачкулик, Ситникова, 1981). В статистическом смысле это означает переход к средним величинам, то есть замене истинных флуктуирующих величин статистически средними (Налимов, 1971). Принципиальным моментом использования теории поля для теоретического формализованного описания любых агроэкосистем является, прежде всего, выделение ряда системных параметров, отражающих доступность ресурсов, таких как влага, свет, элементы минерального питания и т.д. (все это присутствует при орошении и соблюдении культуры агротехники).

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

На основе проведенных исследований агроценозов, сформулировано следующее положение: независимо от географического положения района исследований, структуры агроценоза, рядности (их может быть 2-3-8 рядов), породного состава лесных полос, почвы, на которой возделываются культуры (красноземы, черноземы и др.) продуктивность сельскохозяйственных культур (зерновых, кормовых, овощных, масличных, технических и т.д.) при орошении изменяется по экспоненциальному закону (Таблица 1). Лесные полосы являются экологическим фактором, который способствует появлению адаптаций у растений агроценоза на межполосной клетке изменяться строго закономерно от расстояния до лесных полос, по одной зависимости — экспоненциальной: у = a e-bx, где у—изменение биопродуктивности сельскохозяйственных культур; х—расстояние до лесных полос, а и b коэффициенты, определяющие наклон и изгиб экспоненциальной кривой.

Таблица 1 — Изменение биопродуктивности сельскохозяйственных культур (y) от расстояния до лесных полос (x) по экспоненциальной зависимости у = a e-bx Название культуры Параметры уравнения Коэффициент детерминации R2 a b озимая пшеница 1,02 0,20 0,97 яровая пшеница 1,03 0,17 0,96 кукуруза 1,02 0,32 0,86 ячмень 1,10 0,63 0,88 подсолнечник 1,02 0,37 0,88 рис 1,09 0,36 0,96 нут 0,80 0,71 0,66 свекла 0,98 0,16 0,65 картофель 1,0 0,50 0,97 люцерна 1,0 0,88 0,93 хлопок 1,02 0,36 0,93 чайные кусты 1,0 0,28 0,82 Материалы Второй конференции «Математическое моделирование в экологии»

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

Для подтверждения достоверности обобщенных экспоненциальных зависимостей по сельскохозяйственным культурам был проведен сравнительный анализ по параметрам a и b вышеперечисленных культур. Стандартное отклонение во всех случаях, кроме случая с картофелем, ячменем (0,06), и подсолнечником (0,07), не превышает 0,5% уровень вероятности, что подтверждает высокую достоверность полученных зависимостей и применимость данного метода к описываемому явлению.

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

Литература Налимов В.В. Теория эксперимента.- М.: Наука, 1971.- 207 с.

Рачкулик В.И., Ситникова М.В. Отражательные свойства и состояние растительного покрова.– Л.:

Гидрометеоиздат, 1981.– 248 с.

Рулева О.В. Теоретическо-методологическая оценка влияния лесных полос на биологическую продуктивность сельскохозяйственных культур в орошаемых агролесоландшафтах // Вестник Саратовского госагроуниверситета им. Н.И. Вавилова.- 2006 г. № 3.-С.62-67.

РЕАКЦИЯ ЭКОСИСТЕМ БОЛЬШИХ СТРАТИФИЦИРОВАННЫХ ОЗЕР

НА СНИЖЕНИЕ АНТРОПОГЕННОЙ НАГРУЗКИ И УСТОЙЧИВОЕ РАЗВИТИЕ

Руховец Л.А.1, Петрова Н.А.2, Меншуткин В.В.1, Астраханцев Г.П.1, Минина Т.Р.1, Полосков В.Н.1 Учреждение Российской академии наук Санкт-Петербургский экономико-математический институт РАН, Санкт-Петербург, Россия emi@emi.nw.ru Учреждение Российской академии наук Институт озероведения РАН, Санкт-Петербург, Россия lake@spb.org.ru Аннотация: для Ладожского озера лимнологи установили, что снижение фосфорной нагрузки в период 1996-2005 гг. не привело к ожидаемому снижению продуктивности фитопланктона. С помощью математического моделирования в данной работе объяснен механизм этого явления.

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

Введение Ладожское озеро, как и американские Великие озера, подверглось процессу антропогенного эвтрофирования, причиной которого стал рост, начиная с 1962 года, фосфорной нагрузки. К 1982 году нагрузка выросла в три с половиной раза, достигнув величины 8110 т Pобщ/год. При этом озеро перешло из олиготрофного в развитое мезотрофное состояние. Принятые после 1982 года меры и разразившийся в 90-х годах экономический кризис привели к снижению фосфорной нагрузки. В результате процесс деградации озера замедлился. В работах (Петрова и др., 2005; Румянцев и Драбкова, 2006) весь процесс трансформации экосистемы Ладожского озера представлен как Работа выполнялась при финансовой поддержке РФФИ (грант № 10-06-00380а).

Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия последовательность этапов, в течение которых трансформация протекала по-разному. В Таблице из работы (Руховец и др., 2010) представлены в агрегированном виде характеристики этапов, основанные на данных, опубликованных сотрудниками Института озероведения РАН. Согласно данным наблюдений различия этапов связаны, прежде всего, с изменениями фосфорной нагрузки. Этап 1996-2005 гг. отличается тем, что среднегодовая нагрузка существенно ниже среднегодовой нагрузки в период 1975-1995 гг. (Таблица). В период 1996-2005 гг. такие показатели, как концентрации минерального и общего фосфора в воде озера приблизились к показателям, имевшим место в олиготрофный период развития озера до 1962 года. Однако состояние биоты, прежде всего, фитопланктона, во-первых, существенно отличается от ее состояния до 1962 года и, во-вторых, практически не отличается от состояния биоты в период 1984-1995 гг., когда фосфорная нагрузка была существенно выше.

Таблица 1. Средние за периоды фосфорная нагрузка (т Pобщ/год), концентрации общего /минерального фосфора (мкг P/л), биомасса суммарного фитопланктона (мг/л) (данные Института озероведения РАН из работ (Петрова и др.

, 2005; Румянцев и Драбкова, 2006)) Периоды 1959-1962 1975-1983 1984-1995 1996-2005 Фосфорная нагрузка (т Pобщ/год) Концентрации общего и 10 25 20 16 минерального фосфора в воде озера 3 11 8 4.5 (мкг P/л) Сырая биомасса суммарного

2.0 фитопланктона (мг/л) средняя за 0.8 1.6 1.5 летний период в приповерхностном 3.1 слое 0-2 м Одной из основных причин отмеченного явления в работе (Петрова и др., 2005) называется увеличение скорости внутриводоемного круговорота фосфора вследствие увеличения численности бактериопланктона и водных грибов. Подобное явление имеет место не только в Ладожском озере, но и в других великих озерах, например, в озере Эри (Conroy et al., 2005).

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

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

Модель сукцессии фитопланктона Ладожского озера Для воспроизведения функционирования экосистемы Ладожского озера использовалась модель сукцессии фитопланктона. В модели состояние экосистемы описывается трехмерными полями концентраций 14 субстанций: девяти комплексов фитопланктона (Ph1, …, Ph9), зоопланктона (Z), детритного фосфора (DP), растворенного в воде органического вещества (DOP), растворенного в воде кислорода (OX) (Rukhovets et al., 2003). Модель построена с использованием только одного биогена – фосфора, который в Ладожском озере является основным регулятором продуктивности экосистемы. Модель верифицирована и откалибрована на основе данных наблюдений, относящихся, в основном, к периоду 1984-1990 гг. (Петрова и др., 2005; Румянцев и Драбкова, 2006).

Материалы Второй конференции «Математическое моделирование в экологии»

ЭкоМатМод-2011, г. Пущино, Россия Применение модели Вычислительные эксперименты по воспроизведение круглогодичного функционирования экосистемы Ладожского озера с 1962 по 2005 г показали, что для периода 1962-1995гг. модель адекватно воспроизводит трансформацию его экосистемы. Однако, для периода 1996-2005гг. среднелетние значения биомассы суммарного фитопланктона, полученные с помощью модели, не соответствуют данным в Таблице.

Для того, чтобы проверить гипотезу лимнологов об ускорении круговорота фосфора в экосистеме Ладожского озера, авторы изменили в модели блок деструкции детрита и растворенного в воде органического вещества. Затем были выполнены основные вычислительные эксперименты по модели со «старым» и «новым» блоками деструкции при двух фиксированных нагрузках 2430 т Р/год и 4000 т Р/год. Счет проводился до получения периодических решений. Первая из нагрузок соответствует периоду олиготрофного статуса озера, вторая – это нагрузка, при которой в рамках модели со «старым» блоком деструкции Ладожское озеро остается в слабо мезотрофном состоянии.



Pages:     | 1 |   ...   | 3 | 4 || 6 | 7 |



Похожие работы:

«Министерство образования и науки Российской Федерации Федеральное государственное автономное образовательное учреждение высшего образования "НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ТОМСКИЙ ПОЛИТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ"_ Институт неразрушающего контроля Направление подготовки: 20.04.01 "Техносферная безопасность" Кафедра экологии...»

«ПОЧИНКОВСКИЙ МУНИЦИПАЛЬНЫЙ РАЙОН МУНИЦИПАЛЬНОЕ БЮДЖЕТНОЕ ОБЩЕОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ РИЗОВАТОВСКАЯ СРЕДНЯЯ ШКОЛА ПРИНЯТО УТВЕРЖДЕНО на заседании педагогического совета приказом директора Протокол от 29.08.201...»

«НАЦИОНАЛЬНАЯ АКАДЕМИЯ НАУК АРМЕНИИ NATIONAL ACADEMY OF SCIENCES OF ARMENIA ДОКЛАДЫ REPORTS №2 Том Volume БИОФИЗИКА УДК 577.391;547.963.3 Ц. М. Авакян, Н. И. Мкртчян, Н. В. Симонян,...»

«Министерство образования и науки Российской Федерации Федеральное государственное бюджетное образовательное учреждение высшего образования "Нижневартовский государственный университет" Факультет экологии и инжиниринга Рабочая...»

«Приложение 3 к приказу Минздрава Кыргызской Республики от 28.04.2015 года №212 ИНСТРУКЦИЯ О МЕТОДАХ МИКРОБИОЛОГИЧЕСКОЙ ДИАГНОСТИКИ МЕНИНГОКОККОВОЙ ИНФЕКЦИИ И БАКТЕРИАЛЬНЫХ МЕНИНГИТОВ I. МИКРОБИОЛОГИЧЕСКАЯ ДИАГНОСТИКА МЕНИНГОК...»

«Matrix– руководство администратора Версия 1.0 Номер документа: GradSoft-Matrix-AG-03.03.2005-1 Содержание Введение Основные функциональные возможности Интеграция Matrix с системами сторонних производителей Интеграция Matrix с АСР “Град” Процесс установки. Требования к окружению Запуск инсталлятора Выбор па...»

«Экологическое занятие Подготовила и провела воспитатель высшей категории Васильева С.А "Лесное царство" Цель: закрепить знания о лесе, как о экосистеме.Задачи: — дать детям представление о том, что лес – это живой орг...»

«Раздел 5. "Химические технологии. Безопасность жизнедеятельности" УДК 734.35 ГУТОРКА А.Д., 1НУГАИЕВА Е.В., 1ЧЕРНЫШЕВА А.А., 1СОЛУЯНОВА Ю.М. (Карагандинский государственный индустриальный университет, г. Темиртау, Казахстан) ЭКОЛОГИЧЕСКИЕ АСПЕКТЫ ВОЗДЕЙСТВИ...»

«Образовательное учреждение высшего образования Тверской институт экологии и права Кафедра Гражданскоправовых дисциплин РАБОЧАЯ ПРОГРАММА УЧЕБНОЙ ДИСЦИПЛИНЫ (МОДУЛЯ) ПРАВО СОЦИАЛЬНОГО ОБЕСПЕЧЕНИЯ Направление подготовки 030900.62 "Юриспруденция" Профиль подготовки "общий" Квалификация (степень) выпускника Бакалавр Т...»

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

«Изумрудная сеть: инструмент охраны естественной среды обитания в Европе Европа и биологическое разнообразие Во всём мире продолжается сокращение биологического разнообразия. Фрагментация местообитаний, загрязнение, чрезмерная эксплуатация территорий и создание искусственны...»

«Лизунов Антон Юрьевич УЧЕТ ВНУТРИЛИГАНДНЫХ ВЗАИМОДЕЙСТВИЙ ПРИ ДОКИНГЕ С ОЦЕНОЧНОЙ ФУНКЦИЕЙ НА ОСНОВЕ УСРЕДНЕННЫХ ПОТЕНЦИАЛОВ МЕЖАТОМНОГО ВЗАИМОДЕЙСТВИЯ Специальность 03.01.02 — "Биофизика" Автореферат диссертации на соискание учёной степени кандидата физи...»

«2011 БЕЛКИ ПЕПТИДЫ РОССИЙСКАЯ АКАДЕМИЯ НАУК Научный совет по биоорганической химии Учреждения Российской академии наук: Институт биологии Карельского научного центра РАН Институт биоорганической химии им. академиков М.М. Шемякина и Ю.А. Овчинникова РАН Петрозаводский госуда...»

«ISSN 2308-6874 Научно-издательский центр Априори ТЕОРИЯ И ПРАКТИКА АКТУАЛЬНЫХ ИССЛЕДОВАНИЙ Материалы IX Международной научно-практической конференции (24 июня 2015 г.) Сборник научных трудов Краснодар УДК 082 ББК 72я431 Т 11 Редакционная коллегия: Бисалиев Р...»

«Протокол № 2015-ВСМН-19/Д от_03.03.2015 стр. 1 из 6 УТВЕРЖДАЮ Председатель конкурсной комиссии _ С.В. Яковлев "03" марта 2015 года ПРОТОКОЛ № 2015-ВСМН-19/Д заседания конкурсной комиссии ОАО "АК "Тран...»

«Ученые записки Крымского федерального университета имени В. И. Вернадского Серия "Биология, химия". Том 1 (67). 2015. № 3. С. 49–55. УДК 579.64:581.14(653.63) ВЛИЯНИЕ МИКРОБНОГО ПРЕПАРАТА "ЭМБИКО" НА ПРОРАСТАНИЕ СЕМЯН ОГУРЦОВ (CUCUMIS SATIVUS L.) СОРТОВ КОНКУРЕНТ И ФЕНИКС Татаренко Я. И., Отурина И. П., Тепл...»

«ПРАВИЛА ПРОЖИВАНИЯ В КВАРТИРЕ БЛОКИРОВАННОГО ЖИЛОГО ДОМА В ЖИЛОМ КОМПЛЕКСЕ "ПАВЛОВО-2"1. СФЕРА ДЕЙСТВИЯ ПРАВИЛ ПРОЖИВАНИЯ. 1.1. Пользование Квартирой в блокированном жилом доме, приква...»

«Экосистемы. 2016. Вып. 6. С. 100–106 СЕЛЕКЦИЯ И ДЕКОРАТИВНОЕ РАСТЕНИЕВОДСТВО УДК: 582.973:631.526.32 ХОЗЯЙСТВЕННО-БИОЛОГИЧЕСКАЯ ХАРАКТЕРИСТИКА ДЕКОРАТИВНЫХ СОРТОВ И ФОРМ ЖИМОЛОСТИ (LONICERA L.) В РОССИИ Сорокопудов В. Н.1, Куклина А. Г.2 Всероссийский селекционно-технологический институ...»

«СОВЕТ ДЕПУТАТОВ ЗДВИНСКОГО РАЙОНА НОВОСИБИРСКОЙ ОБЛАСТИ второго созыва РЕШЕНИЕ шестой сессии 17.12.2010 года с. Здвинск №3 Об утверждении актуализированной Комплексной программы социаль...»

«РАСПРОСТРАНЕНИЕ ЖИВОРОДЯЩЕЙ ЯЩЕРИЦЫ НА ЮГЕ АРЕАЛА В ПОВОЛЖЬЕ Г.В. Епланова Институт экологии Волжского бассейна РАН, г. Тольятти eplanova_ievb@mail.ru Живородящая ящерица Zootoca vivipara (Lichtenstein, 1823) – вид с обширным транспалеарктическим ареалом. Она...»

«79014_729310 ВЕРХОВНЫЙ СУД РОССИЙСКОЙ ФЕДЕРАЦИИ № 302-ЭС15-12604 ОПРЕДЕЛЕНИЕ г. Москва 21.10.2015 Судья Верховного Суда Российской Федерации Чучунова Н.С., рассмотрев жалобы (заявления) общества с ограниченной ответственностью "Тайга" и Агентства лесного хозяйства Ир...»

«4 истории в картинках Я сделал свою первую мультимедийную презентацию в 17 лет. Она предназначалась для большого проекта на школьном уроке биологии и была посвящена загрязнению окружающей среды. Мое слайд-шоу должно было показать всему...»








 
2017 www.kn.lib-i.ru - «Бесплатная электронная библиотека - различные ресурсы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.