Особенности естественной радиоактивности почв и пород Рязанского региона

УДК 631.438.1

ОСОБЕННОСТИ ЕСТЕСТВЕННОЙ РАДИОАКТИВНОСТИ ПОЧВ И ПОРОД РЯЗАНСКОГО РЕГИОНА

Ю.А. Мажайский

Мещерский филиал ГНУВНИИГиМ им. А.Н. Костякова РАСХН

С.А. Тобратов

ФГОУ ВПО «Рязанский государственный университет им. С.А. Есенина»

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

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

Естественный радиационный фон обусловлен космическим излучением, а также наличием радиоактивных элементов в земной коре, водоемах и атмосфере. При этом на долю эндогенных радионуклидов (без учета радона) приходится свыше 25% облучения населения РФ природными источниками, что составляет около 20% годовой дозы облучения от всех источников, включая антропогенные - 0,77 мЗв/год из суммарных 4,06 мЗв/год (Старков, Мигунов, 2003). Вклад естественных радионуклидов (далее - ЕРН) наиболее распространенных в горных породах (232ТЪ, 226Яа и 40К) в 2 раза значительнее, чем космических лучей, и уступает только радоновым эманациям и медицинскому облучению. Радиоактивные элементы являются необходимой составляющей питательных сред, без которых невозможно развитие не только растительных, но и животных организмов (Белоусова, Штуккенберг, 1961).

Материалы и методы. Район исследований площадью около 15000 км2 охватывал различные по ландшафтным особенностям территории - от Мещерской низменности до Среднерусской возвышенности и Окско-Донского водораздела, включая г. Рязань и зону воздействия Рязанской ГРЭС. Исследовано свыше 100 образцов почв различных генетических типов, а также почво-образующих пород и природных материалов, используемых в строительстве, энергетике, медицине. В пробах выполнено гамма-спектрометрическое определение удельной активности основных ЕРН - 232ТИ, 226Яа и 40К на спектрометрическом комплексе «Гамма Плюс». При обработке и интерпретации данных использовали картографический метод и методы математической статистики (дисперсионный, регрессионный, корреляционный, вариационный, кластерный анализ).

Результаты и обсуждение. Исследованные радиоизотопы - 232ТИ, 226Яа и 40К - принадлежат к различным радиоактивным семействам и характеризуются определенными различиями в ландшафтно-геохимических особенностях. Так, торий присутствует в почвах в форме изоморфных примесей в минералах, в первую очередь, в первичных кислых силикатах, в кристаллическую ре-

шетку которых данный элемент встраивается легче. Калий, в т.ч. радиоактивный 40К, также может быть изоморфной примесью и, кроме того, интенсивно сорбируется глинами - гидрослюдами и монтмориллонитом. Для радия, кроме изоморфного замещения, большое значение имеют радиоактивный распад 238и в породах (радий - продукт данного распада), а также избирательное поглощение высшими растениями. Данное свойство растительных организмов, очевидно, унаследовано ими от древних форм жизни, для которых имела значение собственно радиоактивность радия, как дополнительный источник энергии [1]. По указанным причинам следует ожидать роста содержания ЕРН в почвах тяжелого состава, а для радия также должна играть роль гумусиро-ванность почв и наличие в среде родоначальника радиоактивного семейства - урана.

По данным В.Д. Старкова и В.И. Мигунова (2003), среднее содержание ЕРН в зональных почвах возрастает с севера на юг, от подзолистых почв тайги до сероземов субтропиков в 4-5 раз в связи с улучшением условий для накопления элементов. Таким образом, зональный фон ЕРН находится в обратной зависимости от коэффициента увлажнения. В Рязанском регионе, в условиях сравнительно однородного климата, основную роль играет гранулометрический состав почв. Данные таблицы 1 показывают, что различия в содержании ЕРН в почвах разных типов почвообразования статистически достоверны при р>99%, причем наиболее принципиальные различия наблюдаются между почвами легкого и тяжелого состава. Активности тория и калия в песчаных почвах (Мещера, Раново-Пронский долинный зандр) оказались значительно меньше фонового уровня для зоны подтайги, практически на «северотаежном» уровне. Это связано с многократным перемывом и переотложением почвооб-разующих пород данной территории в ледниковые эпохи плейстоцена. Наибольшая активность всех ЕРН отмечена в почвах зоны воздействия Рязанской ГРЭС, что, видимо, обусловлено атмотехногенезом. Примечательно, что максимумы содержания радионуклидов фиксируются к северо-востоку от предприятия по направлению ос-

