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

Министерство высшего образования

Сыктывкарский государственный университет

Курсовая Работа

по теме: Изучение экологического состояния территории Большеземельской тундры с использованием методов дистанционного мониторинга

Работу выполнил:

студент III курса 234 группы

Шулепов Константин Алексеевич

Научные руководители:

Братцев Андрей Адольфович, к.г.н., зав.отделом геоинформационных и учебных технологий СыктГУ.

Елсаков В.В., к.б.н., зав.отделом экосистемного анализа и ГИС технологий Института биологии Коми НЦ УрО РАН

Сыктывкар 2004

Содержание: стр.

Введение 3

Краткая характеристика физико-географических

условий изучаемой территории 6

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

Обсуждение результатов 13

Основные выводы 16

Литература 17

Введение

Географические информационные системы (ГИС) появились в 1960-х годах как инструмент, позволяющий проводить изучение структурных и функциональных особенностей природных объектов с учетом их пространственной приуроченности. Существует две общераспространенные версии возникновения первых ГИС [8]. Согласно первой, наиболее ранние ГИС были созданы Гарвардском университете и Массачусетском технологическом институте США с целью автоматизированной обработки географической информации. Согласно альтернативной версии – первые ГИС были создана в Канаде и имели цель картирования природных ресурсов (CGIS).

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

Одним из важных источников для ГИС разных уровней (локальных и региональных) являются методы дистанционного зондирования (ДЗ) природных объектов, основанные на использовании электромагнитных излучений, исходящих от предмета исследований и путем их регистрации без непосредственного контакта с ним. Рассматривая ДЗ с позиций системного подхода, необходимо определить входные и выходные элементы системы, ее внутреннюю структуру, границы и окружающую среду. Входными элементами системы являются физические поля, образуемые отражением и/или излучением земной поверхности и естественными процессами в недрах Земли, а также поля техногенного происхождения. Входными элементами системы являются также эмпирические и теоретические закономерности связи физических полей с объектами природной среды. Выходными элементами системы ДЗ следует считать компоненты дистанционной основы карт природоресурсного содержания.

дистанционная основа (ДО) карт определена как оптимальная совокупность материалов ДЗ, результатов их обработок и интерпретации, представленной в цифровом и аналоговом виде. Она состоит из фактографической и интерпретационной частей. Компонентами фактографической части ДО являются нормализованные материалы ДЗ в цифровой и аналоговой формах, а также результаты формализованных преобразований этих материалов. Интерпретационная часть ДО (схемы дешифрирования и интерпретации результатов дешифрирования) создается по результатам экспертного интерактивного анализа изображений и другой информации (Перцов и др., 2000). Наиболее естественным способом система ДЗ подразделяется на следующие три подсистемы: сбор материалов ДЗ, обработка материалов ДЗ и их тематическая интерпретация (Рис.1).

Рис.1. Входные и выходные элементы системы ДЗ.

По способу получения первичных данных дистанционные методы исследования подразделяются на пассивные, т.е. основанные на улавливании излучений от естественных источников (солнца, Луны, звезд, земной поверхности и самих изучаемых объектов), и активные, т.е. предполагающие использование искусственных источников излучения (ламп накаливания, газоразрядных ламп, лазеров). В общеупотребительном смысле термин ДЗ обычно включает в себя регистрацию (запись) электромагнитных излучений посредством различных камер, сканеров, микроволновых приемников, радиолокаторов и других приборов такого рода. Наибольшее применение среди пассивных дистанционных методов получили исследования в оптической области электромагнитного спектра (фотографирование), в том числе различных диапазонах. Достоинство этого метода состоит в том, что фотографические материалы доступны для непосредственного зрительного восприятия и анализа средств используемых в ГИС системах. Космические и аэрофотоснимки обеспечивают территориально полное и непрерывное изучение больших площадей, состояние которых зафиксировано на единый момент времени (Востокова и др, 1988).

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

Наибольшее распространение в настоящее время получили методы ДЗ, основанные на анализе особенностей спектральных характеристик космических снимков разного пространственного разрешения. В России в качестве носителей аппаратуры используются космические аппараты гидрометеорологического (типа «Метеор» и «Электро»), оперативного (типа «Ресурс-О1» и «Океан-О1») и фотографического (типа «Ресурс-Ф») наблюдения [9]. Основными задачами мониторинга, осуществляемого с применением данных средств, являются:

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

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

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

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

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

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

