THE USE OF MATHEMATICAL METHODS TO DEVELOP A GENERALISED SCALE OF ENVIRONMENTAL MODIFICATIONS

Research article
DOI:
https://doi.org/10.23670/IRJ.2023.132.47
Issue: № 6 (132), 2023
Submitted :
05.04.2023
Accepted:
18.05.2023
Published:
16.06.2023
888
2
XML
PDF

Abstract

Water ecosystems respond to anthropogenic impact as a whole by adaptive changes – ecological modifications. There are several stages of changes in the invariant state of ecosystems: ecological progress, ecological regression, ecological modulation. To determine the level of invariant state, based on the statistical analysis and modelling of water quality by biological indicators to the appropriate class, intervals of the most significant hydrochemical predictors were determined. A working scale for evaluating the invariant state of water ecosystems by hydrochemical and hydrobiological parameters has been developed. A statistical model linking changes in the chemical composition of water and changes in the qualitative and quantitative characteristics of hydrobionts as an ecosystem response to changes in the external environment has been obtained.

1. Введение

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

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

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

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

2. Методы и принципы исследования

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

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

.

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

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

.

3. Основные результаты

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Таблица 1 - Оценочная шкала инвариантного состояния водных экосистем по гидрохимическим показателям

Состояние экосистемы

Степень загрязнения

Гидрохимические параметры

БПК5,

мгО/л

ХПК,

мгО/л

NH4,

мг/л

Cu2+,

мкг/л

Фоновое

Условно чистая→загрязненная

<1,2

<20

<0,12

<1,4

Антропогенное напряжение

Очень загрязненная

1,2<C<2,3

20<C<22

0,12<C<0,3

1,4<C<2,0

Элементы Экологического регресса

Грязная

2,3<C<3,0

22,5<C<26

0,3<C<0,6

2,0<C<3,4

Экологическийрегресс

Метаболический регресс

Очень грязная – Экстемально грязная

C>3,0

>26

>0,6

>3,4

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

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

Таблица 2 - Оценочная шкала инвариантного состояния водных экосистем по биологическим показателям

Состояние экосистемы

Фоновое

Антропогенное напряжение

Элементы экологического регресса

Элементы экологического регресса и метаболический регресс

Степень загрязнения

Условно чистая→

загрязненная

Очень загрязненная

Грязная

Очень грязная – Экстремально грязная

Цилиопланктон

Число видов

2-5

>5

0-2

0

Численность млн.экз./м3

0,01-0,5

>0,5

0,1-0,2

0

Hn

>1,0

>2,2

<1,0

0

Фитопланктон

Численность млн.экз./м3

>1,4

1,4<Nf<2,3

2,3<Nf<8,0

Nf>8,0

Биомасса мг/м3

<1,2

1,2<Bf<2,3

2,3<Bf<4,3

4,3<Bf

Зоопланктон

Число видов

0<Qzp ≤3

4 ≤Qzp ≤ 10 и выше

9 < Qzp ≤ 11 и выше

3≤ Qzp<7

Численность тыс.экз./м3

Nzp <1

5 ≤Nzp ≤ 14

10≤Nzp ≤18

1≤ Nzp≤7

Биомасса, мг./м3

Bzp<2.5

2,5<Bzp< 20

20<Bzp< 28

0<Bzp<8

Зообентос

Число видов

9-17 и выше

7-9

3-6

0

Численность, тыс.экз./м2

0,9-4,0 и выше

0,8-1,5

0,1-3,0

0

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

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

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

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

g2=ln 21) = α2 + β2X

g3=ln 31) = α3 + β3X

g4=ln 41) = α4 + β4X

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

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

Таблица 3 - Коэффициенты логит модели при независимой переменной БПК5

Категория 

Коэффициенты

p-значения на основе теста Вальда

αj(Intercept)

βj

αj(Intercept)

βj

Численность зоопланктона

Q2

0,17213536

0,007348934

0,263529700

0,0121306081

Q3

-0,58519142

0,008172713

0,001451651

0,0061024539

Q4

-0,09033194

0,010454177

0,574739059

0,0003138126

 

Число видов зоопланктона

Q2

-0,3442828

0,09167911

2,081418e-02

2,553238e-03