1. Содержание естественных радионуклидов в гумусовом горизонте

Типы почв 232Th 226 Ra 4С K

Х V Х V Х V

Серые лесные 38,61 12,0% 36,28 24,6% 515,9 10,5%

Дерново-подзолистые 8,65 71,0% 14,58 40,6% 151,5 42,0%

Черноземы 35,73 15,4% 34,81 18,3% 451,0 11,1%

Аллювиальные суглинистые 36,51 24,8% 35,53 46,9% 485,7 26,8%

Аллювиальные песчаные 14,90 47,9% 15,85 39,7% 199,6 31,1%

Почвы западин* 32,67 17,4% 37,37 27,1% 514,7 12,2%

Серые лесные зоны воздействия Рязанской ГРЭС (0-12 км) 42,22 10,7% 39,97 22,0% 569,0 8,4%

Болотные торфяные (Мещера) 3,60 23,6% 21,10 66,4% 24,1 64,0%

0,851 0,561 0,873

а 71(Г22 4.1Р9 9.10г24

Примечание. Х - среднее значение, Бк/кг; V - коэффициент вариации, %; rw - внутриклассовый коэф( »ициент корреляции, харак-

теризующий влияние фактора (в данном случае - группировка по типам почв) на результативный признак и рассчитываемый по

формуле: = стА2 - 2 , где стА2 и ст22 - соответственно факториальная и случайная дисперсии значений концентраций, пс -

^а + (nc-1) ctz усредненный объем градаций изучаемого фактора; а - уровень значимости.

*согласно «Классификации почв СССР» (1977) - дерново-подзолистые поверхностно-глеево-элювиальные почвы.

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

Для торфяников Мещеры свойственны минимумы тория и калия - в 2-4 раза более контрастные, чем по средним данным для Европейской части РФ (Старков, Мигу-нов, 2003). Это обусловлено тем, что торфонакопление протекает в окружении озерно-аллювиальных песков, исходно бедных радионуклидами. Таким образом, Мещера - региональный минимум активности большинства радиоизотопов. Однако радий является характерным исключением. Снижение его содержания в песчаных почвах по сравнению с суглинистыми не столь существенное (в 2,5 раза против 4 раз для 232ТЪ и 40К). В данном случае дополнительным источником Яа становится радиоактивный распад урана. Если учесть концентрирование урана в кислых изверженных породах (Белоусова, 1961; Перельман, Касимов, 1999) и то, что четвертичные пески являются продуктом выветривания пород, подвергшихся ледниковой экзарации на Карельском щите, становится объяснимым относительно повышенное - по сравнению с другими радионуклидами - содержание 226Яа в песчаных отложениях и почвах, развитых на них. Кроме того, в мещерских торфах радий, в отличие «ли-тофильных» тория и калия, накапливается за счет биоконцентрирования: его больше, чем в минеральных почвах Мещеры в 1,5 раза (табл. 1). Интенсивное накопление радия растениями на фоне явного блокирования биопоглощения такого техногенного элемента, как 137Сб проявляется во всем Рязанском регионе. Вклад указанных факторов приводит к тому, что тип почвы определяет пространственное варьирование 226Яа лишь на 56%, тогда как других ЕРН - на 85-87% (табл. 1).

Информативным статистическим показателем, отражающим многофакторность внешних воздействий на процесс, служит коэффициент вариации. Как видно из

таблицы 1, максимальное варьирование свойственно содержанию радионуклидов в торфах (из-за различий между верховыми и низинными торфяниками), а также в дерново-подзолистых почвах. В последнем случае это фиксирует тот факт, что такие почвы развиваются на породах, обладающих сходным механическим составом, но значительно отличающимися по генезису и возрасту (от днепровского оледенения до голоцена). На втором месте по вариабельности - накопление ЕРН в условиях нестационарного аллювиального процесса (У=25-50%).

В соответствии с алгоритмом, изложенным в НРБ-99 [2], была рассчитана эффективная удельная активность природных радионуклидов, содержащихся в верхних 020 см горизонтах почв (Аэфф): Аэфф = АЯа + 1,3АТЬ + 0,09АК. Также определяли и долевое участие каждого элемента в суммарной Аэфф. Можно заключить, что в почвенно-геохимических условиях Рязанского региона отсутствуют механизмы, способные обеспечить превышение показателей Аэфф для материалов I и II классов, используемых в строительстве в пределах населенных пунктов (соответственно 370 и 740 Бк/кг). Атмотехно-генное воздействие Рязанской ГРЭС увеличивает Аэфф почвенного покрова, но незначительно: в 1,1 раза по сравнению с зональным фоном. Хорошо заметно снижение удельной активности ЕРН в почвах легкого состава, менее значительное в поймах вследствие частичной аккумуляции мигрирующих элементов. Минимальной Аэфф отличаются органические почвы. Соотношение радиоизотопов, несмотря на различия в генезисе почв и колебания Аэфф, остаются в целом практически неизменными, даже вблизи Рязанской ГРЭС: реализуется ряд торий > калий > радий. Исключениями являются дерново-подзолистые почвы и особенно торфяники, где радий выходит на первое место, обеспечивая соответственно 37 и 75% суммарной эффективной активности почв.

Почвенно-геохимическую неоднородность выборки фиксирует и вариационный анализ [3]. Для 232ТИ

и 40К

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

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

На радиогеохимических картосхемах распределения ЕРН в почвах района исследований хорошо заметна зависимость активности радионуклидов от механического состава почв: минимумы отмечаются в почвах Мещеры, Раново-Пронского зандра, в аллювиальных почвах легкого состава (средние течения многих рек), на внутри-пойменных песчаных останцах (например, в районе сел Коростово, Шумашь, Поляны). В Спасско-Шиловской пойме Оки, в условиях тектонического погружения территории и активной аккумуляции тонкодисперсного суглинистого аллювия, фиксируется устойчивый региональный максимум ЕРН; аналогичный, но менее контрастный максимум обусловлен, по-видимому, техногенным фактором и приурочен к зоне основного переноса выбросов от Рязанской ГРЭС (к северо-востоку от г. Но-вомичуринск). Отличия радия от тория и калия связаны с гораздо меньшей зависимостью радия от степени «опесчаненности»: пространственный градиент активности тория между Мещерой и окской поймой в несколько раз выше. Радий, кроме того, накапливается в верховых и переходных болотах, а торий - в почвах юго-восточной Мещеры, являясь тем самым маркером трансгрессий окской поймы, происходивших здесь в эпоху последнего, валдайского оледенения. Характерно, что и современные процессы развития долины Оки идут по тому же сценарию.

Природные материалы, разрабатываемые в регионе и используемые в качестве топлива, удобрений, строительного сырья и лечебных ресурсов (торфяные лечебные грязи), по величине Аэфф располагаются в следующий ряд (по возрастанию активности ЕРН): кварцевые пески - низинные торфа Мещеры - верховые торфа Мещеры - карбонатные известняки (во всех случаях абсолютно преобладает эффективная активность 22&Яа) - низинные торфа юга области (примерное равенство «тер-ригенного» 232ТИ и «биогенного» радия) - глины юрского возраста (преобладает 40К) - каменноугольные глины (абсолютно преобладает 40К и наблюдается превышение по Аэфф для стройматериалов I класса: 430 Бк/кг против максимально возможных 370 Бк/кг). Таким образом, карбоновые глины - наиболее яркий пример сверхнормативной аккумуляции ЕРН. Зональные условия седиментации (главным образом влажный экваториальный климат) способствовали их избыточному обогащению радиоактивным калием, что ограничивает возможности их использования в качестве стройматериалов. Примечательно при этом, что образец карбоновых глин был отобран в карьере Рязанского кирпичного завода. Другие исследованные материалы можно использовать без ограничений.

Перераспределение ЕРН под влиянием почвообразовательных процессов изучали на примере двух природных объектов. Первый - активно развивающаяся суффо-зионная западина на плоскоместном водоразделе Хупты и Лесного Воронежа (юг Ухоловского района), основным почвообразовательным процессом в которой является кислотный гидролиз [4]. Как показано на рисунке,

наблюдается слабый вынос ЕРН из гумусового горизонта и слабая аккумуляция в элювиальном горизонте А2 -под влиянием типичных для данных почв процессов: соосаждения элементов с новообразованными гидрокси-дами железа и гидрослюдизации. Поглощение гидрослюдами наиболее значительно для калия, поэтому распределение 40К по изученным горизонтам наиболее контрастно (рост в горизонте А2 на 25% по сравнению с содержанием в Аь против 10-15% для других ЕРН). Между тем тяжелые металлы (ТМ) в данных почвах гораздо динамичнее: в частности, свинец в горизонте А1 под влиянием кислотного гидролиза в основном переведен в водорастворимое состояние и активно выносится. Био-фильные элементы - медь и цинк - накапливаются в гумусовом горизонте западины, а токсиканты - свинец и кадмий - выносятся из него, частично аккумулируясь в А2. Таким образом, данные две группы ТМ принципиально различны по миграционным формам, что уже отмечалось нами ранее [5]. Характерно, что радионуклиды по своему поведению в почвах сходны с токсичными ТМ и включены в те же пути миграции. Впрочем, очевидно и то, что формы нахождения ТМ и ЕРН существенно различаются: последние закреплены в почвенной матрице многократно прочнее.

Второй объект исследования - долинный торфяник в эрозионной системе ручья Менек (север Ухоловского района, в 2 км к востоку от с. Соловачево, балка Сизый Менек) - представляет собой «геохимический полигон» изучения сернокислотного гидролиза минералов и гидрогенного ожелезнения. Торфяная залежь эксплуатируется как месторождение лечебных грязей Менек, торф вывозится и в Московский регион. Возраст торфяника оценен нами в 10-12 тыс. лет. В настоящее время он вскрыт попятной эрозией в долине ручья и имеет 2 горизонта ожелезнения - нижний, современный, охристый, ферригидрит-гетитовый и верхний, древний, многослойный, туфовидный, сидерит-ярозит-керченитовый. Причина ожелезнения - влияние грунтовых вод, содержащих продукты выветривания пиритоносных юрских глин, а также меловых фосфоритов, которые располагаются здесь близко к дневной поверхности. Наиболее

ТИ-232 Рэ-226 К-40

Рис. 1. Влияние экстремальной глеегенной деградации почв на удельную активность ЕРН

2. Удельная активность ЕРН в генетических горизонтах торфяника Менек и сопряженных с ним минеральных отложений, Бк/кг

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

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

Торфяная толща в долине ручья Сизый Менек подразделяется на горизонты существенно различной мощности, причем влияние грунтовых вод начинает сказываться с глубины свыше 70 см, к верхним слоям приурочен многослойный панцирь «болотного туфа», а в низах толщи развито ожелезнение. Многослойность туфовид-ных отложений свидетельствует о стадийности процессов увлажнения - иссушения торфяной толщи в прошлом. Новообразования самородной серы развиты только в узкой центральной части толщи, которая соответствует по высоте надпойменной террасе, по поверхности которой осуществляется сток минерализованных грунтовых вод. Нижние слои торфа, по-видимому, формировались по типу сапропелей и содержат минеральные примеси. Подстилающие минеральные отложения интенсивно оглеены.

Как свидетельствуют данные таблицы 2, содержание радионуклидов в органогенных отложениях значительно меньше, чем в почвообразующих породах и аллювии ручья, но минимальных значений активность достигает в современных охристых прослоях (наиболее значительный - на глубине 155-172 см) и особенно в современном (106-110 см) и древнем (60-72 см) горизонтах сернокислотного гидролиза, где происходит (или происходило в недавнем прошлом, до эрозионного вреза ручья) выклинивание грунтовых вод. Можно предположить, что терригенные частицы в агрессивной среде торфяника подвергаются растворению, и радионуклиды высвобождаются в раствор, мигрируя затем в поверхностные воды.

Тем не менее, возможна сорбция ЕРН органическим веществом торфяника, для оценки интенсивности которой рассчитаны соответствующие регрессионные зависимости (табл. 3), представительные только для органогенных горизонтов с зольностью не более 40%. В простейшем случае зависимость имеет вид прямой, а регрессионное уравнение - следующую структуру: У = а + Ьх,

где У - зависимая переменная; х - принятый в расчет фактор в соответствующих единицах измерения; а - У-пересечение, т.е. минимально возможное значение переменной У при нулевом значении фактора; Ь - регрессионный коэффициент, знак и величина которого определяет характер и влияние фактора на зависимую переменную. Положительный коэффициент говорит о росте активности ЕРН под влиянием данного фактора, отрицательный - об уменьшении.

Результаты расчетов свидетельствуют об обратной зависимости содержания ЕРН и органики: соответствующие регрессионные коэффициенты отрицательны. Следовательно, основным фактором аккумуляции радиоизотопов является не сорбция, а привнос частиц твердого стока на этапах роста торфяника. Характерно, что для «биофильного» 22&Яа зафиксирована лишь тенденция отрицательной связи с органикой, тогда как для 232ТЪ и 40К содержание минеральных примесей определяет их варьирование на 44-47%. Характерно и то, что величины У-пересечений совершенно не соответствуют средним активностям ЕРН в торфе (даже с поправкой на зольность: табл. 3), но хорошо соотносятся с содержанием радиоизотопов в частицах твердого стока (см. табл. 2: русловой аллювий ручья, глубина от 230 см). В статистически не значимой зависимости для радия прослеживается соответствие с его содержанием в фоновых почвах (слой 0-25 см в табл. 2).

3. Связь удельной активности ЕРН с содержанием органического вещества (в горизонтах торфяной залежи, не подвергшихся гидрогенной минерализации) _по результатам регрессионного анализа__

Элемент Средняя удельная активность в органогенных горизонтах, Бк/кг Уравнение зависимости R2 а

232Th 8,55 Y = 35,09 - 0,439ОВ 0,470 0,089

40K 151,6 Y = 386,9 - 3,894ОВ 0,435 0,100

226Ra* 16,36 Y = 21,8 - 0,090ОВ 0,023 0,748

Примечание. У - зависимая переменная (удельная активность соответствующего радионуклида, Бк/кг); ОВ - содержание неминерализованного органического вещества, %; Я2 - коэффициент детерминации (статистический параметр, характеризующий адекватность модели экспериментальным данным, в долях единицы); а - уровень значимости. * Для радия значимая зависимость отсутствует.

Глубина, см от 232Th 226Ra 40K

поверхности

0-25 25,4 21,2 400,2

25-30 0,2 7,2 121

30-38 3 7,4 101,2

38-51 0,8 1,6 18,5

51-60 3,7 24,7 97,9

60-72 0,1 7,14 0

72-90 7,2 18,6 149

90-106 12,24 10,31 202,2

106-110 0 8,3 45,3

110-136 9,2 20,7 135,4

136-155 4 11,8 96,5

155-172 0,1 1,39 16,6

172-180 3,9 4,9 45

180-220 20,8 17,4 257,8

220-230 3,7 15,4 117,6

свыше 230 30,3 31,9 411,3

4. Радиологическая характеристика кластеров

Второй уровень иерархии Первый уровень иерархии

Кластер Удельная активность, Бк/кг Кластер Удельная активность, Бк/кг

232Th 226Ra 40K 232Th 226Ra 40K

А] 2,47 9,34 80,18 А 1,76 7,24 62,41

А2 0,25 4,61 20,10

Б] 5,70 21,30 127,43 Б 10,63 18,34 168,46

