Бесплатный автореферат и диссертация по наукам о земле на тему

Метод расчета и климатологического анализа характеристик атмосферы, определяющих перенос и рассеяние примеси

ВАК РФ 25.00.30, Метеорология, климатология, агрометеорология

Автореферат диссертации по теме "Метод расчета и климатологического анализа характеристик атмосферы, определяющих перенос и рассеяние примеси"

Министерство образования и науки РФ Федеральное агентство по образованию

ГОУ ВПО РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ (РГГМУ)

На правах рукописи ' - УДК 551.511.072

Мханна Ааед Исмаил Назир

МЕТОД РАСЧЕТА И КЛИМАТОЛОГИЧЕСКОГО АНАЛИЗА ХАРАКТЕРИСТИК АТМОСФЕРЫ, ОПРЕДЕЛЯЮЩИХ ПЕРЕНОС И РАССЕЯНИЕ ПРИМЕСИ

Специальность 25.00.30 - метеорология, климатология, агрометеорология

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

,, уд Т5®

Санкт-Петербург 2009

003471411

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

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

Официальные оппоненты:

Профессор, доктор физико-математических наук, А.С.Гавршое

Доктор физико-математических наук, профессор А.Д. Егоров

Кандидат физико-математических наук, доцент, А.Г.Попов

Ведущая организация:

Санкт-Петербургский научно-исследовательский центр экологической безопасности РАН (НИЦЭБ РАН)

Защита диссертации состоится "18" июня 2009 г. в ¿5 часов 30 минут на заседании диссертационного совета Д212.197.01 в Российском государственном гидрометеорологическом университете по адресу:

195196, г. Санкт-Петербург, Малоохтинский пр., 98, тел. (812) 444-41-63.

С диссертацией можно ознакомиться в библиотеке Российского государственного гидрометеорологического университета по адресу: 195196, г. Санкт-Петербург, Малоохтинский пр., 98.

Автореферат разослан "18" мая 2009 г.

Ученый секретарь диссертационного совета,

доктор физ.-мат наук, профессор, А.Д. Кузнецов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы. В развивающихся странах за последние десятилетия достигнут существенный прогресс в развитии промышленности и энергетики. Это ставит их перед теми же проблемами в области охраны окружающей среды, что и развитые в промышленном отношении страны Европы и Северной Америки.

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

Целы« работы является: разработка модели и проведение комплекса расчетов основных параметров, определяющих рассеивающие свойства атмосферы применительно к такому интенсивно развивающемуся региону, как Восточное Средиземноморье, их статистическое обобщение и построение, в итоге, набора климатических карт этих параметров.

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

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

• адаптировать разработанную ранее численную модель атмосферного пограничного слоя (АПС) к расчету характеристик рассеивания и осуществить ее верификацию;

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

• провести комплекс расчетов и построить климатические карты параметров, определяющих рассеивающие способности атмосферы.

Методы исследования. Основные методы исследования - компьютерное моделирование с использованием численной модели АПС, а также статистический и пространственный анализ результатов.

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

Основные положения и результаты, выносимые на защиту:

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

2. Результаты статистического анализа временной изменчивости характеристик рассеивающей способности атмосферы на различных территориях.

3. Климатические карты характеристик рассеивающей способности атмосферы для Восточного Средиземноморья.

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

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

Апробация диссертационной работы

Основные результаты исследований, изложенные в диссертации, докладывались и обсуждались:

• на научных семинарах кафедры метеорологи, климатологии и охраны атмосферы Российского государственного гидрометеорологического университета (2003-2008 гг.).

• на заседании Итоговых сессий Ученого совета РГГМУ (2005,2006 гг.);

• на заседаниях Международной научной конференции "Экологические и гидрометеорологические проблемы больших городов и промышленных зон", СПб, РГГМУ (2005г.);

• на заседании международной конференции "Изменение климата и окружающая среда", СПб, РГГМУ, 06-09.12.2005 г.;

Публикации: Основные результаты диссертации опубликованы в 7 печатных работах.

Структура и объем работы:

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

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе рассмотрены существующие проблемы и применяемые в различных странах методы оценки качества атмосферы. Особое внимание уделено природоохранной практике стран Восточного Средиземноморья (Турция, Сирия, Ливан, Египет, Израиль, Палестина, Иордания). Сделан краткий анализ имеющихся в литературе сведений о суммарных выбросах загрязняющих веществ в атмосферу в этих странах и существующих тенденциях такого рода за последние годы.

Загрязнение атмосферы в городах Восточного Средиземноморья, как и южной Европы, за последние десятилетия значительно возросло. По свидетельству Всемирной Организации Здравоохранения (ВОЗ), к числу основных загрязняющих примесей здесь следует отнести окись углерода (СО), оксиды азота (МОх) и оксиды серы (50х). Установлено, что наибольший вклад в загрязнение атмосферного воздуха на территории упомянутых стран вносит теплоэнергетика, автотранспорт и химическая промышленность.

