Введение
Общеизвестно, что почвы являются основой экологической пирамиды и базисом развития наземных экосистем. В условиях техногенно нагруженных территорий почвы претерпевают значительные качественные изменения вследствие негативного воздействия химического, бактериологического и других видов загрязнения, так и вследствие механического преобразования, а именно уплотнения и разубоживания почвенного покрова. В этой связи в геоэкологию было введено понятие «урбоземы» – почвы селитебных территорий и городских агломераций.
Биотестирование почв методологически является определением их экологического состояния в лабораторных условиях (in vitro). В отечественной и зарубежной научной литературе описываются методические подходы к биотестированию экологического качества почв.
В исследовании О. А. Пестовой и А. И. Чупахиной описаны специфические физиологические группы микроорганизмов, которые встречаются лишь в присутствии определенных загрязняющих веществ и реагируют только на специальные химические соединения, например на содержание тяжелых металлов в почве (Пестова, Чупахина, 2017).
В работе самарских авторов показано, что наиболее информативным показателем качества почвы, загрязненной нефтепродуктами, является соотношение длин надземной и подземной частей фасоли. Загрязненность почвы нефтью проявляется в уменьшении данного показателя: 0.8 в нефтезагрязненной почве, 1.2 – в рекультивированной почве и 1.5 – в контрольной пробе (Баландина, Шабанова, 2007).
Исследование авторского коллектива из Российского государственного университета нефти и газа им. И. М. Губкина посвящено оценке качества городских почв с использованием в качестве тест-объектов червей энхитреид (Enchytraeus crypticus). Тест-реакции, регистрируемые в процессе проведения эксперимента, включают в себя смертность и изменение прироста численности особей (Смирнова и др., 2020).
Учеными из Лотарингии (Франция) в горнодобывающих районах, расположенных на юге Марокко, были отобраны пробы почв для проведения исследования с помощью набора для биотестирования MetPAD™. При этом тест на токсичность проводился одновременно с химическим анализом водных экстрактов хвостового материала и почв с целью оценки потенциальной доступности тяжелых металлов. Показано, что основными поллютантами являются цинк, медь и кадмий. Выявлено, что общая тенденция заключалась в увеличении токсичности металлов, измеряемой биотестом, с увеличением содержания доступных металлов в хвостохранилищах и грунтах. Поэтому испытание MetPAD™ может использоваться в качестве быстрого и чувствительного средства прогнозирования для оценки наличия тяжелых металлов в почвах, сильно загрязненных в результате горнодобывающей деятельности (Boularbah et al., 2005).
В работе исследователей из Монпелье (Франция) приведен сравнительный анализ модели активности свободных ионов (FIAM), модели наземного биотического лиганда (TBLM), модели диффузионных градиентов в тонких пленках (DGT) и биотеста растительного происхождения RHIZOtest по способности предсказывать концентрацию корневой меди в выращенной твердой пшенице (Triticum turgidum durum L.). Предполагалось, что ни один из методов не предусматривал адекватного прогноза концентрации корневой меди, что в основном коррелировало с общим количеством почвенной меди. Результаты измерений DGT и тем более прогноза FIAM отрицательно коррелировали с рН почвы и завышенным содержанием корневой меди в кислых почвах. Внедрение TBLM улучшило численное предсказание FIAM, но все же не смогло предсказать адекватную концентрацию корневой меди, т. к. формализм TBLM не учитывал щелочение ризосферы, как это наблюдалось на месте. Напротив, измерения RHIZOtest учитывали щелочение ризосферы и были в основном соотнесены с общим содержанием меди в почве (Bravin et al., 2010).
Учеными из Италии и Швейцарии в качестве тест-объектов почв, загрязненных нефтепродуктами, были использованы эктомикоризные грибы, произрастающие на саженцах норвежской ели и тополя. Продолжительность времени между загрязнением почвы и посевом влияет как на рост саженцев, так и на микоризальный инфекционный потенциал почвы. Полученные результаты подтверждают важность микоризальных грибов в биоремедиации почв, загрязненных сырой нефтью (Nicolotti, Egli, 1998).
Исследование польских ученых было направлено на определение степени снижения незначительного загрязнения почв бензином и дизельным топливом за счет внесения бентонита, компоста и оксида кальция. При этом в качестве тест-объектов были использованы культуры рапса и посевного овса, в которых проводилось измерение содержания азота. Показано, что добавление бентонита, окиси кальция или компоста в почву, загрязненную нефтепродуктами, обычно снижало отрицательное влияние бентонита и дизельного топлива на рост растений, содержание белкового азота и увеличивало общее содержание азота в растениях. Наиболее эффективным оказался бентонит, окись кальция и компост – чуть менее эффективными. Наиболее положительные результаты были получены для ярового рапса как основной культуры. Добавление в почву компоста, бентонита и окиси кальция оказало более сильное модифицирующее влияние на содержание азота в растениях на загрязненных дизельным топливом почвах, чем бензином (Wyszkowski, Ziólkowska, 2009).
Китайскими специалистами оценивалось влияние тяжелых металлов (Cd, Cu, Pb, Zn) на активность EROD и CYP3A4 с помощью дождевого червя (Eisenia fetida) для выяснения их возможной индукции и потенциала в качестве биомаркеров загрязнения почвы тяжелыми металлами. Показано, что среди четырех тяжелых металлов кадмий был наиболее сильно индуцирован EROD и CYP3A4. В то время как активность EROD и CYP3A4 показала схожую тенденцию, EROD является более чувствительным, чем активность CYP3A4 при использовании E. fetida в качестве биомаркера загрязнения тяжелыми металлами (Cao et al., 2012).
Все вышеперечисленные подходы к определению загрязнения почв объединяет детерминированный аспект, который позволяет оценить лишь вероятность наступления кризисного состояния в развитии растения или популяции особей животных – тест-объектов. Одним из перспективных методов определения загрязнения почвы является метод биотестирования, при котором оценивается динамика изменения параметров структуры тест-объекта во времени. В этом аспекте фракталы оказались чрезвычайно удобным математическим инструментом описания динамики развития природного объекта в изменяющихся условиях среды, в которых природный объект сохраняет свое самоподобие, т. е. ведет себя как одно целое.
Таким образом, одним из методов, позволяющим оценивать размер природного объекта или развитость его структуры в изменяющихся условиях внешней и внутренней среды, является фрактальная геометрия (Mandelbrot, Hudson, 2017). Такой подход широко применяется при оценке воздействия на биологические объекты, поскольку они обладают фрактальными свойствами, или самоподобием (рис. 1) (Розенберг, 2018; Насонов и др., 2018).
Рис. 1. Биоиндикаторы: лист березы и таллом лишайника (а) и моделирующие их фракталы (б)
Fig. 1. Bioindicators: birch leaf and lichen thalla (a) and fractals modeling them (б)
В основе предлагаемого метода лежит оценка фрактальной размерности – меры изменения сложности природного объекта в зависимости от масштаба. Математически это величина описывает статистическую меру сложности используемого шаблона фрактала при его масштабировании (Молчатский и др., 2016).
Материалы
Были выбраны три точки наблюдения, расположенные в различных локациях Москвы: ул. Фестивальная (Северный административный округ), Лосиный остров (Восточный административный округ) и станция метро Лермонтовский проспект (Юго-восточный административный округ).
10–14 января 2020 г. интегральные пробы почв отбирались с глубины от 0 до 20 см около жилого дома на ул. Фестивальной, у станции метро Лермонтовский проспект и в национальном парке «Лосиный остров» методом конверта (10 * 10 м). Масса каждой объединенной пробы составила около 1 кг. Для исключения ситуаций вторичного загрязнения исследуемых материалов все образцы были помещены в отдельные герметичные пакеты, чтобы обеспечить защиту от внешнего воздействия. Эксперимент проходил с 16 по 25 января 2020 г. Для проращивания был взят овес посевной (Avena sativa L.). Заранее пророщенные семена овса по одному были посажены в стаканчики с почвой.
Применительно к нашему эксперименту динамика пророста овса посевного в почвах с различным содержанием поллютантов будет различаться значениями фрактальной размерности D∈(1;2), которые для удобства можно нормировать к показателю Херста H = 2-D, отражающему динамику системного метаболизма (Насонов и др., 2019).
При этом развитость структуры проростков овса, последовательно оцениваемых в ходе эксперимента, связана с фрактальной размерностью степенной зависимостью:
где M(δ) – размер проростков; μδ – масштабы измерения; D – фрактальная размерность.
Чрезмерное содержание загрязняющих веществ может приводить к нарушениям системного метаболизма (Жукова, 2010), математически описываемого фракталом как статистическое нарушение самоподобия обменных процессов. Это означает, что природный объект не способен к своему воспроизводству как единого целого со средой существования: среда и объект существуют как отдельные системные целостности, не связанные друг с другом (Насонов и др., 2019).
При этом самоорганизация развития биотеста заключается в конкурентной динамике двух разнонаправленных процессов: ускорение метаболизма (при увеличении содержания биогенных компонентов – неорганических соединений азота, фосфора, серы и калия) и замедление обмена веществ (при увеличении концентрации, например, нефтепродуктов).
В работе (Корчагина, 2015) представлены результаты исследований особенностей пигментного аппарата, фотосинтетической активности листьев и продуктивной активности растений в условиях нефтяного загрязнения. В качестве объектов исследования были выбраны растения пушицы влагалищной (Eriophorum vaginatum), осоки острой (Carex acuta) и рогоза широколистного (Typha latifolia), произрастающие на территории Самотлорского месторождения. Показано, что нефтяное загрязнение значительно снижает интенсивность фотосинтеза, общее количество пигментов в ассимилирующих органах растений и биологическую продуктивность. То есть исследуется воздействие поллютантов, действующих вкупе с абиотическими факторами, в природных условиях.
В нашем случае при биотестировании загрязнения почв влияние абиотических факторов стремится к нулю. В ходе нашего эксперимента решается задача влияния загрязнения почв на динамику фрактальных параметров проростков овса посевного. При этом чрезмерное загрязнение приводит к хаотическим процессам роста овса, а умеренное – к балансу детерминированных и хаотических процессов в пределах экологической емкости среды, что порождает так называемый предельный экологический цикл воспроизводства овса в условиях загрязненной почвы.
При этом наблюдается устойчивая тенденция: чем выше загрязнение почвы, тем более хаотичной становится динамика роста тест-объекта, и наоборот. Поэтому сравнение наблюдаемых динамик в ходе эксперимента на предмет усиления хаотичности (увеличение амплитуды и частоты колебаний до пределов устойчивости Dd–Dk) может являться критерием загрязнености почв, которую также можно ранжировать по степени опасности. То есть по характеру фрактальной динамики мы можем определить степень загрязненности почвы, получая тот же результат, что и при дорогостоящем количественном химическом анализе почв. Иными словами, в схеме техно-природного объекта пропадает контур обратной связи с внешней средой, и объект далее не может устойчиво воспроизводиться как системная целостность (рис. 2).
Рис. 2. Схема техно-природного объекта с обратной связью (Кульнев и др., 2020)
Fig. 2. Scheme of a techno-natural object with feedback (Kulnev et al., 2020)
В этом смысле фрактальность, как свойство природного объекта, обеспечивающее системно достаточную развитость структуры при изменении параметров внешней среды, объясняет биотическую компенсацию привносимых в него техногенных субстанций вещества и энергии за счет самоорганизации (Манжуров, 2002). Отсюда следует, что устойчивость природного объекта, обеспечиваемая самоорганизацией его развития, является первичным качеством, без которого он неизбежно погибает (Розенберг, 2018). В условиях техногенеза только сложноорганизованные природные объекты способны к самоорганизации за счет масштабно-инвариантной перестройки своей структуры (Насонов, Цветков, 2018).
В соответствии с принципом толерантности нарушение самоорганизации природного объекта происходит в случае, когда значения системных показателей достигают лимитирующих значений, а именно:
– дефицита факторов, при котором объекту не хватает ресурсов для выживания в среде (показатель предельно низкого метаболизма, равносильного отсутствию обратной связи экосистемы с внешней средой);
– избыточности факторов, когда объект не успевает адаптироваться к частым изменениям внешней среды (показатель предельно высокого метаболизма с внешней средой, приводящего к нарушению функциональной целостности экосистемы) (Насонов и др., 2019).
В моделировании и анализе динамики сложноорганизованных техно-природных объектов ключевым атрибутом выступает их самоорганизация как компенсаторный биотический механизм внешних возмущений, который выстраивается на двух типах конкурирующих процессов: одни обеспечивают самосохранение системы, закрепляют ее строение и функциональные особенности за счет действия биотических факторов, другие обеспечивают развитие системы за счет действия техногенных факторов. Благодаря масштабно-инвариантной перестройке структуры экосистема имеет возможность балансировать эти процессы в обновляющейся среде существования (Иудин, Копосов, 2010).
С учетом принципа толерантности экосистемы уравнение самоорганизации техно-природных процессов задается соотношениями:
где H(t) – динамика метаболизма экосистемы в пределах лимитирующих факторов, измеренных в шкале показателя Херста; (Ho) – состояние экологического оптимума; Hd – состояние дефицита факторов; Hk – состояние избыточности факторов;– возмущающие техногенные факторы; – компенсирующие возмущение биотические факторы; – весовые коэффициенты соответствующих факторов.
Уравнение самоорганизации говорит о том, что в устойчивой экосистеме после снятия нагрузки ее состояния всегда обратимы за счет действия биотического фактора. При этом соотношение [2] определяет состояние сбалансированности техно-природных процессов относительно экологического оптимума (Ho), а соотношение [3] определяет область устойчивости экосистемы и направленность процессов (рис. 3).
Рис. 3. Тренды метаболизма экосистемы. Стрелками на оси ординат обозначена направленность техно-природных процессов (Кульнев и др., 2020)
Fig. 3. Trends in ecosystem metabolism. The arrows on the ordinate axis indicate the orientation of techno-natural processes (Kulnev et al., 2020)
Очевидно, что аппроксимация полученной статистики фрактальных показателей дает тренд смещения развития овса от экологического оптимума, вызванного нарушением процессов техно-природного обмена, которое тем больше, чем выше загрязнение почвы. При этом корреляция процессов, определяющая самоорганизацию параметров овса, задается соотношением:
Корреляция процессов [4] происходит в пределах устойчивого развития экосистемы, которая нарушается в точках Hk, Hd, когда они переходят в неустойчивое равновесие со средой, а их рост сменяется деградацией (Насонов и др., 2019).
В статистическом плане принципиальным отличием применяемых нами фрактальных моделей является их большая информативность в сравнении с традиционно используемыми линейными детерминированными и случайными гауссовскими распределениями, являющимися частными случаями фрактальных моделей с показателями Херста H = 1 и H = 1/2 соответственно (рис. 4).
Рис. 4. Распределение вероятности статистических событий при различных значениях постоянной Херста (H)
Fig. 4. Probability distribution of of statistical events at different values of Hearst constant (H)
Поскольку случайные (гауссовские) процессы в шкале показателя Херста имеют параметры (Подлазов, 2005), то такой режим трактуется нами как экологический оптимум, обеспечивающий устойчивость обменных взаимодействий природного объекта с внешней средой.
Динамические смещения статистических значений от экологического оптимума к показателю увеличивают вероятность развития детерминированных (трендоустойчивых) процессов, обуславливающих развитие природного объекта. И наоборот, приближение статистических значений к показателю увеличивает вероятность развития хаотических процессов со знакопеременным трендом, обусловливающих самосохранение природного объекта.
Если в динамике происходит резкое упрощение состояний экосистемы , то наступает истощение (дефицит) ресурса экосистемы. Вероятностное наступление такого события говорит о том, что экосистема перешла в неустойчивое равновесие со средой и уже не может регулировать свое состояние. При этом резко возрастает роль случайностей в развитии кризисных процессов.
Методы
Схема эксперимента и обработка полученных статистических данных для верификации уравнения самоорганизации овса в различных образцах почв состояла в реализации следующих этапов:
1. В различных районах Москвы (Лосиный остров, Фестивальная ул., Лермонтовский пр-т) были отобраны пробы почвы для проращивания заранее подготовленных семян овса (рис. 5). В качестве модельных тест-растений были использованы семена Овса посевного (Avéna satíva L.).
Рис. 5. Места отбора проб почв для проращивания овса
Fig. 5. Soil sampling sites for oat germination
2. Проращивание семян Avéna satíva L. на отобранных образцах почв производилось в течение 10 дней с ежедневной фотофиксацией образцов (рис. 6).
Важным обстоятельством является обеспечение для всех экспериментальных образцов одинаковых абиотических факторов: времени освещения, комнатной температуры 24 °C, ежедневного полива образцов 20 мл. Это необходимо для выделения влияния техногенного фактора, обусловленного наличием поллютантов в почвах.
Рис. 6. Динамика проростков овса на отобранных пробах почв
Fig. 6. Dynamics of oats sprouts on selected soil samples
3. Расчет фрактальных размерностей D по изображениям проростков овса через свободно распространяемую программу «Gwyddion» и построение графиков динамики фрактальных показателей в шкале постоянной Херста (рис. 7).
Рис. 7. Интерфейс программы «Gwyddion»
Fig. 7. Gwyddion Interface
Результаты
В основе способа расчета фрактальной размерности, используемого в программе «Gwyddion», лежит наиболее распространенный клеточный метод: квадратная решетка с постоянной (δ) накладывается на поверхность изображения. Изначально (δ) задается равной X/2 (где X – это длина края поверхности). Тогда N(δ) – число всех квадратов, содержащих хотя бы один пиксель изображения. Постоянная решетки на каждом шаге уменьшается в два раза, и процесс повторяется до тех пор, пока (δ) не станет равной расстоянию между двумя соседними пикселями. Наклон аппроксимирующей прямой, выделяющий область скейлинга (область масштабов, в которых объект сохраняет самоподобие), позволяет определить фрактальную размерность D (Терехов, 2011).
Применительно к изображению проростков овса этот метод основан на подсчете числа квадратов, покрывающих его изображение:
где D – фрактальная размерность; N(δ) – число квадратов, покрывающих изображение проростков овса; (δ) – варьируемый масштаб решетки покрытия.
Для нашего случая уравнение самоорганизации динамики роста овса, определяемое корреляцией процессов действующих факторов [2], задается уравнением тренда:
где – весовые коэффициенты действующих факторов (техногенных и биотических).
В нашем эксперименте мы анализируем прямое влияние загрязнений почв на адаптационные возможности овса посевного (Avéna satíva L.) и по характеру этого влияния делаем вывод о способности почв (при прочих благоприятных условиях) поддерживать сложившееся биоразнообразие. Решение таких задач в настоящее время входит в проблему формирования экологического каркаса территории (Манжуров, 2002).
Обсуждение
4. На основании расчетных статистических динамик геометрических параметров проростков овса строим тренды смещений от экологического оптимума и по характеру их направленности и весовым коэффициентам делаем формальную оценку влияний загрязнения почв на развитие биотеста.
Для почвы, отобранной на Лосином острове, функция фрактальных показателей для овса имеет вид:
Соответствующий график функции [7] изображен на рис. 8.
Характер динамики тренда говорит о том, что здесь мы имеем наименее загрязненную из всех образцов почву, практически не влияющую на динамику развития проростков овса, поскольку обменные процессы имеют слабо выраженную тенденцию приближения к экологическому оптимуму, а адаптационный процесс проходит в пределах устойчивости.
Рис. 8. Динамика фрактальных показателей биоиндикатора для почвы, отобранной на Лосином острове
Fig. 8. Dynamics of fractal parameters of the test object for the soil selected on Losiny Island
Более выраженная реакция тест-объекта на загрязнения почвы наблюдается на Фестивальной улице, для образцов которых уравнение динамики фрактальных показателей имеет вид:
Характер динамики аппроксимирующего тренда уравнения [8] говорит о том, что здесь среда развития овса оказывается менее благоприятна, чем на Лосином острове, что может быть следствием наличия в таких почвах более токсичных поллютантов. Однако образец демонстрирует стабилизацию процесса к концу эксперимента, что не свидетельствует о критичности загрязнений (рис. 9).
Рис. 9. Динамика самоорганизации биоиндикатора для почвы, отобранной на ул. Фестивальной
Fig. 9. Dynamics of self-organization of the test object for the soil selected on Festival Street
Наиболее загрязненная почва, оказывающая слабо выраженный отрицательный эффект на развитие Avéna satíva L., отмечается для образцов почв, отобранных на Лермонтовском проспекте (рис. 10). Здесь уравнение самоорганизации имеет вид
Слабо выраженная направленность тренда к показателю H = Hd говорит о том, что процессы стабилизации тест-объекта незначительно превалируют над процессами его развития. Однако, находясь в области знакопеременных трендов, ситуация, скорее всего, придет к стабилизации за более длительный промежуток времени, чем время эксперимента. Поэтому с точки зрения критичности загрязнения почвы ситуация здесь равнозначна ситуации на ул. Фестивальной.
Рис. 10. Динамика самоорганизации тест-объекта для почвы, отобранной на Лермонтовском проспекте
Fig. 10. Dynamics of self-organization of the test object for the soil selected on Lermontovsky Avenue
Результаты эксперимента наглядно демонстрируют прямое влияние загрязнений почв на адаптационные возможности тест-объекта, по характеру влияния которых можно делать вывод об их способности (при прочих благоприятных условиях) поддерживать необходимое биоразнообразие, что наиболее актуально в условиях техногенной нагруженности территорий.
Заключение
Результаты полностью соотносятся с местами отбора проб почв. Лосиный остров относится к особо охраняемым природным территориям (национальный парк), наименее подверженным техногенному загрязнению. Территории, относящиеся к Фестивальной улице и Лермонтовскому проспекту, демонстрируют более загрязненные почвы, однако ситуация на них далека от критической. Показателен факт, что на почвах, отобранных на границе Лосиного острова, примыкающей к автодороге и жилому сектору, ситуация по загрязнению почвы оказалась равнозначна Фестивальной улице.
Описанный метод может также эффективно использоваться при мониторинге загрязнения почв, а при дополнении – в уточнении границ санитарно-защитных зон через распределения плотности поверхностных загрязнений почв от стационарных и передвижных источников (Манжуров, 2002).
Таким образом, мы получаем возможность оценки влияния и классификации загрязнений почв урбанизированных территорий по характеру самоорганизации динамики развития тест-объекта. Это становится возможным благодаря определению структурных параметров биологических объектов в динамике, которые дают основу для разработки новых малозатратных и более информативных методов биотестирования.
Библиография
Баландина Л. П., Шабанова А. В. Опыт оценки качества рекультивации нефтезагрязненных почв методом биотестирования // Экология и промышленность России. Т. 11. № 11. М.: Калвис, 2007. С. 46–47.
Жукова Л. А. Экологические шкалы и методы анализа экологического разнообразия растений : Монография. Йошкар-Ола: Изд-во Марийского ун-та, 2010. 368 с.
Иудин Д. И., Копосов Е. В. Фракталы: от простого к сложному : Монография. Н. Новгород: Изд-во ННГАСУ, 2012. 200 с.
Кульнев В. В., Насонов А. Н., Жогин И. М., Цветков И. В., Грабарник В. Е., Карелин Н. В. Об опыте проведения управляемой альгоремедиации рекреационного водоема // Экология и промышленность России. Т. 24. № 3. М.: Калвис, 2020. С. 58–64. DOI: 10.18412/1816-0395-2020-3-58-64.
Манжуров И. Л. Фрактальная модель распределения плотности поверхностных загрязнений : Автореф. дис. ... канд. физ.-мат. наук. Екатеринбург, 2002. 24 с.
Молчатский С. Л., Казанцев И. В., Матвеева Т. Б. Применение метода фрактального анализа для биоиндикационной оценки состояния окружающей среды . // Самарский научный вестник. 2016. № 4 (17). С. 28–31.
Насонов А. Н., Цветков И. В. Особенности управления развитием техно-природных процессов на основе мультифрактальной динамики // Управление развитием крупномасштабных систем MLSD'2018: Труды одиннадцатой Международной конференции: В 3 т. / Под общ. ред. С. Н. Васильева, А. Д. Цвиркуна. М.: Изд-во ИПУ РАН им. В. А. Трапезникова, 2018. С. 83–88.
Насонов А. Н., Кульнев В. В., Графкина М. В. Моделирование динамики и прогноза экологических состояний воздушной среды селитебных территорий // Экология и развитие общества. 2019. № 2 (29). С. 56–64.
Насонов А. Н., Цветков И. В., Кизеев А. Н., Кульнев В. В., Мартынов Д. Ю., Сметанин В. И. Применение фрактального анализа в лихеноиндикации загрязнения атмосферного воздуха техногенно нагруженных территорий // Экология и промышленность России. 2019. Т. 23. № 3. С. 34–38. DOI: 10.18412/1816-0395-2019-03-34-38.
Насонов А. Н., Цветков И. В., Жогин И. М., Кульнев В. В., Репина Е. М., Кирносов С. Л., Звягинцева А. В., Базарский О. В. Фракталы в науках о Земле : Учеб. пособие. Воронеж: Ковчег, 2018. 82 с.
Пестова О. А., Чупахина А. И. Биотестирование почв с помощью микроорганизмов // Декада экологии: Материалы XI Международного конкурса. Омск, 2017. С. 39–43.
Подлазов А. В. Будущее прикладной математики: Лекции для молодых исследователей . М.: Эдиториал УРСС, 2005. С. 404–426.
Розенберг Г. С. Фрактальные методы анализа структуры сообществ // Принципы экологии. 2018. № 4. С. 4–43. DOI: 10.15393/j1/art.2018.8406.
Смирнова Т. С., Челознова К. В., Галкина А. А. Применение червей энхитреид в биодиагностике состояния городских почв // Экологические системы и приборы. 2020. № 2. С. 15–22.
Терехов С. В. Фракталы и физика подобия . Донецк, 2011. 255 с.
Benoit B. Mandelbrot, Richard L. Hudson. The (mis) Behavior of Markets. A Fractal View of Risk, Ruin, and Reward. Basic Books, 2017.
Boularbah A., Schwartz C., Bitton G., Morel J. L. Heavy metal contamination from mining sites in South Morocco: 1. Use of a biotest to assess metal toxicity of tailings and soils // Chemosphere. 2006. Vol. 63. Issue 5. P. 802–810. DOI: 1016/j.chemosphere.2005.07.079.
Matthieu N. Bravin, Aurélia M. Michaud, Bourane Larabi. Philippe Hinsinger RHIZOtest: A plant-based biotest to account for rhizosphere processes when assessing copper bioavailability // Environmental Pollution. 2010. Vol. 158. Issue 10. P. 3330–3337. DOI: 1016/j.envpol.2010.07.029.
Nicolotti G., Egli S. Soil contamination by crude oil: impact on the mycorrhizosphere and on the revegetation potential of forest trees // Environmental Pollution. 1998. Vol. 99. Issue 1. P. 37–43. DOI: 10.1016/S0269-7491(97)00179-6.
Wyszkowski M., Ziólkowska A. Role of compost, bentonite and calcium oxide in restricting the effect of soil contamination with petrol and diesel oil on plants // Chemosphere. 2009. Vol. 74. Issue 6. P. 860–865. DOI: 10.1016/j.chemosphere.2008.10.035.
Xiufeng Cao, Yufang Song, Jianrong Kai, Xiaoxia Yang, Puhui Ji. Evaluation of EROD and CYP3A4 activities in earthworm Eisenia fetida as biomarkers for soil heavy metal contamination // Journal of Hazardous Materials. 2012. Vol. 243. P. 146–151. DOI: 10.1016/j.jhazmat.2012.10.012.
Благодарности
Авторский коллектив выражает глубокую признательность Владимиру Николаевичу Большакову за поддержку и стимул к проведению исследований.