CC BY
77
УДК 594 : 632.152; ББК: 28.691 : 20.18
Майшанова М.И., Стойко Т.Г., Демаков Ю.П.
ВЛИЯНИЕ ДЛИТЕЛЬНОГО ИЗВЕСТКОВОГО ЗАГРЯЗНЕНИЯ НА СООБЩЕСТВО МОЛЛЮСКОВ СОСНОВОГО БИОГЕОЦЕНОЗА
Maishanova M.I, Stojko T.G., Demakov Y.P.
INFLUENCE OF LONG-TERMED LIME AIR POLLUTION ON MOLLUSCA COENOSIS IN PINE BIOGEOCENOSIS
Ключевые слова: известковое загрязнение, сосновый биогеоценоз, сообщество наземных моллюсков, подстилка.
Key words: lime air pollution, pine biogeocenosis, terrestrial mollusca coenosis, forest
floor.
Аннотация
Приводятся сведения о влиянии длительного воздействия известковой пыли, выбрасываемой в атмосферу заводом силикатного кирпича, на сообщество наземных моллюсков (Mollusca, Gastropoda) соснового биогеоценоза. Выявлено 11 видов моллюсков, в том числе 2 новых для Республики Марий Эл вида. Установлено, что изменение плотности населения моллюсков и параметров их видовой структуры закономерно изменяются по градиенту загрязнения, достигая наибольших величин при определенных значениях кислотности и зольности лесной подстилки.
Abstract
The information on the influence of long-termed affection of lime dust dispersed into the atmosphere by the silex brick plant on the terrestrial mollusca coenisis of the pine biogeocenosis is presented. 11 species of mollusca, including 2 new for Mari El Republic species, are found. It is stated that the changing of the mollusca population density and their species structure regular reflects the pollution gradient reaching the largest values by the certain indices of acidity and ash-contents of the forest floor.
ВВЕДЕНИЕ. Ухудшение состояния природной среды, вызванное неуклонным возрастанием на нее техногенной нагрузки, привело к необходимости создания системы глобального экологического мониторинга, важное место в котором занимает биоиндикация с использованием различных групп организмов. Однако использовать биоиндикацию можно лишь в том случае, если хорошо изучены ответные реакции организмов на воздействие загрязняющих веществ, спектр которых, как и самих организмов, необычайно велик.
В настоящее время в биоиндикации широко используется почвенная мезофауна [1, 3,
4, 10, 17], которая в большинстве естественных биогеоценозов довольно обильна и разнообразна. Она, являясь одним из наиболее чувствительных компонентов экосистем, позволяет учитывать слабые и локальные воздействия, что важно при нормировании техногенного загрязнения почв. На текущий момент времени наиболее изучено влияние загрязнения природной среды на таких представителей почвенной мезофауны как дождевые черви, многоножки, пауки и насекомые [2, 6, 11, 12, 14, 16], но слабо изучены ответные реакции на нее фауны наземных моллюсков [5, 8, 13], которые во многих типах экосистем являются полноправными членами биоценоза. Сведения же о влиянии на них известкового загрязнения полностью отсутствуют.
ЦЕЛЬ РАБОТЫ - оценка влияния длительного известкового загрязнения среды в зоне воздействия завода силикатного кирпича на сообщество наземных моллюсков соснового биогеоценоза.
ОБЪЕКТЫ И МЕТОДИКА ИССЛЕДОВАНИЯ. Объектом исследования являлся сосновый биогеоценоз, расположенный с северной стороны ЗАО «Марийский завод силикатного кирпича» (квартал 27 Силикатного участкового лесничества Республики Марий Эл). Завод действует с 1953 года. Основным компонентом выбросов является оксид кальция (известь негашеная), пыль известняка и БЮ2. Древостой на объекте чисто сосновый разновозрастный (60+
100 лет) II класса бонитета, произрастающий на песчаной слабоподзолистой почве, подстилаемой древнеаллювиальными песками. Тип условий произрастания А2 (свежий бор).
Для оценки влияния интенсивности загрязнения на сообщество моллюсков в биотопе заложено девять учетных лент, расположенных на разном удалении от источника загрязнения (рисунок 1), на которых проведена количественная оценка состояния древостоя, подлеска, напочвенного покрова и подстилки (таблица 1). Между пробными площадями серий I и II проходит ЛЭП. Учетные ленты Ы и II-1 расположены на опушке леса, примыкающей к полотну железной дороги. В пределах учетной ленты !-2 имеется действующий очаг корневой губки (Heterobasidion аппо^тт (Бг.) Вге£), где в древостое образовалось «окно» диаметром 5 м. Учетные ленты II-1 и П-2 дублируют ленты Ы и !-2 в связи с наличием на !-2 очага корневой губки.
Рисунок 1 Таблица 1 -
Схема расположения пробных площадок на объекте исследования Параметры состояния биогеоценоза на объекте исследования
Показатели Значения показателей по ме ре удаления от источника загрязнения
80 м 100 м 110 м 130 м 190 м 280 м 340 м 390 м 1500 м
Полнота древостоя относительная 0,90 0,78 0,70 0,63 0,68 0,66 0,80 0,69 0,68
Лиственный подлесок, тыс. шт./га 2,06 1,72 1,96 1,46 0,64 0,16 0,20 0,26 0,38
Хвойный подлесок, тыс. шт./га 0,37 1,68 0,53 0,54 1,04 1,32 0,71 0,37 1,57
Травяной покров, % 12,6 6,9 7,6 1,8 4,8 7,3 3,0 2,1 4,5
Моховой покров, % 1,4 5,5 1,8 3,8 19,2 40,9 12,8 26,4 45,5
Мертвый 98,7 94,5 98,2 96,2 80,8 59,1 87,2 73,7 54,5
Показатели Значения показателей по ме ре удаления от источника загрязнения
80 м 100 м 110 м 130 м 190 м 280 м 340 м 390 м 1500 м
покров, %
Масса подстилки, кг/м2 16,0 8,5 18,2 8,9 20,7 20,2 11,4 9,3 4,0
Зольность подстилки, % 76,5 67,3 69,3 55,7 64,2 64,1 61,0 56,8 43,1
Органика, кг/м 3,8 2,8 5,6 3,9 7,4 7,3 4,5 4,0 2,3
Минеральная часть, кг/м2 12,2 5,7 12,6 5,0 13,3 12,9 7,0 5,3 1,7
Загрязнитель, кг/м2 10,5 4,0 10,9 3,2 11,6 11,2 5,2 3,6 0,0
pH солевой 7,73 7,38 7,54 7,39 7,18 7,18 7,18 7,06 5,73
pH водный 7,78 7,57 7,68 7,46 7,39 7,39 7,39 7,25 6,35
Установлено, что количество валежа по мере удаления от завода колеблется от 6,0 до 37,3 м3/га и изменяется по градиенту загрязнения случайным образом. Густота лиственного подлеска, состоящего в основном из ивы розмаринолистной (Salix rosmarinifolia L.) и ракитника русского (Cytisus ruthenicus Fisch. ex Woloszcz.), убывает по градиенту загрязнения, а хвойного, представленного можжевельником (Juniperus communis L.), изменяется спонтанно. Травяной покров на участке довольно разнообразен, но очень редкий. Его густота убывает по мере удаления от завода. Степень покрытия почвы моховым покровом варьирует в очень больших пределах, возрастая по мере удаления от источника загрязнения.
На учетных лентах провели оценку абсолютно-сухой массы и зольности подстилки, а также массы органического и минерального вещества в ней, которые характеризуют как степень загрязнения среды, так и условия существования моллюсков. Пробы подстилки брали в пятикратной повторности вместе с мхом в виде почвенного монолита размером 10х10 см, очищенного от минерального слоя. Массу выпавшего загрязнителя оценивали по формуле: Мз = Мп - Морг. - Ммф., где Мз - масса загрязнителя, кг/м2 ; Мп - абсолютносухая масса подстилки, кг/м2; Морг. - масса органики в подстилке, кг/м2; Ммф. - масса минеральной части подстилки в фоновом участке, кг/м2. Оценивали также показатель pH вытяжки из подстилки. Перед подготовкой к анализу из подстилки были выбраны мхи и другая живая растительность. Замер значений pH определяли с помощью рН-тестера и электрода Cheker by Hanna.
Для учета сообщества моллюсков в пределах каждой учетной ленты брали пробы подстилки на пяти площадках 25х25 см, расстояние между которыми составляло 10 м. Сбор проводили в сентябре 2011 года (исполнитель М.И. Майшанова). Определение малакофауны проведено профессором Пензенского государственного педагогического университета Т.Г. Стойко.
Для количественной оценки структурной организации сообщества моллюсков были использованы индексы биоразнообразия, широко применяемые в биогеоценологии [7, 9]:
1) видового богатства (ВБ), соответствующему общему числу видов;
2) видовой насыщенности (ВН), соответствующей числу видов, приходящихся на одну учетную площадку;
3) бета-разнообразия Уиттекера Р = ВБ/ВН - 1;
4) сложности Симпсона-Гибсона SG= 1/Е(р;)2;
5) выравненности сообщества Е = GS/ВБ;
6) доминирования Бергера-Паркера, равного обратной величине доли доминирующего вида, т.е. ВР = 1/pmax;
7) плотность населения малакофауны на 1 м .
Для сравнительной оценки степени сходства видовой структуры моллюсков на разных учетных лентах использовали коэффициент Жаккара, вычисляемый по формуле Kg = Zmin(A,B)/Zmax(A,B), где min, max (А,В) - соответственно минимум и максимум из пары значений показателей сравниваемых рядов А и В. Произведен также подбор функций, описывающих распределение численности моллюсков по градиенту загрязнения.
Исходные данные обработаны на ПК с использованием пакетов стандартных программ Exsel и Statistiсa (обработка и анализ собранного материала проведены М.И. Майшановой и Ю.П. Демаковым).
РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ. На учетных площадках было обнаружено в общей сложности 357 экземпляров моллюсков, относящихся к 11 видам, 9 родам и 9 семействам (таблица 2). В сборах отмечены находки двух новых для Республики Марий Эл видов моллюсков: Discus ruderatus и Vertigo pusilla. Семейства Valloniidae и Cochlicopidae представлены двумя видами, остальные семейства - по одному виду. Массовым по численности является семейство Valloniidae (68,4 %), в котором доминирует улитка Vallonia costata (66,4 %). Типичные условия обитания вида - умеренно влажные смешанные леса [15]. Наибольшая численность Vallonia costata наблюдается на учетной ленте I-2, удаленной от завода на 130 м. Высокое обилие вида в нехарактерных для него сухих борах в условиях сильного известкового загрязнения среды является интересным фактом, заслуживающим внимания. На фоновом участке этот вид отсутствует, зато отмечены единичные экземпляры трех видов моллюсков: Euconulus fulvus, Cochlicopa lubricella, и Vertigo pusilla, который также отмечен на ленте I-4, удаленной от завода на 280 м. Два вида (Fruticicola fruticum и Succinella oblonga) представлены единственными экземплярами на учетных лентах I-4 и I-2. Моллюск Punctum pygmaeum отмечен только на учетной ленте I-6 в количестве двух особей.
Таблица 2 - Таксономическая структура малакофауны на объекте исследования
Семейство / Вид Численность моллюсков на учетных лентах, экз.
II-1 I-1 II-2 I-2 I-3 I-4 I-5 I-6 Фон I
Valloniidae: 9 40 30 59 40 31 20 15 0 244
Vallonia costata (Muller) 9 40 30 59 40 26 18 15 0 237
Valloniapulchella (Muller) 0 0 0 0 0 5 2 0 0 7
Euconulidae: 0 0 2 8 11 8 9 9 1 48
Euconulus fulvus (Muller) 0 0 2 8 11 8 9 9 1 48
Cochlicopidae: 0 3 0 3 4 6 2 13 1 32
Cochlicopa lubricella (Ziegler in Porro) 0 3 0 3 4 4 2 13 1 30
Cochlicopa lubrica (Muller) 0 0 0 0 0 2 0 0 0 2
Discidae: 1 1 2 0 1 2 5 4 0 16
Discus ruderatus (Ferussac) 1 1 2 0 1 2 5 4 0 16
Zonitidae: 0 0 6 0 0 1 1 3 0 11
Perpolita hammonis (Strom) 0 0 6 0 0 1 1 3 0 11
Punctidae: 0 0 0 0 0 0 0 2 0 2
Punctum pygmaeum (Draparnaud) 0 0 0 0 0 0 0 2 0 2
Vertiginidae: 0 0 0 0 0 1 0 0 1 2
Vertigo pusilla Muller 0 0 0 0 0 1 0 0 1 2
Bradybaenidae: 0 0 0 0 0 1 0 0 0 1
Fruticicola fruticum (Muller) 0 0 0 0 0 1 0 0 0 1
Succeneidae: 0 0 0 1 0 0 0 0 0 1
Семейство / Вид Численность моллюсков на учетных лентах, экз.
ІІ-1 І-1 ІІ-2 І-2 І-3 І-4 І-5 І-6 Фон I
Succinella oblonga (Бгаратаиё) 0 0 0 1 0 0 0 0 0 1
Общее число, экз. 10 44 40 71 56 50 37 46 3 357
Примечание: Семейства и виды распределены в таблице по обилию особей;
семейства и виды с равным обилием представлены в алфавитном порядке.
Дисперсионный анализ показал (таблица 3), что основным источником вариации численности моллюсков и вида Vallonia costata на объекте исследования являются варианты опыта, представленные учетными лентами, а также ошибки, проявляющиеся в неоднородностях пространственного распределения малакофауны, что свидетельствует о необходимости увеличения объема выборки для повышения точности оценки.
Исследования показали, что все параметры структуры малакофауны закономерно изменяются по мере удаления от источника загрязнения (таблица 4). Так, изменение плотности населения моллюсков по градиенту загрязнения описывается куполообразный кривой с максимумом на расстоянии 130 м от завода (рисунок 2). Сходным образом изменяются значения видового богатства, видовой насыщенности и индекса Симпсона-Гибсона (таблица 5). Увеличение видового богатства моллюсков на учетной ленте 1-4 связано с увеличением проективного покрытия мохового покрова, под которым складываются благоприятные для многих видов улиток условий развития. Значение индекса Р-разнообразия и выравненности закономерно повышаются по градиенту загрязнения, а Бергера-Паркера -снижаются (рисунки 3, 4). Данные проведенных расчетов свидетельствуют, таким образом, что загрязнение среды известкой пылью действует на моллюсков положительно, но лишь до определенной степени. Это явление в полной мере вписывается в экологическое правило, известное под названием «доза-эффект».
Таблица 3 - Результаты дисперсионного анализа численности моллюсков
Вариация Сумма квадратов Число степеней свободы Средний квадрат Критерий Фишера Доля влияния, %
Рфакт. Р 0,05
Общая численность моллюсков
Вариантов 729,2 8 91,2 2,58 2,24 37,0
Повторностей 107,2 4 26,8 0,76 2,67 5,4
Ошибки 1132,4 32 35,4 57,5
Общая 1968,8 44 100,0
Численность Vallonia costata
Вариантов 529,2 8 66,2 2,75 2,24 38,5
Повторностей 76,4 4 19,1 0,79 2,67 5,6
Ошибки 769,2 32 24,0 56,0
Общая 1374,8 44 100,0
Таблица 4 - Значение параметров структурной организации фауны моллюсков
Параметр структуры Значение параметра структу ры на учетных лентах
ІІ-1 І-1 ІІ-2 І-2 І-3 І-4 І-5 І-6 фон
Плотность населения, экз./м2 32 141 128 227 179 157 118 147 10
Видовое богатство 2 3 4 4 4 9 6 6 3
Видовая насыщенность 1,2 1,6 2,2 2,4 2,6 4,2 3,0 3,8 0,6
Индекс Р-разнообразия 0,67 0,88 0,82 0,67 0,54 1,14 1,00 0,58 4,00
Индекс Симпсона-Гибсона 1,22 1,21 1,70 1,42 1,81 3,03 3,12 4,20 3,00
Параметр структуры Значение параметра структу ры на учетных лентах
11-1 1-1 11-2 1-2 1-3 1-4 1-5 1-6 фон
Индекс выравненности 0,61 0,40 0,42 0,36 0,45 0,34 0,52 0,70 1,00
Индекс Бергера-Паркера 0,82 0,83 0,59 0,71 0,55 0,33 0,32 0,24 0,33
Расстояние от завода, км
Рисунок 2 - Закономерность изменения плотности населения наземных моллюсков по градиенту известкового загрязнения
Расстояние от завода, км
Рисунок 3 - Закономерность изменения значения индекса Р-разнообразия по
градиенту загрязнения
Расстояние от завода, км
Рисунок 4 - Закономерность изменения значения индекса Бергера-Паркера
по градиенту загрязнения
Таблица 5 - Уравнения регрессии, отражающие зависимость изменения показателей структуры малакофауны по градиенту известкового загрязнения_________________________________
Показатель Уравнение Я2
Плотность населения, экз./м2 У = 100 X / (4,2-Х2 - 1,1 X +0,15) 0,736
Видовое богатство У = 100 X / (27,4-Х2 - 4,8 X +3,5) 0,728
Видовая насыщенность У = 100 X / (92,1^2 - 26,4 X +7,7) 0,881
Индекс Р-разнообразия Уиттекера У = 2,8 X0,819 0,889
Индекс Симпсона-Гибсона У = 100 X / (20,8^2 - 4,1^ +8,7) 0,922
Индекс выравненности У = 0,41 X + 0,39 0,735
Индекс Бергера-Паркера У = 0,55-ехр[- 8,1^ - 0,08 )] + 0,27 0,900
Примечание: X - расстояние от источника загрязнения, км; Я - коэффициент детерминации уравнения.
Расстояние не может служить в полной мере критерием степени загрязненности среды, поскольку количество выпавших загрязняющих веществ и характер их распределения определяются действием множества факторов: мощности производства, применяемых технологий защиты окружающей среды и физико-географических условий территорий. В связи с этим, необходимо подобрать другие показатели степени загрязнения среды. В качестве таковых при известковом загрязнении могут выступать, в частности, кислотность и зольность подстилки, поскольку для построения раковины улиткам необходим кальций и щелочная реакция среды. Расчеты показали, что данное предположение оказалось верным: плотность населения моллюсков тесно связана с данными показателями (таблица 6, рисунки 5, 6), которые, в свою очередь,
взаимозависимы. Характер изменения плотности населения моллюсков по градиенту загрязнения отображается куполообразной кривой с оптимумом, приходящимся на интервал 7,2-7,5 рН и 58 % содержания минеральных веществ в подстилке. Таким образом, как слишком высокое, так и низкое значения показателей, отражающие степень загрязнения окружающей среды известью, приводят к снижению численности моллюсков.
рН солевой
Рисунок 5 - Закономерность изменения плотности населения моллюсков от кислотности подстилки
Зольность подстилки, %
Рисунок 6 - Закономерность изменения плотности населения моллюсков от
зольности подстилки
Таблица 6 - Значение физических показателей подстилки и плотности населения моллюсков на объекте исследования
Показатель Уравнение Я2
Плотность населения, экз./м2 У = 100Х / [14,6-Х2 - 77,6 X +105,0] 0,788
Плотность населения, экз./м2 У = 100 Ъ / [20,4^Ъ2 - 102,9 Ъ +131,4] 0,784
Плотность населения, экз./м2 У =100 / [0,21^2 - 5,3 +58,1] 0,727
Зольность подстилки и pH солевой, % = 0,05 X + 3,9 0,772
Зольность подстилки и pH водный, % = 0,04 Ъ + 4,9 0,828
Примечание: X = 10 - pH солевой вытяжки; Ъ - 10 - pH водной вытяжки подстилки; W = ^ - 40) зольность подстилки, %; Я2 - коэффициент детерминации уравнения.
Расчеты показали, что наиболее сходен состав фауны моллюсков на учетных лентах 1-3 и 1-1, 1-3 и 1-2 (таблица 7), а малакофауна фона резко отличается от загрязненного участка.
Таблица 7 - Матрица коэффициентов количественного сходства Жаккара видовой структуры моллюсков на объекте исследования____________________________________________
Номер ленты Значения коэффициента Жаккара между учетными лентами
II-1 I-1 II-2 I-2 I-3 I-4 I-5 I-6 фон
II-1 1,00
I-1 0,23 1,00
II-2 0,25 0,58 1,00
I-2 0,13 0,60 0,42 1,00
I-3 0,18 0,79 0,52 0,67 1,00
I-4 0,20 0,47 0,53 0,44 0,58 1,00
I-5 0,27 0,35 0,43 0,35 0,48 0,61 1,00
I-6 0,22 0,27 0,34 0,29 0,40 0,45 0,60 1,00
фон 0,00 0,02 0,02 0,03 0,04 0,06 0,05 0,04 1,00
ВЫВОДЫ
1. Сообщество наземных моллюсков соснового биогеоценоза, испытывающего длительное загрязнение известковой пылью, состоит из 11 видов, относящихся к 9 родам и 9 семействам. Массовым по численности является семейство Valloniidae (68,4 %), в котором доминирует Vallonia costata (66,4 %), являющийся индикатором известкового загрязнения.
В сборах отмечены находки двух новых для Республики Марий Эл видов моллюсков: Discus ruderatus и Vertigo pusilla.
2. Изменение плотности населения моллюсков и значений параметров видовой структуры малакофауны закономерно изменяются по градиенту загрязнения.
3. Загрязнение среды известкой пылью действует на моллюсков положительно, но лишь до определенной степени, которое вписывается в экологическое правило, известное под названием «доза-эффект».
4. Известковое загрязнение среды наилучшим образом отражают кислотность и зольность подстилки, а не расстояние от источника загрязнения.
Библиографический список
1. Бутовский, Р.О. Устойчивость комплексов почвообитающих членистоногих к антропогенным воздействиям / Р.О. Бутовский - М.: Издательство День серебра, 2001. -322 с.
2. Воробейчик, Е.Л. Население дождевых червей (Lumbricidac) лесов Среднего Урала в условиях загрязнения выбросами медеплавильных комбинатов / Е.Л. Воробейчик // Экология - 1998. - № 2. - С. 102-108.
3. Гонгальский, К.Б. Почвенные беспозвоночные как биоиндикаторы промышленного воздействия в лесных экосистемах Центра Европейской России: автореф. дис. ... канд. биол. наук. - М., 2004. - 16 с.
4. Криволуцкий, Д.А. Почвенная фауна в экологическом контроле / Д.А. Криволуцкий
- М.: Наука, 1994. - 272 с.
5. Куранова, А.П. Информативность популяционных параметров моллюсков в оценке уровня загрязнения среды обитания / А.П. Куранова, Е.М. Романова, О.А. Индирякова // Современное состояние и приоритеты развития фундаментальных наук в регионах: тр. V Всерос. науч. конф. - Краснодар, 2008. - T. I. - С. 25-28.
6. Лапиня, И.М. Интерпретация изменения видового разнообразия почвенных беспозвоночных при загрязнении почвы индустриальными кальцийсодержащими выбросами / И.М. Лапиня // Влияние промышленных предприятий на окружающую среду: всесоюзная школа. - Пущино, 1984. - С. 115-118.
7. Мэгарран, Э. Экологическое разнообразие и его измерение / Э. Мэгарран - М.: Мир, 1992. - 181 с.
8. Насекина, Е.С. Влияние нефтяного загрязнения на жизнедеятельность моллюсков / Е.С. Насекина, Г.А. Петухова // Современные наукоемкие технологии. - 2008. - №5. -С. 72-72.
9. Песенко, Ю.А. Принципы и методы количественного анализа в фаунистических исследованиях / Ю.А. Песенко. - М.: Наука, 1982. - 286 с.
10. Пространная экология почвенных животных / А.Д. Покаржевский, К.Б. Гонгальский, А.С. Зайцев, А.С. Савин.- М.: Товарищество научных изданий КМК, 2007. -174 с.
11. Растительные сообщества и почвенная мезофауна территорий химических предприятий в степной зоне Башкирского Предуралья / Д.В. Зейферт, И.Х. Бикбулатов, К.М. Рудаков, И.Н. Григорьева - Уфа: Изд-во УГНТУ, 2000. - 166 с.
12. Сидоренко, М.В. Беспозвоночные - индикаторы состояния природных комплексов в условиях большого города (на примере г. Нижнего Новгорода) / М.В. Сидоренко // Изв. Самарского научного центра РАН. - 2001. - Т. 3. - № 2. - С. 358-366.
13. Снегин, Э.А. Использование видов наземных моллюсков в качестве индикаторов реликтовых ценозов / Э.А. Снегин // Вестн. Житомирск. пед. ун-та. - Житомир, 2002. -Вып. 10. - С.128-129
14. Соромотин, А.В. Влияние нефтяного загрязнения на почвенных беспозвоночных (мезофауны) в таежных лесах Среднего Приобья / А.В. Соромотин // Сибирский экологический журнал. - 1995. - №6 - С. 549-552.
15. Стойко, Т.Г. Определитель наземных моллюсков лесостепи Правобережного Поволжья / Т.Г. Стойко, О.В. Булавкина - М.: Товарищество научных изданий КМК, 2010.
- 96 с.
16. Суходольская, Р.А. Способы адаптации насекомых к антропогенному воздействию / Р.А. Суходольская, Г.А. Тимофеева // Принципы и способы сохранения биоразнообразия: материалы IV Всероссийской конференции с международным участием, 22-26 сентября, 2010. - Йошкар-Ола, 2010. - С. 387-388.
17. Эйтминавичюте, И.С. Закономерности формирования комплексов почвенных беспозвоночных под влиянием антропогенных воздействий в зоне дерново-подзолистых почв: автореф. дис. ... докт. биол. наук. - М.: ИЭМЭЖ, 1982. - 34 с.
