Лидарные исследования средней атмосферы в Якутии

Николашкин Семен Викторович; Титов Семен Вячеславович

УДК 53.01, 53.05

С. В. Николашкин, С. В. Титов

Лидарные исследования средней атмосферы в Якутии

Институт космофизических исследований и аэрономии им. Ю.Г. Шафера СО РАН,

г. Якутск, Россия

Аннотация. Среди всех известных методов мониторинга атмосферы, включая всевозможные методы прямых контактных измерений, а также методы активного и пассивного дистанционного зондирований, несомненным преимуществом обладают методы активного дистанционного зондирования с использованием лазерных источников излучения. Методы лазерного зондирования обеспечивают получение профилей или полей различных параметров атмосферы с исключительно высоким временным и пространственным разрешением, обладая при этом рекордными концентрационными чувствительностями. В статье представлены описание и основные характеристики якутского лидара, а также полученные на нем результаты по коэффициенту аэрозольного рассеяния R и температуре атмосферного слоя (25-60 км). В частности, в результате анализа коэффициентов аэрозольного рассеяния в ночь с 20 на 21 февраля 2013 г. был обнаружен аэрозольный слой на высоте примерно 39,5 км, который наблюдался на протяжении примерно 1,5 часов, толщина аэрозольного слоя составила 500 метров. Анализ карт траекторий воздушных масс для дат с 15 по 28 февраля, выполненный в Сибирском федеральном университете для уровней от 30 до 42 км с шагом 0,25 км по высоте, показал, что регистрируемый аэрозольный слой был вызван прохождением болида в атмосфере в районе Челябинска 15 февраля 2013 года. При анализе результатов наблюдений коэффициента аэрозольного рассеяния с января по март 2012 г., во время солнечных протонных событий (СПС), обнаружено аэрозольное наполнение на уровне стратосферы после СПС 23-24 января и на уровне мезосферы после события 7-8 марта. В статье исследуется поведение температуры стратосферы во время внезапных стратосферных потеплений (ВСП), наблюдавшихся в зимний период над Якутском. Для анализа и проверки достоверности были привлечены данные со спутника «Aura». Рассматриваются особенности формирования и распада ВСП. Посредством вейвлет анализа вариаций температуры выявлены основные квазимонохроматические волнообразные структуры (КВС) в вариациях температурного профиля на длинах волн 2-4 км, 5-7 км и 10-15 км.

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

DOI 10.25587/SVFU.2018.65.14071

Работа выполнена при финансовой поддержке Проекта Минобрнауки и СО РАН 11.16.2.1. 03752018-0004 (№ госрегистрации АААА-А17-117021450054-8).

НИКОЛАШКИН Семен Викторович - к. ф-м. н., зав. лабораторией оптики атмосферы Института космофизических исследований и аэрономии им. Ю.Г. Шафера СО РАН.

E-mail: nikolashkin@ikfia.ysn.ru

NIKOLASHKIN Semen Viktorovich - head of the Laboratory of Atmospheric Optics at the Yu.G. Shafer Institute of Cosmophysical Research and Aeronomy of SB RAS.

ТИТОВ Семен Вячеславович - м. н. с. лаборатории оптики атмосферы Института космофизических исследований и аэрономии им. Ю.Г. Шафера СО РАН.

E-MAIL: stitov@ikfia.ysn.ru

TITOV Semen Vyacheslavovich - junior research associate of the Atmospheric Optics Laboratory of Yu.G. Shafer Institute of Cosmophysical Research and Aeronomy of SB RAS.

S. V. Nikolashkin, S. V. Titov

Lidar Investigations of the Middle Atmosphere in Yakutia

Yu.G. Shafer Institute of Cosmophysical Research and Aeronomy of SB RAS,

Yakutsk, Russia