Для реализации поставленной цели были поставлены следующие задачи:

1. Дать характеристику физико-географических условий территории для выявления основных факторов, влияющих на особенности пространственного распределения основных компонентов экосистем в ландшафтах Большеземельской тундры;

2. Собрать и подготовить серию снимков спутника LANDSAT и ASTER для выявления закономерностей распределения спектральных характеристик района исследований;

3. Выявить динамику нарушений растительного и почвенного покрова под влиянием объектов добычи и транспортировки нефти на основании использования разновременных снимков 1998-2000 гг. на участках испытывающих влияние нефтегазового комплекса

4. Подготовить серию изображений, для проведения дешифрирования в течении летнего полевого сезона 2004 г.

Краткая характеристика физико-географических условий изучаемой территории

Большеземельская тундра, холмистая моренная равнина, расположенная в пределах междуречья рек Печора и Уса, Уралом и Пай-Хоем, в административном отношении в пределах Ненецкого автономного округа (Архангельской области) и республики Коми. Преобладающие высоты варьируют в пределах 100-150 м, наибольшая 242 м. низменное, сильно заболоченное приморское побережье к югу повышается террасами, сложенными морскими песками и глинами, и переходит в сильно всхолмленную местность с довольно высокими грядами: Вангурей, Еней, Лыммусюр и др. От Хайпудырской губы к устью Цильмы протягивается возвышенность, называемая Большеземельский хребет, служащий водоразделом рек, впадающих в Баренцево море и в реку.

Климат формируется преимущественно под воздействием арктических и атлантических воздушных масс. Частая смена воздушных масс, перемещение фронтов и связанных с ними циклонов обуславливают неустойчивую погоду. Климат субарктический, с продолжительной холодной зимой (средняя температура января от -16°С на Северо-Западе до -20°С на Юго-Востоке) и коротким прохладным летом (средняя температура июля от +8 до +12°С); в летние месяцы возможны заморозки, средне годовая температура отрицательна. Осадков в год от 450 мм на Юге до 250 мм на Севере. минимум осадков наблюдается, как правило, в феврале, максимум – в августе – сентябре. Не менее 30 % осадков выпадает в виде снега. Избыточное увлажнение, обусловленное низким термическим уровнем в сочетании с равнинным рельефом, слабоводопроницаемыми и многолетнемерзлыми грунтами, определяет обилие поверхностных вод, способствует широкому распространению болот. Часты туманы, от 37 до 72 дней в году. Летом и весной преобладают ветры северных направлений, зимой и осенью – южных. Средняя скорость ветра составляет около 48 м/с. Для климата округа характерны метели до 60 дней в году. Мощность многолетнемерзлые породы (ММП) в Большеземельской тундре изменяется в широком диапазоне и достигает 500 м. Температура ММП в подзоне сплошного распространения изменяется от -5 °С до -2 °С; в местах с не сплошным распространением температура пород выше.

Вегетационный период со среднесуточными температурами свыше +5 °С составляет на юге 95 – 110 дней, на севере 72 – 94 дня. Сумма положительных температур колеблется от 400 градусов на севере до 1100 градусов на юге.

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

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

Реки – большей частью притоки Печоры и Усы - в верховьях текут в узких долинах, ниже долины их расширяются и течение становится спокойным. В верховьях главных рек (Шапкина, Колва, Адзьва и др.) много озёр (Вашуткины, Шапкинские и др.) (Атлас Арктики, 1985; Ненецкий автономный…, 2001; Большая Советская…, 1988).

Широко развито оленеводство, пушной промысел и молочное животноводство. В пределах территории находится часть Печорского угольного бассейна, Тимано-Печорской нефтегазоносной провинции; здесь открыты месторождения нефти и газа. недостаточная геологическая изученность территории, слабая разведанность проявлений рудного и нерудного минерального сырья позволяют выделять в качестве наиболее значимых ограниченный круг полезных ископаемых: нефть, некоторые виды стройматериалов. Другие виды полезных ископаемых требуют дополнительной разведки и изучения для оценки их перспективности.