Б2 12,24 10,31 202,20

Бз 20,80 17,40 257,80

В] 30,30 31,90 411,30 В 30,30 31,90 411,30

Регрессионный анализ еще раз подтвердил, что источником ЕРН в низинном торфянике Менек является эрозия в совокупности с русловыми процессами в период его формирования на рубеже плейстоцена и голоцена. При этом последующие процессы гидролиза и минерализации существенно трансформировали формы нахождения ЕРН в торфе, но не повсеместно, а в пределах некоторых генетических горизонтов.

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

Как видно из данных таблицы 4, кластеры А, Б и В объединяют генетически различные горизонты торфяника. Кластер А сформировался под влиянием древнего и современного гидрогенного ожелезнения, а также кислотного гидролиза, и отличается минимальным содержанием ЕРН. Действительно, ламинарные потоки грунтовых вод не могут быть значительным источником радионуклидов, поэтому интенсивное современное обох-ривание торфяника - как и древние аккумуляции сидерита и керченита - наоборот, способствуют «вытеснению» ЕРН соединениями Fe, как подвижного (в глеевых условиях) элемента с гораздо более высоким кларком. Аналогично действует и кислотный гидролиз (кластер А2: слой 106-110 см). При этом кластер А] объединяет «пограничные» горизонты, где ожелезнение только началось. Кластер Б (особенно Б1), напротив, связан с органогенными горизонтами - с относительно пониженной зольностью и с заметной ролью «биофильного» радия. Обособление кластера Б3 в нижней части торфяника обусловлено значительным накоплением терригенных частиц за счет русловых процессов в долине ручья. Но наиболее обособлен кластер В - минеральное ложе торфяника, плейстоцен-голоценовый аллювий с максимальной активностью радионуклидов.

