Необходимость исследования скоростей массообмена между компонентами водных экосистем с учетом антропогенных воздействий, изменяющих естественный ход обменных процессов, выделение антропогенных эффектов на фоне естественных изменений скоростей процессов и исследование ответных реакций экосистем на внешние воздействия, которые не поддаются непосредственному наблюдению, обуславливают необходимость развития новых подходов к изучению водных экосистем. Существующие в настоящее время когнитивные и пространственно-однородные, боксовые (блочные, резервуарные), непрерывные имитационные модели позволяют решать широкий спектр задач: исследовать влияния химического и биологического состава и физических свойств среды на первичную продуктивность экосистем [1, 2], её изменение при различных сценариях формирования внешней и внутренней биогенной нагрузки [3], изучать влияние антропогенного эвтрофирования и воздействие токсикантов на водные экосистемы [4], воспроизводить на моделях биогеохимические циклы химических элементов [5] и др.
Развитие методов имитационного моделирования водных экосистем [4] и методов построения композитных индексов химического и биологического состава; простых и сложных (присущих системе в целом) свойств водных объектов [6], определило возможность сочетания этих подходов для решения ряда научных задач.
В первую очередь к таким задачам можно отнести количественную оценку результата воздействия на водную экосистему и оценку допустимости/недопустимости для экосистемы рассмотренного воздействия или действия многих факторов и их сочетаний [7]. При этом допустимость воздействия оценивается не по воздействию на некоторый компонент или «слабое звено» водной экосистемы, а на основе интегральных показателей состояния системы в целом до и после воздействия. Если система после оказанного воздействия по совокупности параметров и иерархий переходит в более старший класс состояния (качества и токсического загрязнения, трофности, трофосапробности, трофотоксобности и др.), то такое воздействие может рассматриваться как недопустимое при условии выявления его отрицательной значимости на человека и/или на биоту водоемов. Выявление такой возможности и критических концентраций компонентов реализуется в экспериментах с имитационной моделью, а подтверждение перехода системы в новый класс состояния (трофности, качества, токсического загрязнения, экологического благополучия/неблагополучия и др.) подтверждается на основе использования композитных индексов и субиндексов. В сложных случаях, когда временное изменение параметров состояния системы охватывает два и более класса, можно отдельно рассматривать сезонную изменчивость параметров или выполнять оценку только по отдельным субиндексам, а не по вектору состояния системы в целом. Антропоцентризм и биоцентризм исследования обусловливают выбор параметров и приоритетов в задании весовых коэффициентов учитываемых параметров при построении интегральных показателей. Возможно сочетание подходов или введение еще одного уровня свертки, агрегирующего оба подхода при равновесомом их учете или с выделением основного приоритета в прогнозировании сбалансированного, бескризисного (устойчивого) развития.
Во вторую очередь в перечень исследовательских задач можно включить дополнительные возможности проверки адекватности результатов имитационного моделирования на основе построения композитных индексов по результатам натурных наблюдений и имитационного моделирования. Для этого можно рассчитывать интегральные показатели для полного набора компонентов имитационной модели водной экосистемы, или формировать субиндексы оцениваемых свойств водного объекта. В перечень субиндексов целесообразно включить рекомендуемые нами [7–9]: интегральный показатель трофности, интегральный показатель качества и токсического загрязнения воды (и грунтов, если в модели имеется блок донных отложений), интегральные показатели устойчивости (уязвимости, самоочищения) водоема и др. Для более сложных ситуаций (оценка геоэкологического статуса), при моделировании систем «водоем + водосбор» или «водоток + водосбор» необходим учет экологического статуса ландшафтов водосборной территории (и/или локальных береговых экосистем). Для этого необходимо рассчитывать субиндексы интегративных свойств для доминирующих ландшафтов (береговых экосистем).
Степень адекватности модели реальному пространственно-временному изменению компонентов водной экосистемы в этом случае можно характеризовать на основе традиционных для классической статистики индексов Нэша – Сатклифа и его модификаций [10, 11], индекса Тейла (Тила), средней относительной ошибки прогноза, индексам S/σ и ρ [12]. Для этих целей можно также актуализировать применение классических статистических критериев выявления аномальных значений рядов по критериям Ирвина, Смирнова – Граббса, Титьена – Мура, Томпсона (Рошера), Хоглина – Иглевича, Шовене и др. [13, 14]. Выявление адекватности моделей реализуется на основе сочетания нескольких подходов: 1 – сравнения интегральных показателей компонентного состава всех моделируемых и наблюдаемых компонентов. Близость наблюденных и полученных по модели компонентов водной экосистемы устанавливается на основе сравнения рядов по традиционно применяемым индексам; 2 – сравнения интегральных показателей (индексов) и субиндексов неаддитивных свойств исследуемой системы (продуктивности, качества и токсического загрязнения воды и др.) по результатам наблюдений и модельных сценариев. В этом случае изменение класса продуктивности, качества, благополучия системы и др., выявленное на основе применения композитных индексов, потребует акцентированного подтверждения на основе натурных данных.
В третью очередь на основе сравнения модельных сценариев до и после воздействия, применяя классические индексы, упомянутые выше, или их модификации, можно оценить, насколько сильным оказалось воздействие для экосистемы в целом, а также сможет ли система вернуться в прежнее состояние или воздействие оказалось необратимым. На первом этапе важным здесь является установление фактической возможности перехода системы в другой класс по величине композитного показателя или того, что оказанное воздействие на экосистему способно при определенных условиях вывести систему за пределы того класса, в котором она находилась до воздействия.
Целью исследования является доказательство возможности использования композитных индексов для: 1 – проверки адекватности результатов системного моделирования на основе построения композитных индексов по результатам натурных наблюдений и имитационного моделирования; 2 – оценки воздействия на водный объект по результатам имитационного моделирования; 3 – выявление способности системы после оказанного воздействия сохранить класс состояния или интегративного свойства, в котором она пребывала до воздействия, или перейти в другой класс (временно или безвозвратно) на основе статистических критериев оценки близости рядов. Задачи экспериментальных исследований приведены в таблице.
Перечень экспериментов с моделями водных экосистем и их специфика
|
Эксперименты для оценки воздействия на водную экосистему |
М1 СNPXО-модель |
М2 Модель ВОДЭКО |
||||||||||||||||||||||||
|
1 |
2 |
3 |
||||||||||||||||||||||||
|
Э1. Оценка допустимого воздействия на основе изменения системой способности продуцировать органическое вещество |
Задачи: 1 – оценить использование выбранных имитационных моделей для оценки допустимости воздействия на экосистему; 2 – на основе построения композитного индекса по результатам моделирования воздействий выявить изменение способности системы продуцировать органическое вещество; 3 – на основе традиционных статистических индексов выявить изменение компонентов систем и системных свойств после оказанного воздействия |
|||||||||||||||||||||||||
|
Содержание эксперимента Э1М1 |
Содержание эксперимента Э1М2 |
|||||||||||||||||||||||||
|
Сравнивались результаты двух модельных экспериментов с разным фоном мутности воды (2 и 40 мг/л) на годовом этапе развития экосистемы [15]. На основе параметров, используемых для оценки трофического статуса водоема, рассчитывался композитный индекс продуктивности (ИПТ) и выявлялось изменение класса трофности после оказанного воздействия. Для построения композитного индекса продуктивности использовались 6 компонентов: 1. Биомасса фитопланктона, мг/л. 2. Отношение биомассы зоопланктона к биомассе фитопланктона. 3. Концентрация общего фосфора, Робщ, мг/л. 4. Концентрация общего азота, Nобщ, мг/л. 5. Отношение концентраций N:P. 6. Концентрация растворенного кислорода, % насыщения. Индекс ИПТ = ИПТ(q, w) = = ИПТ(q1,…,q6; w1,…,w6) = где qi – нормированное значение характеристики, wi – вес параметра в сводной оценке. Для wi использовалось: 0,0 ≤ wi ≤ 1,0 и Шкала ИПТ получена для 4 классов трофности при равновесомости выбранных индикаторов: О – олиготрофия; М – мезотрофия; Э – эвтрофия; Г – гипер(эв)трофия.
(1) |
Сравнивались результаты модельных экспериментов на годовом этапе развития транзитной экосистемы по изменению продукционной способности речных вод при постоянном содержании условного токсиканта в воде. Внутригодовой ход содержания токсиканта задан на основе анализа внутригодовых динамик различных токсических элементов в речных водах [4]. При этом в модели учитывались типы воздействия токсиканта: 1 тип – стимуляция первичной продукции низкими концентрациями токсиканта или 2 тип – ингибирование продукционных возможностей высокими концентрациями токсиканта. На основе параметров оценки трофического статуса водоема рассчитывался субиндекс продуктивности (ИПТ) и выявлялось изменение класса трофности после оказанного воздействия. Для построения композитного индекса продуктивности водного объекта использовались 9 компонентов: 1. Биомасса фитопланктона, мг/л. 2. Отношение биомассы зоопланктона к биомассе фитопланктона. 3. Концентрация общего фосфора, Робщ, мг/л. 4. Концентрация общего азота, Nобщ, мг/л. 5. Отношение концентраций N:Р. 6. Чистая первичная продукция, мгС/м2·сут. 7. Продукция фитопланктона, мгС/л·сут. 8. Скорость фотосинтеза, мгО/л·сут. 9. БПК5, мгО2/л. Индекс ИПТ = ИПТ(q, w) = = ИПТ(q1,…,q9; w1,…,w9) = Шкала ИПТ получена для 4 классов трофности при равновесомости выбранных индикаторов: О – олиготрофия; М – мезотрофия; Э – эвтрофия; Г – гипер(эв)трофия.
(2) |
|||||||||||||||||||||||||
|
Продолжение таблицы |
||||||||||||||||||||||||||
|
1 |
2 |
3 |
||||||||||||||||||||||||
|
Э2. Выявление сезонных изменений способности системы продуцировать живое вещество после воздействия |
Задачи: 1 – оценить возможность использования имитационных моделей для оценки допустимости воздействия на экосистему; 2 – на основе построения композитного индекса по результатам моделирования воздействия выявить особенности сезонного изменения продуцирования органического вещества в системе |
|||||||||||||||||||||||||
|
Содержание эксперимента Э2М1 |
Содержание эксперимента Э2М2 |
|||||||||||||||||||||||||
|
Рассматривались сезонные изменения способности системы продуцировать живое органическое вещество по результатам двух модельных экспериментов с разным фоном мутности воды (2 и 40 мг/л). На основе параметров, используемых для оценки трофического статуса водоема, рассчитывался композитный индекс продуктивности (трофности) и выявлялось изменение класса трофности после оказанного воздействия. Параметры и шкала для построения композитного индекса продуктивности водного объекта выбраны аналогично эксперименту Э1М1 |
Рассматривались сезонные изменения способности системы продуцировать живое органическое вещество по результатам модельных экспериментов годового цикла развития экосистемы с учетом типа воздействия токсиканта на первичные продуценты (1 тип – стимуляция первичной продукции низкими концентрациями токсиканта или 2 тип – ингибирование продукционных возможностей высокими концентрациями токсиканта). На основе параметров, используемых для оценки трофического статуса водоема, рассчитывался композитный индекс продуктивности (трофности) и выявлялось изменение класса трофности после оказанного воздействия. Параметры и шкала для построения композитного индекса продуктивности водного объекта выбраны аналогично эксперименту Э1М2 |
|||||||||||||||||||||||||
|
Э3. Оценка допустимого воздействия на систему на основе изменения качества и токсического загрязнения воды |
Задачи: 1 – оценить возможность использования имитационных моделей для оценки допустимости воздействия на экосистему; 2 – на основе построения композитного индекса качества и токсического загрязнения по результатам моделирования воздействия выявить изменение качества и токсического загрязнения воды в рассматриваемых экспериментах; 3 – на основе традиционных статистических индексов выявить изменение поведения компонентов систем и системных свойств после оказываемого воздействия |
|||||||||||||||||||||||||
|
Содержание эксперимента Э3М1 |
Содержание эксперимента Э3М2 |
|||||||||||||||||||||||||
|
Сравнивались результаты двух модельных экспериментов с разным фоном мутности воды (2 и 40 мг/л) на годовом этапе развития экосистемы. На основе параметров, используемых для оценки качества и токсического загрязнения воды, рассчитывался композитный индекс качества вод и выявлялось изменение класса качества после оказанного воздействия. Для построения композитного индекса качества воды использовались 7 компонентов: 1. Аммонийный азот NH4, мг/л. 2. Нитратный азот NO3, мг/л. 3. Концентрация растворенного кислорода, % насыщения. 4. Фосфор фосфатный PO4, мг/л. 5. Нитритный азот NO2, мг/л. 6. Концентрация общего фосфора, Робщ, мг P/л. 7. Концентрация общего азота, Nобщ, мгN/л. Индекс ИПК = ИПК(q, w) = = ИПК(q1,…,q7; w1,…,w7) = где qi – нормированное значение характеристики, wi – вес параметра в сводной оценке. Для wi использовалось: 0,0 ≤ wi ≤ 1,0 и |
Сравнивались результаты модельных экспериментов на годовом этапе развития экосистемы по изменению качества и токсического загрязнения речных вод при постоянном содержании условного токсиканта в воде с учетом типа его воздействия на первичные продуценты (1 тип – стимуляция первичной продукции низкими концентрациями токсиканта или 2 тип – ингибирование продукционных возможностей высокими концентрациями токсиканта). На основе параметров, используемых для оценки качества воды, рассчитывался композитный индекс качества вод и выявлялось изменение класса качества после оказанного воздействия. Для построения композитного индекса продуктивности водного объекта использовались 15 компонентов: 1. Аммонийный азот NH4, мг/л. 2. Нитратный азот NO3, мг/л. 3. Растворенный кислород, мг/л. 4. Фосфор фосфатный PO4, мг/л. 5. Нитритный азот NO2, мг/л. 6. Концентрация общего фосфора, Робщ, мг/л. 7. Концентрация общего азота, Nобщ, мг/л. 8. Взвешенные вещества, мг/л. 9. ХПК, мгO/л. 10. БПК5, мгO2/л. 11. БПК20, мгO2/л. 12. Растворенный органический углерод, мг С/л. 13. Биомасса фитопланктона, мг/л. 14. Разряды трофности в предположении однонаправленного изменения трофности и качества воды, балл. 15. Растворенный органический азот, мг N/л. |
|||||||||||||||||||||||||
|
Окончание таблицы |
||||||||||||||||||||||||||
|
1 |
2 |
3 |
||||||||||||||||||||||||
|
Шкала ИПК получена для 6 классов качества воды при равновесомости выбранных индикаторов: I ОЧ – очень чистые; II Ч – чистые; III УЗ – умеренно загрязненные; IV З – загрязненные; V Г – грязные; VI ОГ – очень грязные
(3) |
Индекс ИПК = ИПК(q, w) = = ИПК(q1,…,q15; w1,…,w15) Шкала ИПК получена для 6 классов качества воды при равновесомости выбранных индикаторов: I ОЧ – очень чистые; II Ч – чистые; III УЗ – умеренно загрязненные; IV З – загрязненные; V Г – грязные; VI ОГ – очень грязные.
(4) |
|||||||||||||||||||||||||
|
Э4. Выявление сезонных изменений качества и токсического загрязнения вод после оказанного на систему воздействия |
Задачи: 1 – оценить возможность использования имитационных моделей для оценки допустимости воздействия на экосистему; 2 – на основе построения композитного индекса качества и загрязнения по результатам моделирования воздействия выявить особенности сезонного изменения качества и токсического загрязнения вод |
|||||||||||||||||||||||||
|
Содержание эксперимента Э4М1 |
Содержание эксперимента Э4М2 |
|||||||||||||||||||||||||
|
Рассматривались сезонные изменения качества и токсического загрязнения воды по результатам двух модельных экспериментов с разным фоном мутности воды (2 и 40 мг/л). На основе параметров, используемых для оценки качества и токсического загрязнения вод водоема, рассчитывался композитный индекс качества воды и выявлялось изменение класса качества после оказанного воздействия на разных временных масштабах. Параметры и шкала для построения композитного индекса качества воды выбраны аналогично эксперименту Э3М1 |
Рассматривались сезонные изменения качества и токсического загрязнения речных вод по результатам модельных экспериментов развития экосистемы с учетом типа воздействия токсиканта на первичные продуценты (1 тип – стимуляция первичной продукции низкими концентрациями токсиканта или 2 тип – ингибирование продукционных возможностей высокими концентрациями токсиканта). На основе параметров, используемых для оценки качества вод, рассчитывался композитный индекс качества и выявлялось изменение класса качества после оказанного воздействия на разных временных масштабах. Параметры и шкала для построения композитного индекса качества и токсического загрязнения выбраны аналогично эксперименту Э3М2 |
|||||||||||||||||||||||||
,
.
.Материалы и методы исследования
Для выполнения поставленных целей и задач использовались две авторские имитационные модели пространственно-однородной водной экосистемы: «СNPXО-модель», далее модель «М1» [15] и модель «ВОДЭКО», далее модель «М2» [4]. На моделях были реализованы эксперименты, представленные в таблице.
Выполнение оценочных исследований эмерджентных свойств водных экосистем, а также построение композитных индексов рассматривалось во многих наших публикациях [7, 8, 15]. С учетом рассмотренных ранее особенностей интегральной оценки сложных систем акцентируем основное внимание на этапах формирования композитных индексов и построения композитных (интегральных) показателей по результатам экспериментов, проведенных на моделях М1 и М2. При этом процесс получения интегральных показателей реализуется в ходе выполнения нескольких этапов: 1 – отбор параметров для интегральной оценки, разработка и обобщение авторских классификаций и оценочных шкал для оценивания продукционной способности, качества и токсического загрязнения водной экосистемы; 2 – выбор нормирующей функции; 3 – выбор вида синтезирующей функции для композитного индекса; 4 – выбор весовых множителей для учета приоритетов в интегральной оценке; 5 – выполнение нескольких уровней обобщения информации для субиндексов и разработанных для них классификаций; 6 – получение оценочных шкал для субиндексов и композитных индексов для оценки системных свойств в целом; 7 – интегральная оценка системных свойств и эффектов для ключевых объектов и их временных или пространственных трансформаций.
Результаты исследования и их обсуждение
Ниже приведены результаты анализа экспериментов и основные выводы, полученные на основе анализа.
Эксперимент Э1М1. По мере увеличения фона мутности с 2 мг/л до 40 мг/л поведение компонентов экосистемы претерпевает значительные изменения. Изменение мутности оказывает негативное влияние на водную экосистему – происходит увеличение содержания «невостребованных» биогенных веществ в воде в связи со снижением продукционных способностей экосистемы. Отмечен сдвиг максимума цветения фитопланктона на более поздние сроки и высокое содержание фитопланктона в воде летом. Это приводит к росту абсолютных значений биомассы фитопланктона летом и незначительного увеличения среднегодового значения интегрального показателя трофности (ИПТ, шкала 1 в таблице) в пределах одного класса (мезотрофный, ближе к правой границе класса, первые проценты).
По рассчитанному критерию Нэша – Сатклиффа (классический вид [10]), использованному для выявления отрицательной значимости близости рядов, отмечается «неудовлетворительная» сходимость рядов для всех компонентов водной экосистемы, что свидетельствует о наличии серьезных различий в поведении компонентов модели в двух сравниваемых экспериментах. По рассчитанному критерию Тейла, согласно которому «0» соответствует полному совпадению рядов, а близость к «1» – отрицательной пропорциональности рядов, отмечается хорошая сходимость рядов для таких компонентов, как процентное содержание кислорода (0,035), общий азот (0,058), свидетельствующая о незначительном различии в поведении этих компонентов модели в сравниваемых экспериментах. Для остальных компонентов водной экосистемы индекс Тейла изменяется от (0,250 до 1), что подтверждает наличие различий в поведении этих компонентов.
Эксперимент Э1М2. На основе полученных значений ИПТ (шкала 2 в таблице) для речных вод получено, что при отсутствии поступления токсиканта и отсутствии антропогенного воздействия речные воды изначально относятся к классу высокопродуктивных (ИПТ = 0,378 – эвтрофия, левая граница класса). При постоянном содержании токсиканта в речных водах по 1-му типу воздействия (воздействие на первичные продуценты стимулирует жизненные функции при низких концентрациях и подавляет при высоких концентрациях), ИПТ остается в том же классе, сдвигаясь ближе к середине класса (ИПТ = 0,408) относительно эксперимента без воздействия токсиканта. При постоянном содержании токсиканта в речных водах, воздействие которого приводит только к подавлению жизненных функций (2-й тип воздействия), система также остается в эвтрофном классе, как и в предыдущем эксперименте (ИПТ = 0,408). Это объясняется докритическими значениями концентраций токсиканта в экспериментах, не изменяющими естественный ход процессов массообмена. Таким образом, экспериментально доказан постепенный переход водной экосистемы внутри одного класса от левой границы класса к его середине (отмечается рост ИПТ на 7,9 %).
В паре экспериментов «без имитации токсиканта – 1-й тип воздействия токсиканта» согласно рассчитанному критерию Нэша – Сатклиффа, отмечается хорошая сходимость рядов практически для всех компонентов системы, что подтверждает незначительные различия в поведении компонентов модели в экспериментах. Исключением являются: концентрации общего азота, общего фосфора, а также отношение их концентраций – отмечается неудовлетворительная сходимость рядов, указывающая на наличие изменений в динамике данных компонентов после оказанного воздействия. Согласно критерию Тейла отмечается хорошая сходимость рядов для таких компонентов, как скорость фотосинтеза (0,089), чистая первичная продукция (0,094), общий азот (0,096), свидетельствующая о незначительном различии в поведении этих компонентов модели в сравниваемых экспериментах. Для остальных компонентов водной экосистемы значения индекса Тейла указывают на наличие различий в поведении компонентов модели. В паре экспериментов «без имитации токсиканта – 2-й тип воздействия токсиканта» согласно рассчитанным критериям Нэша – Сатклиффа и Тейла, отмечается ситуация идентичная первой паре экспериментов.
Эксперимент Э2М1. Для выявления сезонных изменений продукционной способности водной экосистемы рассматривалось изменение значений композитных индексов продуктивности при повышении фона мутности. Так, для зимнего периода отмечается практически полное совпадение ИПТ (шкала 1) с 0,392 до 0,397 (мезотрофный, ближе к правой границе класса). Для весеннего периода ИПТ не изменяется (ИПТ = 0,390). В летний период наблюдается изменение композитного показателя ИПТ с 0,374 (мезотрофный, правая граница класса) до 0,403 (эвтрофный, левая граница класса). Для осеннего периода, напротив, отмечается снижение ИПТ в пределах одного класса с 0,425 (эвтрофный, левая граница класса) до 0,387 (мезотрофный, ближе к правой границе класса). Для вегетационного периода (май – октябрь) прослеживается незначительное увеличения значения ИПТ с 0,387 до 0,403 (4,1 %). Для летнего периода отмечается рост ИПТ (до 7,8 %). Основной результат эксперимента – сдвиг максимума цветения фитопланктона и высокое содержание фитопланктона летом приводят к росту ИПТ в весенне-летний период.
Эксперимент Э2М2. При выявлении сезонных изменений значения композитного индекса ИПТ (шкала 2, таблица) речных вод получено, что в зимний период при отсутствии поступления токсиканта и отсутствии антропогенного воздействия, система изначально относится к мезотрофному классу, ближе к его правой границе (ИПТ = 0,315); при постоянном содержании токсиканта в речных водах по 1-му типу воздействия отмечается переход системы в эвтрофный класс (ИПТ = 0,377, ближе к левой границе класса); при 2-м типе воздействия также отмечается переход системы в класс эвтрофных вод (ИПТ = 0,377, ближе к левой границе). Для весеннего периода при отсутствии поступления токсиканта в воду отмечается рост композитного индекса с 0,375 (эвтрофный класс, ближе к левой границе класса) до 0,442 (эвтрофный класс, ближе к левой границе класса) при 1-м типе воздействия и до 0,442 при 2-м типе воздействия. В летний период при отсутствии воздействия наблюдается ИПТ = 0,431 (эвтрофный, ближе к правой границе класса), а затем его рост до 0,456 внутри класса при 1-м типе воздействия и до 0,451 (эвтрофный класс, ближе к правой границе класса) при 2-м типе воздействия. В осенний период отмечается уменьшение значения интегрального показателя трофности с 0,392 (эвтрофный, ближе к правой границе класса) при отсутствии поступления токсиканта в воду до 0,359 (эвтрофный, ближе к правой границе класса) при 1-м типе воздействия токсиканта, и до 0,360 (при 2-м типе воздействия). В целом за вегетационный период (май – октябрь) отмечается рост ИПТ, но в рамках одного класса продуктивных (эвтрофных) вод, значения ИПТ близки для 1 и 2 типов воздействия токсиканта (с 0,420 до 0,436 при 1-м типе и до 0,434 при 2-м типе воздействия).
В итоге в зимний период отмечается переход речных вод из второго класса в третий при воздействии токсиканта при увеличении ИПТ (на 19,7 %). Для весны отмечается значительный рост ИПТ (на 17,9 %) (рис. 1).
Рис. 1. Внутригодовое изменение интегрального показателя продуктивности (трофности) ИПТ в эксперименте Э2М2
Эксперимент Э3М1. По мере увеличения мутности поведение компонентов претерпевает изменения, что приводит к росту (на 31,2 %) среднегодового значения интегрального показателя качества и токсического загрязнения воды (ИПК, шкала 3). Он изменяется от 0,032 (чистые, ближе к левой границе класса) до 0,042 (чистые, ближе к правой границе класса).
По рассчитанному критерию Нэша – Сатклиффа отмечается «неудовлетворительная» сходимость рядов для всех компонентов водной экосистемы, что свидетельствует о наличии серьезных различий в поведении компонентов модели в двух сравниваемых экспериментах. Согласно рассчитанному критерию Тейла для всех компонентов отмечается хорошее совпадение рядов для таких компонентов, как процентное содержание кислорода (0,035), общий азот (0,058), нитратный азот (0,053). Это свидетельствует о незначительном различии в поведении этих компонентов модели в сравниваемых экспериментах. Для остальных компонентов водной экосистемы индекс Тейла изменяется от (0,330 до 1), что подтверждает наличие различий в поведении этих компонентов в сравниваемых модельных экспериментах.
Эксперимент Э3М2. На основе полученных значений композитных показателей ИПК (шкала 4, таблица) речных вод получено, что при отсутствии поступления токсиканта и отсутствии других антропогенного воздействия речные воды изначально относятся к III классу умеренно загрязненных вод (ИПК = 0,214, правая граница класса). При постоянном содержании токсиканта в речных водах по 1-му типу воздействия, система остается в том же классе (ИПК = 0,216) относительно эксперимента без воздействия токсиканта. При постоянном содержании токсиканта в речных водах, воздействие которого приводит только к подавлению жизненных функций («2-й тип воздействия»), система остается в умеренно-загрязненном классе, как и в предыдущем эксперименте (ИПК = 0,216). Таким образом, ИПК практически не изменяется (на 1,2 %).
В паре экспериментов «без имитации токсиканта – 1-й тип воздействия токсиканта» согласно рассчитанному критерию Нэша – Сатклиффа, отмечается хорошая сходимость рядов практически для всех компонентов системы, что также говорит о незначительных различиях в поведении компонентов модели в рассматриваемых экспериментах. Исключением являются концентрации общего азота, общего фосфора, аммонийного азота, азота нитратов, азота нитритов, фосфора фосфатов, растворенного органического углерода. Для них отмечается неудовлетворительная сходимость рядов, указывающая на наличие кардинальных изменений в динамике данных компонентов после оказанного воздействия. Согласно критерию Тейла, для данной пары экспериментов наблюдается хорошая сходимость рядов для таких компонентов системы, как взвешенное органическое вещество (0,0), органический азот (0,05), свидетельствующая о крайне незначительном различии в поведении компонентов модели. Для остальных компонентов отмечается неудовлетворительная или полная несходимость рядов, свидетельствующая о наличии серьезных изменений в поведении компонентов модели. В паре экспериментов «без имитации токсиканта – 2-й тип воздействия токсиканта» согласно рассчитанным критериям Нэша – Сатклиффа (классический) и Тейла, отмечается абсолютно идентичная первой паре экспериментов ситуация.
Эксперимент Э4М1. Для выявления сезонных изменений качества и токсического загрязнения вод рассматривалось изменение значений ИПК (шкала 3, таблица) при повышении фона мутности с 2 мг/л до 40 мг/л. Для зимнего периода отмечается увеличение ИПК в пределах II класса с 0,043 (чистые, правая половина класса) до 0,050. Для весеннего периода характерно увеличение значения ИПК с 0,033 (чистые, ближе к левой границе класса) при фоне мутности = 2 мг/л до 0,041 (чистые, ближе к правой границе класса) при фоне мутности = 40 мг/л. В летний период наблюдаются изменения ИПК в пределах класса с 0,019 (чистые, ближе к левой границе класса) до 0,032. Для осеннего периода отмечается увеличение ИПК с 0,034 (чистые, ближе к левой границе класса) до 0,042 (чистые, ближе к правой границе класса). В целом для вегетационного периода (май – октябрь) прослеживается увеличение значения ИПК на 52 % с 0,023 (чистые, ближе к левой границе класса) до 0,035. Но переход в III класс не выявлен. Таким образом, изменение качества воды происходит в рамках II класса «чистые» воды, при этом эффект сильнее выражен в летний период (увеличение ИПК на 68,4 %) (рис. 2).
Рис. 2. Внутригодовое изменение интегрального показателя качества воды (ИПК) в эксперименте Э4М1
Эксперимент Э4М2. При прослеживании сезонных изменений значения ИПК речных вод не отмечается сильных изменений ИПК (шкала 4, таблица). В экспериментах получено, что в зимний период при отсутствии поступления токсиканта и отсутствии других антропогенных воздействий вóды изначально относятся к умеренно загрязненному классу, ближе к левой границе класса (ИПК = 0,185); при постоянном содержании токсиканта в речных водах по 1-му и 2-му типам воздействия» отмечается рост ИПК на 14 % (до 0,211). Для весеннего (ИПК = 0,216) и летнего (ИПК = 0,221) периодов значения интегрального показателя не изменяются. В осенний период отмечается уменьшение значения интегрального показателя качества с 0,236 (умеренно загрязненные, ближе к правой границе класса) при отсутствии поступления токсиканта в воды до 0,217 (умеренно загрязненные, ближе к левой границе класса) при «1-м типе воздействия» и при «2-м типе воздействия» токсиканта.
Заключение
Оценена реакция водной экосистемы на оказанное воздействие на основе построенных субиндексов продуктивности и качества воды, характеризующих экологический статус водного объекта. По результатам имитационного моделирования воздействий на двух разных имитационных моделях выявлено изменение среднегодовых и сезонных значений интегральных показателей продукционной способности, качества и токсического загрязнения вод. Важным результатом явилось выявление сезонных и пространственных особенностей рассмотренных системных эффектов. Показано увеличение интегральных показателей продукционной способности вод на 17,9 % (модель М2) в весенне-летний период при имитации антропогенного эвтрофирования и поступления токсиканта, а также рост интегрального показателя качества воды в летний период на 68,4 % (модель М1) в эксперименте с увеличением фона мутности до 40 мг/л.
Экспериментально доказан переход системы в другой класс трофности (из мезотрофного в эвтрофный) в летний период при увеличении фона мутности (М1); в эксперименте с имитацией антропогенного эвтрофирования и поступления токсиканта (М2) – в зимний период. Переход системы в другой класс качества вод (из загрязненных в умеренно загрязненные) отмечен в эксперименте с М2 в осенний период. Этими экспериментами доказан сезонный эффект изменения классов продуктивности и качества воды в водоеме и в проточной водной экосистеме.
В результате анализа полученных значений традиционных статистических индексов (индекса Нэша – Сатклиффа, индекса Тейла), которые были использованы для выявления изменения поведения компонентов системы и системных свойств после оказанного воздействия, было получено, что некоторые подсистемы (или компоненты) можно охарактеризовать как более чувствительные или наоборот – инвариантные, т.е. часто остающиеся неизменными по отношению к воздействию определенного типа на водную экосистему. В экспериментах к инвариантным компонентам (в которых отмечены минимальные изменения), согласно критерию Тейла, можно отнести: Nобщ, N-NO3 по модели М1.
При имитации поступления докритических концентраций токсиканта в речные воды (М2), согласно критериям Нэша – Сатклиффа и Тейла, такими компонентами являются: первичное звено трофической цепи – фитопланктон, компоненты РОВ (N-Nорг), взвешенное органическое вещество (детрит). Но во всех случаях это зависит от величины и специфики воздействия на экосистему. По-видимому, некоторые системные параметры являются более чувствительными к внешним воздействиям, чем другие. Показано, что разные системные трансформации проявляются с разной скоростью или требуют разного времени для формирования отклика на воздействие.
Следующими шагами исследования допустимости/недопустимости антропогенного воздействия на водную экосистему и выявления значений параметров, характеризующих необратимость процесса на основе имитационного моделирования и построения композитных индексов системных свойств будет выявление критических значений факторов, при которых системы при заданных климатических условиях и внешних воздействиях на длительное время или безвозвратно переходят в соседние классы состояния (продуктивности, качества и токсического загрязнения воды), а также других эмерджентных свойств системы. Для этого необходима реализация вариантов моделирования водных экосистем, соответствующих выбранным климатическим условиям и сценариям внешнего воздействия.
Исследования выполнялись при поддержке грантов РФФИ № 19-05-00683А; 18-05-60291 «Арктика».