Восточноевропейские тундры на сегодняшний день в значительной мере сохранили естественный облик вследствие слабого освоения и признаны эталоном равнинных тундр Европы. Тем не менее, на большинстве площадей региона отмечены изменения структуры и состава растительного покрова, связанные с влиянием хозяйственной деятельности различной интенсивности. Среди основных факторов влияния выделим – не всегда контролируемое оленеводство, ведущее к перевыпасу и изменению состава доминирующих видов растений и интенсивное развитие топливно-энергетического комплекса (ТЭК), особо возросшее в последние десятилетия.

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

Модельным участком была выбрана территория, представленная естественными и антропогенно трансформированными сообществами – выпас оленей различной интенсивности (ПСК «Ижемский оленевод»), деятельность объектов нефте– и газодобычи (Средне–Харьягинское месторождение) (рис.1).

Рис. 1. Расположение модельного участка (1) в пределах Большеземельской тундры. Показаны основные месторождения углеводородов ТПНГП, границы ПСК «Ижемский оленевод» и пространственно привязанное изображение снимка landsatETM+

В качестве основы для выполнения работы была подобрана серию снимков спутников ASTER (6.18.2001, 16_009) и landsatETM+ за период с 1995 по 2001 года. Характеристика снимков спутника Landsat приведена в табл.1.

Таблица 1. Характеристики использованных снимков landsatETM+.

Aster - усовершенственный космический тепловой эмиссионный и отражательный радиометр (advanced Spaceborne Thermal Emission and Reflector Radiometer). Датчик установлен на искусственном спутнике Terra, который запущен в декабре 1999 года. ASTER используется, чтобы получить детальные карты температуры поверхности земли, коэфициента излучения, отражательной способности и превышения. Радиометр был построен в Японии для Министерства экономики, торговли и промышленности (Ministry of Economy Trade and Industry – METI). Из всех датчиков, установленных на спутнике TERRA, ASTER – единственный радиометр с высоким пространственным разрешением. Он воспринимает излучения в видимом, ближнем, среднем ИК (или тепловом) диапазонах, размер сцены 60х60 км.

Американский спутник Landsat–7 (спроектирован и создан компанией Lockheed Martin Missiles and Space)был успешно выведен на солнечно–синхронную орбиту с временем нисходящего пересечения экватора 10ч. 00 мин. 15 апреля 1999 г. и имеет расчетный срок эксплуатации 5 лет. Спутник продолжает серию природно–ресурсных спутников Landsat (первый аппарат этой серии был запущен в 1972г.). Миссия Landsat-7 является совместным проектом трех крупнейших американских правительственных организаций: NASA, NOAA и USGS, и призвана обеспечивать национальных и зарубежных потребителей спутниковой информацией высокого разрешения. Установленная на спутнике съемочная аппаратура – сканирующий радиометр ETM+ (Enhanced Thematic Mapper Plus), обеспечивает съемку земной поверхности в шести каналах с разрешением 30 м, в одном ИК канале – с разрешением 60 м и одновременную панхроматическую съемку с разрешением 15 м при ширине полосы обзора для всех каналов около 185 км. Основные параметры орбиты:

Номинальная высота: 705 км;

Периодичность повторения трассы: 16 сут.;

Наклонение: 98.2 град.

Сканирующий радиометр ETM+ создан по контракту NASA компанией Hughes Santa Barbara Remote Sensing. Этот прибор является усовершенствованным вариантом хорошо зарекомендовавших себя сканеров TM (Thematic Mapper), которые работали на борту предыдущих спутников серии Landsat. От предшественников его отличают следующие важные свойства:

· наличие панхроматического канала высокого (15м) разрешения

· наличие теплового ИК-канала

· погрешность абсолютной калибровки – 5%

Представленные радиометры различаются количеством каналов, пространственным разрешением и спектральными диапазонами. Характеристики последних представлены в табл. 2.

Таблица 2. Характеристики спектральных диапазонов радиометров ASTER и ETM+.

В районе исследований растительные сообщества представлены преимущественно кустарниковыми тундрами (Salixlanata, S.phylicifolia), крупноерниковыми (Betulanana) тундрами, приуроченными преимущественно к ложбинам стока и плоскобугристыми болотами с кустарничково-мохово-лишайниковыми сообществами на буграх и осоковыми и пушициево-сфагновыми сообществами на плато водораздела.