Кроме происхождения генетических горизонтов, кластерный анализ позволяет оценить и контрастность поч-венно-геохимических границ. Результаты оценки свидетельствуют о чрезвычайно высокой контрастности распределения элементов, которая снижается лишь под

влиянием колебаний высоты капиллярной каймы грунтовых вод (относительная размытость кислородного барьера: слои 135-172 см) и привноса терригенных частиц в верхние горизонты (25-38 см). Следовательно, распределение ЕРН, несмотря на их отмечаемую в литературе [5] стабильность и даже инертность в ландшафте, находится в соответствии с геохимическими процессами в торфяной толще.

Таким образом, естественный радиационный фон контролируется почвенно-геохимическими закономерностями, а для некоторых элементов (226Яа) -также и избирательной активностью живого вещества. Кислотный гидролиз и сопряженные процессы при их высокой интенсивности способны приводить к мобилизации и контрастному перераспределению ЕРН в почвах. Нуждается в рассмотрении вопрос о росте водной миграции радиоактивных элементов и изменении их доступности растениям под влиянием кислотного гидролиза, который, как известно [4], активизируется под влиянием текущих климатических изменений. В современных ландшафтах Рязанской области, тем не менее, нет природных и техногенных механизмов сверхнормативной аккумуляции естественных радионуклидов в почвах. В то же время в глинах каменноугольного возраста такое накопление наблюдается, что свидетельствует об актуальности радиоэкологического мониторинга в нашем регионе.

Литература

1. Вернадский В.И. О концентрации радия живыми организмами. - Докл. АН СССР, 1929, №2. - С. 33-34.

2. СП 2.6.1. 758-99. Нормы радиационной безопасности (НРБ-99).

3. Михеева И.В. Изменение пространственной вариабельности свойств почвы при антропогенном воздействии. // Почвоведение. - 1997, № 1. - С. 102-109.

4. Зайдельман Ф.Р. Процесс глееобразования и его роль в формировании почв. - М.: Изд-во МГУ, 1998. - 316 с.

5. Мажайский Ю.А., Тобратов С.А., Дубенок Н.Н., По-жогин Ю.П. Агроэкология техногенно загрязненных ландшафтов. - Смоленск: «Маджента», 2003. 384 с.

6. Пузаченко Ю.Г. Математические методы в экологических и географических исследованиях. - М.: Академия, 2004. - 416 с.