Abstract. Among all the known methods of atmospheric monitoring, including all possible methods of direct contact measurements, as well as those of active and passive remote sensing, the active remote sensing using laser sources of radiation have an undoubted advantage. Laser-sensing methods provide the production of profiles or fields of various atmospheric parameters with an exceptionally high temporal and spatial resolution, while possessing record concentration sensitivities. The main characteristics and description of the Yakutsk lidar, as well as the data on the aerosol scattering coefficient R, and the atmospheric layer temperature (25-60 km) are presented in this article. As a result of the aerosol scattering coefficients analysis on the night of February 20th to 21st, 2013, an aerosol layer was detected at an altitude of approximately 39.5 km which was observed for approximately 1.5 hours, the thickness of the aerosol layer being 500 meters. The analysis of maps of air mass trajectories for dates from February 15 to 28, performed at the Krasnoyarsk University for levels from 30 to 42 km, in 0.25 km height increments, showed that the detected aerosol layer was caused by the a car passage of in the atmosphere in Chelyabinsk region on February 15, 2013 . When analyzing the results of observations of the aerosol scattering coefficient R from January to March 2012 during solar proton events (SPS), aerosol filling at the level of the stratosphere after solar proton events on January 23-24 and at the level of the mesosphere after solar proton events of March 7th to 8th. The behavior of the stratosphere temperature during sudden stratospheric warming (SSW) observed during winter period over Yakutsk is investigated. For the analysis and validation data from the "Aura" satellite was used. Features of SSW formation and decay are considered, by means of the wavelet transformation of temperature variations, the main quasi-monochromatic wave-like structures (QWS) are revealed in the variations of the temperature profile at wavelengths of 2-4 km, 5-7 km and 10-15 km.

Keywords: lidar, aerosol scattering coefficient, Chelyabinsk bolide, sudden stratospheric warming, internal gravity waves, wavelet transformation.

