Scientific journal
European Journal of Natural History
ISSN 2073-4972
ИФ РИНЦ = 0,301

USING COMPOSITE INDICES TO STUDY THE ADEQUACY OF AQUATIC ECOSYSTEM MODELS AND ASSESS THE IMPACT ON THEM

Sedova S.A. 1 Dmitriev V.V. 1
1 Saint Petersburg State University
At the present stage of research aimed at studying and modeling natural transformations of ecosystems and geosystems, as well as assessing the state of complex systems in nature and society, an important stage is the assessment of complex (emergent, integrative) properties of aquatic ecosystems, analysis of the ability of ecosystems to maintain integrative properties and parameters modes under external influence on them. The article deals with the study of such integrative properties of an aquatic ecosystem as trophicity, quality and toxic pollution of water, as well as ecological (systemic) regulation of external influences on the ecosystem. A methodology for constructing composite indices for studying these integrative properties of the system is considered. The stages and results of the research are given. The assessment was carried out based on the results of the author’s imitation spatially homogeneous model of the cycle of carbon, nitrogen, phosphorus and other chemical elements, and the dynamics of oxygen in a shallow-water ecosystem (CNPXO model). To do this, the model was used to reproduce the annual development cycle of a «clean» aquatic ecosystem (standard version of the model) and an ecosystem under anthropogenic impact. For this, the planned version of the “impact” is implemented on the model. The possibility of investigating and assessing the impact on the system based on the construction of composite indices based on the results of simulation modeling is shown.
composite index
convolution of indicators
model-classification
integral assessment
trophic status
water quality
simulation
impact assessment

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

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

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

Развиваемая авторами статьи, а также зарубежными [4-6] и отечественными исследователями [7, 8] методология изучения сложных природных систем включает в себя 5 основных этапов: 1 – первичное измерение геокомпонентов системы и ее свойств полевыми и дистанционными методами. На 2 этапе выполняются анализ и обобщение собранной информации, на основе полученных данных разрабатываются оценочные шкалы и классификации компонентов и интегративных свойств системы (продуктивность, устойчивость, экологический статус, экологическое благополучие и др.). Этап 3 характеризуется моделированием отдельных свойств и связей в системах – создаются многофакторные зависимости, позволяющие описать совместный эффект влияния факторов среды на процессы переноса, синтеза, трансформации и транслокации вещества в модельной экосистеме; выполняются оценки интегративных (эмерджентных) свойств системы на основе построения моделей-классификаций и композитных индексов, проводятся многокомпонентные и многоуровневые свертки показателей по индексам. На этапе 4 выполняется синтез общей имитационной модели на основе моделей подсистем – осуществляется реализация модели на ПК, с последующей проверкой ее адекватности по отобранным критериям и интегративным свойствам [8]. После этого на модели проводятся тестовые эксперименты, позволяющие определить ее чувствительность к заданию параметров в функциях, описывающих скорости переноса и массообмена компонентами, что в дальнейшем позволяет исследователю проводить эксперименты в рамках экологического нормирования воздействия на водную экосистему [8-10]. На последнем этапе реализуется оптимизация модели. Проводятся эксперименты по оценке воздействия на экосистему, исследуются ее интегративные свойства, рассчитываются нормы допустимого/недопустимого воздействия на систему и др. [8].

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

Для оценки адекватности выбранной имитационной модели, учитывающей влияние факторов среды на скорости массообмена в экосистеме (CNPXО-модель) [10], и определения нормы воздействия на систему предлагается использовать серию композитных индексов, агрегирующих интегративные (эмерджентные) свойств водной экосистемы. При этом решаются две задачи: 1 – доказывается, что результаты моделирования и оценка эмерджентного свойства на основе композитного индекса дают возможность (или не позволяют) отнести исследованную систему к одному классу (трофности, благополучия, экологического статуса и др.); 2–проводится оценка нормы состояния системы и допустимого воздействия на систему по совокупности оцениваемых критериев на основе композитных индексов.

Теоретико-методологические основы исследования

Имитационные модели, отражающие естественную территориальную и хронологическую дифференциацию природных комплексов, в настоящее время представляют наибольший интерес для географических и геоэкологических исследований. Интегральная оценка трофности экосистемы, качества и токсического загрязнения воды в водоеме, оценка нормы воздействия на водную экосистему выполнялись на основе пространственно-однородной (точечной) модели круговорота углерода (C), азота (N), фосфора (Р) и других элементов (X), и динамики кислорода (O2) в мелководной экосистеме (CNPXО-модель) [10]. Полученные на модели решения позволили исследовать влияние факторов среды и биоты на изменение удельных скоростей массообмена в экосистеме; оценить скрытые от наблюдений составляющие баланса скоростей.