Картирование нарушений растительного покрова проводили на основании анализа снимков спутника Landsat (TM и ETM) за 1987, 1988 и 2000 гг., полученные в середине вегетационного периода. Изображения были предварительно обработаны с помощью операции TasseledCap (TC) программного обеспечения ErdasImagine 8.5. Для выявления годичной динамики нарушений растительного покрова модельного участка было составлено комбинированное изображение (рис.2), рассчитанное как:

Красный канал = TC 1 канал 2000 г. – TC 1 канал 1987 г.

Зеленый канал = TC 2 канал 2000 г. – TC 2 канал 1987 г.

Синий канал = TC 3 канал 2000 г. – TC 3 канал 1987 г.

Рис.2. Динамика площадных нарушений на модельной территории

за период 1987 – 2000 гг.

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

Для количественного учета изменений нарушенных площадей изображения каждого года были подвергнуты векторизации (ArcVew 3.2а, UTM, зона 40), для чего была произведена оцифровка антропогенно трансформированных участков. Затем у полигональных объектов измерялась площадь, у линейных – длина, а также толщина линий, соответствующая таковым на снимках (рис.3).

а б

Рис. 3. Изменения площадных нарушений растительного покрова (1:100 000), выявленные на основании анализа космических снимков 1987 г. (а) и 2000 г. (б). Для удобства данные спроецированы на схему района (1:1 000 000), содержащую рельеф, озера,

основные водотоки.

Обсуждение результатов

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

- влиянием вездеходного транспорта;

- влиянием трубопровода;

- влиянием деятельности буровых площадок.

Влияние гусеничного транспорта на растительный и почвенный покров естественных сообществ тундры изучено в работах многих отечественных и зарубежных исследователей (Груздев, Умняхин, 1984; Творогов, Неустроева, 1987; Чалышева 1992). При однократном проезде гусеничной техники удельное давление на грунт составляет порядка 0.47 кгс/см² (Груздев, Умняхин, ), при этом наиболее сильные повреждения испытывает мохово-лишайниковый покров. Устойчивость к повреждениям определяется составом растительности и приуроченностью растительных сообществ к участкам разного рельефа. За три-четыре проезда наиболее легкой машины – ГАЗ–71 происходит полное уничтожение растительного покрова в кустарничково-лишайниковой тундре. При этом моховой слой разрезается гусеницами трактора, отделяется от минерального грунта. Наилучшей устойчивостью обладали ерниковые сообщества. Их уничтожение наблюдали после восьмикратного проезда. На влажной почве нарушения способны проявляться сильнее. Колеи вездеходов на таких участками становятся «резервуарами» для накопления воды и при наличии уклона поверхности они превращаются в водотоки. Усиливается смыв верхнего почвенного горизонта, развивается глубинная эрозия, приводящая к развитию термокарста на участках с мерзлотой. Развитие ускоренной эрозии возможно на площадях с частично или полностью уничтоженной естественной растительностью с уклоном более 1° (Зеликов, 1999). В результате происходит образование делювиальных отложений у подножья склонов или поступление взвешенных частиц почвы, минеральных и органических веществ в водотоки и водоемы. Особо ранима растительность переувлажненных (заболоченных) участков Уже после однократного проезда вездехода образовывается глубокая колея, что также способно привести к развитию термокарстовых процессов.

Влияние трубопроводного транспорта. На представленном участке был использован наземный открытый способ прокладки трубопровода. Нарушения, связанные с трубопроводом связаны с периодом его прокладки, перераспределением снежного покрова вблизи трубопровода и используемого технологического режима при перекачке сырья. Разность температур трубопровода и почвенного слоя может существенно влиять на температурный режим верхних почвенных слоев и мерзлотных слоев. Использование термоизоляции только замедляет процессы образования талых зон или мерзлых ядер: они образуются за три, четыре года при отсутствии изоляции, и через 10–12 лет при ее наличии (Сумина, 1992).

Площадки буровых скважин. Нередко в литературе встречаются указания на то, что нарушения природных экосистем занимают площади гораздо большие, чем это планировалось. По мнению В.П.Гладкова (1989) площадь техногенных нарушений вокруг буровых в тундровой зоне на 9–25% больше, чем в лесотундре. Причины этого следующие.

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

Во-вторых, в проектах порой отсутствуют рекомендации по размещению базы монтажников-строителей. Часто она формируется вне зоны официального отвода, и тем самым площадь нарушений увеличивается на 25–40%. Нередко в процессе монтажа конструкции буровую приходится передвигать на новое место в ходе подготовительных работ.

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

