CC BY
12
УДК 502.3+528.854.2+551.511
ОЦЕНКА ПАРАМЕТРОВ ВЫБРОСОВ НОВОСИБИРСКИХ ТЭЦ С ИСПОЛЬЗОВАНИЕМ СПУТНИКОВОЙ ИНФОРМАЦИИ
Владимир Федотович Рапута
Институт вычислительной математики и математической геофизики СО РАН, б30090, Россия, г. Новосибирск, пр. Академика Лаврентьева, б, доктор физико-математических наук, главный научный сотрудник, тел. (383)330-61-51, e-mail: [email protected]
Анатолий Александрович Леженин
Институт вычислительной математики и математической геофизики СО РАН, 630090, Россия, г. Новосибирск, пр. Академика Лаврентьева, б, кандидат физико-математических наук, старший научный сотрудник, тел. (383)330-64-50, e-mail: [email protected]
Татьяна Владимировна Ярославцева
ФБУН «Новосибирский НИИ гигиены» Федеральной службы по надзору в сфере защиты прав потребителей и благополучия человека, б30108, Россия, г. Новосибирск, ул. Пархоменко, 7, кандидат технических наук, старший научный сотрудник, тел. (383)330-61-51, e-mail: [email protected]
На спутниковых снимках территорий городов и промышленных площадок непеременным атрибутом являются дымовые шлейфы, которые характеризуют источники выбросов и процессы распространения примесей в атмосфере. Цель работы состоит в использовании этой информации для оценки параметров источников и характеристик пограничного слоя атмосферы. Предложенный подход базируется на уравнениях гидротермодинамики. Для восстановления вертикального профиля ветра в пограничном слое используется уравнения Эк-мана и данные аэрологического зондирования. Эффективная высота источника рассчитывается по высоте солнца над горизонтом и положением тени дымовой струи на спутниковом снимке. Контроль точности оценивания дополнительной высоты подъема дымового факела проводится с помощью известных соотношений, включающих в себя динамический импульс выброса и тепловой поток плавучести. По зимним спутниковым снимкам г. Новосибирска проведены оценки эффективной высоты подъема факелов для крупных ТЭЦ города. Представлены результаты численного восстановления вертикальных профилей ветра по дымовым шлейфам. Предлагаемые подходы могут быть использованы для анализа процессов загрязнения атмосферы городов.
Ключевые слова: пограничный слой атмосферы, метеорологические условия, дымовой шлейф, эффективная высота источника, поток плавучести, спутниковые снимки, обратная задача, численное моделирование.
ESTIMATING THE EMISSION PARAMETERS OF NOVOSIBIRSK THERMAL POWER STATIONS WITH THE USE OF SATELLITE INFORMATION
Vladimir F. Raputa
Institute of the Computational Mathematics and Mathematical Geophysics SB RAS, 6, Prospect Аkademik Lavrentiev St., Novosibirsk, 630090, Russia, D. Sc., Chief Researcher, phone: (383)330-61-51, e-mail: [email protected]
Anatoly A. Lezhenin
Institute of the Computational Mathematics and Mathematical Geophysics SB RAS, 6, Prospect Аkademik Lavrentiev St., Novosibirsk, 630090, Russia, Ph. D., Senior Researcher, phone: (383)330-64-50, e-mail: [email protected]
Tatyana V. Yaroslavtseva
FBSI «Novosibirsk scientific research institute of hygiene» of Rospotrebnadzor, 7, Parhomenko St., Novosibirsk, 630108, Russia, Ph. D., Senior Researcher, phone: (383)330-61-51, e-mail: [email protected]
Satellite images of cities and industrial platforms show plumes of smoke which characterize sources of emissions and transport of pollutions in the atmosphere. The purpose of this work is the use of this information for assessment of emission sources and characteristics of the atmospheric boundary layer. An approach based on the equations of atmospheric dynamics is proposed. For the calculation of a vertical wind's profile in the atmospheric boundary layer Ekman's equations and data of aerological sounding are used. The effective source height is calculated using the sun's angle over the horizon and the projection of the shadow of the smoke plume in the satellite image. Assessment of the accuracy of estimation of the additional height of plumes of smoke is carried out by means of some ratios including the dynamic impulse of emissions of pollutions and the buoyancy heat flux. Using winter satellite images of Novosibirsk of effective heights of a plumes of smoke for large combined heat and power plants of the city are estimated. The results of calculations of vertical profiles of wind are presented. This approach can be used for the analysis of the processes of pollution in the urban atmosphere.
Key words: atmospheric boundary layer, meteorological conditions, plumes of smoke, effective height, flux of buoyancy, satellite images, inverse problem, numerical modeling.
Введение
Экспериментальные и теоретические исследования атмосферного загрязнения играют важную роль в решении проблем окружающей среды. Они необходимы для проведения оценок риска здоровья населения, воздействия на биосферу и климатическую систему [1-5].
Дымовые шлейфы природных пожаров, выбросов промышленных предприятий и крупных ТЭЦ хорошо видны на снимках из космоса [6-8]. Траектории дымовых факелов фиксируют направления и скорости ветра на высотах выбросов. Возможности же визуальных наблюдений за распространением дымовых загрязнений с помощью наземных средств весьма ограничены. Использование спутниковых снимков позволяет прослеживать дымовые шлейфы до нескольких десятков километров. Такая картина особенно часто наблюдается в зимний период времени. Протяженность и ширина дымовых шлейфов может изменяться в больших диапазонах [2, 4].
Модели, базирующиеся на численном решении уравнений гидротермодинамики атмосферы, позволяют получить детальное описание процессов переноса примеси в локальном и региональном масштабах [9 - 14]. Для проведения расчетов требуется задание значительных объемов входной информации, включающей сведения о метеорологических условиях и параметрах источников выбросов примеси [15 - 18]. При этом основные характеристики источника такие, как интенсивность выброса и эффективная высота подъема дымового шлейфа требуют большей точности определения. Также важна информация о вертикальном профиле ветра. Для получения этих сведений нами используются как расчетные методы, так и аэрологические, метеорологические наблюдения и спутниковые снимки. Решение поставленных задач требует построения соот-
ветствующих математических моделей оценивания параметров с использованием данных с различных систем наблюдений состояния атмосферы и процессов переноса и диффузии примеси [19]. В практическом плане это приведет к повышению точности расчетов полей загрязнения.
Объекты и методы исследования
Материалами исследований служили зимние спутниковые снимки территории г. Новосибирска. На них хорошо видны траектории дымовых выбросов от труб ТЭЦ-2, ТЭЦ-3, ТЭЦ-4, ТЭЦ-5. Высота основных труб ТЭЦ-2, ТЭЦ-3 составляет 120 м, их диаметр равен 8,2 м. Высота и диаметр устья большой трубы ТЭЦ-5 равны 260 м и 10,8 м соответственно.
На рис. 1 представлен снимок г. Новосибирска с ИСЗ «Ресурс-П» №1 за 13 февраля 2018 г. на 12 часов местного времени, полученный в Сибирском центре ФГБУ «Научно-исследовательский центр космической гидрометеорологии «ПЛАНЕТА» (http://www.rcpod.ru). Предварительный анализ показывает, что направления выносов примеси от разновысоких труб ТЭЦ заметно различаются. Дымовые шлейфы от труб высотных труб ТЭЦ-5, сносятся в западно-северо-западном направлении. Для более низких источников ТЭЦ-2, ТЭЦ-3, ТЭЦ-4 вынос примеси происходит в секторе запад-юго-запад.
Рис. 1. Спутниковый снимок г. Новосибирска от 13 февраля 2018 г. на 12 часов местного времени с ИСЗ «Ресурс-П» № 1. 1 - ТЭЦ-2, 2 - ТЭЦ-3, 3 - ТЭЦ-4, 4 - ТЭЦ-5
Отметим, что при слабом ветре и устойчивой температурной стратификации в нижней атмосфере эффективная высота подъема газоаэрозольной примеси, выбрасываемой из труб ТЭЦ-5, может достигать несколько сотен метров. По сути эти дымовые струи расположены в верхней части зимнего пограничного слоя атмосферы. На таких уровнях наблюдается правый поворот ветра с высотой [20]. Это подтверждается анализом рис. 1 и данными, приведенными в табл. 1, 2. Относительно направления приземного ветра на близкорасположенной к г. Новосибирску метеостанции Огурцово (индекс станции 29638) угол правого поворота дымовых факелов от труб ТЭЦ-5 составляет не менее 35 - 40 градусов. Для анализа метеорологических условий в течение 13 февраля 2018 г. использованы данные аэрологической станции г. Новосибирска. Станция расположена на высоте 143 м. над уровнем моря. В табл. 1, 2 представлены распределения по высоте основных метеорологических величин на рассматриваемый день: давления, температуры, направления и скорости ветра в 07 и 19 часов местного времени.
Таблица 1
Метеорологические параметры в нижней атмосфере по данным аэрологической станции г. Новосибирска 07 часов местного времени 13 февраля 2018 г.
Высота над Давление, Температура, Направление Скорость
уровнем моря, м ИРа °С ветра, град. ветра, м/с
143 1008 -15,3 0 0
203 1000 -15,1 40 2
473 965 -14,9 47 3
737 932 -13,9 54 3
794 925 -14,1 55 3
1436 850 -14,7 70 4
Таблица 2
Метеорологические параметры в нижней атмосфере по данным аэрологической станции г. Новосибирска на 19 часов местного времени 13 февраля 2018 г.
Высота над уровнем Давление, Температура, Направление Скорость
моря, м ИРа °С ветра, град. ветра, м/с
143 1008 -11,5 30 2
208 1000 -11,7 35 3
528 959 -13,1 54 5
747 932 -10,3 67 7
805 925 -10,5 70 7
1152 884 -11.8 80 10
1452 850 -12,9 85 11
Из табл. 1, 2 следует, что в течение дня в пограничном слое атмосферы температурная стратификация была близка к нейтральной и наблюдались слабые ветра северо-восточного направления.
Для расчетов полей концентрации примесей в атмосфере необходимо задание параметров источников выбросов и текущих метеорологических условий [15, 21]. Основными характеристиками источника являются мощность эмиссии, его геометрическая высота, диаметр трубы, температура и скорость исходящих газов. Существенное влияние на процессы распространения примеси оказывают метеорологические условия: распределение с высотой температуры, влажности воздуха и скорости ветра [7, 9, 15-18].
Для вычисления эффективной высоты источника и вертикального профиля ветра существует множество подходов [15, 21]. Фактической высотой источника выброса следует считать его эффективную высоту
Кф = Кр + АН, (1)
где Нтр - геометрическая высота трубы, Ан - добавка к высоте за счет скорости
выброса примеси и теплового потока плавучести.
При стратификации близкой нейтральной расчет Ан производится по формуле [15]
ан = а^д, + ъФг, Ф = (2)
и и3 4Та
где g - ускорение силы тяжести, WQ - скорость выброса примеси в атмосферу, и - скорость ветра, Б0- внутренний диаметр трубы, Ф - поток плавучести, АТ - превышение температуры исходящих газов над температурой атмосферы Та.
При измерении скорости ветра на уровне флюгера: а = 1,9, Ъ = 4,95 [15].
Для расчета вертикального профиля ветра в пограничном слое используем следующую модель [22]:
дк— + ¡(V _ у ) = 0 , ^к— _ ¡{и _ ия) (3)
дг ^ дг дг g W
Здесь и, V - горизонтальные компоненты скорости ветра в направлении осей х, у соответственно, ось г направлена вертикально вверх, I - параметр Кориолиса, к - коэффициент вертикального турбулентного обмена, и , Vg -
горизонтальные составляющие скорости геострофического ветра. Для определения коэффициента вертикального турбулентного обмена применяются различные схемы замыкания [12, 13, 18-20]. Для случая постоянного коэффициента турбулентности в работе [22] предложен метод восстановления компонентов скорости ветра в пограничном слое атмосферы по данным спутниковых наблюдений распространения дымовых струй от разновысоких источников.
Результаты
Расчет дополнительной высоты подъема дымового шлейфа выброса ТЭЦ проводился двумя способами. В первом случае использовался спутниковый снимок (рис. 1). По углу высоты Солнца и измеренным проекциям дымовых шлейфов от труб ТЭЦ-3 и ТЭЦ-5 на поверхность земли были определены эффективные высоты выбросов примеси. В рассматриваемом случае высота
Солнца над горизонтом достигала 19° 03измеренная ширина тени от шлейфа ТЭЦ-3 составила примерно 1000 м, от шлейфа ТЭЦ-5 - 1500 м. Рассчитанные по этим значениям эффективные высоты подъема дымовых струй равны: 344 м для ТЭЦ-3, 516 м для ТЭЦ-5.
По формулам (1), (2) были выполнены расчеты дополнительных высот подъема дымовых шлейфов от труб ТЭЦ-3, ТЭЦ-5. Исходные данные и результаты расчетов для 6 м/с представлены в табл. 3. Скорость ветра на высоте флюгера составляла 3 м/с (метеостанция Огурцово).
Таблица 3
Оценки дополнительных высот подъема дымовых шлейфов ТЭЦ-3 и ТЭЦ-5
ТЭЦ-3 ТЭЦ-5
ат ,°с Ф , м4/с3 Ан, м АТ ,°С Ф , м4/с3 Ан, м
150 571 136 150 848 194
200 760 171 200 1130 245
Численное восстановление вертикального профиля ветра проводилось с использованием модели (3) [22]. На рис. 2 приведены вычисленные компоненты скорости ветра.
01234567
м/с
Рис. 2. Восстановленные с помощью модели (4) составляющие скорости ветра и , V
Геострофический ветер задавался 6 м/с. Угол между направлением приземного ветра и направлением дымового факела от трубы ТЭЦ-5 был равен 40°.
Обсуждение
Многочисленные примеры фиксации на космических снимках дымовых выбросов в атмосферу от антропогенных и природных источников указывают на возможности широкого использования этой информации для изучения процессов распространения примесей. Проведенное исследование позволяет с помощью наблюдений из космоса определить параметры источников.
Оценки эффективных высот подъема дымовых шлейфов ТЭЦ-3, ТЭЦ-5 г. Новосибирска, вычисленных по формулам (1), (2) (табл.3) и по спутниковой информации вполне сопоставимы. Наиболее высокий уровень совпадения при АТ = 200°С. Оценка Нэф по формулам (1), (2) для ТЭЦ-3 составила 291 м, а для
ТЭЦ-5 - 505 м. Значения Нэф, полученные на основе спутниковых данных
(рис. 1), составили 344 м и 516 м соответственно. Следует отметить, что для конкретного момента времени неадекватное определение параметров, таких как скорость ветра, температуры выбросов, может привести к значительным погрешностям при вычислении Нэф . В то же время, оценки Нэф, полученные по
космоснимку, могут оказаться более объективными.
По спутниковой информации о траекториях шлейфов (рис. 1), данных с метеостанции, уравнений (3) рассчитаны компоненты скорости ветра в пограничном слое атмосферы г. Новосибирска (рис. 2). Полученные ветровые характеристики позволяют более точно описывать процессы распространения загрязнения территории города и его окрестностей от высотных источников.
Заключение
Предложен метод косвенного оценивания эффективной высоты подъема шлейфа выбросов от точечного источника, базирующийся на спутниковой информации. Наиболее эффективно применение этого метода при слабых ветрах, устойчивой или нейтральной стратификации пограничного атмосферы.
В таких метеорологических условиях реализуются благоприятные возможности для визуализации траекторий дымовых шлейфов в достаточном диапазоне расстояний от источника выброса. Наличие снегового покрова усиливает контрастность проекции шлейфа (тень) на земную поверхность. Проведена апробация предложенного похода для ряда крупных ТЭЦ г. Новосибирска и выполнено сравнение с результатами расчетов эффективной высоты подъема шлейфа по общепринятым методикам.
При отсутствии достоверной информации об источнике примеси имеется дополнительная возможность оценивать динамические и тепловые характеристики выбросов на основе спутниковых наблюдений. Подобная ситуация может возникнуть при аварийных выбросах. Полученные результаты указывают на необходимость комплексного использования различных видов информации
и расчетных методов для повышения контроля процессов загрязнения атмосферы от природных и антропогенных источников.
Работа выполнена при поддержке проекта РФФИ (грант № 17-47540342), госзадания 0315-2016-0004, проекта ПРАН№ 51 (0315-2018-0016).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Balter B. M., Faminskaya M. V. Irregularly emitting air pollution sources: acute health risk assessment using aermod and the Monte Carlo approach to emission rate // Air Quality, Atmosphere and Health. - 2017. Vol. 10 (4). - P. 401-409.
2. Балтер Б. М., Балтер Д. Б., Егоров В. В., Стальная М. В. Использование данных ИСЗ Landsat для определения концентрации загрязнителей в шлейфах от продувки газовых скважин на основании модели источника // Исследование Земли из космоса. - 2014. - № 2. - С. 55-66.
3. Lezhenin A. A., Raputa V. F., Yaroslavtseva T. V. Numerical Analysis of Atmospheric Circulation and Pollution Transfer in the Environs of Norilsk Industrial Region // Atmospheric and Oceanic Optics. - 2016. Vol. 29 (6). - P. 565-569.
4. Obolkin V. A., Potemkin V. L., Makukhin V. L., Chipanina Y. V., Marinayte I. I. Low-level atmospheric jets as main mechanism of long-range transport of power plant plumes in the Lake Baikal Region // Int. J. Environ. Studies. - 2014. Vol. 71 (3). P. 391-397.
5. Ding A. J., Huang X., Nie W., Sun J. N., Kerminen V.-M., Petäjä T., Su H., Cheng Y. F., Yang X.-Q., Wang M. H., Chi X. G., Wang J. P., Virkkula A., Guo W. D., Yuan J., Wang S. Y., Zhang R. J., Wu Y. F., Song Y., Zhu T., Zilitinkevich S., Kulmala M., Fu C. B., 2016: Enhanced haze pollution by black carbon in megacities in China // Geophysical Research Letters. - 2016. Vol.43 (6). - P. 2873-2879.
6. Кондратьев К. Я., Григорьев Ал. А., Покровский О. М., Шалина Е. В. Космическое дистанционное зондирование атмосферного аэрозоля. - Л.: Гидрометеоиздат, 1983. - 216 с.
7. Solomos S., Amiridis V., Zanis P., Gerasopoulos E., Sofiou F. I., Herekakis T., Brioude J., Stohl A, Kahn R. A., Kontoes C. Smoke dispersion modeling over complex terrain using high resolution meteorological data and satellite observations - The FireHub platform // Atmos. Environ. -2015. - Vol. 119. - P. 348-361.
8. Tohidi A., Kaye N. B. Highly buoyant bent-over plumes in a boundary layer // Atmos. Environ. - 2016. - Vol. 131. - P. 97-114.
9. Baklanov A., Smith Korsholm U., Nuterman R., Mahura A., Nielsen K. P., Sass B. H., Rasmussen A., Zakey A., Kaas E., Kurganskiy A., S0rensen B., and González-Aparicio, I.: Enviro-HIRLAM online integrated meteorology-chemistry modelling system: strategy, methodology, developments and applications (v7.2) // Geosci. Model Dev., - 2017. - Vol 10. - P. 2971-2999.
10. Siljamo P., Sofiev M., Filatova E., Grewling L., Jäger S., Khoreva E., Linkosalo T., Ortega Jimenez S., Ranta H., Rantio-Lehtimäki A., Svetlov A., Veriankaite L., Yakovleva E., Kukkonen J. A numerical model of birch pollen emission and dispersion in the atmosphere. Model evaluation and sensitivity analysis // Int. J. Biometeorol., - 2013. - Vol. 57. - P. 125-136.
11. Sokhi R., Baklanov A., Schluenzen H. (Eds.): Air Pollution and Meteorological Modelling for Atmospheric Research and Policy Applications, Anthem Press, 2017 - 260 p.
12. Gao Z., Bresson R., Qu Y., Milliez M., Munck C., Carissimo B. High resolution unsteady RANS simulation of wind, thermal effects and pollution dispersion for studying urban renewal scenarios in a neighborhood of Toulouse // Urban Climate. - 2018. - Vol. 23. - P. 114-130.
13. Li H., Cui G., Zhang Z. A New Scheme for the Simulation of Microscale Flow and Dispersion in Urban Areas by Coupling Large-Eddy Simulation with Mesoscale Models // Boundary-Layer Meteorol. - 2018. - Vol. 167. - P. 145-170.
14. Fallah-Shorshani M., Shekarrizfard M., Hatzopoulou M. Evaluation of regional and local atmospheric dispersion models for the analysis of traffic-related air pollution in urban areas // Atmos. Environ. - 2017. - Vol. 167. - P. 270-282.
15. Бызова Н. Л., Гаргер Е. К., Иванов В. Н. Экспериментальные исследования атмосферной диффузии и расчеты рассеяния примеси. - Л.: Гидрометеоиздат, 1991. - 278 с.
16. Islam M. A., Roy G. D. A mathematical model in locating an unknown emission source // Water, Air, and Solid Pollut. - 2002. - Vol. 136. - N 1/4. - P. 331-345.
17. Mikhailuta S. V., Lezhenin A. A., Pitt A., Taseiko O. V. Urban wind fields: Phenomena in transformation // Urban Climate. - 2017. - Vol. 19. - P. 122-140.
18. Baas P., van de Wiel B. J. H., van der Linden S. J. A., Bosveld F. C. From Near-Neutral to Strongly Stratified: Adequately Modelling the Clear-Sky Nocturnal Boundary Layer at Cabauw // Boundary-Layer Meteorol. - 2018. - Vol. 166. - P. 217-238.
19. Рапута В. Ф., Шлычков В. А., Леженин А. А., Романов А. Н., Ярославцева Т. В. Численный анализ данных аэрозольных выпадений примесей от высотного источника // Оптика атмосферы и океана. - 2014. - Т. 27. - № 8. - С. 713-718.
20. Зилитинкевич С. С. Атмосферная турбулентность и планетарные пограничные слои. - М.: ФИЗМАТЛИТ, 2013. - 252 с.
21. Методы расчетов рассеивания выбросов вредных (загрязняющих) веществ в атмосферном воздухе // Утверждены приказом Минприроды России от 06.06.2017 № 273.
22. Леженин А. А., Ярославцева Т. В., Рапута В. Ф. Использование спутниковой информации о траекториях дымовых факелов для расчета полей ветра // Интерэкспо ГЕ0-Сибирь-2016. XII Междунар. науч. конгр. : Междунар. науч. конф. «Дистанционные методы зондирования Земли и фотограмметрия, мониторинг окружающей среды, геоэкология» : сб. материалов в 2 т. (Новосибирск, 18-22 апреля 2016 г.). - Новосибирск : СГУГиТ, 2016. Т. 1. - C. 63-67.
REFERENCES
1. Balter, B. M., & Faminskaya, M. V. (2017). Irregularly emitting air pollution sources: acute health risk assessment using aermod and the Monte Carlo approach to emission rate. Air Quality, Atmosphere and Health, 10 (4), 401-409. doi: 10.1007/s11869-016-0428-x
2. Balter, B. M., Balter, D. B., Egorov, V. V., & Stalnaya, M. V. (2014). Fitting the emission model and meteorology to Landsat data on pollutant plumes from gas well pumping. Issledovanie Zemli iz kosmosa [Earth Observation and Remote Sensing], 2, 55-66 [in Russian].
3. Lezhenin, A. A., Raputa, V. F., & Yaroslavtseva, T. V. (2016). Numerical analysis of atmospheric circulation and pollution transfer in the environs of Norilsk industrial region. Atmospheric and Oceanic Optics, 29(6), 565-569.
4. Obolkin, V. A., Potemkin, V. L., Makukhin, V. L., Chipanina, Y. V., & Marinayte, I. I. (2014). Low-level atmospheric jets as main mechanism of long-range transport of power plant plumes in the Lake Baikal Region. Int. J. Environ. Studies, 71(3), 391-397. doi:10.1080/00207233.2014.918396.
5. Ding, A. J., Huang, X., Nie W., Sun, J. N., Kerminen, V.-M., Petaja, T., Su, H., Cheng, Y. F., Yang, X.-Q., Wang, M. H., Chi, X. G., Wang, J. P., Virkkula, A., Guo, W. D., Yuan, J., Wang, S. Y., Zhang, R. J., Wu, Y. F., Song, Y., Zhu, T., Zilitinkevich, S., Kulmala, M., Fu, C. B. (2016). Enhanced haze pollution by black carbon in megacities in China. Geophysical Research Letters, 43, 6, 2873-2879. doi: 10.1002/2016GL067745.
6. Kondrat'yev, K. Ya., Grigor'yev, Al. A., Pokrovskiy, O. M., Shalina, E. V. (1983). Kosmicheskoe distantsionnoe zondirovanie atmosfernogo aerozolya [Space remote sensing of an atmospheric aerosol]. Leningrad: Gidrometeoizdat [in Russian].
7. Solomos, S., Amiridis, V., Zanis, P., Gerasopoulos, E., Sofiou, F. I., Herekakis, T., Brioude, J., Stohl, A, Kahn, R. A., Kontoes, C. (2015). Smoke dispersion modeling over complex terrain using high resolution meteorological data and satellite observations - The FireHub platform. Atmos. Environ., 119, 348-361. https://doi.org/10.1016/j.atmosenv.2015.08.066
8. Tohidi, A., Kaye, N. B. (2016). Highly buoyant bent-over plumes in a boundary layer. Atmos. Environ., 131, 97-114. https://doi.org/10.1016/j.atmosenv.2016.01.046
9. Baklanov, A., Smith Korsholm, U., Nuterman, R., Mahura, A., Nielsen, K. P., Sass, B. H., Rasmussen, A., Zakey, A., Kaas, E., Kurganskiy, A., S0rensen, B., and González-Aparicio, I. (2017). Enviro-HIRLAM online integrated meteorology-chemistry modelling system: strategy, methodology, developments and applications (v7.2). Geosci. Model Dev., 10, 2971-2999. https://doi.org/10.5194/gmd-10-2971-2017
10. Siljamo, P., Sofiev, M., Filatova, E., Grewling, L., Jäger, S., Khoreva, E., Linkosalo, T., Ortega Jimenez, S., Ranta, H., Rantio-Lehtimäki, A., Svetlov, A., Veriankaite, L., Yakovleva, E., and Kukkonen, J. (2013). A numerical model of birch pollen emission and dispersion in the atmosphere. Model evaluation and sensitivity analysis. Int. J. Biometeorol., 57, 125-136. https://doi.org/10.1007/s00484-012-0539-5
11. Sokhi, R., Baklanov, A., & Schluenzen, H. (Eds.) (2017). Air Pollution and Meteorological Modelling for Atmospheric Research and Policy Applications. Anthem Press.
12. Gao, Z., Bresson, R., Qu, Y., Milliez, M., Munck, C., Carissimo, B. (2018). High resolution unsteady RANS simulation of wind, thermal effects and pollution dispersion for studying urban renewal scenarios in a neighborhood of Toulouse. Urban Climate, 23, 114-130. https://doi.org/10.1016/j.uclim.2016.11.002
13. Li, H., Cui, G., & Zhang, Z. (2018). A New Scheme for the Simulation of Microscale Flow and Dispersion in Urban Areas by Coupling Large-Eddy Simulation with Mesoscale Models. Boundary-Layer Meteorol.,167, 145-170. https://doi.org/10.1007/s10546-017-0323-5
14. Fallah-Shorshani, M., Shekarrizfard, M, Hatzopoulou, M. (2017). Evaluation of regional and local atmospheric dispersion models for the analysis of traffic-related air pollution in urban areas. Atmos. Environ., 167, 270-282. https://doi.org/10.1016Zj.atmosenv.2017.08.025
15. Byzova, N. L., Garger, E. K., & Ivanov, V. N. (1991). Eksperimental'nye issledovaniya atmosfernoy diffuzii i raschety rasseyaniya primesi [Experimental researches of atmospheric diffusion and calculations of dispersion of pollutions]. Leningrad: Gidrometeoizdat [in Russian].
16. Islam, M. A., Roy, G. D. (2002). A mathematical model in locating an unknown emission source. Water, Air, and Solid Pollut., 136(1/4), 331-345.
17. Mikhailuta, S. V., Lezhenin, A. A., Pitt, A., Taseiko, O. V. (2017). Urban wind fields: Phenomena in transformation. Urban Climate, 19, 122-140. http://dx.doi.org/10.1016/j.uclim.2016.12.005
18. Baas, P., van de Wiel, B. J. H., van der Linden, S. J. A., Bosveld, F. C. (2018). From Near-Neutral to Strongly Stratified: Adequately Modelling the Clear-Sky Nocturnal Boundary Layer at Cabauw. Boundary-Layer Meteorol., 166, 217-238. https://doi.org/10.1007/s10546-017-0304-8
19. Raputa, V. F., Shlychkov, V. A., Lezhenin, A. A., Romanov, A. N., Yaroslavtseva, T. V. (2014). Numerical analysis of aerosol substance fallout from a high-altitude source. Optika Atmosfery i Okeana. [Atmospheric and Oceanic Optics], 27(8), 713-718 [in Russian].
20. Zilitinkevich, S. S. (2013). Atmosfernaya turbulentnost' i planetarnye pogranichnye sloi [Atmospheric turbulence and planetary boundary layers]. Moscow: FIZMATLIT [in Russian].
21. Metody raschetov rasseivaniya vybrosov vrednykh (zagryaznyayushchikh) veshchestv v atmo-sfernom vozdukhe (2017) [Methods of calculations of dispersion of the harmful (polluting) substances in atmospheric air]. Utverzhdeny prikazom Minprirody Rossii ot 06.06.2017 № 273 [in Russian].
22. Lezhenin, A. A., Raputa, V. F. & Yaroslavtseva, T. V. (2016). Use of satellite information on trajectories of the plumes of smoke for the calculation of wind fields. In Sbornik materialov Interekspo Geo-Sibir' - 2016: Mezhdunarodnoy nauchnoy konferentsii: T. 1. Distantsionnye metody zondirovaniya Zemli i fotogrammetriya, monitoring okruzhayushchey sredy, geoekologiya - 2016 [Proceedings of Interexpo GEO-Siberia-2016: International Scientific Conference: Vol. 1. Remote sensing methods of the Earth and photogrammetry, environmental monitoring, geoecology-2016] (pp. 63-67). Novosibirsk: SSUGT [in Russian].
© B. O. Panyma, A. A. Hewenun, T. B. flpocme^ea, 2018
