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

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

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

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

В качестве исходного материала были использованы многолетние ряды данных (2008-2018 гг.) по гидрохимическим и гидробиологическим показателям, полученные в результате мониторинга водных объектов Республики Татарстан. Мониторинг проводился на Куйбышевском водохранилище и его притоках (рр. Вятка, Казанка, Степной Зай, Меша, Свияга, Илеть, Юшут). В анализ были включены 21 гидрохимический компонент и 12 гидробиологических (число видов, численность, биомасса фитопланктона, инфузорий, зоопланктона, зообентоса), объединенных в одну общую базу данных.

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

Для построения модели в пакете R использовали функцию multinom() из пакета nnet, которая осуществляет оценку коэффициентов системы логит-моделей с использованием алгоритмов построения искусственных нейронных сетей.

Статистический анализ был выполнен с использованием программного обеспечения R

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

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

Рисунок 1 - Интервалы содержания органических веществ по ХПК (а), БПК5 (б), азота аммонийного (в), меди (г) в зависимости от инвариантного состояния экосистемы по показателям фитопланктона

Примечание: по оси абсцисс цифрами обозначены уровни инвариантного состояния: 1 – фон, 2 – антропогенное напряжение 3 – элементы экологического регресса+экологический регресс

DOI:10.23670/IRJ.2023.132.47.1

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

Рисунок 2 - Интервалы содержания органических веществ по ХПК (а), БПК5 (б), азота аммонийного (в), меди (г) в зависимости от инвариантного состояния экосистемы по показателям зоопланктона

Примечание: по оси абсцисс цифрами обозначены уровни инвариантного состояния: 1 – фон, 2 – антропогенное напряжение, 3 – элементы экологического регресса, 4 – экологический регресс

DOI:10.23670/IRJ.2023.132.47.2

Зообентос, также как и фитопланктон, слабо реагирует на переход системы из состояния экологического регресса в метаболический регресс, что отражается на разделении медианных значений гидрохимических показателей в зависимости от инвариантного состояния (рис. 3).

Рисунок 3 - Интервалы содержания органических веществ по ХПК (а), БПК5 (а), азота аммонийного (в), меди (г) в зависимости от инвариантного состояния экосистемы по показателям зообентоса

Примечание: по оси абсцисс цифрами обозначены уровни инвариантного состояния: 1 – фон, 2 – антропогенное напряжение, 3 – элементы экологического регресса + экологический регресс

DOI:10.23670/IRJ.2023.132.47.3

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

Рисунок 4 - Интервалы содержания органических веществ по ХПК (а), БПК5 (б), азота аммонийного (в), меди (г) в зависимости от инвариантного состояния экосистемы по показателям цилиопланктона

Примечание: по оси абсцисс цифрами обозначены уровни инвариантного состояния: 1 – фон, 2 – антропогенное напряжение, 3 – элементы экологического регресса + экологический регресс

DOI:10.23670/IRJ.2023.132.47.4

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

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

Аналогичное разделение по типам инвариантного состояния экосистем было сделано для биологических показателей (табл. 2).

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

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

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

Поскольку нами взяты 4 зависимых переменных (БПК5, ХПК, NH4- и Cu2+), которые включают 4 интервала значений, соответствующих состоянию экосистемы (Q1 – фоновое, Q2 – антропогенное напряжение, Q3 – элементы экологического регресса, Q4 – метаболический, экологический регресс), то для определения вероятностей отнесения объекта (по численности N или по чиcлу видов Q) к этим 4 категориям можно сформировать 3 недублированных логита, причём 1 категория (Q1) будет использоваться как эталонная. Тогда совокупность логит-моделей (baseline-category logit) будет иметь вид (модель 1):

g1=ln (π1/π1)=0, π1=1-( π2+ π3+ π4)

g2=ln (π2/π1) = α2 + β2X

g3=ln (π3/π1) = α3 + β3X

g4=ln (π4/π1) = α4 + β4X

где πj вероятность того, насколько предпочтительно присвоение произвольному объекту метки класса j. X – независимая переменная (либо N-суммарная численность, либо Q – число видов, либо В – суммарная биомасса).

Рассматривая в качестве зависимой переменной БПК5, а независимой переменной X численность, биомассу, число видов гидробионтов, получаем следующие значимые коэффициенты модели (1) (табл. 3).

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

Рассматривая в качестве зависимой переменной ХПК, а независимой переменной X численность, биомассу число видов групп гидробионтов, получаем следующие значимые коэффициенты модели (1) (табл. 4).

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

Рассматривая в качестве зависимой переменной NH4-, а независимой переменной X численность, биомассу число видов групп гидробионтов, получаем следующие значимые коэффициенты модели (1) (табл. 5).

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

Рассматривая в качестве зависимой переменной Cu2+, а независимой переменной X численность, биомассу число видов групп гидробионтов, получаем следующие значимые коэффициенты модели (1) (табл. 6).

Анализируя остальные коэффициенты, можно заметить, что число видов зоопланктона и зообентоса снижается при переходе от категории Q1 к категориям Q2, Q3 Q4. Число видов фитопланктона может варьировать, но закономерно снижается при переходе в категорию Q4.

Оценка инвариантного состояния отдельных водных объектов Республики Татарстан по утвержденным методикам, проводилась ранее

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

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

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

В ходе проведенного статистического анализа были выделены 4 гидрохимических показателя (содержание органических веществ по ХПК и БПК5, аммонийного азота и меди), которые оказались наиболее показательными в оценке, как качества воды исследованных водных объектов, так и инвариантного состояния водных сообществ. С 95% вероятностью были выделены интервалы медианных значений вышеуказанных показателей, соответствующие определенному инвариантному состоянию по показателям фито-, зоопланктона, зообентоса, цилиат. Для поиска связей между выделенными главными химическими и биологическими показателями отдельно для фито-, зоопланктона и зообентоса были построены модели мультиномиального логита. Построенная модель показала, что наилучшими предикторами для прогноза отнесения водоема к соответствующему инвариантному состоянию являются число видов зоопланктона (разделение по величине БПК5, ХПК, ионов меди), число видов фитопланктона (разделение по величине БПК5), число видов зообентоса (разделение по ХПК, ионам аммония), число видов, численность и биомасса зообентоса при разделение на инвариантные состояния по БПК5.