Q3

-1,5012390

0,13793464

1,319017e-05

6,795876e-05

Q4

-1,1035093

0,16881321

1,946375e-04

3,978071e-08

Число видов фитопланктона

Q2

-0,8675664

0,04509488

0,0035127250

0,01557266

Q3

-1,0275592

0,05258722

0,0006995803

0,00495290

Q4

-1,1442695

0,06729152

0,0001260606

0,00019671

Биомасса фитопланктона

Q2

-0,5853540

0,1029370

4,036171e-04

7,729969e-04

Q3

-0,8920084

0,1347561

2,493293e-07

7,102055e-06

Q4

-0,9008073

0,1437082

1,846181e-07

1,627621e-06

Биомасса зообентоса

Q2

0,5113953

0,008964565

-0,5113953

0,008964565

Q3

-0,6070980

0,009310865

-0,6070980

0,009310865

Q4

-0,4315975

0,008974466

-0,4315975

0,008974466

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

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

Таблица 4 - Коэффициенты логит модели при независимой переменной ХПК

Категория

Коэффициенты

p-значения на основе теста Вальда

αj (Intercept)

βj

αj (Intercept)

βj

Число видов зоопланктона

Q2

-0,4969706

0,05859616

0,04371748

0,04371748

Q3

-0,5301726

0,05802088

0,04880583

0,032446113

Q4

-0,7122237

0,07609687

0,01256325

 0,01216185

Число видов зообентоса

Q2

0,3272833

-0,01617618

0,19725900

0,44274466

Q3

-0,1197736

-0,01198983

0,67052310

0,60698818

Q4

0,5926468

-0,05800079

0,02131093

0,01216185

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

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

Таблица 5 - Коэффициенты логит-модели при независимой переменной NH4

Категория

Коэффициенты

p-значения на основе теста Вальда

αj (Intercept)

βj

αj (Intercept)

βj

Число видов зоопланктона

Q2

-0,6631983

0,08639962

0,017099886

0,0007313045

Q3

-0,9195684

0,06624728

0,002922628

0,002922628

Q4

0,0715110 -

-0,01324632

0,79276805*

0,6371094290*

Число видов фитопланктона

Q2

-0,6049099

0,034964048

0,03998080

0,03054309

Q3

-0,7752531

-0,002368037

0,03188750

0,91150441*

Q4

0,6304148

-0,039089857

0,02297215

0,02572796

Число видов зообентоса

Q2

-0,06106748

0,02162308

0,813302409*

0,2910704582*

Q3

0,20654497

-0,05338239

0,060102473

0,0361728354

Q4

0,69529834

-0,08189699

0,007705239

0,0009815952

Примечание: *- коэффициенты не значимы

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

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

Таблица 6 - Коэффициенты логит-модели при независимой переменной Cu2+

Категория

Коэффициенты

p-значения на основе теста Вальда

αj (Intercept)

βj

αj (Intercept)

βj

Число видов зоопланктона

Q2

-0,5622712

0,05340661

4,061118e-02

5,853669e-02

Q3

-0,9559745

0,10104604

7,054292e-04

2,360975e-04

Q4

-1,2919422

0,12389183

1,422073e-05

9,865713e-06

Число видов фитопланктона

Q2

-0,7289491

0,0006017818

0,9912885*

0,48755566*

Q3

-1,0031341

0,0309536576

0,9996465*

0,03699848

Q4

-1,0894030

0,0471975543

0,9999292*

0,00189902

Биомасса зообентоса

Q2

0,32200021

-0,045978276

0,0212311

0,04597335

Q3

0,24607096

-0,050993378

0,0354476

0,03415627

Q4

-0,02864323

0,006107941

0,9104870*

0,76673592*

Примечание: *- коэффициенты не значимы

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

4. Обсуждение

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

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

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

. Последняя была применена нами при разработке комплексной шкалы.

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

. В основу модели положен индекс видового разнообразия (Н) инфузорий и 7 наиболее значимых химическими предикторов природного и антропогенного происхождения. Модель позволила успешно прогнозировать отнесение качества воды исследуемого водоема к категории «чисто/грязно».

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

5. Заключение

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

Article metrics

Views:888
Downloads:2
Views
Total:
Views:888