The research was funded by the Ministry of Education and Sciences / SB RAS Project II.16.2.1. 03752018-0004 (state registration # ^^^17-117021450054-8).

Введение

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

Описание лидара

C 2004 г. в ИКФИА СО РАН (г. Якутск) на полигоне ШАЛ (61.66^ 129.37S) запущен

Рис. 1. Принципиальная схема лидара

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

В качестве передатчика используется твердотельный NdYAG лазер LQ-129B (фирма Solar LS, Минск) мощностью 300 мДж на второй гармонике с длиной волны 532 нм, с длительностью импульса 13 нс в моноимпульсном режиме, работающий на частоте 20 Гц. Приемником служит направленный в зенит телескоп системы Ньютона с диаметром главного зеркала 0,6 м и фокусным расстоянием 2 м. Пространственное разрешение системы от 75 до 300 м в зависимости от задачи. Регистрация ведется на компьютере через специальный блок счета фотонов, подсоединяемый через USB 2.0 порт, при помощи интерфейсной программы, написанной с использованием пакета виртуальных приборов LabView 6.1.

За период наблюдений с 2004 по 2016 гг. проведено около 2000 сеансов лидарного зондирования средней атмосферы. Наблюдения проводятся в ночное время в безоблачную погоду.

Аэрозольные исследования

Присутствие аэрозолей в атмосфере играет чрезвычайно важную роль в формировании климата всей Земли. Их наличие существенно влияет на термический режим атмосферы и земной поверхности. Известно, что количество приходящей солнечной радиации зависит от состава атмосферы и прежде всего от наличия аэрозолей. Аэрозоли поглощают и рассеивают коротковолновую солнечную радиацию, уменьшая тем самым величину потока, приходящего на подстилающую поверхность. Это приводит к понижению температуры подстилающей поверхности [1-3]. Актуальной задачей является исследование аэрозольных компонентов, их преобразований в атмосфере, определение

Таблица

Основные характеристики лидара

Тип лазера Nd:YAG Расходимость луча, мрад ~40

Частота, Гц 20 Телескоп Система Ньютона

Длина волны,нм 532 Диаметр зеркала, м 0.6

Выходная энергия, мДж 300 Фокусное расстояние, м 2

Длительность импульса, нс 10-12 Поле зрения, мрад 0.2

Диаметр луча/после расширителя луча, мм 6/60 ФЭУ Hamamatsu Н8259-01

Рис. 2. Профили коэффициента аэрозольного рассеяния 20.03.13 г. по лидарным измерениям в Якутске. Вертикальные линии сетки соответствуют R=1, и каждое измерение смещено на 1 для удобства восприятия

Рис. 3. Карты траекторий воздушных масс с 20 февраля 2013 г., построены суточные отрезки траекторий за эту дату от 0 до 24 часов, цветом показаны высоты траекторий

механизмов диффузии в пространстве и дальности распространения, уровня концентрации и оптического действия [4, 5]. К настоящему времени накоплен большой объем информации о физико-химических свойствах атмосферного аэрозоля, фундаментальных закономерностях его образования и трансформации, пространственно-временной изменчивости полей оптических и микрофизических характеристик [6]. Вместе с тем для совершенствования моделей, учитывающих динамичность процессов, проходящих в атмосфере с участием аэрозоля, изучения особенностей формирования и трансформации аэрозольных слоев и их изменчивости в зависимости от глобальных, региональных и локальных геофизических факторов требуются систематические данные на основе долговременных наблюдений [7, 8].

В ходе обработки данных лидарного зондирования за 20 февраля 2013 г. на высоте 39-40 км был обнаружен аэрозольный слой с коэффициентом аэрозольного рассеяния R~1,5. Слой отчетливо наблюдался на 5 профилях, полученных с временным интервалом 10 минут, начиная с 23:41 до 00:21 ЦТ (Ш+10), его нет на профилях в 20:34, на 03:02 и последующих. Максимальное значение коэффициента аэрозольного рассеяния составляет 1,69 на высоте 39,55 км в 23:51, а минимальное значение коэффициента аэрозольного рассеяния составило 1,36 на высоте 39,25 км в 00:21 (рис. 2 ).

Нашими коллегами из Сибирского федерального университета (г. Красноярск) были произведены расчеты траекторий воздушных масс для интервала дат с 15 по 28 февраля, начатых с момента наблюдения образования метеорного следа, 15 февраля в 3:20 иТ на высотах от 30 до 42 км, с шагом 0,25 км по высоте. Переносимый атмосферной циркуляцией метеорный след достиг в первый раз Якутска 15-16 февраля, следующий раз при втором витке - 20-21 февраля (рис. 3) и на третьем витке - 24-25 февраля.

15-16 февраля лидарное зондирование не проводилось, а 24-25 февраля в данных аэрозольный слой не выявлен. Также нами для проверки гипотезы рассматривались и другие возможные варианты происхождения аэрозольного слоя, такие как вулканическое и солнечное, но до 20 февраля ни одного события солнечной активности и вулканического извержения не было. Следует также отметить тот факт, что на протяжении 10 лет лидарных зондирований в Якутске на высоте 40 км аэрозольные слои не наблюдались. Это подтверждает, что аэрозольный слой был вызван прохождением болида через атмосферу [9].

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

■2012.01.26_21.0В —2D12.01.27_19.09 —2012.01.2S.0633 —201J.i8i6_il.21 — 2012 ,»КП_03.00 —2012i3,W_21,13 — 20u03.l5_0i.13

2012.01.29.06.44 —2012.01.29.20.50 —2012.01.30.19,30 —W12.03.15_00.O8 —201J.03.15_21.29 —2012.03.16JB.3)

Рис. 4. Профили коэффициента аэрозольного Рис. 5. Профили коэффициента аэрозольного рассеяния с 26 по 30 января 2012 г. по рассеяния 6, 7, 14, 15, 16 марта 2012 г. по

лидарным измерениям в Якутске лидарным измерениям в Якутске

стратосферного лидара было обнаружено аэрозольное наполнение на уровне стратосферы после СПС 23-24 января и на уровне мезосферы после события 7-8 марта (рис. 4 и 5). Такая высотная разность в образовании аэрозоля, по-видимому, связана с разностью высоты ионизации атмосферы заряженными частицами, в последующем эти заряженные частицы послужили ядрами конденсации для образования аэрозоля [10-14].

Исследования температуры средней атмосферы во время внезапных зимних стратосферных потеплений

В настоящее время основной и наиболее информативной характеристикой атмосферы является ее температура. Исследование температуры разных слоев атмосферы и последующее построение температурного профиля атмосферы в каждом районе земного шара дает нам представление о химических и физических процессах, которые протекают на различных высотах атмосферы в различных районах Земли. Зная эти процессы, мы можем проследить общую динамику атмосферы, и с той или иной долей вероятности прогнозировать ее поведение [15, 16]. С этой точки зрения Дальний Восток России, в частности Якутия, представляет собой огромный интерес из-за его климатических особенностей.

