CC BY
54
УДК 551.510.42 © 2012: В.Л. Макухин, В.Л. Потемкин; ФНИ «XXI век»
МОДЕЛИРОВАНИЕ ПЕРЕНОСА И ТРАНСФОРМАЦИИ ЗАГРЯЗНЯЮЩИХ ПРИМЕСЕЙ, В ТОМ ЧИСЛЕ РТУТИ, НА АКВАТОРИИ ОЗ. БАЙКАЛ
В.Л. Макухин, В.Л. Потемкин
Лимнологический институт СО РАН, Иркутск, Россия
Эл. почта: [email protected] Статья получена редакцией 05.07.2012, принята к печати 29.08.2012
Выполнен обзор работ по моделированию атмосферных процессов в регионе оз. Байкал. Построена математическая модель распространения и осаждения ртути в регионе Южного Байкала — объекта Всемирного природного наследия ЮНЕСКО. Источниками выбросов ^ в этом регионе являются промышленные объекты городов Усолье-Сибирское, Ангарск, Шелехов, Иркутск, Слюдянка, Байкальск. Рассчитанные величины концентраций ^ оказались соответствующими по порядку величины данным немногочисленных измерений концентраций ртути в этом регионе. С использованием рассчитанных полей концентраций ртути выполнены оценки интенсивности осаждения ртути на подстилающую поверхность. Результаты по порядку величины соответствуют данным по содержанию ртути в снеговом покрове на опорных станциях региона. Выбросы ^ в атмосферу предприятиями Слюдянки и Байкальска составляют 72%, Иркутско-Черемховского промышленного комплекса — 28%. Менее значимый вклад этого комплекса объясняется удаленностью от озера и орографическими неоднородностями, препятствующими переносу примесей. Оценки влияния рельефа местности на перенос ртути в Прибайкалье показали, что при северо-западном ветре со скоростью 2 м/с Приморский хребет и Олхинское плато задерживают примерно третью часть от массы поступившей на озеро. В радиусе 5 км от одиночного источника оседает примерно 10% валового выброса ртути, что соответствует результатам других исследователей.
Ключевые слова: математическое моделирование, атмосферные процессы, загрязнение атмосферы, оседание на поверхность, ртуть, водные бассейны, оз. Байкал.
MODELING OF THE TRANSFER AND TRANSFORMATION OF POLLUTANTS, INCLUDING MERCURY, IN THE LAKE BAYKAL BASIN V.L. Makukhin, V.L. Potiomkin
Limnological Institute of the Siberian Branch of the Russian Academy of Science, Irkutsk, Russia
E-mail: [email protected]
Approaches to the modeling of atmospheric processes in the lake Baykal region are reviewed. A mathematical model of mercury transfer and deposition in South Baykal region, which is rated by UNESCO as a world heritage object, is developed. Hg is discharged in this region by industrial objects located in the cities Usol'ye-Sibirskoye, Angarsk, Shelekhov, Irkutsk, Sliudianka and Baykalsk. The calculated Hg concentrations were found to correspond, by the order of magnitude, to data obtained by direct Hg determinations, which are scarce in this region. The calculated Hg concentration fields were used to estimate the rates of Hg deposition onto underlying surfaces. The results correspond, by the order of magnitude, to data about Hg content in the snow mantle at the reference stations located in the region. Hg discharges into the air by industrial objects related to Sliudianka and Baikalsk are responsible for 72% of the total air pollution with Hg, and by those related to the Irkutst-Cheremkhovask industrial complex, for 28%. The lesser importance of the latter complex is caused by its relative remoteness from the lake and by orographic unevennesses making obstacles to pollutant transfer. The estimates of the effects of relief on Hg transfer in the region suggest that, when northwest windspeed is 2 m/sec, the Primorskiy Ridge and Olkhinskoye Plateau retain about on third of the amount of Hg that reaches the lake basin. About 10% of a total Hg discharge is deposited within 5 km from its source, in accordance with estimates made by other authors.
Keywords: mathematical modeling, atmospheric processes, air pollution, mercury, pollutant sedimentation, water basins, the lake Baykal.
Введение
Развитие энергетики, машиностроения, химии, транспорта в XX в. привело к тому, что человеческая деятельность стала сравнима по масштабам с естественными энергетическими и материальными процессами, происходящими в биосфере. Антропогенные воздействия приводят к нарушению практически всех природных биогеохимических циклов, в том числе включающих в себя тяжелые металлы.
К наиболее опасным загрязнителям относят ртуть и ее соединения. По степени воздействия на организм ртуть относится к первому классу опасности (вещества чрезвычайно опасные). Несмотря на это, ртуть и ее соединения широко применяются в металлургии, электротехнике, медицине, в производстве целлюлозно-бумажной продукции, ацетальдегида, хлорвинила, красок, хлора и каустика, при изготовлении детонаторов, управляемых снарядов и в сельском хозяйстве. За счет
хозяйственной деятельности человека в окружающую среду ежегодно поступает до 10 000 т ртути [22]. Основными техногенными источниками ее эмиссии являются три отрасли промышленности: комбинаты цветной металлургии, получающие данный металл из ртутьсо-держащих пород; предприятия по добыче и переработке руд различных металлов, где ртуть выступает как сопутствующий компонент; предприятия химической и электрохимической промышленности. Загрязнение окружающей среды металлами и другими примесями происходит и при сжигании жидкого и твердого топлива, а также при использовании ртути для экстрагирования и концентрирования золота [47, 51]. Важно установить, как и в каких количествах распространяется ртуть в атмосфере и осаждается на подстилающую поверхность региона оз. Байкал, участка Всемирного природного наследия ЮНЕСКО.
В настоящее время для описания распространения и трансформации атмосферных примесей и выявления причинно-следственных связей широко используются математические модели, позволяющие оценить возможные последствия тех или иных воздействий на окружающую среду. Использовались модели и при исследованиях процессов загрязнения региона оз. Байкал. Предпринимались попытки оценки антропогенных потоков ряда микроэлементов из атмосферы на зеркало озера [5] на основе простых балансных соотношений и аналитических решений [48]. Для оценки вклада 15 основных промышленных источников, расположенных в регионе Байкала, в загрязнение северной части озера пылью и тяжелыми металлами [6] использовалась простая траекторная модель [19, 21, 43, 49].
Трансформация соединений серы и азота над оз. Байкал исследовалась в [23] с помощью аналитических моделей, балансовой с элементами статистики [42, 50] и двумерной стационарной диффузионной [41]. Существенным недостатком этих моделей, не учитывающих влияние рельефа, является то, что вертикальное распределение примесей предполагается равномерным.
С помощью модели эйлерова типа Г.И. Марчуком, К.Я. Кондратьевым, В.В. Пененко, А.Е. Алояном была предпринята попытка оценить влияние антропогенного загрязнения пылью, окислами серы и азота на оз. Байкал и окружающие его территории при различных метеорологических ситуациях [27, 28, 30, 33]. А.Е. Алоян, В.А. Загайнов, А.А. Лушников, С.В. Макаренко решали эту же задачу с учетом эволюции атмосферных аэрозолей за счет коагуляции [3]. Под коагуляцией авторы работы [3] подразумевают процесс роста аэрозольных частиц в результате слияния друг с другом. Для описания изменения концентрации и спектра размеров частиц для пространственно однородного случая использовано уравнение Смолуховского. Отметим, что все численные эксперименты в этих работах проведены с шагом по горизонтальным координатам 25 км, что больше средней ширины хребтов в регионе озера. Г.И. Марчук, А.Е. Алоян, В.Н. Пискунов в работах [2, 4, 29] с помощью этого же комплекса моделей выполнили исследования распространения и трансформации аэрозольных частиц от пяти крупных источников Приан-гарья и Южного Прибайкалья с учетом кинетических процессов конденсации, испарения и коагуляции, причем шаг сетки по горизонтали был равен 10 км.
Расчеты проводились в следующей последовательности. По модели гидротермодинамики были получены
характеристики атмосферной циркуляции в рассматриваемом регионе и вычислены поле течения и турбулентные характеристики. На фоне полученной атмосферной циркуляции решалась задача переноса газовых примесей с учетом фотохимической трансформации, в результате чего образуются пары серной кислоты, а также другие вторичные загрязнители. Далее рассматривались процессы образования нуклеационной моды, после чего была смоделирована динамика формирования аэрозольных частиц с учетом кинетических процессов конденсации, испарения и коагуляции. Модель имеет следующие достоинства. Характеристики атмосферной циркуляции рассчитываются с помощью модели гидротермодинамики. Учтена фотохимическая трансформация примесей, кинетические процессы конденсации, испарения и коагуляции. Однако по расчетам по модели А.Е. Алояна на высоте 250 м наиболее загрязняется средняя часть оз. Байкал, хотя по данным инструментальных измерений наименее загрязнены средняя и северная части озера [38, 39, 40, 45]. Наблюдаются области повышенного загрязнения сульфатами в Южном Забайкалье, занятом достаточно высокими хребтами (высота отдельных вершин достигает 1700-1800 м), а данных измерений, что на вершинах хребтов повышенные концентрации сульфатов, не имеется. Складывается впечатление, что при расчетах в работах А.Е. Алояна шаг сетки выбран большим, и влияние горных хребтов, окружающих оз. Байкал, учтено слабо.
Для оценки вклада конкретных промышленных источников в загрязнение атмосферы и подстилающей поверхности района Южного Байкала А.В. Аргучин-цевой использован стохастический подход описания климатических характеристик в виде многомерных функций плотности вероятности реализации конкретного метеокомплекса, замыкающих дифференциальное уравнение переноса и турбулентной диффузии загрязняющего вещества [12, 13]. В работах [8, 14-16] А.В. Ар -гучинцевой, В.К. Аргучинцевым, С.Ж. Вологжиной рассмотрен более общий подход к стохастическому моделированию. Основа модели базировалась на прямом (втором) уравнении Колмогорова для переходных вероятностей, записанном в фазовых координатах. С помощью модели А.В. Аргучинцевой, С.Ж. Вологжиной оценивалось влияние стационарных источников загрязняющих веществ, расположенных в Байкальской котловине, на особо охраняемые природные территории [17].
Методы исследования
Процессы распространения и трансформации примесей в регионе оз. Байкал изучались с помощью нелинейной нестационарной пространственной модели, основанной на численном решении полуэмпирического уравнения турбулентной диффузии примеси [7, 10]. В качестве исходной используется система уравнений в следующем виде:
3
^+ V■gradS-Ж = У— К — -ЛБ+ F (1)
д ^ 5 £ д х3 д х1 ' д х1
Здесь ( - время; V = (u,v,w) - вектор скорости; u,v -горизонтальные и w - вертикальная компоненты вектора скорости движения воздуха вдоль осей декартовой прямоугольной системы координат (х1, х2, х3); S = - тензор массовых концентраций исследуемых примесей (/ = 1,...,и); К, К , К3 - коэффициенты турбулентной диффузии по осям х , х , х3 соответственно;
А(х,^) = {А к(х.^)} - матричный оператор, описывающий взаимодействие различных су бстанций между собой и их локальные изменения; г = 1,2,3, к = 1,..., L; L=nx(n2+ п +1);
,() - вектор-функция, описывающая источники примесей; W ={W } - тензор скоростей гравитационного оседания субстанций; п — количество субстанций в многокомпонентной среде.
Имея в виду решение задачи о распространении примесей над региональной областью, предполагаем фоновое распределение концентраций субстанций известным. Из-за отсутствия детальной информации из наблюдений в качестве начальных условий приняты Б , равные фоновому распределению, а при его отсутствии Б. = 0.
В качестве краевых условий приняты
д Б.
_I=о
д х1 при х1=0,Х
аз.
д х.
д х-
-=0 при х2=0,Х
■ = 0 при х3=0,Н,
(2)
О
В — + GS = дf
д f
+ о1 F1,
где
G _\0 для хеD, 1 0 для х е D1.
F1 =
№ для хеD, 0 для хе D
1
Задача (1)-(2) решается конечно-разностным методом. Решение ищется в области D, представляющей собой прямоугольный параллелепипед, на сетке с основными узловыми точками х., у zk и неравномер-
ными шагами:
Ах. = х. - х
/-г
I = 1,...,1 + 1,
Ау = у) - У-1, ]' = 1,..^ + 1,
Агк = ^ - ^-р
к = \,..,К + 1.
где х1 = 0, х1 = Х, х2 = 0, х2 = У, х3 = Н - границы области счета. На уровне подстилающей поверхности поставлено граничное условие, учитывающее отражение и поглощение примесей в зависимости от свойств подстилающей поверхности.
Уравнение турбулентной диффузии примеси (1) вместе с начальными и граничными условиями (2) численно интегрируется в декартовой системе координат с применением метода фиктивных областей [25]. Метод фиктивных областей основан на том, что в среде с относительно большими значениями коэффициента диффузии изменение плотности диффундирующего вещества относительно мало. Суть метода состоит в дополнении фактической области интегрирования со сложной геометрией фиктивными областями до более удобной для численного интегрирования формы и в доопределении системы уравнений модели специальными условиями ее продолжения в фиктивных областях. Запишем уравнение (1) в операторной форме
Используются также полуцелые индексы для обозначения вспомогательных точек х(.+1/2), у(.+1/2), z(k+1/2), являющихся серединами основных интервалов. Граничные поверхности будут представлять собой
Х1/2 = -Х, Х1+1/2 = Х, У1/2= У;+1/2 = ^ ^/2 = 0, ^+1/2 = Н. Конечно-разностные аппроксимации производных
по пространственным переменным построены на ос -нове интегро-интерполяционного метода [25]. Дифференциальные операторы первого порядка по пространственным переменным аппроксимировались центральными разностями, а операторы второго порядка - трехточечными разностными операторами.
При решении задачи используется метод расщепле -ния по физическим процессам [26]. На каждом малом интервале времени рассматривается схема,
состоящая из двух этапов. На первом учитываются процессы переноса примесей по траекториям и турбулентной диффузии, на втором - локальные преобразования и влияние источников. Аппроксимация задачи по времени на первом этапе построена с помощью двуциклического полного расщепления [26].
1
т+—
(Е+ 0,5 АЛ 1) Б 7 = ( Е- 0,5 АЛ 1) Бт,
±2 Л
т+— т+—
77
(Е+ 0,5 АЛ 2) Б 7 = ( Е- 0,5 АЛ 2) Б
2
т+ —
(Е+ 0,5 АЛ ) Б 7 = ( Е- 0,5 АЛ ) Б 7,
'3
4 3
т+— т+—
Б 7 = Б 7 + 2Д №"
5
4
т+ —
(Е+ 0,5 АЛ ) Б 7 = ( Е- 0,5 АЛ ) Б 7,
л
т+ — 7
(Е+ 0,5 АЛ 2) Б 7 = ( Е- 0,5 АЛ 2) Б
т+ —
7
(3)
где функция Б определена в основной области Д Дополним область D до прямоугольного параллелепипеда и рассмотрим в нем задачу:
Через D1 обозначено дополнение основной области D до прямоугольного параллелепипеда. Введение таких областей позволяет вести расчеты с произвольным рельефом, обеспечивая более универсальное применение модели.
, т+ —
(Е+ 0,5 АЛ 1) 5т+1 =( Е-0,5 АЛ 1) Б 7,
где Л1, Л2, Л3 - разностные операторы, аппроксимирующие дифференциальный оператор задачи (3); г = 1,...,1, ]' = 1,...Д к = 1,...,К, Е - единичная матрица.
Конечно-разностные уравнения реализуются немонотонной прогонкой. Используемая неявная схема покомпонентного расщепления дает решение для некоммутативных операторов со вторым порядком аппроксимации по времени и координатам. На гладких решениях устойчивость схемы зависит от устойчивости немонотонной прогонки.
На втором этапе при решении системы уравнений химической кинетики применена монотонная полунеявная схема первого порядка аппроксимации [37]. Сравнение численных решений уравнения турбулентной диффузии примеси с известными аналитическими решениями нестационарного [31] и стационарных уравнений распространения примесей [18] показало, что они согласуются с аналитическими решениями с удовлетворительной степенью точно -
3
сти, как на высоте источника, так и у поверхности земли [7].
Для верификации математической модели [7, 10] была проведена серия численных экспериментов по данным экспедиционных наблюдений летом 1992 г. в регионе оз. Байкал [11, 44]. В рассмотрение были включены промышленные объекты городов Иркутск, Ангарск, Усолье-Сибирское, Черемхово, Зима, Шелехов, Байкальск. Получено, что коэффициенты корреляции между расчетными и измеренными концентрациями имеют значения: 0,7 для сульфатов и 0,8 для нитратов. Относительные ошибки не превосходили для сульфатов 40%, а для нитратов 50%, при этом средние квадратические отклонения этих ошибок составляли соответственно 30% и 25%. Наибольшие отклонения рассчитанных значений концентраций от экспериментальных связаны с влиянием удаленных источников.
Выполненные с помощью модели расчеты (с шагом по горизонтали 5 км) позволили оценить средние годовые выпадения сульфатной серы и нитратного азота в южной части Байкала [9, 46], хорошо согласующиеся с данными наблюдений сети контроля снежного покрова [32].
Результаты
Целью данной работы было исследование процессов распространения и осаждения ртути в регионе
Южного Байкала с помощью математической модели распространения и трансформации примесей [7, 10]. Модель успешно применялась ранее, например, при оценке влияния метеорологических процес -сов на загрязнение атмосферы в регионе оз. Байкал [24]. Основными источниками выбросов ^ в этом регионе являются промышленные объекты городов Усолье-Сибирское, Ангарск, Шелехов, Иркутск, Слюдянка, Байкальск. Данные об интенсивности источников были взяты из работ [20, 34, 35].
Моделирование процессов распространения примесей проводилось в области площадью 200x200 км2 и высотой 4 км над поверхностью оз. Байкал. Шаги по времени и горизонтали составляли соответственно 150 с и 1 км; шаг по вертикали задавался следующим образом: до высоты 350 м он равнялся 50 м, далее - 150 м, 500, 1000 и 2000 м. Коэффициенты турбулентной диффузии рассчитывались с использованием соотношений полуэмпирической теории турбулентности [10].
На рис. 1 представлено распределение концентраций ртути при северо-западном ветре со скоростью 2 м/с. Примесь по долине Ангары выносится к оз. Байкал. Проведено сравнение рассчитанных вели -чин концентраций с данными немногочисленных измерений концентраций ртути в этом регионе [1], по порядку величины наблюдается удовлетворительное соответствие (табл. 1).
Рис. 1. Изолинии рассчитанных приземных концентраций ртути при северо-западном ветре в регионе Южного Байкала, в нг/м3.
Табл. 1
Рассчитанные и измеренные концентрации ртути в г. Усолье-Сибирское
Станция измерения Концентрация, нг/м3
измеренная рассчитанная
Усолье-Сибирское 22-300 10-250
С использованием рассчитанных полей концентраций ртути выполнены оценки интенсивности осаждения ^ на подстилающую поверхность региона Южного Байкала (рис. 2). Сравнение с данными по содержанию ртути в снеговом покрове на опорных станциях Южного Прибайкалья [20] также показало их удовлетворительное по порядку величины соответствие (табл. 2).
Табл. 2
Рассчитанные и измеренные интенсивности осаждения ртути в Приангарье
Значительно меньше влияние Иркутско-Черемхов-ского промышленного комплекса - 28% (табл. 3). Менее значимый вклад этого комплекса объясняется удаленностью источников выбросов от озера и наличием орографических неоднородностей, препятствующих переносу примесей.
Табл. 3
Вклад отдельных источников выбросов ртути в загрязнение Южного Байкала
Источник выбросов Вклад, %
Усолье-Сибирское 20
Ангарск 1
Шелехов 1
Иркутск 6
Слюдянка 42
Байкальск 30
Станция измерения Интенсивность осаждения, г/(км2*месяц)
измеренная рассчитанная
Усолье-Сибирское 0,7-8,2 2-15
Ангарск 1-1,2 0,4-2
Оценен вклад предприятий Приангарья и Прибайкалья в загрязнение Южного Байкала ртутью при атмосферных выбросах. Получено, что наибольшее влияние на озеро оказывают выбросы предприятий Слюдянки и Байкальска, их вклад равен 72%.
Выполнены оценки влияния рельефа местности на перенос ртути в Прибайкалье. Для этого были выполнены расчеты с учетом рельефа местности и без него. Получено, что при северо-западном ветре со скоростью 2 м/с Приморский хребет и Олхинское плато задерживают примерно третью часть от величины массы поступившей на озеро.
Представляет интерес оценить, какая доля валового выброса оседает вблизи источника и на удалении от него. Такие оценки были выполнены с помощью модельных расчетов. Получено, например, что в радиусе 5 км от Слюдянки оседает примерно 10% валового выброса ртути, в радиусе 80 км - 40% (табл. 4). Эти
Рис. 2. Изолинии плотности массового расхода ртути в регионе Южного Байкала, в г/(км2 X месяц).
оценки справедливы и для других источников выбро -сов: Усолья-Сибирского, Ангарска, Иркутска, Шеле-хова и Байкальска и удовлетворительно согласуются с результатами других исследователей [20, 36].
Табл. 4
Доля валового выброса Hg, оседающая вокруг организованного источника выбросов (Слюдянка) на площади радиуса R
R, км Доля валового выброса, %
5 10
10 15
20 20
40 30
80 40
Заключение
Сравнение рассчитанных по модели и измеренных характеристик распространения ртути показало их удовлетворительное количественное соответствие.
Наибольшее влияние на оз. Байкал оказывают выбросы предприятий Слюдянки и Байкальска. Значительно меньше влияние Иркутско-Черемховского промышленного комплекса из-за удаленности источников выбросов от озера и наличия орографических неодно -родностей. При северо-западном ветре со скоростью 2 м/с Приморский хребет и Олхинское плато задерживают третью часть ртути от попавшей на озеро.
В радиусе 5 км от одиночного источника оседает примерно 10% валового выброса ртути, что соответствует результатам других исследователей.
Литература
1. Алехин Ю.В., Лапицкий С.А., Мухамади-ярова Р.В., Пухов В.В. Новые результаты исследования отдельных составляющих геохимического цикла ртути // Вестник Отделения наук о Земле РАН. - 2007. - № 1(25). - С. 1-4.
2. Алоян А.Е. Моделирование динамики и кинетики газовых примесей и аэрозолей в атмосфере. - М. : Наука, 2008. - 415 с.
3. Алоян А.Е., Загайнов В.А., Лушни-ков А.А., Макаренко С.В. Перенос трансфор-мирующегося аэрозоля в атмосфере // Изв. АН СССР. Сер. Физика атмосферы и океана. - 1991. - Т. 27. - С. 1232-1240.
4. Алоян А.Е., Пискунов В.Н. Моделирование региональной динамики газовых примесей и аэрозолей // Изв. РАН. Физика атмосферы и океана. - 2005. - Т. 41. - С. 328-340.
5. Анохин Ю.А., Остромогильский А.Х., Пословин А.Л., Хицкая Е.В. Оценка антропогенного потока микроэлементов из атмосфе -ры на зеркало оз. Байкал // Проблемы экологического мониторинга и моделирования экосистем. Т. IV. - 1981. - С. 32-40.
6. Анохин Ю.А., Кокорин А.О., Прохорова Т.А., Анисимов М.П. Аэрозольное загрязнение атмосферы над озером Байкал и влияние на него промышленных источников // Мониторинг состояния озера Байкал. - Л. : Гидрометеоиздат, 1991. - С. 44-50.
7. Аргучинцев В.К. Численное моделирование распространения аэрозолей в пограничном слое атмосферы // Оптика атмосферы и океана. - 1994. - Т. 7. - С. 1106-1111.
8. Аргучинцев В.К., Аргучинцева А.В. Моделирование мезомасштабных гидротермодинамических процессов и переноса антропогенных примесей в атмосфере и гидросфере региона оз. Байкал. - Иркутск : Изд-во Иркутского гос. ун-та, 2007. - 255 с.
9. Аргучинцев В.К., Куценогий К.П., Маку-хин В.Л. и др. Экспериментальное исследование и численное моделирование аэрозолей и газовых примесей в атмосфере Южного Байкала // Оптика атмосферы и океана. -1997. - Т. 10. - С. 598-604.
10. Аргучинцев В.К., Макухин В.Л. Математическое моделирование распространения
аэрозолей и газовых примесей в пограничном слое атмосферы // Оптика атмосферы и океана. - 1996. - Т. 9. - С. 804-814.
11. Аргучинцев В.К., Макухин В.Л., Обол-кин В.А. и др. Исследование распределения соединений серы и азота в приводном слое оз. Байкал // Оптика атмосферы и океана. -1996. - Т. 9. - С. 748-754.
12. Аргучинцева А.В. Математическое моделирование распределения антропогенных аэрозолей // Оптика атмосферы и океана. -1996. - Т. 9. - С. 800-803.
13. Аргучинцева А.В. Климатическое распределение загрязняющих веществ от Се-ленгинского целлюлозно-картонного комбината (СЦКК) // Оптика атмосферы и океана. - 1997. - Т. 10. - С. 605-609.
14. Аргучинцева А.В. О вероятностном подходе к моделям экологического районирования и рационального природопользования // Оптика атмосферы и океана. - 1998. - Т. 11. -С. 606-609.
15. Аргучинцева А.В. Вероятностный подход к моделированию задач рационального природопользования // Оптика атмосферы и океана. - 1999. - Т. 12. - С. 499-502.
16. Аргучинцева А.В., Аргучинцев В.К., Вологжина С.Ж. Тенденции и уровень решения проблемы оценки загрязнения атмосферы // Известия ИГУ. Сер. Науки о Земле. - 2009. - Т. 2(2). - С. 20-36.
17. Аргучинцева А.В., Вологжина С.Ж. Загрязнение атмосферного воздуха Байкальской котловины // Известия ИГУ. Сер. Науки о Земле. - 2011. - Т. 4(2). - С. 28-41.
18. Берлянд М.Е. Современные проблемы атмосферной диффузии и загрязнения атмосферы. - Л. : Гидрометеоиздат, 1975. - 448 с.
19. Гальперин М.В. Модель для расчета дальнего трансграничного переноса соединений серы в атмосфере (выпадения и концентрации) // Тр. ИПГ. - 1988. - Вып. 71. - С. 9-14.
20. О состоянии окружающей природной среды Иркутской области в 2007 году. Государственный доклад. - Иркутск : Мин-во природ. ресурсов и экологии Ирк. обл., 2008. - 354 с.
21. Израэль Ю.А., Назаров И.М., Прес-сман А.Я., Ровинский Ф.Я., Рябошапко А.Г. Кислотные дожди. - Л. : Гидрометеоиздат, 1989. - 270 с.
22. Краснопеева И.Ю. Распространение ртути и ее соединений в окружающей среде и влияние на организм человека // Сибирский медицинский журнал. - 2005. -Т. 54(5). - С. 7-12.
23. Кудрявцева Л.В., Устинова С.И. Оценка вклада дальнего переноса соединений серы и азота в их поступление в оз. Байкал // Мониторинг и оценка состояния Байкала и Прибайкалья. - Л. : Гидрометеоиздат,
1991. - С. 86-92.
24. Латышева И.В., Макухин В.Л., Потемкин В. Л. Исследование характеристик Азиатского максимума и его влияния на загрязнение атмосферы в регионе оз. Байкал // Оптика атмосферы и океана. - 2005. -Т. 18. - С. 466-469.
25. Марчук Г.И. Методы вычислительной математики. - М. : Наука, 1980. - 534 с.
26. Марчук Г.И. Математическое моделирование в проблеме окружающей среды. -М. : Наука, 1982. - 320 с.
27. Марчук Г.И., Алоян А.Е. Математическое моделирование в задачах экологии. - М. : ОВМ АН СССР. Препринт № 234, 1989. - 36 с.
28. Марчук Г.И., Алоян А.Е. Математическое моделирование в задачах окружающей среды // Проблемы механики и некоторые современные аспекты науки. - М. : Наука, 1993. - С. 12-25.
29. Марчук Г.И., Алоян А.Е. Динамика и кинетика газовых примесей и аэрозолей в атмосфере и их значение для биосферы // Биосфера. - 2009. - Т. 1. - С. 48-57.
30. Марчук Г.И., Кондратьев К.Я. Приоритеты глобальной экологии. - М. : Наука,
1992. - 264 с.
31. Монин А.С., Яглом А.М. Статистическая гидромеханика. Ч. 1. - М. : Наука, 1965. - 639 с.
32. Обзор фонового состояния окружающей природной среды в СССР за 1987 г. -М. : Гидрометеоиздат, 1988. - 97 с.
33. Пененко В.В., Алоян А.Е. Математические модели взаимосвязей между термодинамическими и химическими процессами в атмосфере промышленных регионов // Изв. АН. Сер. ФАО. - 1995. - Т. 31. - С. 372-384.
34. Петросян В.С. Глобальное загрязнение окружающей среды ртутью и ее соединениями // Россия в окружающем мире: 2006 (Аналитический ежегодник) / Отв. ред. Н.Н. Марфенин. - М. : МНЭПУ Авант., 2007. - С. 149-163.
35. Ревич Б.А. «Горячие точки» химического загрязнения окружающей среды и здоро -вье населения в городах России // Россия в окружающем мире: 2006 (Аналитический ежегодник) / Отв. ред. Н.Н. Марфенин. - М. : МНЭПУ Авант, 2007. - С. 108-148.
36. Ртуть: обзор экотоксикологических свойств и некоторых промышленно-экологи-ческих проблем. - М. : Эколайн, 1997. - 137 с.
37. Тимофеева С.С., Седов Д.В., Макухин В.Л., Малыхин А.В. Загрязнение атмосферы залповыми выбросами при пожарах и авариях на промышленных объектах Байкальского региона: оценка, прогноз, превентивные меры. - Иркутск : Изд-во ИрГТУ, 2009. - 192 с.
38. Ходжер Т.В., Потемкин В.Л., Оболкин
B.А. Химический состав аэрозоля и малые газовые примеси в атмосфере над Байкалом // Оптика атмосферы и океана. - 1994. -Т. 7. - С. 1059-1066.
39. Ходжер Т.В., Оболкин В.А., Потемкин В.Л. О роли атмосферы в формировании химического состава вод оз. Байкал // Оптика атмосферы и океана. - 1999. - Т. 12. -
C. 512-515.
40. Bashurova V.S., Dreiling V., Hodger T.V. et al. Measurements of atmospheric condensation nuclei size distributions in Siberia // J. Aerosol Sci. - 1992. - Vol. 23. - P. 191-199.
41. Fay J.A., Rosenzweig J.J. An analytical diffusion model for long distance transport of air pollutions // Atmosph. Environ. - 1980. -Vol. 14. - P. 355-365.
42. Fisher B.E.A. A review of the processes and models of long-rang transport of air pollutants // Atmosph. Environ. - 1983. - Vol. 17. -P. 1865-1880.
43. Giorgi F. A particle dry-deposition parameterization scheme for use in the trace transport models // J. Geophys. Res. D. - 1986. -Vol. 91. - P. 9794-9806.
44. Khodger T.V., Obolkin V.A., Arguchint-sev V.K. et al. Studies of aerosols spreading in the Southern Pribaikalye atmosphere // J. Aerosol Sci. - 1995. - Vol. 26 (Suppl. 1). - P. S461.
45. Koutsenogii P.K., Bufetov N.S., Drozdo-va V.I. et al. Ion composition of atmospheric aerosol near lake Baikal // Atmosph. Environ. -1993. - Vol. 27 A. - P. 1629.
46. Makukhin V.L., Potemkin V.L. Numerical simulation of sulfur and nitrogen compounds distribution near the lake Baikal // J. Aerosol Sci. - 2000. - Vol. 31. - Suppl. 1. -P. S424-S425.
47. Pai P., Heisler S., Joshi A. An emissions inventory for regional atmospheric modeling of mercury // Water Air, Soil Pollut. - 1998. -Vol. 101. - P. 289-308.
48. Rodhe H. A study of sulphur budget for the atmosphere over Northern Europe // Te-llus. - 1972. - Vol. 24(2). - P. 128-137.
49. Sievering H. Small particles dry deposition on natural waters: How large the incertainity? // Atmosph. Environ. - 1984. -Vol. 18. - P. 2271-2272.
50. Venkatram A. Statistical long-range transport models // Atmosph. Environ. - 1986. -Vol. 20. - P. 1317-1324.
51. Vieira J.L.F., Hassarelli M.M. Determi-nagao de mercurio total em amostras de agua, sedimento e solidos em suspensao de corpos aguaticos por espectrofotometria de absorgao atomica com gerador de Vapor a frio // Rev. saude publ. - 1996. - Vol. 30(3). - P. 260.
