ОПРЕДЕЛЕНИЕ НАИБОЛЕЕ ИНФОРМАТИВНЫХ ПАРАМЕТРОВ ВОДНЫХ ОБЪЕКТОВ ДЛЯ ИНТЕГРАЛЬНОЙ ОЦЕНКИ ИХ
ЭКОЛОГИЧЕСКОГО СТАТУСА НА ОСНОВАНИИ
ИМИТАЦИОННОГО МОДЕЛИРОВАНИЯ
ВВЕДЕНИЕ 6
1 ОБЗОР ЛИТЕРАТУРЫ 8
1.1 Государственный мониторинг водных объектов 9
1.1.1 Задачи и нормативно-правовая база мониторинга водных объектов 9
1.1.2 Информативность показателей интегральной оценки состояния
водных объектов 11
1.2 Примеры математических моделей водных объектов 14
1.2.1 Математические модели озёр 14
1.2.2 Математические модели проточных водоёмов 19
1.2.3 Математические модели нагрузки на водосбор 20
1.2.4 Модели экономической оценки ассимиляционного потенциала 22
2 ОБЪЕКТ И МЕТОДЫ ИССЛЕДОВАНИЯ 23
2.1 Модель «ВОДЭКО» 23
2.2 Информационная база для модели 25
2.3 Критерий Нэша-Сатклиффа 27
3 РЕЗУЛЬТАТЫ 31
3.1 Влияние «олиготрофного» и «эвтрофного» уровней содержания азота и
фосфора в речном стоке на поведение модели 31
3.2 Влияние начальных условий на поведение модели 34
3.3 Влияние залпового поступления биогенов 37
3.4 Моделирование воздействия токсиканта 39
3.4.1 Сравнение воздействия разных типов токсиканта 39
3.4.2 Влияние поступления токсикантов в разных концентрациях 48
3.5 Оценка чувствительности отдельных параметров 51
3.5.1 Оценка чувствительности параметров в олиготрофном и
эвтрофном экспериментах 52
3.5.2 Оценка чувствительности параметров в экспериментах с
поступлением токсиканта 54
ВЫВОДЫ 61
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ 63
Актуальность работы. Методы исследований геоэкологии как междисциплинарного научного направления основаны на системном подходе, т.е. на широком применении различных способов моделирования, включая математическое имитационное моделирование. Системы окружающей среды рассматриваются как кибернетические системы, в которых происходит обмен веществом, энергией и информацией между самими блоками этих систем, а также с внешними по отношению к ним системами.
Системный анализ - это совокупность специфических методов и средств, используемых для подготовки и обоснования решений по сложным проблемам. Он предполагает формализацию современных представлений о происходящих в экосистемах процессах в формы математических зависимостей и уравнений, и их исследование. Основа системного анализа - общая теория систем и системный подход. Последний представляет собой методологическое направление в науке, отвечающее за разработку средств и методов исследования систем как сложноорганизованных объектов.
Имитационное моделирование обеспечивает возможность исследования интенсивностей массобмена между компонентами модели экосистемы, что невозможно установить непосредственным наблюдением за реальной экосистемой. Актуальность работы объясняется необходимостью определения параметров водных экосистем, наиболее информативных для оценки их состояния.
Цель заключается в определении на основании результатов численных экспериментов с имитационной моделью экосистемы проточного водоёма комплекса наиболее информативных параметров для многокритериальной оценки экологического состояния водных объектов. При этом эти параметры должны быть достаточно просто определяемы при мониторинге.
Для достижения поставленной цели были решены следующие задачи. Выполнялся анализ показателей, характеризующих экологическое состояние водных объектов, и определяемых при государственном мониторинге поверхностных вод, и оценена их информативность.
Изучались современные программы мониторинга водных объектов и сгруппированы по их отличительным признакам. Рассматривались параметры, применяемые для оценки качества водных объектов и их экологического состояния.
Выполнялся обзор современных имитационных моделей функционирования экосистем водных объектов.
По результатам численных экспериментов с имитационной моделью функционирования экосистемы выполнялся анализ информативности параметров функционирования водной экосистемы для оценки её экологического состояния.
Объект исследования. В данной работе использованы результаты численных экспериментов с имитационной моделью функционирования экосистемы проточного водоёма.
Обработка численных экспериментов с компьютерной моделью «Программа для моделирования антропогенного эвтрофирования и отравления токсикантами водных экосистем (EVTOX)» (автор - В.Ю. Третьяков, свидетельство о государственной регистрации программы для ЭВМ № 2018661228, дата регистрации в Реестре программ для ЭВМ 04.09.2018). Исследования степени различий между результатами разных численных экспериментов осуществлены с помощью расчётов значений модифицированного В.Ю. Третьяковым критерия Нэша-Сатклиффа (А.Н. Гельфан, 2017; С. А. Седова, В. В. Дмитриев, В. Ю. Третьяков [и др.], 2021). Расчёты модифицированного критерия Нэша-Сатклиффа выполнены с помощью разработанной В.Ю. Третьяковым компьютерной программы «Программа для сравнения результатов двух модельных расчётов с помощью 7
модификации критерия Нэша-Сатклиффа (МК Нэша-Сатклиффа)», свидетельство о государственной регистрации программы для ЭВМ № 2020612788, дата регистрации в Реестре программ для ЭВМ 03.03.2020.
1. При имитации изменчивости содержания субстанций в поступающем речном стоке для «олиготрофного» и «эвтрофного» экспериментов, соответствующей их внутригодовых динамикам в реках 1-ой и 4-й групп, показали колоссальные значения модифицированного критерия Нэша-Сатклиффа, рассчитанного сравнением динамик биогенных элементов. Это свидетельствует о большой чувствительности модели экосистемы к антропогенному эвтрофированию.
В то же время в экспериментах с имитацией динамики биогенов и органики в речном стоке в соответствии с их временной изменчивостью в реках 3-ей группы обнаружено незначительное различие динамик биомасс зоопланктона.
2. В экспериментах с «олиготрофным» и «эвтрофным» уровнями поступления азота и фосфора с речным стоком сильнее всего изменение начальных значений компонентов модели влияет на динамики биомасс рыб-прошлолеток, хищных рыб-сеголеток и бентофагов, олигохет и моллюсков, а также концентраций биогенов, растворённого кислорода, углекислого газа и органических углерода и азота. В «олиготрофных» экспериментах существенны различия динамик биомасс хищного зоопланктона.
3. При имитации залпового поступления азота и фосфора существенные изменения динамик наблюдались у таких компонентов модели, как биомасса рыб-прошлолеток, олигохет, а также концентрации растворённого углекислого газа.
4. В численных экспериментах с разными типами токсикантов как с имитацией их постоянного поступления, так и залпового значительные величины модифицированного критерия Нэша-Сатклиффа получены при сравнении динамик интенсивностей поступления токсиканта с пищей для бентосоядных рыб-прошлолеток и его поступления в диатомовые и синезелёные водоросли при биосинтезе, а также с водой для хищных прошлолеток и сеголеток и для моллюсков. У остальных динамик поступления токсиканта в организмы не обнаружено существенных различий.
5. При имитации поступления токсикантов в разных концентрациях в случае токсиканта 1-го типа наблюдается значительное различие динамик содержания токсиканта в органическом веществе взмученных донных илов; в случае токсиканта 2-го типа - значительное различие динамик содержания токсиканта в фитопланктоне и во взмученном органическом иле, а также значительный отклик даёт изменение динамик при биосинтезе фитопланктона и с поступлением токсиканта с пищей для растительноядных прошлолеток. Однако динамика содержания токсиканта в зоопланктоне, ихтиофауне, моллюсках и олигохетах практически не изменяется при значительном повышении концентрации токсиканта в речных водах.
6. Расчёты чувствительности отдельных параметров к антропогенному эвтрофировнию показывают значимый отклик на это воздействие динамик концентраций органического азота, а также азота в аммонийной, нитритной и нитратной формах. Чувствительность динамики концентрации органического фосфора существенно ниже. Чувствительность динамики биомассы олигохет к имитации антропогенного эвтрофирования значима только в экспериментах с имитацией внутригодовой динамики субстанций, характерной для рек 3-ей и 4-ой групп.
7. Расчёты чувствительности отдельных параметров к поступлению токсикантов 1-го или 2-го типов показывают значимый отклик на это воздействие динамик биомасс растительноядных рыб сеголеток, бактериопланктона и детрита, а для вегетационного периода - динамики биомассы растительноядных рыб сеголеток и бактериопланктона.
1. Айзель, Г.В. Расчеты речного стока на основе модели SWAP для водосборов с недостаточным информационным обеспечением. 2. Использование методов физико-географического подобия и пространственной геостатистики / Г.В. Айзель, Е.М. Гусев, О.Н. Насонова // ВОДНЫЕ РЕСУРСЫ, 2017. - Т. 44. - № 4. - С.419-431.
2. Астраханцев, Г.П. Использование моделей экосистем больших озер для получения оценок ассимиляционного потенциала / Г.П. Астраханцев, В.В. Меншуткин, Т.Р. Минина // УБС, 2015. — № 55. — С. 17-34.
3. Брюханов, А.Ю. Методика определения биогенной нагрузки сельскохозяйственного производства на водные объекты / А.Ю. Брюханов, С.А. Кондратьев, Н.С. Обломкова, А.С. Оглуздин, И.А. Субботин // Технологии и технические средства механизированного производства продукции растениеводства и животноводства, 2016. — № 89. — С. 175-182.
4. Владимирова, О.М. Модельные оценки составляющих баланса азота и фосфора в экосистеме Финского залива / О.М. Владимирская, Т.Р. Еремина, А.В. Исаев, В.А. Рябченко, О.П. Савчук // Ученые записки РГГМУ, 2018. — № 53. — С. 72-82.
5. Гельфан, А.Н. О проблеме валидации гидрологической модели для диагностических задач // Водные ресурсы: новые вызовы и пути решения, 2017. — С. 143-149.
6. Дмитриев, В.В. Моделирование влияния факторов среды на скорости массообмена в водной экосистеме / В.В. Дмитриев, Ю.Н. Сергеев, В.П. Кулеш, В.Ю. Третьяков // Труды II Всероссийской конференции «Гидрометеорология и экология: достижения и перспективы развития». - СПб.: ХИМИЗДАТ, 2018. - 753 с. С. 228-232.
7. Дмитриев, В.В. Мониторинг и моделирование Продукционно-деструкционных отношений в водных экосистемах / В.В. Дмитриев, О.Н. Боброва, И.В. Грачева, П.А. Колодкин, Е.А. Примак, С.А. Седова, А.А. Четверова // УСПЕХИ СОВРЕМЕННОГО ЕСТЕСТВОЗНАНИЯ, 2019. - № 1. - С. 82-87
8. Дмитриев, В.В. Оценка эмерджентных свойств водных объектов: трофический статус, устойчивость, экологическое благополучие /В.В. Дмитриев, Г.В. Пряхина, А.Н. Огурцов, Е.А. Примак, Д.Р. Амаро Медина // ТРЕТЬИ ВИНОГРАДОВСКИЕ ЧТЕНИЯ. ГРАНИ ГИДРОЛОГИИ, 2018. - С. 347-354.
9. Дмитриев, В.В. Подходы к интегральной оценке и ГИС-картографированию устойчивости и экологического благополучия геосистем. Часть IV. Интегральная оценка экологического благополучия наземных и водных геосистем /В.В. Дмитриев, И.В. Федорова, А.С. Бирюкова. // Вестник Санкт-Петербургского университета, 2016. — № 2. - С. 37-53.
10. Дмитриев, В.В. Региональные особенности влияния факторов среды на продуцирование органического вещества в водных экосистемах /В.В. Дмитриев, Д.Р. Амаро Медина, А.А. Булатова, М.С. Шабан, А.С. Саганович, А.В. Валичев, А.М. Ильев, Ю.А. Крамчанина // Международный журнал прикладных и фундаментальных исследований, 2016. — № 9-2. — С. 205-209.
11. Дмитриев, В.В. Экологический статус водоема: аксиология, аксиометрия, оценка / В.В. Дмитриев, А.А. Четверова, А.Н. Огурцов, Д.Р. Амаро Медина, С.А. Седова, Е.Н. Вакула, М.Р. Кузнецова, А.Д. Колупаева // ОЗЕРА ЕВРАЗИИ: ПРОБЛЕМЫ И ПУТИ ИХ РЕШЕНИЯ, 2019. — С. 66-70.
12. Дрововозова, Т.И. Интегральный показатель качества сточных вод, отводимых в водный объект / Т.И. Дрововозова, Н.Н. Паненко, А.В. Лещенко // ИВД, 2019. - №3 (54). - С. 27.
13. Клейнер, Г.Б. Производственные функции : Теория, методы, применение / Г. Б. Клейнер. — М. : Финансы и статистика, 1986. — 238 С.
14. Кондратьев, С.А. Биогенная нагрузка на Онежское озеро от рассеянных источников по результатам математического моделирования / С.А. Кондратьев, М.В. Шмакова, Е.Г. Маркова, П.А. Лозовик, А.Ю. Брюханов, Е.Ф. Чичкова. — Из-во РГО, 2016. —Т. 148. — № 5. — С. 53-64.
15. Кондратьев, С.А. Внешняя и внутренняя фосфорная нагрузка на водоем (на примере водохранилища Сестрорецкий разлив) / С.А. Кондратьев, Н.В. Игнатьева, С.Г. Каретникова // Региональная экология, 2016. — № 4 (46), — С. 59-70...37