Полярная и субполярная зимняя стратосфера северного полушария характеризуется весьма динамичным поведением, связанным с взаимодействием распространяющихся в тропосфере планетарных волн со стоячей волной орографического происхождения, также средним зональным потоком. При определенных условиях это взаимодействие приводит к явлению т.н. внезапного стратосферного потепления (ВСП). ВСП характеризуются резким повышением температуры в стратосфере и наблюдаются в зимний период в полярных и субполярных зонах. В зависимости от величины потепления и наличия обращения направления среднего зонального ветра с западного на восточное его подразделяют на типы «major» и «minor» [16-18].

В период с 2004 по 2016 гг. нами было рассмотрено 10 ВСП. Их длительность составила от 3 дней до 7 дней с повышением температуры в очаге потепления до 50 К. Во время всех ВСП очаг потепления опускался со скоростью примерно 1 км в сутки и, как следствие, стратопауза опускалась до 35 км. Во время ВСП в январе 2010 и 2012 гг. наблюдается сдвоение стратопаузы.

Для сравнения данных нами были построен высотно-временной ход температуры в зимний период 2011-2012 гг. по данным лидара и Microwave Limb Sounder (MLS) на

Рис. 6. Высотно-временной ход температуры по данным лидара (слева) и прибора MLS (справа)

спутнике «Aura» (рис. 6). Из рисунка 6 видно, что измерения обоих приборов в основном схожи, а небольшое разногласие вызвано некоторым несовпадением места измерения спутником с местом зондирования лидара и самим принципом измерения температуры слоя атмосферы спутником.

В настоящее время хорошо известна значимость внутренних гравитационных волн (ВГВ) в глобальной циркуляции, климатологии и структуре атмосферы [19]. Основными источниками ВГВ являются орография, конвекция и дисбаланс в системах струйных течений и фронтов [5]. Существующие в широком диапазоне горизонтальных масштабов и собственных частот [15] ВГВ и их влияние на средний поток представляют собой серьезную проблему для моделей общей циркуляции атмосферы. Из-за недостаточного пространственного разрешения эти модели не могут явно описывать полный спектр волн и, следовательно, нуждаются в параметризации для представления их эффектов. Недавнее исследование [20] выявило большие различия между существующими моделями и наблюдениями и указало на необходимость проведения дополнительных исследований.

Для исследования волновых процессов нами было использовано непрерывное вейвлет-преобразование с помощью материнского вейвлета Гауссиана 8-го порядка [20]. Флуктуации профиля температуры, связанные с внутренними гравитационными волнами, были получены путем сглаживания температурного профиля атмосферы скользящим окном 2 км для ВЧ фильтрации и полиномом 8-го порядка для НЧ фильтрации (рис. 7) Численные расчеты выполнялись с использованием пакета прикладных программ MATLAB (wavelet toolbox).

О 100 И» ЗСО «О 500 бОО ТОО BOO

_П2_

О 100 200 300 <*СЮ 500 600 700 900

Рис. 7. Пример сглаживания полиномом седьмого порядка температурных данных лидара, полученные 25.01.10 г. 21:15

Рис. 8. Графики непрерывного вейвлет-преобразования с материнским вейвлетом «Гауссиан 8-го порядка»: 16.10.2014 г. 22:36, 23.11.2014 г. 21:50

Результаты вейвлет-преобразования представлены в виде поверхности в декартовой системе координат, где по оси Х представлена высота в км, по оси Y - длина волны в км и по оси Ъ - амплитуда вейвлет - коэффициента (рис. 8). Из графиков можно выделить основные квазимонохроматические волнообразные структуры (КВС) с вертикальными длинами волн 2-4 км, 5-7 км и 10-15 км [21]. Предполагается, что эти флуктуации связаны с вертикальным распространением внутренних гравитационных волн. Амплитуды этих КВС, как и ожидалось, увеличиваются с повышением высоты из-за уменьшения плотности атмосферы. Во время ВСП амплитуды КВС резко возрастают, возможно, это связано с обращением среднего зонального ветра. Основное изменение длины волны квазимонохроматических волнообразных структур происходит на высоте стратопаузы. По-видимому, это связано с изменением градиента температуры в стратопаузе, т.е. на уровне стратопаузы образуется критический уровень, на котором в зависимости от направления среднего зонального ветра планетарные волны затухают с выделением энергии в виде тепла. Стоит отметить, что выявленные нами основные КВС на длинах волн 2-4 км, 5-7 км и 10-15 км хорошо согласуются с подобными лидарными исследованиями, проводимыми за рубежом [22].