Модель водной экосистемы записывается в виде системы уравнений баланса скоростей массообмена между компонентами водной экосистемы [10]:

1) sedov01.wmf

2) sedov02.wmf

3) sedov03.wmf

4) sedov04.wmf

5) sedov05.wmf

6) sedov06.wmf

7) sedov07.wmf

8) sedov08.wmf

9) sedovVV01.wmf

10) sedovVV02.wmf

11) sedovVV03.wmf

12) sedovVV04.wmf

13) sedovVV05.wmf

14) sedovVV06.wmf

где БF – скорость первичного биосинтеза фитопланктона; sedov15.wmf – скорости выедания, B, D зоопланктоном; SF, SZ, SB – скорости естественного отмирания F, Z, B; RF, RZ, RB – скорости трат на обмен F, Z, B; AF, AB, АD – скорости ассимиляции пищи зоопланктоном; BD – скорость внеклеточного бактериального гидролиза детрита бактериями (валовая бактериальная продукция органического вещества); ДD – скорость деструкции детрита (D) бактериями; ГВС, ГBN, ГВР, ГВХ – скорости поступления в среду С, N, Р, Х в процессе внеклеточного бактериального гидролиза, ФZD – скорость поступления в детрит неусвоенной пищи (фекалий) при питании зоопланктона; OZD – скорость поступления в среду D в процессе питания зоопланктона («твердые отбросы пищи»); ОZC, ОZN, OZP, OZX – скорости поступления в среду компонентов РОВ: С, N, Р, Х в процессе питания зоопланктона («жидкие отбросы пищи» или псевдофекалии); ЭВС, ЭBN, ЭВР, ЭВХ, ЭFС, ЭFN, ЭFP, ЭFX, ЭZС, ЭZN, ЭZP, ЭZX – скорости поступления в среду С, N, P, X с тратами на обмен F, Z и В; sedov16.wmf – скорости поступления в среду NH4, и PO4 c экскретируемыми метаболитами Z; MC, MN, MP – скорости минерализации растворенных органических C, N, P; sedov17.wmf – скорости нитрификации I и II ступени; sedov18a.wmf sedov18b.wmf sedov19.wmf – скорости изъятия из воды NH4, NO2, NO3 , PO4,CO2, X в процессе биосинтеза водорослей; sedov20.wmf – скорость поступления в воду кислорода в процессе биосинтеза F; sedov21.wmf – скорости выделения в вводу СО2 в процессе дыхания F, Z, B; sedov22.wmf – скорости изъятия из воды кислорода в процессе дыхания F, B, Z; sedov23.wmf – скорости газообмена СО2 и О2 через поверхность раздела «вода – воздух». В каждом уравнении присутствуют скорости Qi, описывающие внешнее поступление (изъятие) компонентов в системе на шаге по времени.

Материалы и методы исследования

Примеры внутригодового изменения компонентов водной экосистемы на основе СNPXО-модели приведены на рис. 1. и рис. 2. Именно эти результаты имитационного моделирования использовались авторами для: 1 – определения по результатам моделирования интегративных свойств водной экосистемы, чтобы ответить на вопрос, какую трофность и какое качество воды формирует сочетание климатических условий, условий среды (и биоты) для каждого месяца и сезона года (или вегетационного периода в целом); 2 – системного (экологического, экосистемного) нормирования воздействий и количественного отражения эффекта внешнего воздействия на водную экосистему.

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

Sed1.tif

Рис. 1. Реализация климатического варианта точечной СNPXО-модели («КЧ»). Горизонтальными линиями на графиках подписаны средние значения компонентов по месяцам года. Цифрами подписаны начальные условия для компонентов

Sed2.tif

Рис. 2. Влияние мутности воды на внутригодовое изменение компонентов в водной экосистеме на фоне варианта «КЧ». Вариант КЧ, сплошная линия – М = 2 мг/л; пунктир – М = 6 мг/л; точки – М = 40 мг/л ежесуточно в течение года

