CC BY
60
Влияние природных локальных факторов на загрязнение снежного покрова (на примере Печенгского района)
H.Е. Раткин1'1, В.Э. Асминг3, В.В. Кошкин1
1 Институт проблем промышленной экологии Севера КНЦ РАН 1 Естественно-технический факультет МГТУ, кафедра экологии и природопользования
Региональный сейсмический центр КНЦ РАН
Аннотация. Целью работы являлась разработка метода расчета накопления сульфатов, никеля и меди в снеге зоны локального загрязнения комбината "Печенганикель". Решалась основная задача выявления характера формирования и загрязнения снежного покрова на территории, имеющей сложно-пересеченный рельеф. На основе пятилетнего эксперимента была разработана модель аэротехногенного загрязнения снежного покрова в ландшафтах локальной зоны. По данным о годовых выбросах комбината, распределении ветра и сумме осадков по метеостанции "Никель" за зимний период модель обеспечивает:
• представление материала в форме изолинейного отображения современного, ретроспективного и перспективного накопления запаса влаги, сульфатов, никеля и меди в снеге за любой зимний период, а также интегрального накопления загрязняющих веществ и запаса влаги с нарастающим итогом с момента возникновения источника выбросов;
• расчет влагозапаса и накопления загрязнителей на площади водосборных бассейнов и других природно-территориальных комплексов;
• решение обратной задачи: по замерам загрязнения снега в контрольных точках определение выброса сульфатов, никеля и меди от предприятия.
Модель может быть применена в других северных районах Российской Федерации со схожими физико-географическими и климатическими условиями. В основу модели положены зависимости суммы осадков от влагозапасов в однотипных ландшафтах, а также зависимости содержания в снеге загрязнителей от рельефа и расстояния до источника выбросов.
Abstract. The purpose of the work is the method of account of accumulation of sulphate, nickel and copper in the snow of the local pollution zone of the industrial complex "Pechenganikel". The basic task of studying character of formation and pollution of snow cover in the territory having crossed relief has been solved. On the basis of five years' experiment the model of aerotechnogenic pollution of the snow cover in landscapes of the local zone has been developed. According to the data of annual emissions of the industrial complex, distribution of wind and amount of deposits on the meteostation "Nikel" for a winter period the model provides:
• representation of the material in the form of isoline display of modern, retrospective and perspective accumulation of water stock, sulphates, nickel and copper in snow for any winter period, and also integrated accumulation of polluting substances and water stock with an increasing result from the moment of emission source occurrence;
• calculation of water stock and polluter accumulation on the area of drainage pools and other natural and territorial complexes;
• the decision of a reverse task: measuring snow pollution in control points to define emission of sulphate, nickel and copper from the enterprise.
The model can be applied in other northern areas of the Russian Federation with similar physico-geographical and climatic conditions. The basis of the model are dependences of the amount of deposits from water stock in the same landscapes and also dependences of snow polluter contents from a relief and distance up to a source of emissions.
I. Введение
В силу очагового характера размещения промышленных предприятий на территории Кольского полуострова ряд районов, к которым относятся горно-металлургические центры, испытывает исключительно высокую антропогенную нагрузку на природные ландшафты (Лопатина, 1983).
С точки зрения экологических требований к качеству природной среды основными загрязнителями территории Мурманской области являются никель, медь, сернистый газ (Крючков, Макарова, 1989). В этой связи вызывают интерес исследования, направленные на определение количества никеля, меди и серы в природной среде и особенностей их пространственного распределения. В настоящее время при оценке степени загрязнения окружающей среды наибольший интерес представляют природные индикаторы. Одним из современных методов индикации процессов воздушной миграции природных и техногенных веществ является гидрохимическое исследование снежного покрова. Снежный покров обладает рядом свойств, делающих его удобным индикатором загрязнения не только самих атмосферных осадков, но и атмосферного воздуха, а также последующего загрязнения вод и почв (Назаров и др., 1976; Воеводова, 1979). Установлено, что снежный покров может быть использован для решения и более сложных задач: выявления закономерностей пространственного распределения загрязняющих веществ в природной среде, определения вещественного состава и мощности выбросов предприятий, доли вещества, увлекаемого в дальний и локальный перенос (Глазовский и др., 1983; Андрианов и др., 1975; Василенко и др., 1985). Практическая значимость сведений о запасах и динамике накопления снега в ландшафтах в плане контроля загрязнения территорий огромна (Василенко и др., 1985; Израэлъ и др., 1987).
В Мурманской области под промышленными узлами занято около 0.1% ее площади (Крючков, 1979). В этом плане Кольский полуостров может рассматриваться как своеобразная модель для изучения антропогенной нагрузки на весь Крайний Север в целом. С другой стороны, особенности рельефа во многом усложняют решение этой задачи. Проблема оценки характеристик снежности и загрязнения снежного покрова в горных условиях является сложной и многоплановой, поскольку требует учета множества природных факторов: абсолютных отметок, ориентации и крутизны склонов, особенностей ветрового режима, типа подстилающей поверхности и растительности, микрорельефа. (Иверонова, 1956; Соседов, 1967; Северский и др., 1987; Раткин, Макарова, 1992).
Вопросы влияния природных факторов на распределение, режим и загрязнение снежного покрова в горах изучены недостаточно. Это связано с недостатком фактической информации из-за сложности постановки многолетних и весьма трудоемких натурных наблюдений. Поэтому особую актуальность приобретает разработка методов расчета показателей снежности и загрязнения снежной толщи.
2. Материалы и методы
Район исследований охватывает небольшую часть территории на северо-западе Мурманской области и часть территории Норвегии пограничной с Россией по реке Пасвик.
В работе использованы результаты гидрохимического опробования снежного покрова в импактной зоне комбината "Печенганикель" за 1973-1994 гг., данные метеостанции "Никель" за 19401994 гг., статистические данные о выбросах загрязняющих веществ комбинатом "Печенганикель" за 1979, 1983, 1990-1994 гг. Пробы снега отбирались в конце марта - начале апреля на постоянной сети мониторинга (рис.1), отражающей типичные ландшафты исследуемой территории: равнины, склоны и вершины, с разным типом растительного покрова: хвойный лес, лиственный лес и безлесные кустарничково-мохово-лишайниковые ландшафты. Работы по отбору, подготовке и химическому анализу проб проведены по методикам, применяемым в Институте проблем промышленной экологии Севера (ИППЭС) КНЦ РАН и многими другими исследователями (Василенко и др., 1985; Глазовский и др., 1983; Крючков и др., 1989; Раткин и др., 1992). В марте 1994 года ИППЭС провел контрольную проверку качества аналитической информации, участвуя в совместном эксперименте с экологическими службами Норвегии и Финляндии. Получено 97, 97 и 100% правильных результатов соответственно. Обе зарубежные лаборатории аккредитованы в соответствии с общеевропейскими стандартами обеспечения качества измерений серии БМ 45000, и присланные ими результаты могут рассматриваться как арбитражные. Для изучения влагозапаса в точках мониторинга в 1990-1992 гг. одновременно с отбором проб снега производились промеры высоты снежного покрова и определение его плотности на площадках размером 150x150 м в узлах равномерной сети с шагом 25 м (рис.2). Отбор проб снега на гидрохимический анализ осуществлялся в вершинах равносторонних треугольников, расположенных в центре этой площадки.
Расчеты произведены на основе созданной базы данных ЭГИС, в которую введены следующие показатели: рельеф, тип растительности, абсолютные отметки, углы наклонов склонов, полярные координаты точек мониторинга (азимут, расстояние), результаты гидрохимического опробования снежного покрова, высота и плотность снега, данные об осадках, выбросах серы и количестве суток за зимний сезон, в течение которых ось факела была направлена от источника выброса загрязняющих
веществ на снегопункт. Программное обеспечение работы произведено программистами В.В.Кошкиным иВ.Э. Асмингом.
Рис.1. Расположение снегопунктов.
точки отбора проб снега на гидрохимический анализ
точки определения высоты и плотности снежного покрова для определения влагозапаса
Рис.2. Схема снегопункта.
3. Результаты и обсуждение
В результате исследований было изучено влияние природных локальных факторов на формирование влагозапасов и загрязнение снежного покрова на территории, характеризующейся сложно-пересеченным рельефом. Установлено, что для района исследований характерны однородные климатические условия во времени с "мягким" климатом и умеренным ветровым режимом. Осадки на Кольском полуострове главным образом связаны с прохождением циклонов и имеют фронтальное происхождение. Ширина фронтов достигает 200-250 км (Яковлев, 1961). Поэтому данные по осадкам на метеостанции "Никель" являются репрезентативными для всех точек мониторинга. Изменение запаса влаги по сезонам в реперных точках происходит согласованно с изменением сумм осадков, выпавших за зимний период на метеостанции "Никель": коэффициенты вариации запасов влаги в ландшафтах и сумм осадков близки (21 и 23% соответственно). На основании этого для каждого типа ландшафта за каждый сезон были рассчитаны ряды коэффициентов, полученные как отношение запаса влаги в реперных точках к сумме осадков на метеостанции "Никель". Произведен их сравнительный анализ с помощью ЭВМ в программном пакете 8ТАТОЯАР. Различие коэффициентов оценивалось по нуль-гипотезе через ^критерий Стьюдента при уровне значимости 0.05.
В результате произведенного анализа было выявлено, что во времени значения коэффициентов статистически не различаются. Влияние форм рельефа проявляется не значимо, а видового состава растительности - значимо. На основе этого получены средние достоверные при 95% доверительном уровне значения коэффициентов, статистически значимо различающиеся по растительному признаку (табл.1).
Таблица 1. Средние статистически достоверные коэффициенты отношения запаса влаги к сумме осадков за зиму по метеостанции "Никель"
Тип рельефа Тип растительного покрова
Без леса Хвойный лес Лиственный лес
Равнина 0.66 0.80 1.03
Склон 0.69 0.83 1.15
Вершина 0.66 0.84 1.08
Средние 0.67 0.82 1.09
Полученные коэффициенты позволяют рассчитывать запас влаги в точках мониторинга по формуле:
^ = БКг, (1)
где: - запас влаги /-ой точки, мм;
£ - сумма осадков на метеостанции "Никель", мм;
К - эмпирический коэффициент, выбираемый в зависимости от типа растительного покрова в /-ой точке.
Расчеты содержания загрязняющих веществ в снежном покрове проводились с учетом особенностей примеси и источника, направления ветрового переноса, рельефа местности в следующей последовательности:
1. Определение на основе данных о повторяемости направлений ветра и штилей на метеостанции "Никель" общей продолжительности распределения ветра за ноябрь - март в сутках по восьми основным румбам и во время штилей (табл.2).
Таблица 2. Распределение ветра и штилей за зимний период, сутки
Румб/год С СВ В ЮВ Ю ЮЗ 3 СЗ Штиль
1979 4.7 7.6 8.6 11.5 50.0 22.8 11.8 5.6 28.4
1983 3.6 5.6 5.1 8.7 55.5 30.3 13.6 5.9 22.7
1990 4.7 3.9 8.1 32.2 35.8 25.4 13.9 6.8 20.2
1991 3.3 4.5 7.0 30.6 52.6 24.5 10.0 6.7 11.8
1992 5.6 4.8 3.9 31.1 52.4 27.7 8.1 6.4 10.9
1993 2.9 3.5 5.3 30.2 51.3 21.6 10.7 8.0 17.5
1994 5.3 5.8 5.8 19.8 49.8 22.2 9.3 5.9 27.0
2. Разбиение точек мониторинга по трем морфологическим классам ландшафта: вершина, склон, равнина, и составление для этих классов рядов времени экспозиции факела за сезоны исследований по направлению от источника на точку мониторинга.
3. Расчет среднесуточной концентрации веществ за сезон в точках мониторинга по формуле:
Чсут.1 = 4/ / Рг (2)
где: 4сут.1 - среднесуточная концентрация за сезон в /-ой точке, мг/л. сутки; 4 - измеренная концентрация вещества в /-ой точке, мг/л;
Рг - продолжительность направленности факела от источника на /-ую точку, сутки.
В результате за каждый сезон получены ряды пространственного распределения среднесуточной концентрации загрязнителей при выбросах, показанных в табл.3.
Таблица 3. Выбросы вредных веществ комбинатом "Печенганикель" (тыс. т/год)
1979 1983 1990 1991 1992 1993 1994
802 399 293 273 258 251 230 200
N1 0.44 0.38 0.51 0.40 0.38 0.35 0.39
Си 0.35 0.29 0.31 0.27 0.24 0.23 0.25
4. Для учета выбросов и приведения рядов к одной единице измерения произведен расчет удельной среднесуточной концентрации веществ в точках мониторинга по формуле:
Цудл. Цсутл / V (3)
где: 4удЛ - удельная среднесуточная концентрация в /-ой точке, мг/л. сутки на 1 тыс. т. выброса в год; 4 сутл - среднесуточная концентрация за сезон в /-ой точке, мг/л. сутки; V - годовой выброс веществ, тыс. т.
Произведенный расчет по формуле (3) позволил объединить пять статистических рядов данных за отдельные сезоны в один ряд, отдельно для каждого элемента по трем формам рельефа. На основе полученных рядов определено выборочное уравнение связи и расстояния от источника. Уравнения определялись методом простой регрессии в программном пакете 8ТАТвКАР. Определен общий вид связи:
4 удл = а хЬ , (4)
где: 4УдЛ - удельная среднесуточная концентрация в /-ой точке, мг/л. сутки на 1 тыс. т выброса в год; а== \ogai - свободный член уравнения регрессии /-того элемента; х/ - расстояние от источника /-ой точки, км; Ь - коэффициент регрессии /-того элемента.
Уравнения связи и критерии оценки линий регрессии приведены в табл.4-6.
Таблица 4. Уравнения связи ду¿. и расстояния от источника для Б042- в разных формах рельефа и критерии оценки линий регрессии
Тип рельефа Функция Критерии оценки
Число п Б-критерий К кор. Р я2
Вершина дш = -5.4Х°9 78 82 -0.73 0.000 0.54
Склон = -5.68Х-063 88 87 -0.72 0.000 0.52
Равнина дш = -6.18Х/0'59 80 60 -0.69 0.000 0.48
Таблица 5. Уравнения связи и расстояния от источника для М в разных формах рельефа
и критерии оценки линий регрессии
Тип рельефа Функция Критерии оценки
Число п Б-критерий К кор. Р я2
Вершина Яш = -1.51ХГ1'17 78 117 -0.80 0.000 0.64
Склон Цуд. / = -1.52ХГ133 88 137 -0.79 0.000 0.62
Равнина дш = -1.39ХГ145 80 98 -0.76 0.000 0.57
Таблица 6. Уравнения связи дуд. и расстояния от источника для Си в разных формах рельефа
и критерии оценки линий регрессии
Тип
Функция
К
ритерии оценки
рельефа Число п Б-критерий К кор. Р я2
Вершина дш = -1.35ХГ1'18 78 127 -0.80 0.000 0.64
Склон дш = -1.14Х11'44 88 138 -0.80 0.000 0.64
Равнина дуд. г = -1.31ХГ1."2 80 103 -0.77 0.000 0.59
Из формул (2), (3) и (4) выводится общая формула расчета полной концентрации загрязнителей в точке мониторинга на момент максимального влагозапаса:
д, = аХЫ- Р, •V , (5)
где: д, - концентрация за зиму в /-ой точке, мг/л;
аХгЫ = дуьл - удельная концентрация в /-ой точке, мг/л. сутки на 1 тыс. т выброса в год; Р, - время экспозиции факела на /-ую точку, сутки; V- годовой выброс, тыс. т.
Сопоставление расчетов по формуле (5) с экспериментальными данными показало, что рассчитанные концентрации сульфатов, никеля и меди достаточно хорошо согласуются с величинами, полученными в результате примененной методики отбора и обработки образцов снега. В среднем за 1993-1994 гг. 61% всех точек по сульфатам (по никелю и меди - 55%) имеют различия расчетной и измеренной концентрации в 10^25%. Количество точек, имеющих различия >100% невелико (табл.7).
Таблица 7. Число точек от общего количества (%) по градациям различия (%) измеренной и расчетной концентрации
Вещество Градации
10-25 30-45 50-90 >100
бо42- 61 6 14 19
N 55 13 14 18
Си 55 14 18 13
Известно, что
Я, = д, • Ж, (6)
где: Qi - содержание веществ в ,'-ой точке, мг/м2;
дг - концентрация веществ на момент максимального снегонакопления в ,'-ой точке, мг/л;
Ж, - запас влаги в ,'-ой точке, мм.
Подставляя в (6) выражение д (5) и Ж (1), получаем общую формулу перехода от концентрации к накоплению веществ в снежном покрове:
Я, = Р, -V- дудГ 8 К,, (7)
где: Qi - накопление веществ в снежном покрове за зиму в ,'-ой точке, мг/м2;
Р, - продолжительность направленности факела от источника на ,'-ую точку, сутки;
V - годовой выброс вещества, тыс. т;
д^ = аХ!" ,-ой точки, мг/л. сутки при V = 1 тыс. т в год;
8 - сумма осадков за зиму, мм;
К, - поправка на тип растительного покрова в ,'-ой точке.
Полученные уравнения (5) и (7) могут быть использованы для оценки аэротехногенного загрязнения территории исследований и контроля за выбросами через определение концентраций веществ в снежном покрове на контрольных пунктах мониторинга, расположенных в зоне локального загрязнения.
Разработан блок создания цифровой карты нагрузок в изолинейном отображении и масштабировании при выводе на печать, картографического отображения ареалов водосборных бассейнов и расчета накопления веществ на их площади. На рис.3 показано суммарное накопление сульфатов, никеля и меди в снежном покрове за период 1973-1994 гг. На рис.5 показаны ареалы водосборных бассейнов главных рек Печенгского района и их первых притоков. Интегральные нагрузки запаса влаги, сульфатов, никеля и меди на площади водосборных бассейнов показаны в табл.8.
Известно, что распространение примесей в приземном слое воздуха существенно зависит от скорости ветра и коэффициента турбулентного обмена (Соломатина, 1965). В свою очередь коэффициент турбулентности зависит от характера подстилающей поверхности, степени её увлажнения, температурных условий и радиационного баланса территории. Естественно, что в районе исследований коэффициенты турбулентного обмена холодного и тёплого периодов будут разными. Различается и вымывающая
способность жидких и твердых осадков. Это даёт основание предположить, что в тёплый период года траектории полёта частиц будут отличаться от тех траекторий, которые описываются функциями зимнего периода. Этот вопрос требует дальнейшего изучения. Однако, учитывая, что средние многолетние скорости ветра за год, а также холодного и тёплого периодов, практически одинаковы, соответственно 3.8, 3.9 и 3.8 м/с (Научно-прикладной справочник, 1988), то с определённой степенью допущения функции распределения , приведённые в таблицах 4-6, на современной стадии изученности можно использовать для приблизительной оценки годового аэротехногенного загрязнения рассматриваемой территории. Тогда расчет годовой концентрации можно произвести по формуле:
5._5_1
Рис.3. Накопление загрязнителей в снеге, мг/м2 . Рис.4. Выпадение загрязнителей за период, г/м2 .
Вверху - 804 (1973-94 гг.), в середине - N1 (1975-94 гг.), внизу - Си (1979-94 гг.).
Рпс.5. Ареалы водосборных бассейнов.
Таблица 8. Интегральные нагрузки запаса влаги (ш, млн. т), сульфатов, никеля и меди на площади водосборных бассейнов, т
N 8, км2 Накопление(период)
8042-(1973-94) N1 (1975-94) Си (1979-94)
1 182.7 548.1 1670 6.9 5.9
2 50.8 132.1 706 3.7 2.3
3 200.3 520.8 1902 1 3.8 6.4
4 30.8 80.1 715 4.6 2.2
5 98.8 321.1 1535 14.9 9.2
6 1530.1 4590.3 51864 573.4 413.0
7 57.8 187.8 1601 11.4 8.4
8 134.0 435.5 2492.4 20.1 10.0
9 75.8 227.4 2657 39.7 27.4
10 189.6 663.6 2301 41.7 29.9
11 237.7 832.0 3390 54.7 30.9
12 120.6 392.0 2290 25.9 13.3
13 114.8 321.4 1505 40.6 14.4
14 174.0 487.2 6947 200.9 152.7
15 136.1 442.3 8670 347.8 283.8
16 278.5 905.1 5759 105.5 55.2
17 771.1 2159.1 4548 72.5 44.0
18 448.8 1256.6 7840 73.3 47.9
19 73.1 255.8 1031 8.0 3.4
20 228.6 617.2 1626 11.3 5.9
21 210.9 569.4 722 4.7 2.1
22 65.6 200.1 549 3.8 1.5
23 143.6 524.1 659 5.5 2.3
24 261.9 654.7 771 4.6 2.1
25 85.3 255.9 462 1.2 0.8
26 230.1 690.3 783 3.5 2.1
год = Цудл*Р, год*У , (8)
где: ЪудЛгод - годовая концентрация веществ в /-ой точке, мг/л;
= аХ!" - удельная концентрация веществ зимнего периода в /-ой точке, мг/л. сутки при выбросе 1 тыс. т в год;
Р, год - время экспозиции факела за год на ,-ую точку, сутки; V- годовой выброс, тыс. т.
Тогда уравнение для расчета накопления примет следующий вид:
Яг год = Р, год^^.П^Ш^], (9)
где: Qiгос) - годовое накопление в /-ой точке, мг/м2;
Р, год - время экспозиции факела за год на ,-ую точку, сутки; V- годовой выброс, тыс. т;
Чудл - удельная концентрация зимнего периода в ,-ой точке, мг/л при выбросе 1 тыс. т в год; £ - сумма осадков за зиму, мм;
К, - поправка на тип растительного покрова в ,-ой точке; £ - сумма жидких осадков за теплый период (апрель - октябрь), мм.
По формуле (9) рассчитаны годовые интегральные нагрузки сульфатов, N1 и Си на территорию за период 1973-1994 гг., показанные в изолинейном отображении на рис.4.
4. Заключение
В результате проведенных исследований было изучено влияние природных локальных факторов на формирование влагозапасов и загрязнение снежного покрова на территории, характеризующейся сложно-пересеченным рельефом. Установлено, что району характерны однородные климатические условия во времени. В условиях спокойного и умеренного ветрового режима значимое влияние на режим и распределение снежного покрова оказывает тип растительности. Влияние форм рельефа второстепенно. Изменение запаса влаги в ландшафтах по сезонам происходит согласованно с изменением сумм осадков на метеостанции "Никель". На основании этого получены постоянные поправочные коэффициенты на тип растительного покрова, позволяющие определять запас влаги в ландшафтах через сумму осадков за ноябрь - март на метеостанции "Никель".
Получены выборочные уравнения связи удельной концентрации сульфатов, никеля и меди в снежном покрове на момент максимального снегонакопления и расстояния от источника для трех форм рельефа: вершина, склон, равнина. Общий вид связи выражен степенной функцией. Все расчеты
произведены на основе созданной базы данных ЭГИС, в которую вошли основные параметры, влияющие на режим, распределение и загрязнение снежного покрова. Показан подход к расчёту годового поступления загрязняющих веществ на подстилающую поверхность. Разработано программное обеспечение, позволяющее картографически моделировать аэротехногенное загрязнение территории в изолинейном изображении и рассчитывать нагрузку на любых природно-территориальных выделах.
Авторы благодарят Советника губернатора Финнмарка по вопросам охраны окружающей среды Пера-Енара Фескебека за оказание финансовой и технической помощи при проведении полевых работ. Авторы также выражают благодарность Л.М. Терещенко и Е.А. Леоновой за выполнение анализов проб снега, В.П. Манжуле, С.А. Рязанцевой, Л.К. Врублевской, Е.Р. Крепсу, Д.Ю. Смирнову за помощь при проведении снегосъемки.
Литература
Андрианов А.Н., Дроздова В.М. Исследования химического состава снега вокруг г. Ленинграда. В кн.:
Атмосферная диффузия и загрязнение воздуха. Л., Гидрометеоиздат, с.208-212, 1975. Василенко В.Н., Назаров И.М., Фридман Ш.Д. Мониторинг загрязнения снежного покрова. Л.,
Гидрометеоиздат, 181с., 1985. Воеводова З.И. Возможность определения влияния загрязнения атмосферы на водные ресурсы путем отбора проб снега. В кн.: Влияние хозяйственной деятельности человека на водные ресурсы Коми АССР. Сыктывкар, с.80-88, 1979. Глазовский Н.Ф., Злобина А.И., Учватов В.П. Химический состав снежного покрова некоторых районов Верхнеобского бассейна. В кн.: Региональный экологический мониторинг. М., Наука, с.67-83, 1983.
Иверонова М.И. Закономерности распределения снежного покрова на склонах и лесолуговом поясе хребта Терскей Алатау. М., Труды Института географии АН СССР. Т.67. Работы Тянь-Шаньской физико-географической станции, вып.5, c.127-137, 1956. Израэль Ю.А., Назаров И.М., Фридман Ш.Д. Мониторинг трансграничного переноса загрязняющих
воздух веществ. Л., Гидрометеоиздат, 303с., 1987. Крючков В.В. Север: природа и человек. Перспективы освоения. М., "Наука", 127с., 1979. Крючков В.В., Макарова Т.Д. Аэротехногенное воздействие на экосистемы Кольского Севера.
Апатиты, изд. КНЦАНСССР, 96с., 1989. Лопатина Е.Б. Пространственный анализ интегральной антропогенной нагрузки. Региональный
экологический мониторинг. М., Гидрометеоиздат, с.24-28, 1983. Назаров И.М., Ренне О.С., Фридман Ш.Д., Шаповалова Л.Г., Махонько Э.П. Содержание примесей в атмосферных осадках, атмосферные аэрозоли. В кн.: Защита атмосферы от загрязнений. Вильнюс, АН Лит ССР, Институт физики и математики, вып.3, с.7-11, 1976. Научно-прикладной справочник по климату СССР. Серия 3, многолетние данные 4.1-6, вып.2,
Мурманская область. Л., Гидрометеоиздат, 316с., 1988. Обзоры состояния загрязнения атмосферного воздуха и выброс вредных веществ в атмосферу на
территории деятельности Мурманского УГКС за 1979-1994 гг. Мурманск, МУГКС, с.185, 1994. Раткин Н.Е., Макарова Т.Д. Роль снежного покрова в загрязнении ландшафтов Мурманской области. В кн.: Эколого-географические проблемы Кольского Севера. Апатиты, изд-во КНЦ РАН, с.20-35, 1992.
Северский И.В., Северский С.И. Влияние локальных природных факторов на распределение снежного
покрова в горах. В кн.: Снежный покров в горах и лавины. М., Наука, с.16-23, 1987. Соломатина И.И. О влиянии рельефа на метеорологические характеристики в приземном слое воздуха. Труды ГГО: Вопросы атмосферной диффузии и загрязнения воздуха. Л., Гидрометеоиздат, вып.172, с.58-69, 1965.
Соседов И.С. Исследование баланса снеговой влаги на горных склонах. Алма-Ата, Наука, 197с., 1967. Яковлев Б.А. Климат Мурманской области. Мурманск, 180с., 1961.