Заключение

В данной работе мы рассмотрели результаты лидарного исследования средней атмосферы в Якутии. К ним относятся исследования аэрозольных слоев и ВСП в средней атмосфере. Как показали аэрозольные исследования, 20 февраля 2013 г. был обнаружен аэрозольный слой на высоте 39-40 км с коэффициентом аэрозольного рассеяния R~1,5, который образовался в результате пролета через толщу атмосферы Челябинского метеорита. Показано образование аэрозольного наполнения на уровне мезосферы и стратосферы после солнечных протонных событий 23-24 января и 7-8 марта 2012 года. Установлено, что высота аэрозольного наполнения зависит от энергетического спектра солнечных протонов. Исследования температурного профиля средней атмосферы в зимний период показали, что каждую зиму над Якутском наблюдаются ВСП, их

длительность составила от 3 дней до 7 дней с повышением среднемесячной температуры в очаге потепления до 50 К. Во время всех ВСП очаг потепления опускался со скоростью примерно 1 км в сутки. При помощи вейвлет-анализа вариаций температурного профиля атмосферы вычислены основные КВС с длиной волны 2-4 км, 5-7 км и 10-15 км. Показано, что основное затухание КВС происходит на уровне стратопаузы.

Л и т е р а т у р а

1. Розенберг Г. В., Горчаков Г. И., Георгиевский Ю. С. и др. Оптические параметры атмосферного аэрозоля // Физика атмосферы и проблемы климата / Под ред. Голицына Г. С. - М.: Наука, 1980. - С. 216-256.

2. Розенберг Г. В., Сандомирский А. Б. Оптическая стратификация атмосферного аэрозоля // Изв. АН СССР. Физика атмосферы и океана. - 1971. - Т. 7, № 7. - С. 737-749.

3. Shipley Scott T., Tracy D. H., Eloranta Edwin W., Trauger John T., Sroga J. T., Roesler F. L. and Weinman James A. "High spectral resolution lidar to measure optical scattering properties of atmospheric aerosols. 1: Theory and instrumentation." Applied optics 22, no. 23. - 1983: 3716-3724. 1.

4. McNeill V. Faye. "Atmospheric Aerosols: Clouds, Chemistry, and Climate." Annual review of chemical and biomolecular engineering 8. - 2017: 427-444.

5. Plougonven R. and Zhang F., Internal gravity waves from atmospheric jets and fronts, Rev. Geophys. 52. - 2014, doi:10.1002/2012RG000419.

6. Sun, Tianze, Huizheng Che, Bing Qi, Yaqiang Wang, Yunsheng Dong, Xiangao Xia, Hong Wang et al. "Aerosol optical characteristics and their vertical distributions under enhanced haze pollution events: effect of the regional transport of different aerosol types over eastern China." Atmos. Chem. Phys. Discuss. - 2017, https://doi. org/10.5194/acp-2017-805.

7. Di Biagio, Claudia, Jacques Pelon, Gérard Ancellet, Ariane Bazureau, and Vincent Mariage. "Sources, Load, Vertical Distribution, and Fate of Wintertime Aerosols North of Svalbard From Combined V4 CALIOP Data, Ground-Based IAOOS Lidar Observations and Trajectory Analysis." Journal of Geophysical Research: Atmospheres 2018, 123, no. 2: 1363-1383.

8. Verichev Konstantin, Manuel Carpio, and Cristian Salazar. "Effects of El Niño/La Niña on the climatic values of aerosol optical thickness for the Valdivia city, Chile." In Geoscience and Remote Sensing (GRSS-CHILE), 2017 First IEEE International Symposium of, pp. 1-4. IEEE, 2017.