Так, например, в Турции, за период 1980-2005 гг. суммарный годовой выброс оксидов азота вырос почти в 3 раза, а оксидов серы - примерно в 1.8 раза. Интенсивный рост выбросов в атмосферу отмечен в Ливане: лишь за 3 года с 1997 по 2000гг количество выбросов ЫОх увеличилось почти в 2.5 раза, а БОх - почти в 3 раза. За последние годы за счет интенсивного расширения автопарка столицы Египта г. Каира, концентрация КОх выросла в несколько раз, что совместно с эффектом непрерывного поступления из прилегающей пустыни взвешенных аэрозолей привело к образованию почти непрерывно висящего над городом облака загрязнения.

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

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

За последние десятилетия неоднократно осуществлялись попытки установления некоторых интегральных показателей, характеризующих климатические условия загрязнения атмосферы. Так, Э.Ю. Безуглая (1980) предложила методологию использования такого показателя, как потенциал загрязнения атмосферы (ПЗА). Кроме того, для прогноза неблагоприятных метеорологических условий (НМУ) были сформулированы так называемые «индексы НМУ» для шести групп источников (высокие горячие, высокие холодные, средние горячие, средние холодные, низкие и

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

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

- значения скорости ветра в приземном слое для использования при расчете индекса НМУ;

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

- вертикальные градиенты скорости ветра и температуры для расчета категорий устойчивости Пэскуилла-Гиффорда.

Известно, что данные вертикального зондирования атмосферы оказываются недостаточно точными для нижнего слоя тропосферы, поскольку радиозонд не дает возможности корректно оценить вертикальный профиль температуры в нижнем 100200 м слое и высоту слоя перемешивания, да и производятся они на очень небольшом числе станций зондирования. Другие методы (содарное и лидарное зондирование, высотные градиентные мачты) предоставляют вообще лишь эпизодический материал и непригодны для получения надежной статистики и последующего пространственного анализа.

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

- организация специальных масштабных экспериментальных исследований для каждой территории;

- путем расчета с использованием математических моделей атмосферы на основе стандартной гидрометеорологической информации.

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

Во второй главе излагается метод расчета характеристик рассеивающей способности атмосферы на основе численной модели атмосферного пограничного слоя (АПС) и некоторые его уточнения, реализованные в настоящей работе.

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

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

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

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

В разделе 2.2 данной главы сформулированы исходные уравнения, а в пункте 2.3 - описана модель суточных колебаний. Для АПС над достаточно однородной подстилающей поверхностью систему уравнений динамики для отклонений 4 = можно записать:

дид„ди ,,

ся СИ, ох,

-к0^--2 0)ги (2)

ди _ 3 Эу 3/ Эх, и Эдг,

Здесь ось дг, декартовой системы координат направлена на восток, ось х2 - на север, ось х, - вертикально вверх, и,,С/, - горизонтальные компоненты скорости ветра, а ис, Ус — соответствующие компоненты геострофического ветра, со. = й«т <р {со - угловая скорость вращения Земли, <р - широта).

Особенностью используемого подхода является то, что уравнение для потенциальной температуры 0 записывается здесь следующим образом:

3* Э эв (3)

Эг аху дху

где функция ^в(х1,х2,х3,0 - совокупность всех остальных членов уравнения притока тепла, не учитываемых в явной форме в (3). В данной интерпретации это уравнение не представляют собой ничего нового по сравнению с обычным до тех пор, пока не будет указан способ определения этой функции.

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

специальные исследования (A.C. Гаврилов, Д.Л. Лайхтман, 1973 и др.), ее вертикальный профиль имеет целый ряд характерных особенностей, причем максимальные по модулю значения достигаются в непосредственной близости к подстилающей поверхности (за счет интенсивного длинноволнового излучения вследствие значительных локальных градиентов температуры). Если выбрать нижнюю границу области расчета совпадающей с высотой измерения температуры на уровне ZK-2 м (высота метсобудки), то влияние радиационных притоков тепла на формирование вертикальной структуры поля температуры уже не будет столь существенным.

В отличие от проведенных ранее исследований (И.Н. Липовицкая, 2008г.), где для расчета Ч^(г) в уравнении (3) привлекались лишь данные стандартных гидрометеорологических измерений на метеостанциях, в нашей работе предпринята попытка использования для этих целей трехмерных синоптических полей температуры (данные анализа температуры на уровнях 1000, 925 и 850 гПа), синхронизированные по времени с данными метеостанций.

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

Qm(t) = B™+B™t + B2lm)t\ (4)

где т определяет здесь номер поверхности (т=2 для 925 гПа, т=3 - для 850 гГТа ). На нижнем уровне ст=1 используются данные обработки аналогичным способом результатов измерений на метеостанциях (4).

Определение коэффициентов В■ осуществляется здесь методом наименьших квадратов из условия минимальной среднеквадратической ошибки представления:

j>je(tK)- в0- BttK - B2tl f =тт(В,). (5)

/Ы