Для исследования системного нормирования возможных воздействий на экосистему на модели был воспроизведен внутригодовой цикл развития: «чистой» водной экосистемы (климатический или стандартный вариант модели, обозначим его «КЧ») и экосистемы, находящейся под антропогенным воздействием (климатический с воздействием – «КВ»). При этом на модели нами реализован некоторый запланированный вариант «воздействия», а именно – повышение фона мутности воды c 2 мг/л до 40 мг/л. При построении композитных индексов для «чистого» варианта и варианта с воздействием использовались одни и те же веса выбранных для свертки показателей. В этом случае важно было оценить, изменятся ли трофический статус (трофность) и качество воды на временном интервале исследования, и, в случае, если система перейдет в другой класс, сделать выводы по допустимости/недопустимости воздействия или о влиянии воздействия на водную экосистему. Допустимым воздействием при этом считалось смещение экосистемы внутри класса, без выхода за его границы. Если же воздействие приводит к выходу системы по величине интегрального показателя за пределы класса, в котором она находилась до воздействия, то такое воздействие должно быть охарактеризовано как недопустимое. На основе этого мы можем также оценить устойчивость экосистемы к оказанному воздействию: 1 – изменит ли система свои свойства и параметры режимов и перейдет ли (частично, полностью) в другой класс (трофности, качества и др.); 2 – сможет ли система вернуться к прежним свойствам и параметрам режимов после устранения оказанного воздействия на интервале моделирования.

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

Начальные условия при моделировании уточнялись на основе 3-х летнего непрерывного решения задачи без дополнительного поступления/изъятия компонентов так, чтобы в начале следующего года задавались значения компонентов, полученные в конце предыдущего года. Так продолжалось до тех пор, пока разница между концентрациями компонентов в начале и в конце года не совпадала (с точностью 5 %). Поэтому начальные условия на рис. 1 подписаны с данной точностью.

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

В результате осреднения оценочных шкал разных авторов, рассмотренного в [7] и др. публикациях, нами были получены шкалы последнего уровня свертки для дальнейшей интегральной оценки трофности, качества и токсического загрязнения экосистемы, которые представлены в табл. 1 и табл. 2.

Затем, в результате проведенной оценки интегративных свойств системы – трофности и качества воды, были получены таблицы (табл. 3 и табл. 4) для двух вариантов расчета – «чистого» варианта «КЧ» и варианта с воздействием «КВ» (увеличение мутности с 2 мг/л до 40 мг/л).

Было получено, что в естественных условиях («КЧ» – вариант модели) вода в водоеме относится к «чистой» (II класс). Данный класс отмечается во все месяцы, кроме мая: в переходный сезон от весны к лету воды относились к I классу («очень чистые»), ближе к правой границе качества воды. Сезонная изменчивость качества воды не выявлена – в каждый сезон воды относились ко II классу («чистые»). Однако, можно отметить смещение ИПК внутри класса: весна и лето – ближе к левой границе класса, осень и зима – ближе к правой границе класса.

Полученные результаты расчета эксперимента с воздействием позволяют сделать вывод о том, что воздействие на систему привело к изменению класса качества воды в отдельные месяцы года. Так, качество воды в мае изменилось с «очень чистых» (I класс) на «чистые» (II класс); в июне с «чистых» (II класс) на «умеренно-загрязненные» (III класс), т.е. качество перешло в более старший класс. По сезонам года вода, по-прежнему характеризуется II классом, при этом выявлено смещение значений интегрального показателя (ИПК) в сторону правой границы класса.

Таблица 1

Итоговая шкала интегрального показателя качества и токсического загрязнения воды.

Показатель

Классы качества воды и токсического загрязнения

Очень чистые

Чистые

Умеренно загряз.

Загрязнен.

Грязные

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

I

II

III

IV

V

VI

Интегральный показатель качества и токс. загрязнения воды (ИПК)

1-0,960

0,960-0,917

0,917-0,849

0,849-0,668

0,668-0,264

0,294-0

 

Таблица 2

Итоговая шкала интегрального показателя трофического статуса

 

Классы трофического статуса (трофности)

Олиготрофный

I

Мезотрофный

II

Эвтрофный

III

Гипертрофный

IV

Интегральный показатель трофического статуса водоема (ИПТ)

1-0,794

0,794-0,468

0,468-0,185

0,185-0

 

Таблица 3

Сравнение интегральных показателей качества воды и токсического загрязнения для «КЧ» и «КВ» вариантов

 

Месяц

Янв.

Фев.

Мар.

Апр.

Май

Июн.

Июл.

Авг.

Сен.

Окт.

Ноя.

Дек.

«КЧ» вариант

ИПК

месяц

0,922

II (п)

0,921

II (п)

0,919

II (п)

0,938

II (п)

0,971

I (п)

0,937

II (п)

0,939

II

0,940

II (л)

0,936

II (п)

0,925

II (п)

0,925

II (п)

0,919

II (п)

ИПК сезон

0,921

II (п)

0,942

II (л)

0,939

II (л)

0,929

II (п)

 

«КВ» вариант