Полученные нами статистические показатели позволили рассчитать площадную структуру выявленных нарушений растительного и почвенного покрова на модельных участках (таблица 4).

Таблица 3. Роль формы рельефа и глубины разрабатываемых скважин на площадь нарушений растительного и почвенного покрова (по В.П.Гладкову, 1989)

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

Из полученных данных видно, что к 2000 году произошло значительное увеличение общей площади нарушенных земель: более чем в два раза увеличились площади, отведенные под буровые скважины, доля линейных объектов (дороги, нефтепровод) возросла примерно в три раза. Для пяти буровых площадок, представленных на изображениях 1987, 1988 и 2000 гг., однозначного роста площадных характеристик выявлено не было. Уменьшение площадей в период 1987–2000 гг. отмечено для площадок № 2, 3, 4, что, вероятно, связано с прекращением или снижением интенсивности работ на данных участках. Дальнейшее сопоставление полученных данных с данными проектной документации может способствовать выявлению несанкционированных вездеходных дорог.

Основные выводы:

1. Использование дистанционных методов в сочетании с полевыми методами исследований в полной мере применимо для экологического картографирования масштаба 1:25 000, 1:100 000 и мельче в целях оценки объема и степени нарушений растительного и почвенного покрова.

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

3. Подготовленные изображения станут основой пространственного анализа состояния естественных и нарушенных экосистем (полевое дешифрирование) которое планируется провести в составе полевого отряда Института биологии Коми НЦ УрО РАН летом 2004.

ЛИТЕРАТУРА

1. Антипов В.С., Астахов В.И., брусничкина н.а., бычкова И.А., Викторов С.В., Вострокнутов Е.П., Гальперов Г.В., Карпузов А.Ф., Кильдюшевский Е.И., Кирсанов А.А., Липияйнен К.Л., Перцов А.В., Рукояткин А.А., Русанова А.А., Смирнова И.О., Старостин В.А., Стрельников С.И., Сухачева Л.Л., Турченко С.И. Аэрокосмические методы геологических исследований. СПб картфабрика ВСЕГЕИ, 2000. С. 15–18.

2. атлас арктики. М.: Главное управление геодезии и картографии при Совете министров СССР, 1985 г. 204 с.

3. Большая советская энциклопедия: В 30-ти т. Т.3. М.: Советская энциклопедия, Востокова Е.А., Сущеня В.А., Шевченко Л.А. М., 1988.

4. Гл. ред. А.М.Прохоров.3–е изд.,1970. 640 с.

5. Гладков В.П. Проектирование и охрана окружающей среды в районах проведения буровых разведочных работ // Труды Коми НЦ УрО АН СССР. № 104.Сыктывкар, 1989. C.6–17.

6. Груздев Б.И., Умняхин А.С. Влияние вездеходного транспорта на растительность Большеземельской тундры // Устойчивость растительности к антропогенным факторам и биорекультивация в условиях Севера. Матер, всесоюз. совещ. Т. 2. Сыктывкар, 1984. С. 19–22.

7. Зеликов В.Д. Почвоведение с основами геологии. Москва, 1999. С.38.

8. Интернет публикация 1: www.gasu.ru, cnit.pgu.serpukhov.su, www.iworld.ru.

9. Интернет публикация 2 www.rosaviakosmos.ru.

10. ненецкий автономный округ. Энциклопедический словарь. М.: Дом Книги «Аванта+», 2001. 304 с.

11. Новаковская Т.В., Акульшина Н.П. Использование геоботанических показателей экологической шкалы для картирования нарушенных земель на Харьягинском нефтегазовом месторождении // Экология. 1997. № 4. С. 256–262.

12. Сахаев В.Г., Щербицкий Б.В. Справочник по охране окружающей среды. Киев: Будiвульник, 1986. С. 33–34.

13. Творогов В.А. Естественное зарастание нарушенных участков тундры в районе Ямбургского газоконденсатного месторождения (полуостров Тазовский) // Бот. журнал, 1988. Т, 73.№11.C. I577–1583.

14. Чалышева Л.В. Особенности формирования растительного покрова техногенных ландшафтов районов нефтедобычи на Европейском Северо-Востоке // Препринт Коми науч. Центра УрО РАН. Сыктывкар, 1992. Вып. 299. 20 с.