Дифференцируя это выражение последовательно по B(r ßp В2 и приравнивая результаты нулю, получаем в итоге систему трех линейных уравнений для вычисления необходимых коэффициентов.

Вычисление искомых значений функции Ve(xl,x2,x3,t) осуществляется после линейной интерполяции значений (4) в расчетную сетку модели АПС с координатами сеточных узлов:

Х(=Х0 + Д,(|-1); Yj=Y0+byU-iy, Z, = Z, ехр[Д,(* -1)], (б)

где для горизонтальных координат Х,У приняты постоянные шаги сетки Д х, Д „ а для вертикальной координаты 2 - так называемая логарифмическая сетка, в которой

Гг„/' 'г,

постоянным является шаг логарифма высоты Ah = —— 1 rt

причем Z|,ZH -

высоты, соответственно нижнего и верхнего уровней расчетной области, а число расчетных уровней по вертикали. Величины Хо,1о здесь - координаты нижнего левого угла области расчета.

При этом расчет функции 4*,,О,,г) для каждого сеточного узла (ij,k) здесь может быть произведен следующим образом:

*?<0 = (f) ■ (7)

Над поверхностью суши высота нижнего расчетного уровня принимается равной Z\- 2 м и совпадает с высотой измерений в стандартной метеорологической будке. Это дает возможность задавать нижние граничные условия для функции 6(t).

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

~ = Mw(ei-Tw), (8)

дхъ

где Mw = у^ ! ^ ^, а для расчета величины Zw как функции скорости ветра использовалась известная формула Чарнока.

Схема замыкания уравнений (1)~(3) опиралась на подход, основанный на дополнительном привлечении уравнений для вторых и, частично, третьих моментов одноточечных моментов турбулентных пульсаций, хорошо апробированный и верифицированный применительно к расчету структуры АПС (A.C. Гаврилов, 1986, 1992).

Систему уравнений (1) - (3) следует интегрировать по времени, начиная с некоторого момента t0, с заданием фиктивных начальных условий, поскольку фактические сведения о реальной вертикальной структуре АПС отсутствуют.

Для компонент скорости ветра в этом качестве привлекались функции uST,vST{xy), удовлетворяющие стационарному решению уравнений для горизонтальных компонент скорости, а для температуры - линейный профиль потенциальной температуры, получаемый на основе расчета ее среднего градиента ув по данным реанализа на поверхностях 925, 850 и 700 гПа: 0(х,,х2,х,,го) = в(/0) + /ff(х,,х2)х3.

В реальных условиях не во всех случаях суточный ход температуры носит «правильный» характер это происходит за счет, прежде всего, резкой смены синоптических условий при прохождении атмосферных фронтов и влияния облачности, которая может существенно искажать суточный ход инсоляции. Понятно, что в этих условиях аппроксимация (4) оказывается неприемлемой и вместо нее следует ограничиться получением путем той или иной формой интерполяции зависимости температуры от времени между сроками наблюдений. Отметим, что подобная ситуация имеет место не более чем в 10 % случаев.

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

~ = Ми(и + ис),

Ъх-.

(7)

где V,,, =д1(и+исУ~ +{у + Ус)2 - модуль скорости ветра, Ми =

ге1п{2Б/гау а 2<

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

На верхней границе расчетной области Х„, заведомо превосходящей высоту распространения суточных колебаний (около 2 км), все отклонения температуры и компонент скорости ветра от своих фоновых значений.

Для верификации предложенного метода применительно к расчету характеристик термической устойчивости АПС привлекались уникальные данные метеорологических наблюдений, осуществляемых в настоящее время на Останкинской телевизионной башне (ОТБ) в г. Москва (температура, влажность, скорость и направление ветра до высоты 500м).

На основе измерений на двух нижних уровнях (82 и 2 м) за период летних месяцев 2006г. рассчитывалось так называемое интегральное число Ричардсона

падения, Т - температура в шкале Кельвина), АИ = 80м - толщина слоя, а Ав,АУ -разности значений потенциальной температуры и скорости в этом слое, соответственно. Результаты сопоставления измеренных и расчетных величин числа Ш представлены на рисунке 1.

г, где р = ■ - параметр плавучести (g -ускорение свободного

в

-4

-2

0

по измерениям

2

4

6

Рисунок 1 - Корреляционный график рассчитанных и измеренных на ОТБ интегральных чисел Ричардсона в слое 2 - 82 м за летний период 2006 г.

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

На рисунке 2 представлена рассчитанная за летний период 2006 г диаграмма соответствия категорий устойчивости Пэскуилла-Гиффорда, рассчитанных с использованием предложенного нами метода и непосредственно из данных градиентных наблюдений скорости ветра и температуры в слое 2 - 82 м на ОТБ. Из представленных результатов следует, что вероятность согласования расчета и измерений с точностью до категории оказывается около 80%, а если включать в число успешных расчетов и по одной соседней категории, то эта цифра достигнет 90 %. Этот результат можно признать вполне удовлетворительным с учетом приближенного характера собственно понятия «категории устойчивости».

Рисунок 2 - Диаграмма соответствия категорий устойчивости Пэскуилла-Гиффорда, рассчитанных с использованием модели АПС; 1 - вероятность совпадения с точностью до категории, 2-е точностью до соседней категории.

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

В третьей главе представлен климатологический анализ результатов расчетов применительно ко всем странам Восточного Средиземноморья. Исходная информация была получена из общедоступного архива метеоданных для метеостанций, имеющих пятизначные номера Всемирной Метеорологической Организации (ВМО). Для получения необходимой расчетов функции (7) привлекались накопленные за период 2000 - 2006 гг. архивы данных реанализа для Северного полушария в

формате GRIB с сервера национальной Службы Погоды США (National Weather Service, http://nomad3.ncep.noaa.gov/cgi-bin/ftp2u_6p_r2.sh).

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

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

В качестве примера на рисунке 3 представлена общая повторяемость высоты приземной инверсии в целом для всего региона, а на рисунке 4 - ее суточная изменчивость.

Рисунок 3 - Повторяемость высоты верхней границы инверсии зимой в целом для всего региона Восточного Средиземноморья.

Час

—•—Январь "—Июнь —*—Среднегодовая

Рисунок 4 - Средняя высота инверсии в течение суток по региону Восточного Средиземноморья

Как можно видеть из приведенных примеров расчетов, кривая среднегодового хода средней высоты инверсии достигает максимума в 08:00, и составляет 55 м. Рост высоты инверсии с 03:00 до 09:00 (с 25 м до 55 м) можно объяснять тем, что в утренние часы после восхода солнца, происходит интенсивный нагрев поверхности, приводящий к развитию конвекции и быстрому разрушению инверсионного слоя. Затем наблюдается некоторое снижение высоты инверсионного слоя до 20 метров.

Еще один пример результатов расчетов высоты приземной инверсии приведен на рисунке 5, из которого можно выявить различия в высотах средних ночных при-, земных инверсий в северной области рассматриваемого региона (пункт Анталья), в центральной (пункт Хайфа) и южной его областях (пункт Порт-Саид).

М

80-100 Ш—

70-80

60-70 Iтттщш щ и I .,....

50-60 "|ЩЩЩШШЩшя

40-50 тттшыщшт*

30-40 щжр^

2о-зо штшшштттштш

]0_20 ЩЩЩЩЩШЖ''»»'»

о-ю И™

0 10 20 30

Повторяемость, %

Рисунок 5- Повторяемость высот приземных инверсий ночью в среднем за год по городам Хайфа, Порт-Саид и Анталья.

Как следует из представленных результатов, наибольшую повторяемость на севере и юге региона имеют высоты 20 - 30 м, а в центральной части - около 60 -70 м. Связано это, вероятно, с проявлением в синоптических процессах данного региона характерного для умеренной зоны Северного полушария западно-восточного переноса, в процессе которого на территорию таких стран, как Израиль и Сирия происходит заток в зимний период теплого морского воздуха со Средиземного моря. Этот эффект, вследствие конфигурации береговой черты, не является характерным для территорий Турции (п. Анталья) и Египта (п. Порт-Саид).

ЕЭ Анталья В Порт-Саид □ Хайфа

В качестве еще одного примера исследования временной изменчивости, приведем здесь рисунок 6, на котором представлен суточный ход повторяемостей различных категорий устойчивости атмосферы Пэскуилла-Гиффорда в целом для всего региона Восточного Средиземноморья за год.

1 3 5 7 9 11 13 15 17 19 21 23

Час

Рисунок 6 - Пример расчета повторяемостей различных категорий устойчивости атмосферы в различные часы суток в среднем за год.

Категории «А» и «В» в этой классификации соответствуют свободной конвекции, «С» - вынужденной, «О» - нейтральной стратификации, «Е», «Б» и «в» соответственно, слабой, умеренной и сильной устойчивости атмосферы. Как можно видеть из рисунка, в среднем, примерно в 60% случаев стратификация атмосферы близка к нейтральной, в дневные сроки наблюдений доминируют первые три категории устойчивости (А,В,С), а в ночные - последние (Е, Е и в).

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

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

привлекались также и данные о температуре подстилающей поверхности, задаваемой в градусной сетке 1° х 1°, также имеющиеся в нашем распоряжении.

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

На начальном этапе производится совмещение всех исходных данных в точках проведения расчетов. Для поверхности суши в качестве таких точек выбрано местоположение метеостанций, куда из градусной сетки реанализа (2,5° х 2,5°) переносились синхронизированные по времени сведения о значениях температуры на стандартных изобарических поверхностях.

Далее, для выбранных точек расчета с использованием описанной выше методологии восстанавливались фоновые характеристики атмосферы и для каждых суток за изучаемый период путем интегрирования уравнений модели на каждый час восстанавливалась вертикальная структура АПС, и рассчитывались искомые характеристики наблюдаемости.

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

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

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

В итоге сопряжения всех этих источников информации восстанавливалась единая для всей территории карта тех или иных статистических характеристик, которые подготавливались непосредственно в обменных форматах (*.5НР) ГИС Агсу1еш 3.2, графические изображения которых и анализируются в диссертации.

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

В Заключении сформулированы основные результаты и выводы:

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

2. Результаты верификации численной модели путем сопоставления с данными измерений на Останкинской телебашне позволили оценить

точность расчета характеристик вертикальной структуры поля температуры и скорости ветра, а также характеристик устойчивости атмосферы. Показано, что коэффициент корреляции между расчетными и измеренными значениями интегрального числа Ричардсона составляет величину 0.9, а вероятность точной оценки категории устойчивости -80 %.

Рисунок 7 - Расчетная повторяемость образования в течение года приподнятых инверсий.

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

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

Основные публикации по теме диссертации.

Статьи в изданиях из списка ВАК

1. Прогноз и климатологический анализ характеристик атмосферы, определяющих рассеяние антропогенных загрязнений // Естественные и технические науки. М., 2008. № 6 с. 221-225 (соавторы Гаврилов А.С, Алталули Р.)

Статьи в научно-техиических сборниках

2. Моделирование атмосферного пограничного слоя применительно к проблемам климатологии. // Математическое моделирование, численные методы и комплексы программ: Межвуз.темат.сб.тр., вып.12. / СПбГАСУ. - СПб.,2006, - с. 23-31 (соавторы Липовицкая И.Н., Василенко C.B., Гаврилов А.С).

3. Метод климатологического анализа вертикальной структуры атмосферного пограничного слоя с использованием численной модели. // Ученые записки № 2. Научно-теоретический журнал. - СПб.: изд. РГГМУ, 2006.- с.74-87 (соавторы Липовицкая И.Н., Василенко C.B., Гаврилов А.С).

4. Некоторые результаты климатический восточного средиземного моря. // Математическое моделирование, численные методы и комплексы программ: Межвуз.темат.сб.тр., вып.12. / СПбГАСУ. - СПб.,2006, с. 31- 35.

5. Численное моделирование атмосферного пограничного слоя применительно к прогнозу неблагоприятных метеорологических условий. // Математическое моделирование, численные методы и комплексы программ: Межвуз.темат.сб.тр., вып. 10. / СПбГАСУ. - СПб.,2004. с. 99-108 (соавторы Алталули Р., Гаврилов А.С).

Публикации в материалах конференций

6. Метод расчета карт потенциала загрязнения атмосферы города Москвы и его окрестностей // Экологические и гидрометеорологические проблемы больших городов и промышленных зон. Международная научная конференция 25-27 мая 2005. Материалы конференции / СПб.: изд. РГГМУ, 2005. с. 13-14. (соавторы Гаврилов A.C., Баранова М.Е.).

7. Метод расчета метеорологического потенциала загрязнения атмосферы для крупного города. Международной конференции "Изменение климата и окружающая среда", СПб, РГГМУ, 06-09.12.2005 г.; с. 107-108. (соавторы Баранова М.Е, Алталули Р.)

Отпечатано с готового оригинал-макета в ЦНИТ «АСТЕРИОН» Заказ № 142. Подписано в печать 14.05.2009 г. Бумага офсетная. Формат 60х84!/16. Объем 1,25 п. л. Тираж 100 экз. Санкт-Петербург, 191015, а/я 83, тел. /факс (812) 275-73-00, 970-35-70 asterion@asterion.ru

Содержание диссертации, кандидата физико-математических наук, Мханна Ааед Исмаил Назир

СОКРАЩЕНИЯ.

ВВЕДЕНИЕ.

1 ОСНОВНЫЕ ПРОБЛЕМЫ И МЕТОДЫ ОЦЕНКИ КАЧЕСТВА АТМОСФЕРЫ.!.!.

1.1 Особенности загрязнения атмосферы Восточного Средиземноморья.

1.2 Роль метеорологических факторов в переносе и рассеивании примесей в атмосфере.

1.3 Основные закономерности рассеивания примесей в атмосфере.

1.4 Категории устойчивости атмосферы Пэскуилла-Гиффорда.

1.5 Общие понятия о неблагоприятных метеорологических условиях.

1.6 Потенциал загрязнения атмосферы.

2 МЕТОД РАСЧЕТА ВЕРТИКАЛЬНОЙ СТРУКТУРЫ АТМОСФЕРНОГО ПОГРАНИЧНОГО СЛОЯ ПРИМЕНИТЕЛЬНО К РАСЧЕТУ ХАРАКТЕРИСТИК РАССЕИВАЮЩИХ СПОСОБНОСТЕЙ АТМОСФЕРЫ.

2.1 Постановка вопроса.

2.2 Исходные уравнения.

2.3 Модель суточных колебаний.

2.4 Апробация метода для различных регионов.

2.5 Метод расчета категорий устойчивости атмосферы по данным стандартных гидрометеорологических измерений и его верификация.

2.6 Обобщение метода на случай учета синоптического фона и температуры морской поверхности.

2.7 Экспериментальная проверка метода.

3 ИССЛЕДОВАНИЕ ПРОСТРАНСТВЕННО-ВРЕМЕННОЙ ИЗМЕНЧИВОСТИ РАССЕИВАЮЩИХ СПОСОБНОСТЕЙ АТМОСФЕРЫ ПРИМЕНИТЕЛЬНО К ВОСТОЧНОМУ СРЕДИЗЕМНОМОРЬЮ.

3.1 Постановка проблемы.

3.2 Анализ исходных данных.

3.3 Данные реанализа атмосферных процессов.

3.4 Программное обеспечение.

3.5 Общий статистический анализ результатов расчета климатических характеристик.

3.5.1 Характеристики инверсий температуры воздуха.

3.5.2 Категории устойчивости.

3.6 Пространственный анализ.

Введение Диссертация по наукам о земле, на тему "Метод расчета и климатологического анализа характеристик атмосферы, определяющих перенос и рассеяние примеси"

Актуальность работы. В развивающихся странах за последние десятилетия достигнут существенный прогресс в развитии промышленности и энергетики. Это ставит их перед теми же проблемами в области охраны окружающей среды, что и развитые в промышленном отношении страны Европы и Северной Америки.

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

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

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

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

• адаптировать разработанную ранее численную модель атмосферного пограничного слоя (АПС) к расчету характеристик рассеивания и осуществить ее верификацию;

• провести анализ доступных информационных ресурсов и сформировать архивы исходных данных для расчета в достаточном для статистического ^ и пространственного) анализа количестве;

• провести комплекс расчетов и построить климатические карты параметров, определяющих рассеивающие способности" атмосферы.

Методы исследования

Основные методы исследования — компьютерное моделирование с использованием численной модели АПС, а также статистический и пространственный анализ результатов.

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

Основные положения и результаты, выносимые на защиту:

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

2. Результаты статистического анализа временной изменчивости характеристик рассеивающей способности атмосферы на различных территориях.

3. Климатические карты характеристик рассеивающей способности атмосферы для Восточного Средиземноморья.

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

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

Основные результаты исследования, изложенные в диссертации, докладывались и обсуждались:

• на научных семинарах кафедры метеорологи, климатологии и охраны атмосферы Российского государственного гидрометеорологического университета (2003-2008 гг.);

• итоговых сессиях Ученого совета РГГМУ (2005,2006 гг.);

• на заседаниях Международной научной конференции "Экологические и гидрометеорологические проблемы больших городов и промышленных зон", СПб, РГГМУ (2005г.);

• на заседании международной конференции "Изменение климата и окружающая среда", СПб, РГГМУ, 06-09.12.2005 г.

Заключение Диссертация по теме "Метеорология, климатология, агрометеорология", Мханна Ааед Исмаил Назир

В заключение сформулируем основные выводы из проделанной работы.

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

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

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

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

Библиография Диссертация по наукам о земле, кандидата физико-математических наук, Мханна Ааед Исмаил Назир, Санкт-Петербург

1. Арабская Республика Египет. Справочник Текст. / Отв. ред. A.M. Васильев-М.: Наука, 1990.-355 с.

2. Основные климатические данные по странам ближнего- и- среднего: востока и' южной Азии? (1966 — 1975) Текст.: в 2 ч. 411 / Под; ред. Ю.С.Чернова;и^АЖ Староверовой.-М.:Гидрометеоиздат, 1981-133 с.

3. Детри, Ж.П1Атмосфера'Должна^быть,чистошТекст.;/Ж.Ш Детри-М.: Прогресс, 1973 .-379 с.

4. Сидоренко, Г.И. Гигиенические критерии максимально допустимой нагрузки Текст. / Г.И. Сидоренко, М.А. Пинигин // Всесторонний анализ окружающей природной среды: сб. работ.-Л.: Гидрометеоиздат,1975.-С. 119 — 128. '

5. Израэль, Ю:А. Допустимая антропогенная, нагрузка на окружающую природную среду Текст. / Ю.А. Израэль // Всесторонний анализ окружающейприродной среды: труды IIСоветско-американского симпозиума, Гонолулу,

6. Гавайи, 20 26 октября 1975 г.-Л.: Гидрометеоиздат, 1976.-С.12 - 19.

7. Лаклэр; Е. Экономические и социальные аспекты загрязнения атмосферы^ Текст. / Е. Лаклэр // Загрязнение атмосферного воздуха.-Женева: ВОЗ, 1963.

8. Air quality criteria for sulphur oxides Text.-Washington, D.C.:U.S. Department of Health, Education and Welfare, Jan. 1969.-P.178.

9. Air quality guidelines for particulate matter, ozone, nitrogen dioxide and sulfur dioxide: Global update 2005 Text.-Geneva: WHO, 2005.

10. Environmental; statistics in the Mediterranean countries Compendium 2005 Text.;—Istanbul: Turkish Ministry of Nature Protection, 2005.

11. Transportation emissions in Lebanon Text.: Extent and mitigation 2005—

12. Beirut: Lebanon Ministry of Nature Protection, 2005.

13. Качество воздуха // Офиц. сайт Палестинской' национальной автономии Электронный ресурс.; / Палестинский: стат. центр.-Рамалла-Режим доступа: www.pcbs.gov.ps.

14. Environmental: balance // World lunks Arab Reagion Electronic resource. / Ministry of education Sirya.-Damascus.-Point of access: www.pnic.gov.ps.

15. Статистические данные о выбросах загрязняющих веществ в атмосферу автомобилями // Офиц. сайт Палестинского ин-та стандарт. Электронный-ресурс. / PSI -Режим-доступа: www.pshgov.ps.

16. Monthly report от. .the.-, air/ quality in Egypt // The Environmental Information and Monitoring Program Electronic resource.?; 7 EIMP- Point of access: http://www.eeaa.gov.eg/eimp/index:html;

17. Эл Абси, Самир. Исследование влияния метеорологических условий на формирование режима загрязнения* большого города;: и его окрестностей Текст.: Магистерская диссертация / Самир Эл Абси.-СПб.: РП'МИ, 1995;102 с. ;.;■" ■.'. . . .

18. Schrenk, Н.Н. Air pollution in Donora. Epidemiology of unusual; smog episode of October 1948 Text. / H;H. Schrenk- USA: By Health Security. Agency, 1949.

19. Берлянд, М.Е. О расчете загрязнения атмосферы выбросами из дымовых труб электростанций Текст. / М.Е. Берлянд- Е.ЖЕенихович, Р.И Оникул // Труды ГГО -1964- Вып. 158.-С.З 21.

20. Оникул, Р.И:, Методика расчета загрязнения атмосферы выбросами промышленных предприятий и тепловых электростанций Текст. / Р.И: Оникул // Метеорологические аспекты загрязнения: атмосферы: сб. науч. работ.-Л.: Гидрометиздат, 1971.-С.70 81.

21. Берлянд, М.Е. К обобщению теории рассеивания промышленных выбросов в атмосферу Текст. / М.Е. Берлянд, Р.И. Оникул // Труды ГГО.— 1971.-Вып. 254.-С.З -38.

22. Holzworth, G.C. Estimates of mean maximum, mixing depths in1 the contiguous United-States Text. / G.C. Holzworth // Mon. Wea. Rev.-1964-Vol.92, № 5.-P.17 19.

23. Fensterstoork, J. Thanksgiving , 1966. Air Pollution episode in the Eastern United States Text. / J: Fensterstoork-U.S. Department of Health, Education and Welfare;.Consumer Protection and Environmental Health Service, 1968.

24. Wolzworth, G.C. Glimatological aspects of the composition and pollution of the atmosphere Text. / G.C. Wolzworth // WMO Techn. Note-1974-№ 139.-P.43.

25. Берлянд, М.Е. Численное исследование атмосферной диффузшг при» нормальных и аномальных условиях стратификации Текст. / М.Е. Берлянд. Е.Л. Гениховвич. В.П. Ложкина, Р.И. Оникул // Тр. ГГО-1964.-Вып. 158-С.22-32.

26. Безуглая, Э.Ю. К определению потенциала загрязнения воздуха Текст. / Э.Ю. Безуглая // Тр. ГГО.-1968.-Вып. 234.-С.69 79.

27. Шевчук, И.А. Повторяемость метеорологических условий, способствующих увеличению загрязнения приземного слоя атмосферы в Новосибирске Текст. / И.А. Шевчук, Л.И. Введенская, Т.Г. Володкевич // Тр. Новосиб. Регион. ГМЦ.-1969.-Вып. 2.

28. Bryson, A. Air pollution Text. / A. Bryson, Е. Kutzbach // Assoc. Am. Geograph. Resource Pap.-1968.-№ 2.-P.62.

29. Holzworth, G.C. A study of air potential for the Western United States Text. / G.C. Holzworth // J. Appl. Met.-1962.-Vol. 1, № 3.

30. Сонькин, Л.Р. Синоптические условия формирования периодов высокого загрязнения воздуха в различных районах СССР Текст. / Л.Р. Сонькин //ТрудыГТО.-1979-Вып. 436.-С.49

31. Берлянд, М.Е. Современные проблемы атмосферной диффузии и загрязнения атмосферы Текст. / М.Е. Берлянд-Л.: Гидрометеоиздат,1975-448 с.

32. Царев, A.M. К вопросу о загрязнении воздуха при туманах Текст. / A.M. Царев // Тр. ГТО.-19721-Вып. 352.- С. 113 118:

33. Hanna, S.R. A simple dispersion model for the analysis of chemically reactive pollutants Text. / S.R. Hanna // Atmos. Environ.-1973.-Vor. 7, № 8 — P.803 817.

34. Hoagen-Smit, A.J. Ozone formation in photochemical oxidation of organic substances Text. / A.J. Hoagen-Smit, M.M. Fox // Ind. Eng. Chem-1956-Vol. 4.

35. Дмитриев, M.T. Механизм фотохимического загрязнения атмосферы городов Текст. / М.Т. Дмитриев, Н.А. Китросский, В. А. Попов // Метеорологические аспекты загрязнения атмосферы: сб. науч.работ—Л.: Гидрометеоиздат, 1971- С.295 309.

36. Дмитриев, М.Т., Соловьева Т.В., Китросский Н.А. О Физико-химическом механизме образования озона и фотооксидантов в атмосферном воздухе Текст. / Гигиена и санитария, 1971.-№ 10 — С.6 10.

37. Безуглая, Э.Ю. Метеорологический потенциал и климатические особенности загрязнения городов; Текст. / Э.Ю.-. Безуглая.-Л.: Гидрометеоиздат, 1980.-184 с.

38. Безуглая Э.Ю. Статистический метод оценки влияния! метеорологических условий на содержание примесей в атмосфере Текст. / Э:Ю. Безуглая , В.В. Клинго // Труды ГГО-1974.-Вып. 314 С.81 - 96.

39. Климатические характеристики условий распространения примесей в атмосфере. Справочное пособие Текст. /.-JI.: Гидрометеоиздат, 1983.-326 с.

40. Лайхтман, Д.Л. Динамическая: метеорология Текст. / Д.Л. Лайхтман.-Л.: Гидрометеоиздат, 1976. •42 0рленко, Л.Р. Строение планетарного пограничного слоя атмосферы Текст.;/ .ШР; Орленко:-Л ^Гидрометеоиздат, 1979::

41. Гаврилов, А.С. ZONE: следующий шаг Текст. / А.С. Гаврилов-СПб.:Изд. АО «ДЕЙТА», 1995.-32 с.-( Сер. "Библиотека пользователя ЭПК "ZONE").

42. Гаврилов, А.С. Интеллектуальная; геоинформационная система Текст. / А.С. Гаврилов // Бюллетень ГИС ассоциации.-1998.-Вып. 1(13).-С.58-59.

43. Гаврилов, А.С. Методика- расчета выбросов; в атмосферу загрязняющих; веществ автотранспортом на городских магистралях Текст. / А.С. Гаврилов.-М., 1996.-22 с.

44. Лайхтман, Д.Л. Физика пограничного слоя атмосферы Текст. / Д.Л. Лайхтман-Л.: Гидрометеоиздат, 1970

45. Матвеев, Л.Т. Курс общей метеорологии. Физика атмосферы Текст. / Л.Т. Матвеев.-Л.: Гидрометеоиздат, 1984.-751 с

46. Атмосферная турбулентность и моделирование распространения примесей Текст. / Под ред. Ф.Т.М. Ньистадт, Х.Ван Доп.-Л.: Гидрометеоиздат, 1985.-350.-(Пер. с англ.)

47. Лайхтман, Д.Л. Физика пограничного слоя атмосферы Текст. / Д.Л. Лайхтман-Л.: Гидрометеоиздат, 1970.-341 с.

48. Бин, Б.Р: Радиометеорология Текст. / Б.Р. Бин, Е.Дж. Даттон.-Л.: Гидрометеоиздат, 1971.-361 с.-(Пер.с англ.)

49. Гутман, Л.Н. Введение в нелинейную- теорию мезометеорологических процессов Текст. / Л.Н. Гутман—Л.: Гидрометеоиздат, 1969.-293 с.

50. Гаврилов, А.С. Математическое моделирование мезометеорологических процессов Текст. / А.С.Гаврилов.-Л.: ЛПИ, 1988.— 96 с.

51. Гаврилов, А.С. О влиянии радиации на режим приземного слоя атмосферы Текст. / А.С. Гаврилов, Д.Л. Лайхтман // Изв. АН СССР. Сер. ФАО.-1973—Т. 9.-С.27 33.

52. Гаврилов, А.С. Нестационарная задача о строении пограничного слоя атмосферы с учетом радиационного теплообмена Текст. / А.С. Гаврилов, Л.Н. Гутман, В.Н. Лыкосов // Труды Зап.-Сиб. РНИГМИ.-1974.-Вып.П.-С.35 48.

53. Экологический программный комплекс для ПЭВМ Текст. / Под ред.

54. A.С.Гаврилов.-СПб.: Гидрометеоиздат, 1992.-300 с.

55. Вызова, H.JI. Экспериментальные исследования атмосферной диффузии и расчеты рассеяния примеси Текст. / H.JI. Вызова, Е.К. Гаргер,

56. B.Н. Иванов—Л.:Гидрометеоиздат, 1991.-274 с

57. Гаврилов, А.С. Прогноз и климатологический анализ характеристик атмосферы, определяющих рассеяние антропогенных загрязнений Текст. / А.С. Гаврилов, А. Мханна, Р. Алтлули // Естественные и технические науки — М., 2008.-№ 6 .-С.221 -225.

58. Вагер, Б.Г. Сплайны при решении прикладных задач метеорологии и гидрологии Текст. / Б.Г. Вагер, Н.К. Серков.-Л.: Гидрометеоиздат, 1987159 с.

59. Bilgen S. Spline boundary element solution of dielectric scattering problems Text. / S. Bilgen, A Wexder- Ankara, Turkey: Department of Management, Middle East Technical University.