ИПК

месяц

0,922

II (п)

0,921

II (п)

0,935

II (п)

0,937

II (п)

0,922

II (п)

0,917

III – II

0,933

II (п)

0,950

II (п)

0,932

II (п)

0,926

II (п)

0,920

II (п)

0,918

II (п)

ИПК сезон

0,920

II (п)

0,931

II (п)

0,933

II (п)

0,926

II (п)

 

Примечание: I, II, III – классы качества; буквами «л», «с», «п» обозначена близость результата к левой, правой границам класса или к его середине.

Таблица 4

Сравнение интегральных показателей трофности для «КЧ» и «КВ» вариантов

 

Месяц

Янв.

Фев.

Мар.

Апр.

Май

Июн.

Июл.

Авг.

Сен.

Окт.

Ноя.

Дек.

«КЧ» вариант

ИПТ

месяц

0,453

III (л)

0,454

III (л)

0,453

III (л)

0,482

II

(п)

0,548

II

(п)

0,616

II

(л)

0,492

II

(п)

0,429

III (л)

0,403

III (л)

0,425

III (л)

0,447

III (л)

0,441

III (л)

ИПТ сезон

0,449

III (л)

0,494

II (п)

0,512

II (п)

0,425

III (л)

 

ИПТ вегет. период

 

0,486

II (п)

 

«КВ» вариант

ИПТ

месяц

0,452

III (л)

0,453

III (л)

0,485

II

(п)

0,449

III (л)

0,566

II

(п)

0,549

II

(п)

0,367

III

(л)

0,353

III (л)

0,470

II

(п)

0,612

II

(л)

0,451

III (л)

0,438

III

(л)

ИПТ сезон

0,448

III (л)

0,500

II (п)

0,423

III (л)

0,511

II (п)

 

ИПТ вегет. период

 

0,483

II (п)

 

Примечание: I, II, III – классы качества; буквами «л», «с», «п» обозначена близость результата к левой, правой границам класса или к его середине.

Трофность системы в «чистом» варианте, остается постоянной в течение года – III класс (эвтрофный водоем), ближе к левой границе, за исключением переходного весенне-летнего периода (апрель – июль), когда по величине ИПТ получен II класс (мезотрофный), ближе к правой границе. Отмечается сезонная изменчивость продуктивности: весна – лето – II класс (мезотрофный), ближе к правой границе: осень – зима – III класс (эвтрофный), ближе к левой границе. В целом за вегетационный период (март-октябрь) водоем относится ко II классу (мезотрофный) (ИПТ = 0,486).

В случае оказываемого на систему воздействия получено, что продуктивность заметно изменяется по месяцам и сезонам года, но при этом интегральный показатель за вегетационный период остался практически неизменным: 0,486-0,483 – II класс – мезотрофный, ближе к правой границе. Стоит отметить, что в летний сезон продуктивность растет, ИПТ перешел из II класса (мезотрофный) в III класс (эвтрофный), в тоже время в осенний период система вернулась из эвтрофного в мезотрофный класс.

Таким образом, полученные результаты свидетельствуют о том, что при увеличении мутности воды в водоеме (до 40 мг/л) поведение компонентов претерпевает значительные изменения. Отмечается переход водоема в эвтрофный класс в летний период по сравнению с «чистым» экспериментом. Качество воды остается во II классе, однако смещается ближе к его правой границе. Результатом увеличения мутности стало увеличение содержания биогенных веществ и, предположительно, сохранение нагрузки приведет далее к ухудшению качества воды и возможному переходу водоема в более старший класс качества, а также к полному переходу экосистемы в более продуктивный эвтрофный класс.

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

Заключение

Рассмотрены важные аспекты исследования сложных природных систем с применением методов имитационного моделирования на примере пространственно-однородной (точечной) модели круговорота углерода, азота, фосфора и других элементов, и динамики кислорода в мелководной экосистеме (CNPXО-модель), и методологии построения композитных индексов для исследования интегративных свойств системы – трофности, качества и токсического загрязнения вод на основе разработанных, и/или существующих оценочных шкал и классификаций. Показана возможность исследования и оценки допустимого/недопустимого воздействия на систему. Для этого были рассмотрены 2 эксперимента, реализованные с помощью CNPXО-модели: 1 – «КЧ»-вариант (климатический или стандартный вариант модели), цикл развития «чистой» водной экосистемы; 2 – «КВ»-вариант (климатический с воздействием), цикл развития экосистемы, находящейся под антропогенным воздействием (увеличением мутности воды с 2 мг/л до 40 мг/л).

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

Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 19-05-00683 А.