9. Иванов В. Н., Зубачев Д. С., Коршунов В. А., Лапшин В. Б., Иванов М. С., Галкин К. А., Губко П. А., Антонов Д. Л., Тулинов Г. Ф., Черемисин А. А., Новиков П. В., Николашкин С. В., Титов С. В., Маричев В. Н. Лидарные наблюдения стратосферных аэрозольных следов от челябинского метеорита. // Оптика атмосферы и океана. - 2014. - Т. 27, № 2 (301). - С. 117-122.

10. Arnold F. Atmospheric aerosol and cloud condensation nuclei formation: A possible influence of cosmic rays? Space Sci. Rev. - 2006, 125, 169-186.

11. Svensmark H., Pedersen J. O. P., Marsh N. D., Enghoff M. B. and Uggerhj U. I. Experimental evidence for the role of ions in particle nucleation under atmospheric conditions, Proc. R. Soc. London, Ser. A, 463. - 2007,385-396, doi:10.1098/rspa.2006.1773.

12. Tinsley B. A. Correlations of atmospheric dynamics with solar wind-induced changes of air-earth current density into cloud tops, J. Geophys. Res., 101, 29. - 1996, 701- 29,714.

13. Tinsley B. A. Influence of the solar wind on the global electric circuit, and inferred effects on cloud microphysics, temperature, and dynamics of the troposphere, Space Science Reviews, 94. - 2000, 231-258.

14. Tinsley B. A., Yu F. Atmospheric ionization and clouds as links between solar activity and climate, pp. 321-340 in Solar Variability and its Effects on Climate, AGU Monograph 141, ed. J. Pap et al., 2004, AGU press, Washington, DC.

15. Holton J. R., An Introduction to Dynamic Meteorology. - 1992, 3rd ed., Academic Press, London.

16. Labitzke K. Stratospheric-mesospheric midwinter disturbances: A summary of observed characteristics // J. Geophys. Res. - 1981. - Vol. 86. - P. 9665-9678. 1.

17. Matsuno T., A dynamical model of the stratospheric sudden warming. // Journal of Atmospheric Sciences. - 1971 ,28(6), 1479-1493.

18. VonZahn U., Fiedler J., Naujokat B., Jangematz U., Kruger K. A note on record-hight temperatures at the northern polar stratopause in winter 1997/98 // Geophys. Res. Lett. - 1998. - V. 25, No. 22. - P. 4169-4172.

19. Alexander M. J., Geller M., McLandress C., Polavarapu S., Preusse, P., Sassi F., Sato K., Eckermann S., Ern M., Hertzog A., Kawatani Y., Pulido M., Shaw T. A., Sigmond M., Vincent R. and Watanabe S., Recent developments in gravity-wave effects in climate models and the global distribution of gravity-wave momentum flux from observations and models. Q.J.R. Meteorol. Soc., - 2010, 136: 1103-1124. doi: 10.1002/ qj.637

20. Geller M. et al., A comparison between gravity wave momentum fluxes in observations and climate models, J. Clim. - 2013, 26, 6383-6405, doi:10.1175/JCLI-D-12-00545.1.

21. Werner R. et al. Application of wavelet transformation to determine wavelengths and phase velocities of gravity waves observed by lidar measurements //Journal of Atmospheric and Solar-Terrestrial Physics. - 2007. - T. 69, №. 17. - C. 2249-2256.

22. Chane-Ming F., Molinaro F., Leveau J., Keckhut P. and Hauchecorne A. Analysis of gravity waves in the tropical middle atmosphere with lidar using wavelet techniques, Ann. Geophys. - 2000, 18, 485-498.

R e f e r e n c e s

1. Rozenberg G. V., Gorchakov G. I., Georgievskij YU. S. i dr. Opticheskie parametry atmosfernogo aehrozolya // Fizika atmosfery i problemy klimata / Pod red. Golicyna G. S. - M.: Nauka, 1980. - S. 216-256.

2. Rozenberg G. V., Sandomirskij A. B. Opticheskaya stratifikaciya atmosfernogo aehrozolya // Izv. AN SSSR. Fizika atmosfery i okeana. - 1971. - T. 7, № 7. - S. 737-749.